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

Fix Mesh::GetSubMesh for 2d meshes

Browse Source
This commit is contained in:
Matthias Hochsteger 2024-09-05 14:46:57 +02:00
parent bb3c3ff565
commit a504372f82
1 changed files with 20 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
libsrc/meshing/meshclass.cpp
View File

@ -7284,17 +7284,27 @@ namespace netgen
regex regex_faces(faces);
regex regex_domains(domains);
for(auto dom : Range(ndomains))
{
if(regex_match(mesh.GetMaterial(dom+1), regex_domains))
keep_domain.SetBit(dom);
}
if(dimension == 3) {
for(auto dom : Range(ndomains))
if(regex_match(mesh.GetMaterial(dom+1), regex_domains))
keep_domain.SetBit(dom);
for(auto fi : Range(nfaces))
{
auto & fd = mesh.FaceDescriptors()[fi];
if (keep_domain[fd.DomainIn()] || keep_domain[fd.DomainOut()] || regex_match(fd.GetBCName(), regex_faces))
keep_face.SetBit(fd.BCProperty());
for(auto fi : Range(nfaces))
{
auto & fd = mesh.FaceDescriptors()[fi];
if (regex_match(fd.GetBCName(), regex_faces)
|| keep_domain[fd.DomainIn()] || keep_domain[fd.DomainOut()])
keep_face.SetBit(fd.BCProperty());
}
}
else {
for(auto fi : Range(nfaces))
{
auto & fd = mesh.FaceDescriptors()[fi];
auto mat = GetMaterial(fd.BCProperty());
if (regex_match(mat, regex_faces))
keep_face.SetBit(fd.BCProperty());
}
}
auto filter_elements = [&keep_point](auto & elements, auto & keep_region)