ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-01-13 14:40:35 +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:
Hochsteger, Matthias 2024-09-27 16:29:47 +02:00 committed by Schöberl, Joachim
parent 13d962acdd
commit 61bed581ec
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);
Array<bool, PointIndex> non_tet_points(points.Size());
non_tet_points = 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())
non_tet_points[pi] = true;
}
});
elementsonpoint = ngcore::CreateSortedTable<ElementIndex, PointIndex>( elements.Range(),
[&](auto & table, ElementIndex ei)
{
@ -344,6 +358,7 @@ namespace netgen
return;
for (PointIndex pi : el.PNums())
if(!non_tet_points[pi])
table.Add (pi, ei);
}, points.Size());
}