ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-12-24 21:10:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

points on 2D boundary

Browse Source
This commit is contained in:
Joachim Schoeberl 2010-03-22 21:41:07 +00:00
parent b5126a479b
commit 7862ac9a12
2 changed files with 6 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
libsrc/geom2d/genmesh2d.cpp
View File

@ -39,6 +39,7 @@ namespace netgen
geometry.PartitionBoundary (h, *mesh); geometry.PartitionBoundary (h, *mesh);
// marks mesh points for hp-refinement // marks mesh points for hp-refinement
for (int i = 0; i < geometry.GetNP(); i++) for (int i = 0; i < geometry.GetNP(); i++)
if (geometry.GetPoint(i).hpref) if (geometry.GetPoint(i).hpref)

6
libsrc/geom2d/spline.hpp
View File

@ -345,8 +345,12 @@ void SplineSeg<D> :: Partition (double h, double elto0,
while (j < curvepoints.Size() && (l >= curvepoints[j] || i == n)) while (j < curvepoints.Size() && (l >= curvepoints[j] || i == n))
{ {
frac = (curvepoints[j]-lold) / (l-lold); frac = (curvepoints[j]-lold) / (l-lold);
mark = pold + frac * (p-pold);
edgelength = i*dt + (frac-1)*dt; edgelength = i*dt + (frac-1)*dt;
// mark = pold + frac * (p-pold);
mark = GetPoint (edgelength);
// cout << "mark = " << mark << " =?= " << GetPoint (edgelength) << endl;
{ {
PointIndex pi1 = -1, pi2 = -1; PointIndex pi1 = -1, pi2 = -1;