netgen/libsrc/csg/extrusion.hpp

156 lines
4.1 KiB
C++
Raw Normal View History

2009-01-13 04:40:13 +05:00
#ifndef _EXTRUSION_HPP
#define _EXTRUSION_HPP
2009-09-07 17:50:13 +06:00
namespace netgen
{
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
class Extrusion;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
class ExtrusionFace : public Surface
{
private:
const SplineSeg<2> * profile;
const SplineGeometry<3> * path;
Vec<3> glob_z_direction;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
bool deletable;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
Array< const SplineSeg3<3> * > spline3_path;
Array< const LineSeg<3> * > line_path;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
mutable Array < Vec<3> > x_dir, y_dir, z_dir, loc_z_dir;
mutable Array < Point<3> > p0;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
mutable Vec<3> profile_tangent;
mutable double profile_par;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
mutable Vector profile_spline_coeff;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
mutable int latest_seg;
mutable double latest_t;
mutable Point<2> latest_point2d;
mutable Point<3> latest_point3d;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
private:
void Orthogonalize(const Vec<3> & v1, Vec<3> & v2) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
void Init(void);
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
public:
double CalcProj(const Point<3> & point3d, Point<2> & point2d,
int seg) const;
void CalcProj(const Point<3> & point3d, Point<2> & point2d,
int & seg, double & t) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
public:
ExtrusionFace(const SplineSeg<2> * profile_in,
const SplineGeometry<3> * path_in,
const Vec<3> & z_direction);
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
ExtrusionFace(const Array<double> & raw_data);
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
~ExtrusionFace();
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
virtual int IsIdentic (const Surface & s2, int & inv, double eps) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
virtual double CalcFunctionValue (const Point<3> & point) const;
virtual void CalcGradient (const Point<3> & point, Vec<3> & grad) const;
virtual void CalcHesse (const Point<3> & point, Mat<3> & hesse) const;
virtual double HesseNorm () const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
virtual double MaxCurvature () const;
//virtual double MaxCurvatureLoc (const Point<3> & /* c */ ,
// double /* rad */) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
virtual void Project (Point<3> & p) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
virtual Point<3> GetSurfacePoint () const;
virtual void Print (ostream & str) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
virtual void GetTriangleApproximation (TriangleApproximation & tas,
const Box<3> & boundingbox,
double facets) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
const SplineGeometry<3> & GetPath(void) const {return *path;}
const SplineSeg<2> & GetProfile(void) const {return *profile;}
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
bool BoxIntersectsFace(const Box<3> & box) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
void LineIntersections ( const Point<3> & p,
const Vec<3> & v,
const double eps,
int & before,
int & after,
bool & intersecting ) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
INSOLID_TYPE VecInFace ( const Point<3> & p,
const Vec<3> & v,
const double eps ) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
const Vec<3> & GetYDir ( void ) const {return y_dir[latest_seg];}
const Vec<3> & GetProfileTangent (void) const {return profile_tangent;}
double GetProfilePar(void) const {return profile_par;}
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
void GetRawData(Array<double> & data) const;
2009-04-20 03:15:26 +06:00
2009-09-07 17:50:13 +06:00
void CalcLocalCoordinates (int seg, double t,
Vec<3> & ex, Vec<3> & ey, Vec<3> & ez) const;
2009-04-20 03:15:26 +06:00
2009-09-07 17:50:13 +06:00
void CalcLocalCoordinatesDeriv (int seg, double t,
Vec<3> & ex, Vec<3> & ey, Vec<3> & ez,
Vec<3> & dex, Vec<3> & dey, Vec<3> & dez) const;
2009-04-20 03:15:26 +06:00
2009-09-07 17:50:13 +06:00
};
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
class Extrusion : public Primitive
{
private:
const SplineGeometry<3> & path;
const SplineGeometry<2> & profile;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
const Vec<3> & z_direction;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
Array<ExtrusionFace*> faces;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
mutable int latestfacenum;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
public:
Extrusion(const SplineGeometry<3> & path_in,
const SplineGeometry<2> & profile_in,
const Vec<3> & z_dir);
~Extrusion();
virtual INSOLID_TYPE BoxInSolid (const BoxSphere<3> & box) const;
virtual INSOLID_TYPE PointInSolid (const Point<3> & p,
double eps) const;
INSOLID_TYPE PointInSolid (const Point<3> & p,
double eps,
Array<int> * const facenums) const;
virtual INSOLID_TYPE VecInSolid (const Point<3> & p,
const Vec<3> & v,
2009-01-13 04:40:13 +05:00
double eps) const;
2009-09-07 17:50:13 +06:00
// checks if lim s->0 lim t->0 p + t(v1 + s v2) in solid
virtual INSOLID_TYPE VecInSolid2 (const Point<3> & p,
const Vec<3> & v1,
const Vec<3> & v2,
double eps) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
virtual int GetNSurfaces() const;
virtual Surface & GetSurface (int i = 0);
virtual const Surface & GetSurface (int i = 0) const;
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
virtual void Reduce (const BoxSphere<3> & box);
virtual void UnReduce ();
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
};
2009-01-13 04:40:13 +05:00
2009-09-07 17:50:13 +06:00
}
2009-01-13 04:40:13 +05:00
#endif //_EXTRUSION_HPP