From c0d9f3d95b5f6fcc247165ec806cc7113d0c51bc Mon Sep 17 00:00:00 2001 From: Michael Neunteufel Date: Sun, 25 Jul 2021 07:47:25 +0000 Subject: [PATCH] Meshpnt surfacemesh --- libsrc/meshing/meshclass.cpp | 16 +++++++++++++++- 1 file changed, 15 insertions(+), 1 deletion(-) diff --git a/libsrc/meshing/meshclass.cpp b/libsrc/meshing/meshclass.cpp index 6d49d55f..4b446db5 100644 --- a/libsrc/meshing/meshclass.cpp +++ b/libsrc/meshing/meshclass.cpp @@ -4863,7 +4863,21 @@ namespace netgen Box<3> box (Box<3>::EMPTY_BOX); for (auto pi : surfelements[sei].PNums()) 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); } }