ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-02-04 17:10:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

fix for getsurfpoint for 3d surfacegeometries if no volume elements exist

Browse Source
This commit is contained in:
Michael Neunteufel 2018-08-07 11:46:00 +02:00
parent 0251b4ad0c
commit 73292485eb
1 changed files with 8 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
libsrc/meshing/meshclass.cpp
View File

@ -4327,10 +4327,12 @@ namespace netgen
elementsearchtree = NULL; elementsearchtree = NULL;
int ne = (dimension == 2) ? GetNSE() : GetNE(); int ne = (dimension == 2) ? GetNSE() : GetNE();
if (dimension == 3 && !GetNE() && GetNSE())
ne = GetNSE();
if (ne) if (ne)
{ {
if (dimension == 2) if (dimension == 2 || (dimension == 3 && !GetNE()) )
{ {
Box<3> box (Box<3>::EMPTY_BOX); Box<3> box (Box<3>::EMPTY_BOX);
for (SurfaceElementIndex sei = 0; sei < ne; sei++) for (SurfaceElementIndex sei = 0; sei < ne; sei++)
@ -4834,7 +4836,7 @@ namespace netgen
const double pointtol = 1e-12; const double pointtol = 1e-12;
netgen::Point<3> pmin = p - Vec<3> (pointtol, pointtol, pointtol); netgen::Point<3> pmin = p - Vec<3> (pointtol, pointtol, pointtol);
netgen::Point<3> pmax = p + Vec<3> (pointtol, pointtol, pointtol); netgen::Point<3> pmax = p + Vec<3> (pointtol, pointtol, pointtol);
if (dimension == 2) if (dimension == 2 || (dimension==3 && !GetNE() && GetNSE()))
{ {
int ne; int ne;
int ps_startelement = 0; // disable global buffering int ps_startelement = 0; // disable global buffering
@ -4988,6 +4990,9 @@ namespace netgen
//(*testout) << "p " << p << endl; //(*testout) << "p " << p << endl;
//(*testout) << "velement " << velement << endl; //(*testout) << "velement " << velement << endl;
if (!GetNE() && GetNSE() )
return velement;
Array<int> faces; Array<int> faces;
topology.GetElementFaces(velement,faces); topology.GetElementFaces(velement,faces);