64 |
|
for (i = 0; i < 3; i++) |
65 |
|
nr.rdir[i] = r->rdir[i] + 2.*r->rod*r->ron[i]; |
66 |
|
/* source we're aiming for next */ |
67 |
< |
nr.rsrc = source[r->rsrc].sa.svnext; |
67 |
> |
nr.rsrc = source[r->rsrc].sa.sv.sn; |
68 |
|
} else { /* ordinary reflection */ |
69 |
|
FVECT pnorm; |
70 |
|
double pdot; |
84 |
|
mir_proj(pm, o, s, n) /* compute a mirror's projection */ |
85 |
|
MAT4 pm; |
86 |
|
register OBJREC *o; |
87 |
< |
register SRCREC *s; |
87 |
> |
SRCREC *s; |
88 |
|
int n; |
89 |
|
{ |
90 |
|
FVECT nv; |
91 |
|
double od; |
92 |
|
/* get surface normal and offset */ |
93 |
< |
if (sfun[o->otype].of->getpleq == NULL) |
94 |
< |
return(0); /* reject non-planar case */ |
95 |
< |
od = (*sfun[o->otype].of->getpleq)(nv, o); |
93 |
> |
od = getplaneq(nv, o); |
94 |
|
/* check for behind */ |
95 |
|
if (DOT(s->sloc, nv) <= (s->sflags & SDISTANT ? FTINY : od+FTINY)) |
96 |
|
return(0); |
108 |
|
register int i, j; |
109 |
|
/* assign matrix */ |
110 |
|
setident4(m); |
111 |
< |
for (i = 0; i < 3; i++) |
112 |
< |
for (j = 0; j < 3; j++) |
111 |
> |
for (j = 0; j < 3; j++) { |
112 |
> |
for (i = 0; i < 3; i++) |
113 |
|
m[i][j] -= 2.*nv[i]*nv[j]; |
116 |
– |
for (j = 0; j < 3; j++) |
114 |
|
m[3][j] = 2.*offs*nv[j]; |
115 |
+ |
} |
116 |
|
} |