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

boundarylayers - don't project on not-in-surface-direction

Browse Source
This commit is contained in:
Christopher Lackner 2022-02-25 12:08:24 +01:00
parent bf02481299
commit 2d70263cda
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
libsrc/meshing/boundarylayer.cpp
View File

@ -360,8 +360,11 @@ namespace netgen
auto v2 = (mesh[prev] - mesh[pi]).Normalize(); auto v2 = (mesh[prev] - mesh[pi]).Normalize();
auto v3 = growthvectors[pi]; auto v3 = growthvectors[pi];
v3.Normalize(); v3.Normalize();
if((v1 * v3 > 1e-12) || (v2 * v3 > 1e-12)) auto tol = v1.Length() * 1e-12;
if((v1 * v3 > -tol) && (v2 * v3 > -tol))
in_surface_direction.SetBit(sel.GetIndex()); in_surface_direction.SetBit(sel.GetIndex());
else
continue;
auto& g = growthvectors[pi]; auto& g = growthvectors[pi];
auto ng = n * g; auto ng = n * g;