Changeset 329

Timestamp:

03/26/08 15:14:16 (8 months ago)

Author:

powell

Message:

Little debugging to eliminate the former segfault.

Files:

• trunk/matml/src/ternary/qhull.c (modified) (3 diffs)

trunk/matml/src/ternary/qhull.c

@@ -50 +50 @@
       qpoints[3*i]   = (coordT) points[i].C2;
       qpoints[3*i+1] = (coordT) points[i].C3;
+      /* TODO: find corners and clip at that plane, not just G=0 */
       qpoints[3*i+2] = (coordT) fmin (points[i].G, 0.);
     }
@@ -63 +64 @@
   printf ("%d total facets\n", *numfacets);
   printf ("facetverts: pre-realloc 0x%lx, ", *facetverts);
-  if (!(*facetverts = realloc (*facetverts, *numfacets * sizeof (int))))
+  if (!(*facetverts = realloc (*facetverts, *numfacets * 3 * sizeof (int))))
     { printf ("qhullCalcHull: could not reallocate memory for vertices\n");
       return -1; }
@@ -83 +84 @@
 
           if (j<3)
-            facetverts [3*i+j] = thepoint;
+            (*facetverts) [3*i+j] = thepoint;
 #ifdef DEBUG
           printf ("  vertex %d, coords %g %g %g\n", thepoint,