--- ray/src/hd/rholo3.c 1997/12/31 20:24:11 3.17 +++ ray/src/hd/rholo3.c 1998/01/06 15:08:51 3.19 @@ -1,4 +1,4 @@ -/* Copyright (c) 1997 Silicon Graphics, Inc. */ +/* Copyright (c) 1998 Silicon Graphics, Inc. */ #ifndef lint static char SCCSid[] = "$SunId$ SGI"; @@ -212,9 +212,11 @@ init_global() /* initialize global ray computation * double frac; int i; register int j, k; - /* free old list */ - if (complen > 0) + /* free old list and empty queue */ + if (complen > 0) { free((char *)complist); + done_packets(flush_queue()); + } /* allocate beam list */ complen = 0; for (j = 0; hdlist[j] != NULL; j++) @@ -225,13 +227,14 @@ init_global() /* initialize global ray computation * /* compute beam weights */ k = 0; for (j = 0; hdlist[j] != NULL; j++) { - frac = 512. * hdlist[j]->wg[0] * - hdlist[j]->wg[1] * hdlist[j]->wg[2]; - if (frac < 0.) frac = -frac; + frac = 512. * VLEN(hdlist[j]->wg[0]) * + VLEN(hdlist[j]->wg[1]) * + VLEN(hdlist[j]->wg[2]); for (i = nbeams(hdlist[j]); i > 0; i--) { complist[k].hd = j; complist[k].bi = i; complist[k].nr = frac*beamvolume(hdlist[j], i) + 0.5; + complist[k].nc = bnrays(hdlist[j], i); wtotal += complist[k++].nr; } }