| 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); |