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

fixes for surface element point search

Browse Source
This commit is contained in:
Joachim Schöberl 2016-05-09 05:45:09 +02:00
parent f07a7a95d0
commit 1337786e73
1 changed files with 12 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
libsrc/interface/nginterface_v2.cpp
View File

@ -739,9 +739,18 @@ namespace netgen
for (int i = 0; i < numind; i++) dummy[i] = indices[i]+1; for (int i = 0; i < numind; i++) dummy[i] = indices[i]+1;
double lam3[3]; double lam3[3];
int ind;
if (mesh->GetDimension() == 2)
{
Point<3> p2d(p[0], p[1], 0); Point<3> p2d(p[0], p[1], 0);
int ind = ind = mesh->GetElementOfPoint(p2d, lam3, &dummy, build_searchtree);
mesh->GetElementOfPoint(p2d, lam3, &dummy, build_searchtree); }
else
{
Point3d p3d(p[0], p[1], p[2]);
ind = mesh->GetSurfaceElementOfPoint(p3d, lam3, &dummy, build_searchtree);
}
if (ind > 0) if (ind > 0)
{ {