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 meshing of close surfaces (use identifications of master mesh)

Browse Source
This commit is contained in:
Matthias Hochsteger 2021-06-14 09:41:57 +02:00
parent 1cf3d2a21c
commit b83d73e919
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
libsrc/meshing/meshfunc.cpp
View File

@ -10,7 +10,7 @@ namespace netgen
extern const char * pyramidrules2[]; extern const char * pyramidrules2[];
extern const char * hexrules[]; extern const char * hexrules[];
void MeshDomain(Mesh & mesh3d, const MeshingParameters & c_mp, int k) void MeshDomain(Mesh & mesh3d, const MeshingParameters & c_mp, int k, const Identifications & identifications)
{ {
MeshingParameters mp = c_mp; // copy mp to change them here MeshingParameters mp = c_mp; // copy mp to change them here
NgArray<INDEX_2> connectednodes; NgArray<INDEX_2> connectednodes;
@ -103,10 +103,10 @@ namespace netgen
for (int i = 1; i <= connectednodes.Size(); i++) for (int i = 1; i <= connectednodes.Size(); i++)
meshing.AddConnectedPair (connectednodes.Get(i)); meshing.AddConnectedPair (connectednodes.Get(i));
*/ */
for (int nr = 1; nr <= mesh3d.GetIdentifications().GetMaxNr(); nr++) for (int nr = 1; nr <= identifications.GetMaxNr(); nr++)
if (mesh3d.GetIdentifications().GetType(nr) != Identifications::PERIODIC) if (identifications.GetType(nr) != Identifications::PERIODIC)
{ {
mesh3d.GetIdentifications().GetPairs (nr, connectednodes); identifications.GetPairs (nr, connectednodes);
for (auto pair : connectednodes) for (auto pair : connectednodes)
meshing.AddConnectedPair (pair); meshing.AddConnectedPair (pair);
} }
@ -327,9 +327,9 @@ namespace netgen
ParallelFor(Range(1, mesh3d.GetNDomains()+1), [&](int k) ParallelFor(Range(1, mesh3d.GetNDomains()+1), [&](int k)
{ {
if(k==1) if(k==1)
MeshDomain(mesh3d, mp, k); MeshDomain(mesh3d, mp, k, mesh3d.GetIdentifications());
else else
MeshDomain(meshes[k-2], mp, k); MeshDomain(meshes[k-2], mp, k, mesh3d.GetIdentifications());
}); });
MergeMeshes(mesh3d, meshes, first_new_pi); MergeMeshes(mesh3d, meshes, first_new_pi);