From b37a3e6cf62bc34e82ae2f985d9b70fefd63b6b3 Mon Sep 17 00:00:00 2001 From: Christopher Lackner Date: Fri, 17 Jul 2020 18:00:38 +0200 Subject: [PATCH] comment code for non orthogonal boundarylayers (not working if multiple surfaces come together) --- libsrc/meshing/boundarylayer.cpp | 36 ++++++++++++++++---------------- 1 file changed, 18 insertions(+), 18 deletions(-) diff --git a/libsrc/meshing/boundarylayer.cpp b/libsrc/meshing/boundarylayer.cpp index 953a96c7..710251a6 100644 --- a/libsrc/meshing/boundarylayer.cpp +++ b/libsrc/meshing/boundarylayer.cpp @@ -238,24 +238,24 @@ namespace netgen } // project growthvector on surface for inner angles - for(const auto& sel : mesh.SurfaceElements()) - if(!blp.surfid.Contains(sel.GetIndex())) - { - auto n = GetSurfaceNormal(mesh, sel); - for(auto pi : sel.PNums()) - { - if(growthvectors[pi].Length2() == 0.) - continue; - auto& g = growthvectors[pi]; - auto ng = n * g; - auto gg = g * g; - auto nn = n * n; - if(fabs(ng*ng-nn*gg) < 1e-12 || fabs(ng) < 1e-12) continue; - auto a = -ng*ng/(ng*ng-nn * gg); - auto b = ng*gg/(ng*ng-nn*gg); - g += a*g + b*n; - } - } + // for(const auto& sel : mesh.SurfaceElements()) + // if(!blp.surfid.Contains(sel.GetIndex())) + // { + // auto n = GetSurfaceNormal(mesh, sel); + // for(auto pi : sel.PNums()) + // { + // if(growthvectors[pi].Length2() == 0.) + // continue; + // auto& g = growthvectors[pi]; + // auto ng = n * g; + // auto gg = g * g; + // auto nn = n * n; + // if(fabs(ng*ng-nn*gg) < 1e-12 || fabs(ng) < 1e-12) continue; + // auto a = -ng*ng/(ng*ng-nn * gg); + // auto b = ng*gg/(ng*ng-nn*gg); + // g += a*g + b*n; + // } + // } if (!blp.grow_edges) {