ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-12-26 13:50:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

fix splinesurface when points are not given exactly

Browse Source
This commit is contained in:
Christopher Lackner 2017-08-02 14:45:03 +02:00
parent 0fe8f9bfcf
commit e910eb5613
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
libsrc/csg/splinesurface.cpp
View File

@ -5,7 +5,8 @@ namespace netgen
{ {
void SplineSurface :: AppendPoint(const Point<3> & p, const double reffac, const bool hpref) void SplineSurface :: AppendPoint(const Point<3> & p, const double reffac, const bool hpref)
{ {
auto pp = Point<3>(p); auto pp = p;
Project(pp);
geompoints.Append(GeomPoint<3>(pp,reffac)); geompoints.Append(GeomPoint<3>(pp,reffac));
geompoints.Last().hpref = hpref; geompoints.Last().hpref = hpref;
} }
@ -20,13 +21,13 @@ void SplineSurface :: AppendPoint(const Point<3> & p, const double reffac, const
string SplineSurface :: GetBCNameOf (Point<3> p1, Point<3> p2) const string SplineSurface :: GetBCNameOf (Point<3> p1, Point<3> p2) const
{ {
double eps = 1e-5;
for(int i=0; i<splines.Size(); i++) for(int i=0; i<splines.Size(); i++)
{ {
auto pp1 = Point<3>(splines[i]->GetPoint(0)); auto pp1 = Point<3>(splines[i]->GetPoint(0));
Project(pp1); Project(pp1);
auto pp2 = Point<3>(splines[i]->GetPoint(1)); auto pp2 = Point<3>(splines[i]->GetPoint(1));
Project(pp2); Project(pp2);
double eps = (p1-p2).Length() * 1e-4;
if (((pp1-p1).Length()<eps && (pp2-p2).Length() < eps) || ((pp1-p2).Length() < eps && (pp2-p1).Length() < eps)) if (((pp1-p1).Length()<eps && (pp2-p2).Length() < eps) || ((pp1-p2).Length() < eps && (pp2-p1).Length() < eps))
{ {
return bcnames[i]; return bcnames[i];