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

fix boundarylayer 2d bug

Browse Source
This commit is contained in:
Matthias Hochsteger 2022-10-03 17:22:44 +02:00
parent 19453911c6
commit 92fb557314
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
libsrc/meshing/boundarylayer.cpp
View File

@ -1382,7 +1382,7 @@ namespace netgen
{ {
Box<3> box{Box<3>::EMPTY_BOX}; Box<3> box{Box<3>::EMPTY_BOX};
for(const auto & seg : mesh.LineSegments()) for(const auto & seg : mesh.LineSegments())
if (seg.si == domain) if (seg.domin == domain || seg.domout == domain)
for (auto pi : {seg[0], seg[1]}) for (auto pi : {seg[0], seg[1]})
box.Add(mesh[pi]); box.Add(mesh[pi]);
@ -1395,11 +1395,13 @@ namespace netgen
PointIndex cnt = PointIndex::BASE; PointIndex cnt = PointIndex::BASE;
auto p2sel = mesh.CreatePoint2SurfaceElementTable(); auto p2sel = mesh.CreatePoint2SurfaceElementTable();
Array<Segment> domain_segments;
PointGeomInfo gi; PointGeomInfo gi;
gi.u = 0.0;
gi.v = 0.0;
gi.trignum = domain; gi.trignum = domain;
for(auto seg : mesh.LineSegments()) for(auto seg : mesh.LineSegments())
{ {
if(seg.domin == domain || seg.domout == domain)
for (auto pi : {seg[0], seg[1]}) for (auto pi : {seg[0], seg[1]})
if (compress[pi]==PointIndex{PointIndex::INVALID}) if (compress[pi]==PointIndex{PointIndex::INVALID})
{ {