netgen/libsrc/occ/occmeshsurf.hpp

204 lines
4.5 KiB
C++
Raw Normal View History

2009-01-13 23:32:10 +05:00
#ifdef OCCGEOMETRY
#ifndef FILE_OCCMESHSURF
#define FILE_OCCMESHSURF
#include "occgeom.hpp"
#define PARAMETERSPACE -1
#define PLANESPACE 1
class OCCGeometry;
class SingularMatrixException
{};
class UVBoundsException
{};
class OCCSurface
{
public:
TopoDS_Face topods_face;
Handle(Geom_Surface) occface;
TopAbs_Orientation orient;
int projecttype;
protected:
Point<3> p1;
Point<3> p2;
/// in plane, directed p1->p2
Vec<3> ex;
/// in plane
Vec<3> ey;
/// outer normal direction
Vec<3> ez;
/// normal vector in p2
Vec<3> n2;
/// average normal vector
Vec<3> nmid;
// for transformation to parameter space
Point<2> psp1;
Point<2> psp2;
Vec<2> psex;
Vec<2> psey;
Mat<2,2> Amat, Amatinv;
// UV Bounds
double umin, umax, vmin, vmax;
public:
OCCSurface (const TopoDS_Face & aface, int aprojecttype)
{
topods_face = aface;
occface = BRep_Tool::Surface(topods_face);
orient = topods_face.Orientation();
projecttype = aprojecttype;
ShapeAnalysis::GetFaceUVBounds (topods_face, umin, umax, vmin, vmax);
umin -= fabs(umax-umin)/100.0;
vmin -= fabs(vmax-vmin)/100.0;
umax += fabs(umax-umin)/100.0;
vmax += fabs(vmax-vmin)/100.0;
// projecttype = PLANESPACE;
/*
TopExp_Explorer exp1;
exp1.Init (topods_face, TopAbs_WIRE);
orient = TopAbs::Compose (orient, exp1.Current().Orientation());
*/
};
~OCCSurface()
{};
void Project (Point<3> & p, PointGeomInfo & gi);
void GetNormalVector (const Point<3> & p,
const PointGeomInfo & geominfo,
Vec<3> & n) const;
/**
Defines tangential plane in ap1.
The local x-coordinate axis point to the direction of ap2 */
void DefineTangentialPlane (const Point<3> & ap1,
const PointGeomInfo & geominfo1,
const Point<3> & ap2,
const PointGeomInfo & geominfo2);
/// Transforms 3d point p3d to local coordinates pplane
void ToPlane (const Point<3> & p3d, const PointGeomInfo & geominfo,
Point<2> & pplane, double h, int & zone) const;
/// Transforms point pplane in local coordinates to 3d point
void FromPlane (const Point<2> & pplane,
Point<3> & p3d,
PointGeomInfo & gi,
double h);
};
///
class Meshing2OCCSurfaces : public Meshing2
{
///
OCCSurface surface;
public:
///
Meshing2OCCSurfaces (const TopoDS_Shape & asurf, const Box<3> & aboundingbox, int aprojecttype);
///
int GetProjectionType ()
{ return surface.projecttype; }
protected:
///
virtual void DefineTransformation (const Point3d & p1, const Point3d & p2,
const PointGeomInfo * geominfo1,
const PointGeomInfo * geominfo2);
///
virtual void TransformToPlain (const Point3d & locpoint,
const MultiPointGeomInfo & geominfo,
Point2d & plainpoint,
double h, int & zone);
///
virtual int TransformFromPlain (Point2d & plainpoint,
Point3d & locpoint,
PointGeomInfo & gi,
double h);
///
virtual double CalcLocalH (const Point3d & p, double gh) const;
};
///
class MeshOptimize2dOCCSurfaces : public MeshOptimize2d
{
///
const OCCGeometry & geometry;
public:
///
MeshOptimize2dOCCSurfaces (const OCCGeometry & ageometry);
///
virtual void ProjectPoint (INDEX surfind, Point<3> & p) const;
///
virtual void ProjectPoint2 (INDEX surfind, INDEX surfind2, Point<3> & p) const;
///
virtual int ProjectPointGI (INDEX surfind, Point<3> & p, PointGeomInfo & gi) const;
///
virtual void GetNormalVector(INDEX surfind, const Point<3> & p, Vec<3> & n) const;
///
virtual void GetNormalVector(INDEX surfind, const Point<3> & p, PointGeomInfo & gi, Vec<3> & n) const;
virtual int CalcPointGeomInfo(int surfind, PointGeomInfo& gi, const Point<3> & p3) const;
};
class OCCGeometry;
class OCCRefinementSurfaces : public Refinement
{
const OCCGeometry & geometry;
public:
OCCRefinementSurfaces (const OCCGeometry & ageometry);
virtual ~OCCRefinementSurfaces ();
virtual void PointBetween (const Point<3> & p1, const Point<3> & p2, double secpoint,
int surfi,
const PointGeomInfo & gi1,
const PointGeomInfo & gi2,
Point<3> & newp, PointGeomInfo & newgi);
virtual void PointBetween (const Point<3> & p1, const Point<3> & p2, double secpoint,
int surfi1, int surfi2,
const EdgePointGeomInfo & ap1,
const EdgePointGeomInfo & ap2,
Point<3> & newp, EdgePointGeomInfo & newgi);
virtual void ProjectToSurface (Point<3> & p, int surfi);
virtual void ProjectToSurface (Point<3> & p, int surfi, PointGeomInfo & gi);
};
#endif
#endif