netgen-mesher 6.2.2506.post35.dev0__cp314-cp314-win_amd64.whl

netgen/include/csg/explicitcurve2d.hpp

@@ -0,0 +1,113 @@
+#ifndef FILE_EXPLICITCURVE2D
+#define FILE_EXPLICITCURVE2D
+
+/**************************************************************************/
+/* File:   explicitcurve2d.hh                                             */
+/* Author: Joachim Schoeberl                                              */
+/* Date:   14. Oct. 96                                                    */
+/**************************************************************************/
+
+
+namespace netgen
+{
+
+  /*
+
+  Explicit 2D Curve representation
+
+  */
+
+
+
+  ///
+  class ExplicitCurve2d : public Curve2d
+  {
+  public:
+    ///
+    ExplicitCurve2d ();
+
+    ///
+    virtual void Project (Point<2> & p) const;
+    ///
+    virtual double ProjectParam (const Point<2> & p) const = 0;
+    ///
+    virtual double NumericalProjectParam (const Point<2> & p, double lb, double ub) const;
+    ///
+    virtual double MinParam () const = 0;
+    ///
+    virtual double MaxParam () const = 0;
+    ///
+    virtual Point<2> Eval (double t) const = 0;
+    ///
+    virtual Vec<2> EvalPrime (double t) const = 0;
+    ///
+    virtual Vec<2> Normal (double t) const;
+    ///
+    virtual void NormalVector (const Point<2> & p, Vec<2> & n) const;
+    ///
+    virtual Vec<2> EvalPrimePrime (double t) const = 0;
+
+    ///
+    virtual double MaxCurvature () const;
+    ///
+    virtual double MaxCurvatureLoc (const Point<2> & p, double rad) const;
+
+    ///
+    virtual Point<2> CurvCircle (double t) const;
+    ///
+    virtual void Print (ostream & /* str */) const { };
+  
+    ///
+    virtual int SectionUsed (double /* t */) const { return 1; }
+    ///
+    virtual void Reduce (const Point<2> & /* p */, double /* rad */) { };
+    ///
+    virtual void UnReduce () { };
+  }; 
+  
+  
+  ///
+  class BSplineCurve2d : public ExplicitCurve2d
+  {
+    ///
+    NgArray<Point<2> > points;
+    ///
+    NgArray<int> intervallused;
+    ///
+    int redlevel;
+  
+  public:
+    ///
+    BSplineCurve2d ();
+    ///
+    void AddPoint (const Point<2> & apoint);
+
+    bool Inside (const Point<2> & p, double & dist) const;
+  
+    ///
+    virtual double ProjectParam (const Point<2> & p) const;
+    ///
+    virtual double MinParam () const { return 0; }
+    ///
+    virtual double MaxParam () const { return points.Size(); }
+    ///
+    virtual Point<2> Eval (double t) const;
+    ///
+    virtual Vec<2> EvalPrime (double t) const;  
+    ///
+    virtual Vec<2> EvalPrimePrime (double t) const;
+    ///
+    virtual void Print (ostream & str) const;
+
+    ///
+    virtual int SectionUsed (double t) const;
+    ///
+    virtual void Reduce (const Point<2> & p, double rad);
+    ///
+    virtual void UnReduce ();
+  };  
+
+}
+
+
+#endif

netgen/include/csg/extrusion.hpp

