Changeset 10211 for branches/HeuristicLab.Analysis.AlgorithmBehavior/qhull-2012.1/src/HeuristicLab.qhull/HeuristicLab_qhull.c

Timestamp:

12/09/13 22:42:34 (10 years ago)

Author:

ascheibe

Message:

#1886 added a wrapper for calculating the convex hull using qhull

File:

• branches/HeuristicLab.Analysis.AlgorithmBehavior/qhull-2012.1/src/HeuristicLab.qhull/HeuristicLab_qhull.c (modified) (3 diffs)

branches/HeuristicLab.Analysis.AlgorithmBehavior/qhull-2012.1/src/HeuristicLab.qhull/HeuristicLab_qhull.c

@@ -31 +31 @@
 #endif
   __declspec(dllexport) double qhull_volume(int dim, int numpoints, double *data);
-  void print_summary (void);
+  __declspec(dllexport) int* qhull_convex_hull(int dim, int numpoints, double *data, int *retCode, int* nrOfFacets); 
+  __declspec(dllexport) void qhull_free(int* data);
 #ifdef __cplusplus
 }
@@ -73 +74 @@
   if (!exitcode) {                      
     print_summary();  
+    qh_vertexneighbors();
     volume = qh totvol;
   } else {
-   return -1.0;
+    return -1.0;
   }
   qh_freeqhull(!qh_ALL);            
@@ -85 +87 @@
 }
 
+
+int* qhull_convex_hull(int dim, int numpoints, double *data, int *retCode, int* nrOfFacets) {
+  boolT ismalloc= False;   
+  char flags[250];       
+  FILE *outfile= stdout;   
+  FILE *errfile= stdout;   
+  facetT *facet;           
+  int curlong, totlong;    
+  int i,j, exitcode;
+  int* result = NULL; 
+
+#if qh_QHpointer  
+  if (qh_qh){
+    printf ("QH6233: Qhull link error.  The global variable qh_qh was not initialized\n\
+            to NULL by global.c.  Please compile user_eg.c with -Dqh_QHpointer_dllimport\n\
+            as well as -Dqh_QHpointer, or use libqhullstatic, or use a different tool chain.\n\n");
+    *retCode = -1;
+    return NULL;
+  }
+#endif
+
+  sprintf (flags, "qhull s Tv Qt FA");
+
+  exitcode= qh_new_qhull (dim, numpoints, data, ismalloc, flags, outfile, errfile);
+  if (!exitcode) {      
+    facetT *facet; 
+    vertexT *vertex, **vertexp;
+    unsigned int n = qh num_facets;
+    *nrOfFacets = n;
+
+    print_summary(); 
+    qh_vertexneighbors();
+
+    result = (int *) malloc(n*dim*sizeof(int));     
+
+    i=0;
+    FORALLfacets {
+      j=0;      
+      FOREACHvertex_ (facet->vertices) {        
+        if (j >= dim)
+          fprintf (errfile,"extra vertex %d of facet %d = %d",j++,i,/*1+*/qh_pointid(vertex->point));    //TODO: handle in return value
+        else
+          result[i+n*j++] = qh_pointid(vertex->point);
+      }
+      if (j < dim) fprintf (errfile,"facet %d only has %d vertices",i,j);
+      i++;
+    }
+  } else { 
+    *retCode = -1;
+    return NULL;
+  }
+  qh_freeqhull(!qh_ALL);            
+  qh_memfreeshort (&curlong, &totlong); 
+  if (curlong || totlong)
+    fprintf (errfile, "qhull internal warning (HeuristicLab_qhull, #1): did not free %d bytes of long memory (%d pieces)\n", totlong, curlong);
+
+  return result;
+}
+
+
+void qhull_free(int* data) {
+  free(data);
+}
+
+