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

use relative tolerance in identifypoints

Browse Source
This commit is contained in:
Christopher Lackner 2020-12-18 15:58:51 +01:00
parent 1e8715dc34
commit 7bd454e385
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
libsrc/csg/identify.cpp
View File

@ -318,6 +318,10 @@ GetIdentifiedPoint (class Mesh & mesh, int pi)
void PeriodicIdentification :: IdentifyPoints (class Mesh & mesh) void PeriodicIdentification :: IdentifyPoints (class Mesh & mesh)
{ {
Point3d p1, p2;
mesh.GetBox(p1, p2);
auto eps = 1e-6 * (p2-p1).Length();
for (int i = 1; i <= mesh.GetNP(); i++) for (int i = 1; i <= mesh.GetNP(); i++)
{ {
Point<3> p = mesh.Point(i); Point<3> p = mesh.Point(i);
@ -327,7 +331,7 @@ void PeriodicIdentification :: IdentifyPoints (class Mesh & mesh)
pp = trafo(pp); pp = trafo(pp);
s2->Project (pp); s2->Project (pp);
for (int j = 1; j <= mesh.GetNP(); j++) for (int j = 1; j <= mesh.GetNP(); j++)
if (Dist2(mesh.Point(j), pp) < 1e-6) if (Dist2(mesh.Point(j), pp) < eps)
{ {
mesh.GetIdentifications().Add (i, j, nr); mesh.GetIdentifications().Add (i, j, nr);
/* /*