--- ray/src/hd/rholo3.c 1999/03/04 10:29:23 3.31 +++ ray/src/hd/rholo3.c 2003/07/07 17:21:51 3.39 @@ -1,16 +1,14 @@ -/* Copyright (c) 1998 Silicon Graphics, Inc. */ - #ifndef lint -static char SCCSid[] = "$SunId$ SGI"; +static const char RCSid[] = "$Id: rholo3.c,v 3.39 2003/07/07 17:21:51 greg Exp $"; #endif - /* * Routines for tracking beam compuatations */ +#include + #include "rholo.h" #include "view.h" -#include #ifndef NFRAG2CHUNK #define NFRAG2CHUNK 4096 /* number of fragments to start chunking */ @@ -79,7 +77,7 @@ register PACKHEAD *b0, *b1; } -int +void dispbeam(b, hb) /* display a holodeck beam */ register BEAM *b; register HDBEAMI *hb; @@ -92,11 +90,11 @@ register HDBEAMI *hb; if (b->nrm > n) { /* (re)allocate packet holder */ n = b->nrm; if (p == NULL) p = (PACKHEAD *)malloc(packsiz(n)); - else p = (PACKHEAD *)realloc((char *)p, packsiz(n)); + else p = (PACKHEAD *)realloc((void *)p, packsiz(n)); CHECK(p==NULL, SYSTEM, "out of memory in dispbeam"); } /* assign packet fields */ - bcopy((char *)hdbray(b), (char *)packra(p), b->nrm*sizeof(RAYVAL)); + memcpy((void *)packra(p), (void *)hdbray(b), b->nrm*sizeof(RAYVAL)); p->nr = p->nc = b->nrm; for (p->hd = 0; hdlist[p->hd] != hb->h; p->hd++) if (hdlist[p->hd] == NULL) @@ -104,7 +102,7 @@ register HDBEAMI *hb; p->bi = hb->b; disp_packet(p); /* display it */ if (n >= 1024) { /* free ridiculous packets */ - free((char *)p); + free((void *)p); p = NULL; n = 0; } } @@ -122,9 +120,9 @@ int nents; /* search for common members */ for (csm = clist+nents; csm-- > clist; ) csm->nc = -1; - qsort((char *)clist, nents, sizeof(PACKHEAD), beamidcmp); + qsort((void *)clist, nents, sizeof(PACKHEAD), beamidcmp); for (i = 0; i < complen; i++) { - csm = (PACKHEAD *)bsearch((char *)(complist+i), (char *)clist, + csm = (PACKHEAD *)bsearch((void *)(complist+i), (void *)clist, nents, sizeof(PACKHEAD), beamidcmp); if (csm == NULL) continue; @@ -163,14 +161,14 @@ int nents; case BS_NEW: /* new computation set */ listpos = 0; lastin = -1; if (complen) /* free old list */ - free((char *)complist); + free((void *)complist); complist = NULL; if (!(complen = nents)) return; complist = (PACKHEAD *)malloc(nents*sizeof(PACKHEAD)); if (complist == NULL) goto memerr; - bcopy((char *)clist, (char *)complist, nents*sizeof(PACKHEAD)); + memcpy((void *)complist, (void *)clist, nents*sizeof(PACKHEAD)); break; case BS_ADD: /* add to computation set */ case BS_MAX: /* maximum of quantities */ @@ -178,7 +176,7 @@ int nents; if (nents <= 0) return; sortcomplist(); /* sort updated list & new entries */ - qsort((char *)clist, nents, sizeof(PACKHEAD), beamcmp); + qsort((void *)clist, nents, sizeof(PACKHEAD), beamcmp); /* what can't we satisfy? */ for (i = nents, csm = clist; i-- && csm->nr > csm->nc; csm++) ; @@ -197,7 +195,7 @@ int nents; /* merge lists */ mergeclists(newlist, clist, n, complist, complen); if (complen) - free((char *)complist); + free((void *)complist); complist = newlist; complen += n; } @@ -218,7 +216,7 @@ int nents; hbarr[i].b = clist[i].bi; } hdloadbeams(hbarr, nents, dispbeam); - free((char *)hbarr); + free((void *)hbarr); if (hdfragflags&FF_READ) { listpos = 0; lastin = -1; /* need to re-sort list */ @@ -264,7 +262,7 @@ int bi; ambient_list() /* compute ambient beam list */ { - long wtotal = 0; + int32 wtotal, minrt; double frac; int i; register int j, k; @@ -275,7 +273,7 @@ ambient_list() /* compute ambient beam list */ complist = (PACKHEAD *)malloc(complen*sizeof(PACKHEAD)); CHECK(complist==NULL, SYSTEM, "out of memory in ambient_list"); /* compute beam weights */ - k = 0; + k = 0; wtotal = 0; for (j = 0; hdlist[j] != NULL; j++) { frac = 512. * VLEN(hdlist[j]->wg[0]) * VLEN(hdlist[j]->wg[1]) * @@ -290,13 +288,15 @@ ambient_list() /* compute ambient beam list */ } /* adjust sample weights */ if (vdef(DISKSPACE)) - frac = 1024.*1024.*vflt(DISKSPACE) / - (wtotal*sizeof(RAYVAL)); + frac = 1024.*1024.*vflt(DISKSPACE) / (wtotal*sizeof(RAYVAL)); else - frac = 1024.*1024.*16384. / (wtotal*sizeof(RAYVAL)); - if (frac > 1.11 || frac < 0.9) - for (k = complen; k--; ) - complist[k].nr = frac*complist[k].nr + 0.5; + frac = 1024.*1024.*2048. / (wtotal*sizeof(RAYVAL)); + minrt = .02*frac*wtotal/complen + .5; /* heuristic mimimum */ + if (minrt > RPACKSIZ) + minrt = RPACKSIZ; + for (k = complen; k--; ) + if ((complist[k].nr = frac*complist[k].nr + 0.5) < minrt) + complist[k].nr = minrt; listpos = 0; lastin = -1; /* flag initial sort */ } @@ -320,7 +320,7 @@ FILE *fp; normaspect(viewaspect(&curview), &pa, &xr, &yr); viewbeams(&curview, xr, yr, &blist); bundle_set(BS_MAX, blist.bl, blist.nb); - free((char *)blist.bl); + free((void *)blist.bl); } } @@ -330,7 +330,7 @@ init_global() /* initialize global ray computation * register int k; /* free old list and empty queue */ if (complen > 0) { - free((char *)complist); + free((void *)complist); done_packets(flush_queue()); } /* reseed random number generator */ @@ -389,25 +389,25 @@ sortcomplist() /* fix our list order */ #endif } if (lastin < 0 || listpos*4 >= complen*3) - qsort((char *)complist, complen, sizeof(PACKHEAD), beamcmp); + qsort((void *)complist, complen, sizeof(PACKHEAD), beamcmp); else if (listpos) { /* else sort and merge sublist */ list2 = (PACKHEAD *)malloc(listpos*sizeof(PACKHEAD)); CHECK(list2==NULL, SYSTEM, "out of memory in sortcomplist"); - bcopy((char *)complist,(char *)list2,listpos*sizeof(PACKHEAD)); - qsort((char *)list2, listpos, sizeof(PACKHEAD), beamcmp); + memcpy((void *)list2,(void *)complist,listpos*sizeof(PACKHEAD)); + qsort((void *)list2, listpos, sizeof(PACKHEAD), beamcmp); mergeclists(complist, list2, listpos, complist+listpos, complen-listpos); - free((char *)list2); + free((void *)list2); } /* drop satisfied requests */ for (i = complen; i-- && complist[i].nr <= complist[i].nc; ) ; if (i < 0) { - free((char *)complist); + free((void *)complist); complist = NULL; complen = 0; } else if (i < complen-1) { - list2 = (PACKHEAD *)realloc((char *)complist, + list2 = (PACKHEAD *)realloc((void *)complist, (i+1)*sizeof(PACKHEAD)); if (list2 != NULL) complist = list2;