792 |
|
alt = (row+samp3[1])*RAH; |
793 |
|
azi = (2.*PI)*(col+samp3[2]-.5)/rnaz(row); |
794 |
|
duvw[2] = cos(alt); /* measured from horizon */ |
795 |
< |
duvw[0] = tcos(azi)*duvw[2]; |
796 |
< |
duvw[1] = tsin(azi)*duvw[2]; |
795 |
> |
duvw[0] = tsin(azi)*duvw[2]; |
796 |
> |
duvw[1] = tcos(azi)*duvw[2]; |
797 |
|
duvw[2] = sqrt(1. - duvw[2]*duvw[2]); |
798 |
|
for (i = 3; i--; ) |
799 |
|
orig_dir[1][i] = -duvw[0]*p->udir[i] - |
844 |
|
|
845 |
|
while (n--) { /* stratified sampling */ |
846 |
|
SDmultiSamp(samp2, 2, (n+frandom())/sampcnt); |
847 |
< |
if (!bo_getvec(duvw, b+samp2[1], kbasis[bi])) |
847 |
> |
if (!bi_getvec(duvw, b+samp2[1], kbasis[bi])) |
848 |
|
return(0); |
849 |
|
for (i = 3; i--; ) |
850 |
|
orig_dir[1][i] = duvw[0]*p->udir[i] + |