mirror of
https://github.com/NGSolve/netgen.git
synced 2024-09-23 02:02:49 +05:00
d08e2daa06
If the mesh contains quads, the edge swapping algorithm switches to generic improve, which introduces quads everywhere. This is not intended if one domain contains a tensor product mesh. Thus, call the optimizer for each face if mesh contains quads but mp.quad is not set.
111 lines
2.7 KiB
C++
111 lines
2.7 KiB
C++
#include <mystdlib.h>
|
|
#include "meshing.hpp"
|
|
|
|
namespace netgen
|
|
{
|
|
|
|
DLL_HEADER void Optimize2d (Mesh & mesh, MeshingParameters & mp)
|
|
{
|
|
static Timer timer("optimize2d"); RegionTimer reg(timer);
|
|
|
|
mesh.CalcSurfacesOfNode();
|
|
|
|
bool secondorder = mesh.GetNP() > mesh.GetNV();
|
|
|
|
|
|
if (secondorder)
|
|
{
|
|
for (SurfaceElementIndex ei = 0; ei < mesh.GetNSE(); ei++)
|
|
mesh[ei].SetType(TRIG);
|
|
}
|
|
mesh.Compress();
|
|
|
|
bool optimize_swap_separate_faces = false;
|
|
if(!mp.quad)
|
|
{
|
|
bool mixed = false;
|
|
ParallelFor( Range(mesh.GetNSE()), [&] (auto i) NETGEN_LAMBDA_INLINE
|
|
{
|
|
if (mesh[SurfaceElementIndex(i)].GetNP() != 3)
|
|
mixed = true;
|
|
});
|
|
if(mixed)
|
|
optimize_swap_separate_faces = true;
|
|
}
|
|
|
|
const char * optstr = mp.optimize2d.c_str();
|
|
int optsteps = mp.optsteps2d;
|
|
|
|
for (int i = 1; i <= optsteps; i++)
|
|
for (size_t j = 1; j <= strlen(optstr); j++)
|
|
{
|
|
if (multithread.terminate) break;
|
|
switch (optstr[j-1])
|
|
{
|
|
case 's':
|
|
{ // topological swap
|
|
MeshOptimize2d meshopt(mesh);
|
|
meshopt.SetMetricWeight (mp.elsizeweight);
|
|
|
|
if(optimize_swap_separate_faces)
|
|
{
|
|
for(auto i : Range(1, mesh.GetNFD()+1))
|
|
{
|
|
meshopt.SetFaceIndex(i);
|
|
meshopt.EdgeSwapping (0);
|
|
}
|
|
}
|
|
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))
|
|
{
|
|
meshopt.SetFaceIndex(i);
|
|
meshopt.EdgeSwapping (1);
|
|
}
|
|
}
|
|
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;
|
|
}
|
|
default:
|
|
cerr << "Optimization code " << optstr[j-1] << " not defined" << endl;
|
|
}
|
|
}
|
|
mesh.Compress(); // better: compress in individual steps, if necessary
|
|
if (secondorder)
|
|
{
|
|
mesh.GetGeometry()->GetRefinement().MakeSecondOrder(mesh);
|
|
}
|
|
}
|
|
|
|
}
|