special point 0-point fixes

This commit is contained in:
Joachim Schoeberl 2011-09-01 07:35:06 +00:00
parent d60abeaa77
commit ba0319d388

View File

@ -16,7 +16,6 @@
*/ */
// #undef DEVELOP
// #define DEVELOP // #define DEVELOP
@ -279,7 +278,7 @@ namespace netgen
dynamic_cast<const Plane*> (geometry->GetSurface(locsurf[k2])), dynamic_cast<const Plane*> (geometry->GetSurface(locsurf[k2])),
dynamic_cast<const Plane*> (geometry->GetSurface(locsurf[k3])), dynamic_cast<const Plane*> (geometry->GetSurface(locsurf[k3])),
pts); pts);
for (int j = 0; j < pts.Size(); j++) for (int j = 0; j < pts.Size(); j++)
if (Dist (pts[j], box.Center()) < box.Diam()/2) if (Dist (pts[j], box.Center()) < box.Diam()/2)
{ {
@ -645,7 +644,7 @@ namespace netgen
{ {
box.GetSubBox (i, sbox); box.GetSubBox (i, sbox);
sbox.Increase (1e-4 * sbox.Diam()); sbox.Increase (1e-4 * sbox.Diam());
sbox.CalcDiamCenter();
Solid * redsol = sol -> GetReducedSolid (sbox); Solid * redsol = sol -> GetReducedSolid (sbox);
if (redsol) if (redsol)