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

Merge branch 'meshpnt_surfacemesh' into 'master'

Browse Source 
Meshpnt surfacemesh

See merge request jschoeberl/netgen!400
This commit is contained in:
Joachim Schöberl 2021-07-25 07:47:25 +00:00
commit 62217b7db6
1 changed files with 15 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
libsrc/meshing/meshclass.cpp
View File

@ -4863,7 +4863,21 @@ namespace netgen
Box<3> box (Box<3>::EMPTY_BOX); Box<3> box (Box<3>::EMPTY_BOX);
for (auto pi : surfelements[sei].PNums()) for (auto pi : surfelements[sei].PNums())
box.Add (points[pi]); box.Add (points[pi]);
auto & el = surfelements[sei];
if(el.IsCurved() && curvedelems->IsSurfaceElementCurved(sei))
{
netgen::Point<2> lami [4] = {netgen::Point<2>(0.5,0), netgen::Point<2>(0,0.5), netgen::Point<2>(0.5,0.5), netgen::Point<2>(1./3,1./3)};
for (auto lam : lami)
{
netgen::Point<3> x;
Mat<3,2> Jac;
curvedelems->CalcSurfaceTransformation(lam,sei,x,Jac);
box.Add (x);
}
box.Scale(1.2);
}
elementsearchtree -> Insert (box, sei+1); elementsearchtree -> Insert (box, sei+1);
} }
} }