760 |
|
} |
761 |
|
while (n--) { /* stratified sampling */ |
762 |
|
SDmultiSamp(samp3, 3, (n+frandom())/sampcnt); |
763 |
+ |
if (row >= RowMax-1) /* avoid crowding at zenith */ |
764 |
+ |
samp3[1] *= samp3[1]; |
765 |
|
alt = (row+samp3[1])*RAH; |
766 |
|
azi = (2.*PI)*(col+samp3[2]-.5)/rnaz(row); |
767 |
|
duvw[2] = cos(alt); /* measured from horizon */ |
768 |
|
duvw[0] = tsin(azi)*duvw[2]; |
769 |
< |
duvw[1] = tcos(azi)*duvw[2]; |
769 |
> |
duvw[1] = -tcos(azi)*duvw[2]; |
770 |
|
duvw[2] = sqrt(1. - duvw[2]*duvw[2]); |
771 |
|
for (i = 3; i--; ) |
772 |
|
orig_dir[1][i] = -duvw[0]*p->udir[i] - |