| 58 |
|
FVECT dv |
| 59 |
|
) |
| 60 |
|
{ |
| 61 |
< |
const double cos_thresh = 0.9999995; /* about 3.44 arcminutes */ |
| 62 |
< |
int ii, jj; |
| 63 |
< |
|
| 61 |
> |
double cos_thresh; |
| 62 |
> |
int ii, jj; |
| 63 |
> |
/* min. spacing = 1/4th division */ |
| 64 |
> |
cos_thresh = (PI/4.)/(double)hp->ns; |
| 65 |
> |
cos_thresh = 1. - .5*cos_thresh*cos_thresh; |
| 66 |
> |
/* check existing neighbors */ |
| 67 |
|
for (ii = i-1; ii <= i+1; ii++) { |
| 68 |
|
if (ii < 0) continue; |
| 69 |
|
if (ii >= hp->ns) break; |
| 75 |
|
if (jj >= hp->ns) break; |
| 76 |
|
if ((ii==i) & (jj==j)) continue; |
| 77 |
|
ap = &ambsam(hp,ii,jj); |
| 78 |
< |
if (ap->d <= .5/FHUGE) continue; |
| 78 |
> |
if (ap->d <= .5/FHUGE) |
| 79 |
> |
continue; /* no one home */ |
| 80 |
|
VSUB(avec, ap->p, hp->rp->rop); |
| 81 |
|
dprod = DOT(avec, dv); |
| 82 |
|
if (dprod >= cos_thresh*VLEN(avec)) |
| 83 |
|
return(1); /* collision */ |
| 84 |
|
} |
| 85 |
|
} |
| 86 |
< |
return(0); |
| 86 |
> |
return(0); /* nothing to worry about */ |
| 87 |
|
} |
| 88 |
|
|
| 89 |
|
|
| 126 |
|
/* avoid coincident samples */ |
| 127 |
|
if (!n && ambcollision(hp, i, j, ar.rdir)) { |
| 128 |
|
spt[0] = frandom(); spt[1] = frandom(); |
| 129 |
< |
goto resample; |
| 129 |
> |
goto resample; /* reject this sample */ |
| 130 |
|
} |
| 131 |
|
dimlist[ndims++] = AI(hp,i,j) + 90171; |
| 132 |
|
rayvalue(&ar); /* evaluate ray */ |
