| 199 |
|
} |
| 200 |
|
|
| 201 |
|
|
| 202 |
< |
hdlseg(lseg, hp, i) /* compute line segment for beam */ |
| 202 |
> |
hdlseg(lseg, hp, gc) /* compute line segment for beam */ |
| 203 |
|
register int lseg[2][3]; |
| 204 |
|
register HOLO *hp; |
| 205 |
< |
int i; |
| 205 |
> |
GCOORD gc[2]; |
| 206 |
|
{ |
| 207 |
– |
GCOORD gc[2]; |
| 207 |
|
register int k; |
| 208 |
|
|
| 210 |
– |
if (!hdbcoord(gc, hp, i)) /* compute grid coordinates */ |
| 211 |
– |
return(0); |
| 209 |
|
for (k = 0; k < 2; k++) { /* compute end points */ |
| 210 |
|
lseg[k][gc[k].w>>1] = gc[k].w&1 ? hp->grid[gc[k].w>>1]-1 : 0 ; |
| 211 |
|
lseg[k][wg0[gc[k].w]] = gc[k].i[0]; |
| 250 |
|
} |
| 251 |
|
|
| 252 |
|
|
| 253 |
+ |
hdworld(wp, hp, gp) /* compute world coordinates */ |
| 254 |
+ |
register FVECT wp; |
| 255 |
+ |
register HOLO *hp; |
| 256 |
+ |
FVECT gp; |
| 257 |
+ |
{ |
| 258 |
+ |
register double d; |
| 259 |
+ |
|
| 260 |
+ |
d = gp[0]/hp->grid[0]; |
| 261 |
+ |
VSUM(wp, hp->orig, hp->xv[0], d); |
| 262 |
+ |
|
| 263 |
+ |
d = gp[1]/hp->grid[1]; |
| 264 |
+ |
VSUM(wp, wp, hp->xv[1], d); |
| 265 |
+ |
|
| 266 |
+ |
d = gp[2]/hp->grid[2]; |
| 267 |
+ |
VSUM(wp, wp, hp->xv[2], d); |
| 268 |
+ |
} |
| 269 |
+ |
|
| 270 |
+ |
|
| 271 |
|
double |
| 272 |
|
hdray(ro, rd, hp, gc, r) /* compute ray within a beam */ |
| 273 |
|
FVECT ro, rd; /* returned */ |
| 296 |
|
|
| 297 |
|
|
| 298 |
|
double |
| 299 |
< |
hdinter(gc, r, hp, ro, rd) /* compute ray intersection with section */ |
| 299 |
> |
hdinter(gc, r, ed, hp, ro, rd) /* compute ray intersection with section */ |
| 300 |
|
register GCOORD gc[2]; /* returned */ |
| 301 |
|
BYTE r[2][2]; /* returned */ |
| 302 |
+ |
double *ed; /* returned (optional) */ |
| 303 |
|
register HOLO *hp; |
| 304 |
|
FVECT ro, rd; /* rd should be normalized */ |
| 305 |
|
{ |
| 360 |
|
return(FHUGE); /* outside wall */ |
| 361 |
|
r[i][1] = 256. * (d - gc[i].i[1]); |
| 362 |
|
} |
| 363 |
< |
/* return distance from entry point */ |
| 364 |
< |
vt[0] = ro[0] - p[0][0]; |
| 365 |
< |
vt[1] = ro[1] - p[0][1]; |
| 366 |
< |
vt[2] = ro[2] - p[0][2]; |
| 351 |
< |
return(DOT(vt,rd)); |
| 363 |
> |
|
| 364 |
> |
if (ed != NULL) /* assign distance to exit point */ |
| 365 |
> |
*ed = t1; |
| 366 |
> |
return(t0); /* return distance to entry point */ |
| 367 |
|
} |