383 |
|
curparams.vup[1] = 1; |
384 |
|
} |
385 |
|
/* determine sample type/bin */ |
386 |
< |
if (tolower(curparams.hemis[0]) == 'u' | curparams.hemis[0] == '1') { |
386 |
> |
if ((tolower(curparams.hemis[0]) == 'u') | (curparams.hemis[0] == '1')) { |
387 |
|
sprintf(sbuf, "if(-Dx*%g-Dy*%g-Dz*%g,0,-1)", |
388 |
|
curparams.nrm[0], curparams.nrm[1], curparams.nrm[2]); |
389 |
|
binv = savqstr(sbuf); |
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] - |
873 |
|
curparams.udir[1] *= -1.; |
874 |
|
curparams.udir[2] *= -1.; |
875 |
|
} |
876 |
< |
if (tolower(curparams.hemis[0]) == 'u' | curparams.hemis[0] == '1') |
876 |
> |
if ((tolower(curparams.hemis[0]) == 'u') | (curparams.hemis[0] == '1')) |
877 |
|
curparams.sample_basis = sample_uniform; |
878 |
|
else if (tolower(curparams.hemis[0]) == 's' && |
879 |
|
tolower(curparams.hemis[1]) == 'c') |