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

fix for close-surface identification (fixes edge xxx-yyy multiple times in surface mesh

Browse Source
This commit is contained in:
Joachim Schöberl 2016-04-05 09:51:12 +02:00
parent 10d0f42ae6
commit d1ebdf35d8
1 changed files with 9 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
libsrc/csg/identify.cpp
View File

@ -1245,13 +1245,20 @@ BuildSurfaceElements (Array<Segment> & segs,
(mesh.GetIdentifications().Get (s2[0], s1[1], nr) && (mesh.GetIdentifications().Get (s2[0], s1[1], nr) &&
mesh.GetIdentifications().Get (s2[1], s1[0], nr))) mesh.GetIdentifications().Get (s2[1], s1[0], nr)))
{ {
Vec<3> ns = surf->GetNormalVector (mesh[s1[0]]);
Vec<3> t1 = mesh[s1[1]] - mesh[s1[0]];
Vec<3> t2 = mesh[s2[1]] - mesh[s2[0]];
Vec<3> nst1 = Cross(t1, ns);
Vec<3> nst2 = Cross(t2, ns);
Vec<3> dvec = Center(mesh[s1[0]], mesh[s1[1]])-Center(mesh[s2[0]], mesh[s2[1]]);
if (nst1 * dvec < 0) continue;
Element2d el(s1[0], s1[1], s2[0], s2[1]); Element2d el(s1[0], s1[1], s2[0], s2[1]);
Vec<3> n = Cross (mesh[el[1]] - mesh[el[0]], Vec<3> n = Cross (mesh[el[1]] - mesh[el[0]],
mesh[el[3]] - mesh[el[0]]); mesh[el[3]] - mesh[el[0]]);
Vec<3> ns = surf->GetNormalVector (mesh[el[0]]);
if (n * ns < 0) if (n * ns < 0)
{ {
Swap (el.PNum(1), el.PNum(2)); Swap (el.PNum(1), el.PNum(2));