fix tolerance in specpoints check (compare Dist2() with sqr(1e-8*geomsize) )

This commit is contained in:
Matthias Hochsteger 2021-06-08 16:08:52 +02:00
parent d6ca80d50e
commit 0b90d24d81

View File

@ -2090,7 +2090,7 @@ namespace netgen
for (int m = 0; m < locsearch.Size(); m++) for (int m = 0; m < locsearch.Size(); m++)
{ {
if (Dist2 (specpoints[locsearch[m]].p, apoints[i]) < 1e-10*geomsize if (Dist2 (specpoints[locsearch[m]].p, apoints[i]) < sqr(1e-8*geomsize)
&& Abs2(specpoints[locsearch[m]].v - t) < 1e-8) && Abs2(specpoints[locsearch[m]].v - t) < 1e-8)
{ {
spi = locsearch[m]; spi = locsearch[m];