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

fix setting boundaries of neighbouring domains in create boundarylayer

Browse Source
This commit is contained in:
Christopher Lackner 2022-07-11 13:46:21 +02:00
parent 68a1a225b6
commit cf992b04da
1 changed files with 17 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
libsrc/meshing/python_mesh.cpp
View File

@ -1199,11 +1199,14 @@ DLL_HEADER void ExportNetgenMeshing(py::module &m)
bool grow_edges, bool limit_growth_vectors) bool grow_edges, bool limit_growth_vectors)
{ {
BoundaryLayerParameters blp; BoundaryLayerParameters blp;
BitArray boundaries(self.GetNFD()+1);
boundaries.Set();
if(int* bc = get_if<int>(&boundary); bc) if(int* bc = get_if<int>(&boundary); bc)
{ {
boundaries.Clear();
for (int i = 1; i <= self.GetNFD(); i++) for (int i = 1; i <= self.GetNFD(); i++)
if(self.GetFaceDescriptor(i).BCProperty() == *bc) if(self.GetFaceDescriptor(i).BCProperty() == *bc)
blp.surfid.Append (i); boundaries.SetBit(i);
} }
else else
{ {
@ -1218,14 +1221,23 @@ DLL_HEADER void ExportNetgenMeshing(py::module &m)
if(dom_pattern) if(dom_pattern)
{ {
regex pattern(*dom_pattern); regex pattern(*dom_pattern);
if((fd.DomainIn() > 0 && regex_match(self.GetMaterial(fd.DomainIn()), pattern)) || (fd.DomainOut() > 0 && regex_match(self.GetMaterial(fd.DomainOut()), pattern))) bool mat1_match = fd.DomainIn() > 0 && regex_match(self.GetMaterial(fd.DomainIn()), pattern);
blp.surfid.Append(i); bool mat2_match = fd.DomainOut() > 0 && regex_match(self.GetMaterial(fd.DomainOut()), pattern);
// if boundary is inner or outer remove from list
if(mat1_match == mat2_match)
boundaries.Clear(i);
// if((fd.DomainIn() > 0 && regex_match(self.GetMaterial(fd.DomainIn()), pattern)) || (fd.DomainOut() > 0 && regex_match(self.GetMaterial(fd.DomainOut()), pattern)))
// boundaries.Clear(i);
// blp.surfid.Append(i);
} }
else // else
blp.surfid.Append(i); // blp.surfid.Append(i);
} }
} }
} }
for(int i = 1; i<=self.GetNFD(); i++)
if(boundaries.Test(i))
blp.surfid.Append(i);
blp.new_mat = material; blp.new_mat = material;
if(project_boundaries.has_value()) if(project_boundaries.has_value())