ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-11-11 16:49:16 +05:00
Code Issues Packages Projects Releases Wiki Activity

Reset MeshTopology before 2d optimize

Browse Source
This commit is contained in:
Matthias Hochsteger 2023-09-18 13:58:43 +02:00
parent 8efc3a2fe2
commit 960fe8a342
3 changed files with 14 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
libsrc/geom2d/genmesh2d.cpp
View File

@ -5,7 +5,7 @@ namespace netgen
{
// extern DLL_HEADER MeshingParameters mparam;
extern void Optimize2d (Mesh & mesh, MeshingParameters & mp);
extern void Optimize2d (Mesh & mesh, MeshingParameters & mp, int faceindex=0);

21
libsrc/meshing/meshfunc2d.cpp
View File

@ -4,7 +4,7 @@
namespace netgen
{
DLL_HEADER void Optimize2d (Mesh & mesh, MeshingParameters & mp)
DLL_HEADER void Optimize2d (Mesh & mesh, MeshingParameters & mp, int faceindex)
{
static Timer timer("optimize2d"); RegionTimer reg(timer);
@ -33,19 +33,25 @@ namespace netgen
optimize_swap_separate_faces = true;
}
if(faceindex)
optimize_swap_separate_faces = false;
const char * optstr = mp.optimize2d.c_str();
int optsteps = mp.optsteps2d;
// reset topology
mesh.GetTopology() = MeshTopology(mesh);
for (int i = 1; i <= optsteps; i++)
for (size_t j = 1; j <= strlen(optstr); j++)
{
if (multithread.terminate) break;
MeshOptimize2d meshopt(mesh);
meshopt.SetMetricWeight (mp.elsizeweight);
meshopt.SetFaceIndex(faceindex);
switch (optstr[j-1])
{
case 's':
{ // topological swap
MeshOptimize2d meshopt(mesh);
meshopt.SetMetricWeight (mp.elsizeweight);
if(optimize_swap_separate_faces)
{
@ -57,16 +63,12 @@ namespace netgen
}
else
{
meshopt.SetFaceIndex(0);
meshopt.EdgeSwapping (0);
}
break;
}
case 'S':
{ // metric swap
MeshOptimize2d meshopt(mesh);
meshopt.SetMetricWeight (mp.elsizeweight);
if(optimize_swap_separate_faces)
{
for(auto i : Range(1, mesh.GetNFD()+1))
@ -77,22 +79,17 @@ namespace netgen
}
else
{
meshopt.SetFaceIndex(0);
meshopt.EdgeSwapping (1);
}
break;
}
case 'm':
{
MeshOptimize2d meshopt(mesh);
meshopt.SetMetricWeight (mp.elsizeweight);
meshopt.ImproveMesh(mp);
break;
}
case 'c':
{
MeshOptimize2d meshopt(mesh);
meshopt.SetMetricWeight (mp.elsizeweight);
meshopt.CombineImprove();
break;
}

8
libsrc/meshing/python_mesh.cpp
View File

@ -78,7 +78,7 @@ namespace netgen
{
extern bool netgen_executable_started;
extern shared_ptr<NetgenGeometry> ng_geometry;
extern void Optimize2d (Mesh & mesh, MeshingParameters & mp);
extern void Optimize2d (Mesh & mesh, MeshingParameters & mp, int faceindex=0);
}
@ -1283,7 +1283,7 @@ DLL_HEADER void ExportNetgenMeshing(py::module &m)
OptimizeVolume (mp, self);
}, py::arg("mp"), py::call_guard<py::gil_scoped_release>())
.def ("OptimizeMesh2d", [](Mesh & self, MeshingParameters* pars)
.def ("OptimizeMesh2d", [](Mesh & self, MeshingParameters* pars, int faceindex)
{
self.CalcLocalH(0.5);
MeshingParameters mp;
@ -1291,8 +1291,8 @@ DLL_HEADER void ExportNetgenMeshing(py::module &m)
else mp.optsteps2d = 5;
if(!self.GetGeometry())
throw Exception("Cannot optimize surface mesh without geometry!");
Optimize2d (self, mp);
}, py::arg("mp")=nullptr, py::call_guard<py::gil_scoped_release>())
Optimize2d (self, mp, faceindex);
}, py::arg("mp")=nullptr, py::arg("faceindex")=0, py::call_guard<py::gil_scoped_release>())
.def ("Refine", FunctionPointer
([](Mesh & self, bool adaptive)