ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-01-19 09:20:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

Don't do mesh smoothing at non-tet elements

Browse Source
This commit is contained in:
Matthias Hochsteger 2024-01-19 09:21:22 +01:00
parent 4e30c0dd64
commit 10377be7d7
1 changed files with 15 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
libsrc/meshing/smoothing3.cpp
View File

@ -335,6 +335,20 @@ namespace netgen
{ {
static Timer tim("PointFunction - build elementsonpoint table"); RegionTimer reg(tim); static Timer tim("PointFunction - build elementsonpoint table"); RegionTimer reg(tim);
Array<bool, PointIndex> bad_point(points.Size());
bad_point = false;
// Don't optimize if point is adjacent to a non-tet element
ParallelForRange(elements.Range(), [&] (auto myrange)
{
for(auto ei : myrange)
{
const auto & el = elements[ei];
if(el.NP()!=4)
for(auto pi : el.PNums())
bad_point[pi] = true;
}
});
elementsonpoint = ngcore::CreateSortedTable<ElementIndex, PointIndex>( elements.Range(), elementsonpoint = ngcore::CreateSortedTable<ElementIndex, PointIndex>( elements.Range(),
[&](auto & table, ElementIndex ei) [&](auto & table, ElementIndex ei)
{ {
@ -344,6 +358,7 @@ namespace netgen
return; return;
for (PointIndex pi : el.PNums()) for (PointIndex pi : el.PNums())
if(!bad_point[pi])
table.Add (pi, ei); table.Add (pi, ei);
}, points.Size()); }, points.Size());
} }