@@ -0,0 +1,185 @@
+#ifndef _EXTRUSION_HPP
+#define _EXTRUSION_HPP
+
+namespace netgen
+{
+
+  class Extrusion;
+
+  class ExtrusionFace : public Surface
+  {
+  private:
+    const SplineSeg<2> * profile;
+    const SplineGeometry<3> * path;
+    Vec<3> glob_z_direction;
+    Array<double> angles;
+
+    bool deletable;
+    int tangential_plane_seg;
+  
+    NgArray< const SplineSeg3<3> * > spline3_path;
+    NgArray< const LineSeg<3> * > line_path;
+  
+    mutable NgArray < Vec<3> > x_dir, y_dir, z_dir, loc_z_dir;
+    mutable NgArray < Point<3> > p0;
+
+    mutable Vec<3> profile_tangent;
+    mutable double profile_par;
+  
+    mutable Vector profile_spline_coeff;
+
+    mutable int latest_seg;
+    mutable double latest_t;
+    mutable Point<2> latest_point2d;
+    mutable Point<3> latest_point3d;
+
+
+  private:
+    void Orthogonalize(const Vec<3> & v1, Vec<3> & v2) const;
+
+    void Init(void);
+
+  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;
+
+  public:
+    ExtrusionFace(const SplineSeg<2> * profile_in,
+		  const SplineGeometry<3> * path_in,
+		  const Vec<3> & z_direction);
+
+    ExtrusionFace(const NgArray<double> & raw_data);
+    // default constructor for archive
+    ExtrusionFace() {}
+
+    ~ExtrusionFace();
+
+    void DoArchive(Archive& ar) override
+    {
+      Surface::DoArchive(ar);
+      ar & profile & path & glob_z_direction & deletable & spline3_path & line_path &
+        x_dir & y_dir & z_dir & loc_z_dir & p0 & profile_tangent & profile_par &
+        profile_spline_coeff & latest_seg & latest_t & latest_point2d & latest_point3d;
+    }
+  
+    int IsIdentic (const Surface & s2, int & inv, double eps) const override;
+  
+    double CalcFunctionValue (const Point<3> & point) const override;
+    void CalcGradient (const Point<3> & point, Vec<3> & grad) const override;
+    void CalcHesse (const Point<3> & point, Mat<3> & hesse) const override;
+    double HesseNorm () const override;
+
+    double MaxCurvature () const override;
+    //virtual double MaxCurvatureLoc (const Point<3> & /* c */ , 
+    //				  double /* rad */) const;
+
+    void Project (Point<3> & p) const override;
+
+    Point<3> GetSurfacePoint () const override;
+    void Print (ostream & str) const override;
+  
+    void GetTriangleApproximation (TriangleApproximation & tas, 
+					   const Box<3> & boundingbox, 
+					   double facets) const override;
+
+    const SplineGeometry<3> & GetPath(void) const {return *path;}
+    const SplineSeg<2> & GetProfile(void) const {return *profile;}
+
+    bool BoxIntersectsFace(const Box<3> & box) const;
+
+    void LineIntersections ( const Point<3> & p,
+			     const Vec<3> & v,
+			     const double eps,
+			     int & before,
+			     int & after,
+			     bool & intersecting ) const;
+
+
+    bool PointInFace (const Point<3> & p, const double eps) const;
+    
+    INSOLID_TYPE VecInFace ( const Point<3> & p,
+			     const Vec<3> & v,
+			     const double eps ) const;
+
+    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;}
+
+    void GetRawData(NgArray<double> & data) const;
+
+    void CalcLocalCoordinates (int seg, double t, 
+			       Vec<3> & ex, Vec<3> & ey, Vec<3> & ez) const;
+
+    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;
+
+    void DefineTangentialPlane(const Point<3>& ap1,
+                               const Point<3>& ap2) override;
+    void ToPlane(const Point<3>& p3d, Point<2>& p2d,
+                 double h, int& zone) const override;
+
+  };
+
+
+
+  class Extrusion : public Primitive
+  {
+  private:
+    shared_ptr<SplineGeometry<3>> path;
+    shared_ptr<SplineGeometry<2>> profile; // closed, clockwise oriented curve
+
+    Vec<3> z_direction;
+
+    NgArray<ExtrusionFace*> faces;
+
+    mutable int latestfacenum;
+
+  public:
+    Extrusion(shared_ptr<SplineGeometry<3>> path_in,
+	      shared_ptr<SplineGeometry<2>> profile_in,
+	      const Vec<3> & z_dir);
+    // default constructor for archive
+    Extrusion() {}
+    ~Extrusion();
+
+    void DoArchive(Archive& ar) override
+    {
+      Primitive::DoArchive(ar);
+      ar & path & profile & z_direction & faces & latestfacenum;
+    }
+    INSOLID_TYPE BoxInSolid (const BoxSphere<3> & box) const override;
+    INSOLID_TYPE PointInSolid (const Point<3> & p,
+                               double eps) const override;
+    INSOLID_TYPE PointInSolid (const Point<3> & p,
+			       double eps,
+			       NgArray<int> * const facenums) const;
+
+    void GetTangentialSurfaceIndices (const Point<3> & p, 
+                                      NgArray<int> & surfind, double eps) const override;
+
+    INSOLID_TYPE VecInSolid (const Point<3> & p,
+                             const Vec<3> & v,
+                             double eps) const override;
+
+    // checks if lim s->0 lim t->0  p + t(v1 + s v2) in solid
+    INSOLID_TYPE VecInSolid2 (const Point<3> & p,
+                              const Vec<3> & v1,
+                              const Vec<3> & v2,
+                              double eps) const override;
+
+  
+    int GetNSurfaces() const override;
+    Surface & GetSurface (int i = 0) override;
+    const Surface & GetSurface (int i = 0) const override;
+
+
+    void Reduce (const BoxSphere<3> & box) override;
+    void UnReduce () override;
+  };
+
+}
+
+#endif //_EXTRUSION_HPP

