fix GetSurfaceElementOfPoint

only call GetElementOfPoint, if there are volume elements (otherwise we get an endless loop of mutual function calls here...)
This commit is contained in:
mhochsteger@cerbsim.com 2021-09-23 09:14:58 +02:00
parent 559a9f2beb
commit 00855acfec

View File

@ -5757,7 +5757,10 @@ namespace netgen
else else
{ {
double vlam[3]; double vlam[3];
int velement = GetElementOfPoint(p,vlam,NULL,build_searchtree,allowindex); int velement = 0;
if(GetNE())
GetElementOfPoint(p,vlam,NULL,build_searchtree,allowindex);
//(*testout) << "p " << p << endl; //(*testout) << "p " << p << endl;
//(*testout) << "velement " << velement << endl; //(*testout) << "velement " << velement << endl;