ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-02-23 20:45:38 +05:00
Code Issues Packages Projects Releases Wiki Activity

fix Powell-Sabin split

Browse Source
This commit is contained in:
Joachim Schoeberl 2024-02-06 21:41:15 +01:00
parent 1c48e552e9
commit adbdf194e0
1 changed files with 9 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
libsrc/meshing/hprefinement.cpp
View File

@ -699,6 +699,8 @@ namespace netgen
{ {
INDEX_2 i2(el.pnums[hprs->splitedges[j][0]-1], INDEX_2 i2(el.pnums[hprs->splitedges[j][0]-1],
el.pnums[hprs->splitedges[j][1]-1]); el.pnums[hprs->splitedges[j][1]-1]);
if (fac1 == 0.5) i2.Sort();
if (!newpts.Used (i2)) if (!newpts.Used (i2))
{ {
Point<3> np; Point<3> np;
@ -787,7 +789,8 @@ namespace netgen
{ {
INDEX_2 i2(el.pnums[hprs->splitedges[j][0]-1], INDEX_2 i2(el.pnums[hprs->splitedges[j][0]-1],
el.pnums[hprs->splitedges[j][1]-1]); el.pnums[hprs->splitedges[j][1]-1]);
if (fac1 == 0.5) i2.Sort();
int npi = newpts.Get(i2); int npi = newpts.Get(i2);
newpnums[hprs->splitedges[j][2]-1] = npi; newpnums[hprs->splitedges[j][2]-1] = npi;
@ -1990,8 +1993,11 @@ bool CheckSingularities(Mesh & mesh, INDEX_2_HASHTABLE<int> & edges, INDEX_2_HAS
} }
case HP_SEGM: case HP_SEGM:
{ {
hpel.type = ClassifySegm(hpel, edges, edgepoint_dom, cornerpoint, edgepoint, if (split == SPLIT_HP)
faces, face_edges, surf_edges, facepoint); hpel.type = ClassifySegm(hpel, edges, edgepoint_dom, cornerpoint, edgepoint,
faces, face_edges, surf_edges, facepoint);
else if (split == SPLIT_POWELL)
hpel.type = HP_SEGM_SINGCORNERL;
dd = 1; dd = 1;
break; break;
} }