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

Fix CheckLimits

Browse Source
This commit is contained in:
Matthias Hochsteger 2024-10-30 21:53:57 +01:00
parent e7e15069f3
commit 9f3ccbb9d0
1 changed files with 17 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
libsrc/meshing/boundarylayer_limiter.hpp
View File

@ -671,14 +671,21 @@ struct GrowthVectorLimiter
} }
} }
void CheckLimits () void CheckLimits (int line)
{ {
for (auto [pi, data] : tool.growth_vector_map) auto check_point = [&] (PointIndex pi) {
if (limits[pi] < 1e-10) if (limits[pi] < 1e-8)
{ {
WriteErrorMesh("error_blayer_intersection_pi" + ToString(pi) + ".vol.gz"); WriteErrorMesh("error_blayer_intersection_pi" + ToString(pi) + ".vol.gz");
throw NgException("Stop meshing in boundary layer thickness limitation: overlapping regions detected at point " + ToString(pi)); throw NgException(__FILE__ + ToString(line) + ": Stop meshing in boundary layer thickness limitation: overlapping regions detected at point " + ToString(pi));
} }
};
for (auto pi : Range(growthvectors))
check_point(pi);
for (auto& [special_pi, special_point] : tool.special_boundary_points)
check_point(special_pi);
} }
void Perform () void Perform ()
@ -695,24 +702,24 @@ struct GrowthVectorLimiter
{ {
PrintMessage(0, "GrowthVectorLimiter pass ", i_pass); PrintMessage(0, "GrowthVectorLimiter pass ", i_pass);
double safety = safeties[i_pass]; double safety = safeties[i_pass];
CheckLimits(); CheckLimits(__LINE__);
// intersect segment with original surface elements // intersect segment with original surface elements
LimitOriginalSurface(2.1); LimitOriginalSurface(2.1);
CheckLimits(); CheckLimits(__LINE__);
// intersect prisms with themself // intersect prisms with themself
LimitSelfIntersection(1.3 * safety); LimitSelfIntersection(1.3 * safety);
CheckLimits(); CheckLimits(__LINE__);
// intesect segment with prism // intesect segment with prism
LimitBoundaryLayer(safety); LimitBoundaryLayer(safety);
CheckLimits(); CheckLimits(__LINE__);
for (auto i : Range(10)) for (auto i : Range(10))
EqualizeLimits(smoothing_factors[i_pass]); EqualizeLimits(smoothing_factors[i_pass]);
CheckLimits(); CheckLimits(__LINE__);
if (i_pass == safeties.size() - 1) if (i_pass == safeties.size() - 1)
FixIntersectingSurfaceTrigs(); FixIntersectingSurfaceTrigs();
CheckLimits(); CheckLimits(__LINE__);
} }
for (auto i : Range(growthvectors)) for (auto i : Range(growthvectors))