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

Also smooth boundary layers at corners if adjacent surface elements have similar normal vectors

Browse Source
This commit is contained in:
Matthias Hochsteger 2025-02-06 18:50:04 +01:00
parent a8a75614c0
commit 7aff94046f
1 changed files with 21 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
libsrc/meshing/boundarylayer_interpolate.cpp
View File

@ -155,8 +155,27 @@ void BoundaryLayerTool ::InterpolateGrowthVectors()
{ {
for (auto* p_seg : edgenr2seg[edgenr]) for (auto* p_seg : edgenr2seg[edgenr])
for (auto pi : p_seg->PNums()) for (auto pi : p_seg->PNums())
if (pi < first_new_pi && point_types[pi] == EDGEPOINT) {
if (pi >= first_new_pi)
continue;
if (point_types[pi] == EDGEPOINT)
point_types[pi] = SURFACEPOINT; point_types[pi] = SURFACEPOINT;
else if (point_types[pi] == FIXEDPOINT)
{
// Check at edge corners if all adjacent surface elements have roughly the same normal.
// If so, also treat this point as surface point for growth vector interpolation
Vec<3> n = 0.0;
for (auto si : p2sel[pi])
n += getNormal(mesh[si]);
n.Normalize();
bool is_corner = false;
for (auto si : p2sel[pi])
if (getNormal(mesh[si]) * n < 0.9)
is_corner = true;
if (!is_corner)
point_types[pi] = SURFACEPOINT;
}
}
continue; continue;
} }
} }