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

Fix abort criteria in boundarylayer thickness limitation

Browse Source
This commit is contained in:
Matthias Hochsteger 2024-10-24 17:10:10 +02:00
parent 7f294a582c
commit 4f927f5469
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
libsrc/meshing/boundarylayer_limiter.hpp
View File

@ -534,16 +534,16 @@ struct GrowthVectorLimiter
ScaleLimit(pi_max_limit, 0.9); ScaleLimit(pi_max_limit, 0.9);
set_points(); set_points();
counter++; counter++;
if (GetLimit(pi_max_limit) < 1e-10)
{
WriteErrorMesh("error_blayer_self_intersection_pi" + ToString(pi_max_limit) + ".vol.gz");
throw NgException("Stop meshing in boundary layer thickness limitation: overlapping regions detected at elements " + ToString(tri) + " and " + ToString(tri2));
}
if (debugparam.debugoutput && counter > 20) if (debugparam.debugoutput && counter > 20)
{ {
cerr << "Limit intersecting surface elements: too many " cerr << "Limit intersecting surface elements: too many "
"limitation steps, sels: " "limitation steps, sels: "
<< Get(sei) << '\t' << Get(sej) << endl; << Get(sei) << '\t' << Get(sej) << endl;
if (GetLimit(pi_max_limit) < 1e-10)
{
WriteErrorMesh("error_blayer_self_intersection_pi" + ToString(pi_max_limit) + ".vol.gz");
throw NgException("Stop meshing in boundary layer thickness limitation: overlapping regions detected at elements " + ToString(tri) + " and " + ToString(tri2));
}
for (auto si : {sei, sej}) for (auto si : {sei, sej})
{ {
auto sel = Get(si); auto sel = Get(si);