netgen/include/csg/gencyl.hpp

@@ -0,0 +1,70 @@
+#ifndef FILE_GENCYL
+#define FILE_GENCYL
+
+/**************************************************************************/
+/* File:   gencyl.hh                                                      */
+/* Author: Joachim Schoeberl                                              */
+/* Date:   14. Oct. 96                                                    */
+/**************************************************************************/
+
+namespace netgen
+{
+
+
+  /*
+  
+  Generalized Cylinder
+  
+  */
+
+
+  ///
+  class GeneralizedCylinder : public Surface
+  {
+    ///
+    ExplicitCurve2d & crosssection;
+    ///
+    Point<3> planep;
+    ///
+    Vec<3> planee1, planee2, planee3;
+  
+    ///  Vec<3> ex, ey, ez;
+    Vec2d e2x, e2y;
+    ///
+    Point<3> cp;
+  
+  public:
+    ///
+    GeneralizedCylinder (ExplicitCurve2d & acrosssection,
+			 Point<3> ap, Vec<3> ae1, Vec<3> ae2);
+  
+    ///
+    virtual void Project (Point<3> & p) const;
+  
+    ///
+    virtual int BoxInSolid (const BoxSphere<3> & box) const;
+    /// 0 .. no, 1 .. yes, 2 .. maybe
+  
+    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;
+    ///
+    virtual double MaxCurvatureLoc (const Point<3> & c, double rad) const;
+    ///
+    virtual Point<3> GetSurfacePoint () const;
+    ///
+    virtual void Print (ostream & str) const;
+  
+    ///
+    virtual void Reduce (const BoxSphere<3> & box);
+    ///
+    virtual void UnReduce ();
+  };  
+
+}
+
+#endif

netgen/include/csg/geoml.hpp

@@ -0,0 +1,16 @@
+#ifndef FILE_GEOML
+#define FILE_GEOML
+
+/* *************************************************************************/
+/* File:   geoml.hh                                                        */
+/* Author: Joachim Schoeberl                                               */
+/* Date:   21. Jun. 98                                                     */
+/* *************************************************************************/
+
+#include <geom/geom.hh>
+
+#include <geom/solid.hh>
+#include <geom/algprim.hh>
+#include <geom/adtree.hh>
+#include <geom/csgeom.hh>
+#endif

netgen/include/csg/identify.hpp

