ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-11-15 18:38:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

Boundarylayers - handle opposing normal vectors correctly

Browse Source
This commit is contained in:
Matthias Hochsteger 2023-10-19 11:49:22 +02:00
parent ec0f848030
commit 7a0d7594c8
1 changed files with 12 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
libsrc/meshing/boundarylayer.cpp
View File

@ -755,7 +755,7 @@ namespace netgen
// combine normal vectors for each face to keep uniform distances // combine normal vectors for each face to keep uniform distances
auto & np = growthvectors[pi]; auto & np = growthvectors[pi];
ArrayMem<Vec<3>, 3> ns; ArrayMem<Vec<3>, 5> ns;
for (auto &[facei, n] : normals) { for (auto &[facei, n] : normals) {
ns.Append(n); ns.Append(n);
} }
@ -784,7 +784,7 @@ namespace netgen
{ {
for (auto j : Range(i + 1, ns.Size())) for (auto j : Range(i + 1, ns.Size()))
{ {
double ip = ns[i] * ns[j]; double ip = fabs(ns[i] * ns[j]);
if(ip > val) if(ip > val)
{ {
val = ip; val = ip;
@ -795,8 +795,17 @@ namespace netgen
} }
removed.Append(ns[maxpos1]); removed.Append(ns[maxpos1]);
removed.Append(ns[maxpos2]); removed.Append(ns[maxpos2]);
ns[maxpos1] = 0.5 * (ns[maxpos1] + ns[maxpos2]); const auto dot = ns[maxpos1] * ns[maxpos2];
auto average = 0.5*(ns[maxpos1] + ns[maxpos2]);
average.Normalize();
ns.DeleteElement(maxpos2); ns.DeleteElement(maxpos2);
// if vectors have nearly opposite directions -> remove both
// otherwise, replace both vectors with their average
if(dot < -1+1e-3)
ns.DeleteElement(maxpos1);
else
ns[maxpos1] = average;
} }
if(ns.Size() == 0) if(ns.Size() == 0)