| 80 | 
  | 
  pointT *interiorPoint; | 
| 81 | 
  | 
  int curlong, totlong; | 
| 82 | 
  | 
   | 
| 83 | 
– | 
  Vector3d boxMax; | 
| 84 | 
– | 
  Vector3d boxMin; | 
| 83 | 
  | 
   | 
| 84 | 
  | 
  vector<double> ptArray(numpoints*dim_); | 
| 85 | 
  | 
   | 
| 456 | 
  | 
  // qh_getarea(qh facet_list); | 
| 457 | 
  | 
  //volume_ = qh totvol; | 
| 458 | 
  | 
  // area_ = qh totarea; | 
| 461 | 
– | 
 | 
| 462 | 
– | 
 | 
| 463 | 
– | 
  int index = 0; | 
| 464 | 
– | 
  FORALLvertices { | 
| 465 | 
– | 
    Vector3d point(vertex->point[0], vertex->point[1], vertex->point[2]); | 
| 466 | 
– | 
    if (index == 0) { | 
| 467 | 
– | 
      boxMax = point; | 
| 468 | 
– | 
      boxMin = point; | 
| 469 | 
– | 
    } else { | 
| 470 | 
– | 
      for (int i = 0; i < 3; i++) { | 
| 471 | 
– | 
        boxMax[i] = max(boxMax[i], point[i]); | 
| 472 | 
– | 
        boxMin[i] = min(boxMin[i], point[i]); | 
| 473 | 
– | 
      } | 
| 474 | 
– | 
    } | 
| 475 | 
– | 
    index++; | 
| 476 | 
– | 
  } | 
| 477 | 
– | 
  boundingBox_ = Mat3x3d(0.0); | 
| 478 | 
– | 
  boundingBox_(0,0) = boxMax[0] - boxMin[0]; | 
| 479 | 
– | 
  boundingBox_(1,1) = boxMax[1] - boxMin[1]; | 
| 480 | 
– | 
  boundingBox_(2,2) = boxMax[2] - boxMin[2]; | 
| 459 | 
  | 
 | 
| 460 | 
  | 
  qh_freeqhull(!qh_ALL); | 
| 461 | 
  | 
  qh_memfreeshort(&curlong, &totlong); |