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