ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-12-28 06:40:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

Cleanup

Browse Source
This commit is contained in:
Matthias Hochsteger 2024-10-07 15:45:19 +02:00
parent 70eb3cc1f1
commit 6160fd4d26
2 changed files with 2 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
libsrc/meshing/improve3.cpp
View File

@ -1382,7 +1382,6 @@ void MeshOptimize3d :: SwapImprove (const NgBitArray * working_elements)
break; break;
auto [pi0, pi1] = edges[i]; auto [pi0, pi1] = edges[i];
try {
double d_badness = SwapImproveEdge (working_elements, elementsonnode, faces, pi0, pi1, true); double d_badness = SwapImproveEdge (working_elements, elementsonnode, faces, pi0, pi1, true);
if(d_badness<0.0) if(d_badness<0.0)
{ {
@ -1390,13 +1389,6 @@ void MeshOptimize3d :: SwapImprove (const NgBitArray * working_elements)
candidate_edges[index] = make_tuple(d_badness, i); candidate_edges[index] = make_tuple(d_badness, i);
} }
} }
catch(const ngcore::Exception &e) {
cerr << "Error in SwapImproveEdge " << pi0 << '-' << pi1 << endl;
if(debugparam.write_mesh_on_error)
mesh.Save("error_swapimproveedge_" +ToString(pi0) + "_" + ToString(pi1) +".vol");
throw;
}
}
}, TasksPerThread (4)); }, TasksPerThread (4));
auto edges_with_improvement = candidate_edges.Part(0, improvement_counter.load()); auto edges_with_improvement = candidate_edges.Part(0, improvement_counter.load());

4
libsrc/meshing/meshfunc.cpp
View File

@ -499,7 +499,6 @@ namespace netgen
PrintMessage (3, "Call remove problem"); PrintMessage (3, "Call remove problem");
RemoveProblem (mesh, domain); RemoveProblem (mesh, domain);
mesh.FindOpenElements(domain); mesh.FindOpenElements(domain);
GetOpenElements( mesh, domain )->Save("open_"+ToString(domain)+"_"+ToString(cntsteps)+".vol");
} }
else else
{ {
@ -553,6 +552,7 @@ namespace netgen
for (auto pi : Range(first_new_pi, m.Points().Range().Next())) for (auto pi : Range(first_new_pi, m.Points().Range().Next()))
pmap[pi] = mesh.AddPoint(m[pi]); pmap[pi] = mesh.AddPoint(m[pi]);
for ( auto el : m.VolumeElements() ) for ( auto el : m.VolumeElements() )
{ {
for (auto i : Range(el.GetNP())) for (auto i : Range(el.GetNP()))
@ -589,7 +589,7 @@ namespace netgen
// extern double teterrpow; // extern double teterrpow;
MESHING3_RESULT MeshVolume (const MeshingParameters & mp, Mesh& mesh3d) MESHING3_RESULT MeshVolume (const MeshingParameters & mp, Mesh& mesh3d)
{ {
static Timer t("MeshVolume1"); RegionTimer reg(t); static Timer t("MeshVolume"); RegionTimer reg(t);
mesh3d.Compress(); mesh3d.Compress();
for (auto bl : mp.boundary_layers) for (auto bl : mp.boundary_layers)