610 |
|
ambcorral(AMBHEMI *hp, FVECT uv[2], const double r0, const double r1) |
611 |
|
{ |
612 |
|
const double max_d = 1.0/(minarad*ambacc + 0.001); |
613 |
< |
const double ang_res = 0.5*PI/(hp->ns-1); |
614 |
< |
const double ang_step = ang_res/((int)(16/PI*ang_res) + (1+FTINY)); |
613 |
> |
const double ang_res = 0.5*PI/hp->ns; |
614 |
> |
const double ang_step = ang_res/((int)(16/PI*ang_res) + 1.01); |
615 |
|
double avg_d = 0; |
616 |
|
uint32 flgs = 0; |
617 |
|
FVECT vec; |
642 |
|
if ((r0*r0*u*u + r1*r1*v*v) * ap->d*ap->d <= u*u + v*v) |
643 |
|
continue; /* occluder outside ellipse */ |
644 |
|
ang = atan2a(v, u); /* else set direction flags */ |
645 |
< |
for (a1 = ang-.5*ang_res; a1 <= ang+.5*ang_res; a1 += ang_step) |
645 |
> |
for (a1 = ang-ang_res; a1 <= ang+ang_res; a1 += ang_step) |
646 |
|
flgs |= 1L<<(int)(16/PI*(a1 + 2.*PI*(a1 < 0))); |
647 |
|
} |
648 |
|
/* add low-angle incident (< 20deg) */ |