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