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

Curve new elements in HPRefinement

Browse Source 
Use correct check for curvature of TRIGs in VSSolution
This commit is contained in:
Matthias Hochsteger 2020-02-11 14:42:16 +01:00
parent bac288ab56
commit dff550a48c
2 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
libsrc/meshing/hprefinement.cpp
View File

@ -1407,6 +1407,7 @@ namespace netgen
el.SetIndex(hpel.index); el.SetIndex(hpel.index);
if(setorders) if(setorders)
el.SetOrder(act_ref+1,act_ref+1,0); el.SetOrder(act_ref+1,act_ref+1,0);
el.SetCurved(true);
mesh.AddSurfaceElement(el); mesh.AddSurfaceElement(el);
break; break;
} }
@ -1422,6 +1423,7 @@ namespace netgen
el.hp_elnr = i; el.hp_elnr = i;
if(setorders) if(setorders)
el.SetOrder(act_ref+1,act_ref+1,act_ref+1); el.SetOrder(act_ref+1,act_ref+1,act_ref+1);
el.SetCurved(true);
mesh.AddVolumeElement(el); mesh.AddVolumeElement(el);
break; break;
} }

2
libsrc/visualization/vssolution.cpp
View File

@ -1651,7 +1651,7 @@ namespace netgen
// NgProfiler::StopTimer(timer1c); // NgProfiler::StopTimer(timer1c);
#else #else
bool curved = curv.IsSurfaceElementCurved(sei); bool curved = (*mesh)[sei].IsCurved();
for (int iy = 0, ii = 0; iy <= n; iy++) for (int iy = 0, ii = 0; iy <= n; iy++)
for (int ix = 0; ix <= n-iy; ix++, ii++) for (int ix = 0; ix <= n-iy; ix++, ii++)