@@ -0,0 +1,213 @@
+
+#ifndef FILE_IDENTIFY
+#define FILE_IDENTIFY
+
+/**************************************************************************/
+/* File:   identify.hh                                                    */
+/* Author: Joachim Schoeberl                                              */
+/* Date:   1. Aug. 99                                                    */
+/**************************************************************************/
+
+
+namespace netgen
+{
+
+  /**
+     Identify surfaces for periodic b.c. or
+     thin domains
+  */
+
+
+  class SpecialPoint;
+  class Identification
+  {
+  protected:
+    const CSGeometry & geom;
+    // identified faces, index sorted
+    INDEX_2_HASHTABLE<int> identfaces;
+    int nr;
+
+  public:
+    DLL_HEADER Identification (int anr, const CSGeometry & ageom);
+    DLL_HEADER virtual ~Identification ();
+    DLL_HEADER virtual void Print (ostream & ost) const = 0;
+    DLL_HEADER virtual void GetData (ostream & ost) const = 0;
+
+    /// obsolete
+    //  virtual void IdentifySpecialPoints (NgArray<class SpecialPoint> & points);
+
+    /// can identify both special points (fixed direction)
+    /// (identified points, same tangent)
+    virtual int Identifiable (const SpecialPoint & sp1, const SpecialPoint & sp2,
+			      const TABLE<int> & specpoint2solid,			  
+			      const TABLE<int> & specpoint2surface) const;
+    ///
+    virtual int Identifiable (const Point<3> & p1, const Point<3> & sp2) const;
+    /// is it possible to identify sp1 with some other ?
+    virtual int IdentifiableCandidate (const SpecialPoint & sp1) const;
+  
+    /// are points (if connected) by a short edge (direction anyhow) ?
+    virtual int ShortEdge (const SpecialPoint & sp1, const SpecialPoint & sp2) const;
+
+    /// add entries in mesh identification tables
+    virtual void IdentifyPoints (class Mesh & mesh);
+
+    /// add entries to identified faces (based on segment infos)
+    virtual void IdentifyFaces (class Mesh & mesh);
+
+    /// get point on other surface, add entry in mesh identifications
+    virtual PointIndex GetIdentifiedPoint (class Mesh & mesh, PointIndex pi1);
+
+    /// copy surfaces, or fill rectangles
+    virtual void BuildSurfaceElements (NgArray<class Segment> & segs,
+				       class Mesh & mesh,
+				       const Surface * surf);
+
+    /// insert volume elements in thin layers
+    virtual void BuildVolumeElements (NgArray<class Element2d> & surfels,
+				      class Mesh & mesh);
+
+    /// get list of identified faces
+    virtual void GetIdentifiedFaces (NgArray<INDEX_2> & idfaces) const;
+
+    friend ostream & operator<< (ostream & ost, Identification & ident);
+  };
+
+
+  class PeriodicIdentification : public Identification
+  {
+    const Surface * s1;
+    const Surface * s2;
+    Transformation<3> trafo; // from s1 to s2
+    Transformation<3> inv_trafo; // from s2 to s1
+  public:
+    PeriodicIdentification (int anr,
+			    const CSGeometry & ageom,
+			    const Surface * as1,
+			    const Surface * as2,
+                            Transformation<3> atrafo = Vec<3>(0,0,0));
+    virtual ~PeriodicIdentification () override;
+    virtual void Print (ostream & ost) const override;
+    virtual void GetData (ostream & ost) const override;
+
+
+    //  virtual void IdentifySpecialPoints (NgArray<class SpecialPoint> & points);
+    virtual int Identifiable (const SpecialPoint & sp1, const SpecialPoint & sp2,
+			      const TABLE<int> & specpoint2solid,
+			      const TABLE<int> & specpoint2surface) const override;
+
+    virtual int Identifiable (const Point<3> & p1, const Point<3> & sp2) const override;
+    virtual PointIndex GetIdentifiedPoint (class Mesh & mesh, PointIndex pi1) override;
+    virtual void IdentifyPoints (class Mesh & mesh) override;
+    virtual void IdentifyFaces (class Mesh & mesh) override;
+    virtual void BuildSurfaceElements (NgArray<class Segment> & segs,
+				       class Mesh & mesh,
+				       const Surface * surf) override;
+  };
+
+
+  ///
+  class TopLevelObject;
+  class CloseSurfaceIdentification : public Identification
+  {
+    const Surface * s1;
+    const Surface * s2;
+    const TopLevelObject * domain;
+    ///
+    int dom_nr;
+    /// number of refinement levels (in Z-refinement)
+    int ref_levels;
+    /// number of refinement levels for layer next to s1 (in Z-refinement)
+    int ref_levels_s1;
+    /// number of refinement levels for layer next to s2 (in Z-refinement)
+    int ref_levels_s2;
+    ///
+    double eps_n;
+    Array<double> slices;
+    /// used only for domain-local identification:
+    NgArray<int> domain_surfaces;
+    ///
+    bool dom_surf_valid;
+
+    ///
+    Vec<3> direction;
+    ///
+    bool usedirection;
+  public:
+    CloseSurfaceIdentification (int anr, 
+				const CSGeometry & ageom,
+				const Surface * as1,
+				const Surface * as2,
+				const TopLevelObject * adomain,
+				const Flags & flags);
+    virtual ~CloseSurfaceIdentification ();
+
+    virtual void Print (ostream & ost) const;
+    virtual void GetData (ostream & ost) const;
+
+
+    //  virtual void IdentifySpecialPoints (NgArray<class SpecialPoint> & points);
+    virtual int Identifiable (const SpecialPoint & sp1, const SpecialPoint & sp2,
+			      const TABLE<int> & specpoint2solid,
+			      const TABLE<int> & specpoint2surface) const;
+    virtual int Identifiable (const Point<3> & p1, const Point<3> & sp2) const;
+    virtual int IdentifiableCandidate (const SpecialPoint & sp1) const;
+    virtual int ShortEdge (const SpecialPoint & sp1, const SpecialPoint & sp2) const;
+    virtual PointIndex GetIdentifiedPoint (class Mesh & mesh, PointIndex pi1);
+    const Array<double> & GetSlices () const { return slices; }
+    virtual void IdentifyPoints (class Mesh & mesh);
+    virtual void IdentifyFaces (class Mesh & mesh);
+    virtual void BuildSurfaceElements (NgArray<class Segment> & segs,
+				       class Mesh & mesh,
+				       const Surface * surf);
+    void BuildSurfaceElements2 (NgArray<class Segment> & segs,
+				class Mesh & mesh,
+				const Surface * surf);
+
+    virtual void BuildVolumeElements (NgArray<class Element2d> & surfels,
+				      class Mesh & mesh);
+
+    int RefLevels () const { return ref_levels; }
+    int RefLevels1 () const { return ref_levels_s1; }
+    int RefLevels2 () const { return ref_levels_s2; }
+
+    bool IsSkewIdentification(void) const {return usedirection;}
+    const Vec<3> & GetDirection(void) const {return direction;}
+
+    const Surface & GetSurface1(void) const
+    { return *s1;}
+    const Surface & GetSurface2(void) const
+    { return *s2;}
+  };
+
+
+  class CloseEdgesIdentification : public Identification
+  {
+    const Surface * facet;
+    const Surface * s1;
+    const Surface * s2;
+  public:
+    CloseEdgesIdentification (int anr,
+			      const CSGeometry & ageom,
+			      const Surface * afacet,
+			      const Surface * as1,
+			      const Surface * as2);
+    virtual ~CloseEdgesIdentification ();
+    virtual void Print (ostream & ost) const;
+    virtual void GetData (ostream & ost) const;
+
+    //  virtual void IdentifySpecialPoints (NgArray<class SpecialPoint> & points);
+    virtual int Identifiable (const SpecialPoint & sp1, const SpecialPoint & sp2,
+			      const TABLE<int> & specpoint2solid,
+			      const TABLE<int> & specpoint2surface) const;
+
+
+    virtual void IdentifyPoints (class Mesh & mesh);
+    virtual void BuildSurfaceElements (NgArray<class Segment> & segs,
+				       class Mesh & mesh,
+				       const Surface * surf);
+  };
+
+}
+
+#endif

