417 |
|
} |
418 |
|
/* determine sample type/bin */ |
419 |
|
if (tolower(curparams.hemis[0]) == 'u' | curparams.hemis[0] == '1') { |
420 |
< |
binv = "0"; /* uniform sampling -- one bin */ |
420 |
> |
sprintf(sbuf, "if(-Dx*%g-Dy*%g-Dz*%g,0,-1)", |
421 |
> |
curparams.nrm[0], curparams.nrm[1], curparams.nrm[2]); |
422 |
> |
binv = savqstr(sbuf); |
423 |
> |
nbins = "1"; /* uniform sampling -- one bin */ |
424 |
|
uniform = 1; |
425 |
|
} else if (tolower(curparams.hemis[0]) == 's' && |
426 |
|
tolower(curparams.hemis[1]) == 'c') { |
519 |
|
{ |
520 |
|
int i; |
521 |
|
|
522 |
< |
uva[1][0] = 0.5 - frandom(); |
520 |
< |
uva[1][1] = 0.5 - frandom(); |
521 |
< |
uva[1][2] = 0.5 - frandom(); |
522 |
< |
for (i = 3; i--; ) |
523 |
< |
if ((-0.6 < nrm[i]) & (nrm[i] < 0.6)) |
524 |
< |
break; |
525 |
< |
if (i < 0) { |
522 |
> |
if (!getperpendicular(uva[0], nrm)) { |
523 |
|
fputs(progname, stderr); |
524 |
|
fputs(": bad surface normal in make_axes!\n", stderr); |
525 |
|
exit(1); |
526 |
|
} |
527 |
< |
uva[1][i] = 1.0; |
531 |
< |
VCROSS(uva[0], uva[1], nrm); |
532 |
< |
normalize(uva[0]); |
533 |
< |
VCROSS(uva[1], nrm, uva[0]); |
527 |
> |
fcross(uva[1], nrm, uva[0]); |
528 |
|
} |
529 |
|
|
530 |
|
/* Illegal sender surfaces end up here */ |
786 |
|
alt = (row+samp3[1])*RAH; |
787 |
|
azi = (2.*PI)*(col+samp3[2]-.5)/rnaz(row); |
788 |
|
duvw[2] = cos(alt); /* measured from horizon */ |
789 |
< |
duvw[0] = tcos(azi)*duvw[2]; |
790 |
< |
duvw[1] = tsin(azi)*duvw[2]; |
789 |
> |
duvw[0] = tsin(azi)*duvw[2]; |
790 |
> |
duvw[1] = tcos(azi)*duvw[2]; |
791 |
|
duvw[2] = sqrt(1. - duvw[2]*duvw[2]); |
792 |
|
for (i = 3; i--; ) |
793 |
|
orig_dir[1][i] = -duvw[0]*p->udir[i] - |