ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-01-27 13:20:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

use convex hulls

Browse Source
This commit is contained in:
Joachim Schöberl 2016-09-07 09:03:29 +02:00
parent 0eaedadf10
commit 103f0e6280
1 changed files with 27 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
libsrc/gprim/spline.hpp
View File

@ -95,6 +95,9 @@ namespace netgen
Array < Point<D> > & points, const double eps) const Array < Point<D> > & points, const double eps) const
{points.SetSize(0);} {points.SetSize(0);}
// is the point in the convex hull (increased by eps) of the spline ?
virtual bool InConvexHull (Point<D> p, double eps) const = 0;
virtual double MaxCurvature(void) const = 0; virtual double MaxCurvature(void) const = 0;
virtual string GetType(void) const {return "splinebase";} virtual string GetType(void) const {return "splinebase";}
@ -140,7 +143,12 @@ namespace netgen
virtual void LineIntersections (const double a, const double b, const double c, virtual void LineIntersections (const double a, const double b, const double c,
Array < Point<D> > & points, const double eps) const; Array < Point<D> > & points, const double eps) const;
virtual bool InConvexHull (Point<D> p, double eps) const
{
return MinDistLP2 (p1, p2, p) < sqr(eps);
}
virtual double MaxCurvature(void) const {return 0;} virtual double MaxCurvature(void) const {return 0;}
virtual void Project (const Point<D> point, Point<D> & point_on_curve, double & t) const; virtual void Project (const Point<D> point, Point<D> & point_on_curve, double & t) const;
@ -186,6 +194,11 @@ namespace netgen
DLL_HEADER virtual void LineIntersections (const double a, const double b, const double c, DLL_HEADER virtual void LineIntersections (const double a, const double b, const double c,
Array < Point<D> > & points, const double eps) const; Array < Point<D> > & points, const double eps) const;
virtual bool InConvexHull (Point<D> p, double eps) const
{
return MinDistTP2 (p1, p2, p3, p) < sqr(eps);
}
DLL_HEADER virtual double MaxCurvature(void) const; DLL_HEADER virtual double MaxCurvature(void) const;
DLL_HEADER virtual void Project (const Point<D> point, Point<D> & point_on_curve, double & t) const; DLL_HEADER virtual void Project (const Point<D> point, Point<D> & point_on_curve, double & t) const;
@ -232,6 +245,11 @@ namespace netgen
virtual void LineIntersections (const double a, const double b, const double c, virtual void LineIntersections (const double a, const double b, const double c,
Array < Point<D> > & points, const double eps) const; Array < Point<D> > & points, const double eps) const;
virtual bool InConvexHull (Point<D> p, double eps) const
{
return (Dist2 (p, pm) < sqr(radius+eps));
}
virtual double MaxCurvature(void) const {return 1./radius;} virtual double MaxCurvature(void) const {return 1./radius;}
}; };
@ -261,6 +279,10 @@ namespace netgen
virtual void GetCoeff (Vector & coeffs) const {;} virtual void GetCoeff (Vector & coeffs) const {;}
virtual double MaxCurvature(void) const {return 1;} virtual double MaxCurvature(void) const {return 1;}
// needs implementation ...
virtual bool InConvexHull (Point<D> p, double eps) const
{ return true; }
}; };
@ -597,6 +619,10 @@ namespace netgen
virtual void GetCoeff (Vector & coeffs) const {;} virtual void GetCoeff (Vector & coeffs) const {;}
virtual double MaxCurvature(void) const {return 1;} virtual double MaxCurvature(void) const {return 1;}
// needs implementation ...
virtual bool InConvexHull (Point<D> p, double eps) const
{ return true; }
}; };
// Constructor // Constructor