netgen/include/csg/manifold.hpp

@@ -0,0 +1,29 @@
+#ifndef FILE_MANIFOLD
+#define FILE_MANIFOLD
+
+/**************************************************************************/
+/* File:   manifold.hh                                                    */
+/* Author: Joachim Schoeberl                                              */
+/* Date:   7. Aug. 96                                                     */
+/**************************************************************************/
+
+namespace netgen
+{
+
+
+  /**
+     Basis class for manifolds in 2d and 3d
+  */
+  class Manifold
+  {
+  public:
+    ///
+    Manifold ();
+    ///
+    virtual ~Manifold ();
+  };
+
+}
+
+
+#endif

netgen/include/csg/meshsurf.hpp

@@ -0,0 +1,46 @@
+#ifndef FILE_MESHSURF
+#define FILE_MESHSURF
+
+namespace netgen
+{
+
+  ///
+  class Meshing2Surfaces : public Meshing2
+  {
+    ///
+    const Surface & surface;
+    
+    /// should be movec to base ... 
+    const MeshingParameters & mparam;
+  public:
+    ///
+    //  Meshing2Surfaces (const Surface & asurf);
+    ///
+    Meshing2Surfaces (const CSGeometry& geo,
+                      const Surface & asurf,
+                      const MeshingParameters & mp,
+		      const Box<3> & aboundingbox);
+
+  protected:
+    ///
+    void DefineTransformation(const Point<3> & p1,
+                              const Point<3> & p2,
+                              const PointGeomInfo * geominfo1,
+                              const PointGeomInfo * geominfo2) override;
+    ///
+    void TransformToPlain(const Point<3> & locpoint, 
+                          const MultiPointGeomInfo & geominfo,
+                          Point<2> & plainpoint, 
+                          double h, int & zone) override;
+    ///
+    int TransformFromPlain(const Point<2>& plainpoint,
+                           Point<3>& locpoint, 
+                           PointGeomInfo & gi,
+                           double h) override;
+    ///
+    double CalcLocalH(const Point<3> & p, double gh) const override;
+  };
+}
+
+#endif
+