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

some thread-safe STLGeometry projection functions

Browse Source
This commit is contained in:
Matthias Hochsteger 2019-10-18 13:40:51 +02:00
parent 17dfd45609
commit 8659e5a3ee
1 changed files with 19 additions and 21 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

36
libsrc/stlgeom/stlgeom.cpp
View File

@ -124,14 +124,18 @@ bool STLGeometry :: CalcPointGeomInfo(int /*surfind*/, PointGeomInfo& gi, const
bool STLGeometry :: ProjectPointGI (int surfind, Point<3> & p, PointGeomInfo & gi) const
{
SelectChartOfTriangle(gi.trignum);
gi.trignum = Project (p);
if (!gi.trignum)
static std::mutex mutex_project_whole_surface;
int meshchart = GetChartNr(gi.trignum);
const STLChart& chart = GetChart(meshchart);
int trignum = chart.ProjectNormal(p);
if(trignum==0)
{
// non-thread-safe implementation
std::lock_guard<std::mutex> guard(mutex_project_whole_surface);
PrintMessage(7,"project failed");
gi.trignum = ProjectOnWholeSurface(p);
if (!gi.trignum)
SelectChartOfTriangle (gi.trignum); // needed because ProjectOnWholeSurface uses meshchartnv (the normal vector of selected chart)
trignum = ProjectOnWholeSurface(p);
if(trignum==0)
{
PrintMessage(7, "project on whole surface failed");
return false;
@ -142,15 +146,7 @@ bool STLGeometry :: ProjectPointGI (int surfind, Point<3> & p, PointGeomInfo & g
void STLGeometry :: ProjectPoint (INDEX surfind, Point<3> & p) const
{
if (!Project (p))
{
PrintMessage(7,"project failed");
if (!ProjectOnWholeSurface(p))
{
PrintMessage(7, "project on whole surface failed");
}
}
throw Exception("ProjectPoint without PointGeomInfo not implemented");
}
void STLGeometry ::
@ -173,11 +169,13 @@ PointBetween (const Point<3> & p1, const Point<3> & p2, double secpoint,
Point<3> np1 = newp;
Point<3> np2 = newp;
SelectChartOfTriangle (gi1.trignum);
int tn1 = Project (np1);
auto ngi1 = gi1;
auto ngi2 = gi2;
// SelectChartOfTriangle (gi1.trignum);
int tn1 = ProjectPointGI (surfi, np1, ngi1);
SelectChartOfTriangle (gi2.trignum);
int tn2 = Project (np2);
// SelectChartOfTriangle (gi2.trignum);
int tn2 = ProjectPointGI (surfi, np2, ngi2);
newgi.trignum = tn1; //urspruengliche version
newp = np1; //urspruengliche version