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

do not restrict mesh size on identified edges

Browse Source
This commit is contained in:
Christopher Lackner 2022-01-18 18:24:09 +01:00 committed by mhochsteger@cerbsim.com
parent fd258a611f
commit 3ab8808fa3
2 changed files with 20 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
libsrc/meshing/basegeom.hpp
View File

@ -196,6 +196,8 @@ namespace netgen
size_t GetNFaces() const { return faces.Size(); } size_t GetNFaces() const { return faces.Size(); }
const GeometryFace & GetFace(int i) const { return *faces[i]; } const GeometryFace & GetFace(int i) const { return *faces[i]; }
const GeometryEdge & GetEdge(int i) const { return *edges[i]; }
const GeometryVertex & GetVertex(int i) const { return *vertices[i]; }
void Clear(); void Clear();

18
libsrc/occ/occgenmesh.cpp
View File

@ -495,6 +495,24 @@ namespace netgen
continue; continue;
} }
bool is_identified_edge = false;
// TODO: change to use hash value
const auto& gedge = geom.GetEdge(geom.edge_map.at(e.TShape()));
auto& v0 = gedge.GetStartVertex();
auto& v1 = gedge.GetEndVertex();
for(auto & ident : v0.identifications)
{
auto other = ident.from == &v0 ? ident.to : ident.from;
if(other->nr == v1.nr && ident.type == Identifications::CLOSESURFACES)
{
is_identified_edge = true;
break;
}
}
if(is_identified_edge)
continue;
double localh = len/mparam.segmentsperedge; double localh = len/mparam.segmentsperedge;
double s0, s1; double s0, s1;