--- ray/src/rt/srcsamp.c 1991/11/25 09:52:49 2.2 +++ ray/src/rt/srcsamp.c 1993/03/08 12:37:39 2.5 @@ -15,6 +15,9 @@ static char SCCSid[] = "$SunId$ LBL"; #include "random.h" +static int cyl_partit(), flt_partit(); + + double nextssamp(r, si) /* compute sample for source, rtn. distance */ register RAY *r; /* origin is read, direction is set */ @@ -82,17 +85,17 @@ nextsample: si->dom = source[si->sn].ss2; if (source[si->sn].sflags & SFLAT) { si->dom *= sflatform(si->sn, r->rdir); - si->dom *= (double)(size[SU]*size[SV])/(MAXSPART*MAXSPART); + si->dom *= size[SU]*size[SV]/(MAXSPART*(double)MAXSPART); } else if (source[si->sn].sflags & SCYL) { si->dom *= scylform(si->sn, r->rdir); - si->dom *= (double)size[SU]/MAXSPART; + si->dom *= size[SU]/(double)MAXSPART; } else { - si->dom *= (double)(size[SU]*size[SV]*size[SW]) / - (MAXSPART*MAXSPART*MAXSPART) ; + si->dom *= size[SU]*size[SV]*(double)size[SW] / + (MAXSPART*MAXSPART*(double)MAXSPART) ; } if (source[si->sn].sflags & SDISTANT) return(FHUGE); - if (si->dom <= FTINY) + if (si->dom <= 1e-4) goto nextsample; /* behind source? */ si->dom /= d*d; return(d); /* sample OK, return distance */