| 105 |
|
theta = PI/180.0/2.0 * so->oargs.farg[3]; |
| 106 |
|
if (theta <= FTINY) |
| 107 |
|
objerror(so, USER, "zero size"); |
| 108 |
< |
src->ss = theta >= PI/4 ? 1.0 : tan(theta); |
| 108 |
> |
src->ss = theta >= PI/4.0 ? 1.0 : tan(theta); |
| 109 |
|
src->ss2 = 2.0*PI * (1.0 - cos(theta)); |
| 110 |
|
} |
| 111 |
|
|
| 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 |
|
} |
| 288 |
|
|
| 289 |
|
d = DOT(sp->aim, nrm); |
| 290 |
|
if (d > FTINY) /* center in front? */ |
| 291 |
< |
return(0); |
| 291 |
> |
return(1); |
| 292 |
|
/* else check horizon */ |
| 293 |
|
d1 = 1. - sp->siz/(2.*PI); |
| 294 |
< |
return(1.-FTINY-d*d > d1*d1); |
| 294 |
> |
return(1.-FTINY-d*d < d1*d1); |
| 295 |
> |
} |
| 296 |
> |
|
| 297 |
> |
|
| 298 |
> |
double |
| 299 |
> |
spotdisk(oc, op, sp, pos) /* intersect spot with object op */ |
| 300 |
> |
FVECT oc; |
| 301 |
> |
OBJREC *op; |
| 302 |
> |
register SPOT *sp; |
| 303 |
> |
FVECT pos; |
| 304 |
> |
{ |
| 305 |
> |
FVECT onorm; |
| 306 |
> |
double offs, d, dist; |
| 307 |
> |
register int i; |
| 308 |
> |
|
| 309 |
> |
offs = getplaneq(onorm, op); |
| 310 |
> |
d = -DOT(onorm, sp->aim); |
| 311 |
> |
if (d >= -FTINY && d <= FTINY) |
| 312 |
> |
return(0.); |
| 313 |
> |
dist = (DOT(pos, onorm) - offs)/d; |
| 314 |
> |
if (dist < 0.) |
| 315 |
> |
return(0.); |
| 316 |
> |
for (i = 0; i < 3; i++) |
| 317 |
> |
oc[i] = pos[i] + dist*sp->aim[i]; |
| 318 |
> |
return(sp->siz*dist*dist/PI/(d*d)); |
| 319 |
> |
} |
| 320 |
> |
|
| 321 |
> |
|
| 322 |
> |
double |
| 323 |
> |
beamdisk(oc, op, sp, dir) /* intersect beam with object op */ |
| 324 |
> |
FVECT oc; |
| 325 |
> |
OBJREC *op; |
| 326 |
> |
register SPOT *sp; |
| 327 |
> |
FVECT dir; |
| 328 |
> |
{ |
| 329 |
> |
FVECT onorm; |
| 330 |
> |
double offs, d, dist; |
| 331 |
> |
register int i; |
| 332 |
> |
|
| 333 |
> |
offs = getplaneq(onorm, op); |
| 334 |
> |
d = -DOT(onorm, dir); |
| 335 |
> |
if (d >= -FTINY && d <= FTINY) |
| 336 |
> |
return(0.); |
| 337 |
> |
dist = (DOT(sp->aim, onorm) - offs)/d; |
| 338 |
> |
for (i = 0; i < 3; i++) |
| 339 |
> |
oc[i] = sp->aim[i] + dist*dir[i]; |
| 340 |
> |
return(sp->siz/PI/(d*d)); |
| 341 |
|
} |
| 342 |
|
|
| 343 |
|
|
