solstice-solver

commit 03c8f1052ed31a29fbd48725f60019db488773ef
parent 32e47a1dec694729df6a2d032840b3be9f9e9416
Author: Christophe Coustet <christophe.coustet@meso-star.com>
Date:   Fri, 13 Jan 2017 17:30:38 +0100

Fix #6: Increase the number of slices in quadric meshes.

Diffstat:
M
src/ssol_shape.c
 | 
4
++--

1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/ssol_shape.c b/src/ssol_shape.c
@@ -955,7 +955,7 @@ ssol_punched_surface_setup
         / (4.0 * psurf->quadric->data.parabol.focal);
       z[1] = (upper[0] * upper[0] + upper[1] * upper[1])
         / (4.0 * psurf->quadric->data.parabol.focal);
-      nslices = MMIN(50, (size_t)(1 + MMAX(z[0], z[1]) * 4));
+      nslices = MMIN(50, (size_t)(3 + sqrt(MMAX(z[0], z[1])) * 6));
       break;
     }
     case SSOL_QUADRIC_PARABOLIC_CYLINDER: {
@@ -964,7 +964,7 @@ ssol_punched_surface_setup
         (4.0 * psurf->quadric->data.parabolic_cylinder.focal);
       z[1] = (upper[1] * upper[1]) /
         (4.0 * psurf->quadric->data.parabolic_cylinder.focal);
-      nslices = MMIN(50, (size_t)(1 + MMAX(z[0], z[1]) * 4));
+      nslices = MMIN(50, (size_t)(3 + sqrt(MMAX(z[0], z[1])) * 6));
       break;
     }
     default: FATAL("Unreachable code\n"); break;