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 point type of geo vertices (FIXEDPOINT) -> locked points

Browse Source
This commit is contained in:
Christopher Lackner 2022-03-01 13:23:06 +01:00
parent aeb6a17255
commit 016b1692e2
2 changed files with 12 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
libsrc/meshing/basegeom.cpp
View File

@ -529,6 +529,8 @@ namespace netgen
tree.Insert(mesh[pi], pi);
vert2meshpt[vert->GetHash()] = pi;
mesh[pi].Singularity(vert->properties.hpref);
mesh[pi].SetType(FIXEDPOINT);
mesh.AddLockedPoint(pi);
if(vert->properties.name)
{

12
libsrc/meshing/boundarylayer.cpp
View File

@ -276,7 +276,11 @@ namespace netgen
bool have_single_segments = HaveSingleSegments(mesh);
Array<Segment> segments;
if(have_single_segments)
{
segments = BuildSegments(mesh);
mesh.SetNextMajorTimeStamp();
mesh.UpdateTopology();
}
else
segments = mesh.LineSegments();
@ -492,25 +496,29 @@ namespace netgen
if(seg.edgenr-1 == edgenr && mesh[seg[0]].Type() == FIXEDPOINT)
{
points.Append(seg[0]);
points.Append(seg[1]);
break;
}
}
while(true)
{
bool point_found = false;
for(auto si : meshtopo.GetVertexSegments(points.Last()))
{
const auto& seg = mesh[si];
if(seg.edgenr-1 != edgenr)
continue;
if(seg[0] == points.Last() &&
(points.Size() < 2 || points[points.Size()-2] !=seg[1]))
if(seg[0] == points.Last() && points[points.Size()-2] !=seg[1])
{
points.Append(seg[1]);
point_found = true;
break;
}
}
if(mesh[points.Last()].Type() == FIXEDPOINT)
break;
if(!point_found)
throw Exception(string("Could not find connected list of line segements for edge ") + edgenr);
}
// tangential part of growth vectors