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

Fix boundary layer thicknesses

Browse Source 
Height values should be interpreted per-layer, not total (accumulated) height
This commit is contained in:
Matthias Hochsteger 2024-07-26 18:09:37 +02:00
parent 433247b09b
commit 24a0f47856
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
libsrc/meshing/boundarylayer.cpp
View File

@ -1201,11 +1201,12 @@ void BoundaryLayerTool ::InsertNewElements(
Array<PointIndex>& new_points) { Array<PointIndex>& new_points) {
Point<3> p = mesh[pi]; Point<3> p = mesh[pi];
PointIndex pi_last = pi; PointIndex pi_last = pi;
double height = 0.0;
for (auto i : Range(params.heights)) { for (auto i : Range(params.heights)) {
// p += params.heights[i] * growth_vector; height += params.heights[i];
auto pi_new = mesh.AddPoint(p); auto pi_new = mesh.AddPoint(p);
new_points.Append(pi_new); new_points.Append(pi_new);
growth_vector_map[pi_new] = {&growth_vector, params.heights[i]}; growth_vector_map[pi_new] = {&growth_vector, height};
if (special_boundary_points.count(pi) > 0) mesh.AddLockedPoint(pi_new); if (special_boundary_points.count(pi) > 0) mesh.AddLockedPoint(pi_new);
pi_last = pi_new; pi_last = pi_new;
} }