| 168 |
|
OBJREC *op; |
| 169 |
|
{ |
| 170 |
|
double maxrad2; |
| 171 |
< |
double d2; |
| 171 |
> |
double d; |
| 172 |
|
register int i, j; |
| 173 |
|
register FACE *f; |
| 174 |
|
|
| 175 |
|
f = getface(op); |
| 176 |
+ |
if (f->area == 0.) |
| 177 |
+ |
return(0.); |
| 178 |
|
for (i = 0; i < 3; i++) { |
| 179 |
|
ocent[i] = 0.; |
| 180 |
|
for (j = 0; j < f->nv; j++) |
| 181 |
|
ocent[i] += VERTEX(f,j)[i]; |
| 182 |
|
ocent[i] /= (double)f->nv; |
| 183 |
|
} |
| 184 |
< |
if (f->area == 0.) |
| 185 |
< |
return(0.); |
| 184 |
> |
d = DOT(ocent,f->norm); |
| 185 |
> |
for (i = 0; i < 3; i++) |
| 186 |
> |
ocent[i] += (f->offset - d)*f->norm[i]; |
| 187 |
|
maxrad2 = 0.; |
| 188 |
|
for (j = 0; j < f->nv; j++) { |
| 189 |
< |
d2 = dist2(VERTEX(f,j), ocent); |
| 190 |
< |
if (d2 > maxrad2) |
| 191 |
< |
maxrad2 = d2; |
| 189 |
> |
d = dist2(VERTEX(f,j), ocent); |
| 190 |
> |
if (d > maxrad2) |
| 191 |
> |
maxrad2 = d; |
| 192 |
|
} |
| 193 |
|
return(maxrad2); |
| 194 |
|
} |