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

Remove search tree in face mapping

Browse Source
This commit is contained in:
Matthias Hochsteger 2024-03-06 11:20:20 +01:00
parent bb7a3fe692
commit 7e4f171b16
2 changed files with 9 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
libsrc/meshing/basegeom.cpp
View File

@ -22,7 +22,7 @@ namespace netgen
ArrayMem<int, 1> points; ArrayMem<int, 1> points;
tree.GetIntersecting(p, p, points); tree.GetIntersecting(p, p, points);
if(points.Size()==0) if(points.Size()==0)
throw Exception("cannot find mapped point"); throw Exception("cannot find mapped point " + ToString(p));
return points[0]; return points[0];
} }
@ -944,11 +944,12 @@ namespace netgen
} }
bool have_identifications = false; bool have_identifications = false;
std::map<tuple<PointIndex, int>, PointIndex> mapto;
for(auto & face : faces) for(auto & face : faces)
if(face->primary != face.get()) if(face->primary != face.get())
{ {
have_identifications = true; have_identifications = true;
MapSurfaceMesh(mesh, *face); MapSurfaceMesh(mesh, *face, mapto);
} }
// identify points on faces // identify points on faces
@ -992,7 +993,8 @@ namespace netgen
if(ident.from == face.get()) if(ident.from == face.get())
for(auto pi : pi_of_face[face->nr]) for(auto pi : pi_of_face[face->nr])
{ {
auto pi_other = tree.Find(ident.trafo(mesh[pi])); auto pi_primary = ident.from->primary->nr == ident.from->nr ? pi : mapto[{pi, ident.to->primary->nr}];
auto pi_other = ident.to->primary->nr == ident.to->nr ? pi_primary : mapto[{pi_primary, ident.to->nr}];
mesh_ident.Add(pi, pi_other, ident.name, ident.type); mesh_ident.Add(pi, pi_other, ident.name, ident.type);
} }
} }
@ -1002,7 +1004,7 @@ namespace netgen
multithread.task = savetask; multithread.task = savetask;
} }
void NetgenGeometry :: MapSurfaceMesh( Mesh & mesh, const GeometryFace & dst ) const void NetgenGeometry :: MapSurfaceMesh( Mesh & mesh, const GeometryFace & dst, std::map<tuple<PointIndex, int>, PointIndex> & mapto ) const
{ {
static Timer timer("MapSurfaceMesh"); static Timer timer("MapSurfaceMesh");
RegionTimer rt(timer); RegionTimer rt(timer);
@ -1083,6 +1085,8 @@ namespace netgen
pmap[pi] = mesh.AddPoint(trafo(mesh[pi]), 1, SURFACEPOINT); pmap[pi] = mesh.AddPoint(trafo(mesh[pi]), 1, SURFACEPOINT);
} }
sel_new[i] = pmap[pi]; sel_new[i] = pmap[pi];
mapto[{pi, dst.nr}] = pmap[pi];
mapto[{pmap[pi], src.nr}] = pi;
} }
if(do_invert.IsTrue()) if(do_invert.IsTrue())
sel_new.Invert(); sel_new.Invert();

2
libsrc/meshing/basegeom.hpp
View File

@ -243,7 +243,7 @@ namespace netgen
virtual void MeshSurface(Mesh& mesh, const MeshingParameters& mparam) const; virtual void MeshSurface(Mesh& mesh, const MeshingParameters& mparam) const;
virtual bool MeshFace(Mesh& mesh, const MeshingParameters& mparam, virtual bool MeshFace(Mesh& mesh, const MeshingParameters& mparam,
int nr, FlatArray<int, PointIndex> glob2loc) const; int nr, FlatArray<int, PointIndex> glob2loc) const;
virtual void MapSurfaceMesh( Mesh & mesh, const GeometryFace & dst ) const; virtual void MapSurfaceMesh( Mesh & mesh, const GeometryFace & dst, std::map<tuple<PointIndex, int>, PointIndex> & mapto) const;
virtual void OptimizeSurface(Mesh& mesh, const MeshingParameters& mparam) const; virtual void OptimizeSurface(Mesh& mesh, const MeshingParameters& mparam) const;
virtual void FinalizeMesh(Mesh& mesh) const; virtual void FinalizeMesh(Mesh& mesh) const;