ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-01-27 13:20:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

rounding num segments per edge

Browse Source
This commit is contained in:
Joachim Schöberl 2019-02-07 23:09:34 +01:00
parent ead955cf04
commit 966c2960fb
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
libsrc/geom2d/genmesh2d.cpp
View File

@ -55,7 +55,8 @@ namespace netgen
sum += dt / fun; sum += dt / fun;
} }
int nel = int (sum+1); int nel = int (sum+0.5);
if (nel == 0) nel = 1;
fperel = sum / nel; fperel = sum / nel;
points.Append (0); points.Append (0);
@ -129,6 +130,7 @@ namespace netgen
Point3d oldmark3(oldmark(0), oldmark(1), 0); Point3d oldmark3(oldmark(0), oldmark(1), 0);
double h = mesh.GetH (Point<3> (oldmark(0), oldmark(1), 0)); double h = mesh.GetH (Point<3> (oldmark(0), oldmark(1), 0));
cout << "edge h = " << h << endl;
Vec<3> v (1e-4*h, 1e-4*h, 1e-4*h); Vec<3> v (1e-4*h, 1e-4*h, 1e-4*h);
searchtree.GetIntersecting (oldmark3 - v, oldmark3 + v, locsearch); searchtree.GetIntersecting (oldmark3 - v, oldmark3 + v, locsearch);