ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/hd/holofile.c
(Generate patch)

Comparing ray/src/hd/holofile.c (file contents):
Revision 3.28 by gwlarson, Mon Nov 9 17:10:53 1998 UTC vs.
Revision 3.33 by gwlarson, Sat Jan 9 09:16:14 1999 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1998 Silicon Graphics, Inc. */
1 > /* Copyright (c) 1999 Silicon Graphics, Inc. */
2  
3   #ifndef lint
4   static char SCCSid[] = "$SunId$ SGI";
# Line 21 | Line 21 | static char SCCSid[] = "$SunId$ SGI";
21   #ifndef PCTFREE
22   #define PCTFREE         15      /* maximum fraction to free (%) */
23   #endif
24 < #ifndef MAXFRAG
25 < #define MAXFRAG         32767   /* maximum fragments/file to track (0==inf) */
24 > #ifndef MAXFRAGB
25 > #define MAXFRAGB        16      /* fragment blocks/file to track (0==inf) */
26   #endif
27 + #ifndef FF_DEFAULT
28 +                                /* when to free a beam fragment */
29 + #define FF_DEFAULT      (FF_ALLOC|FF_WRITE|FF_KILL)
30 + #endif
31 + #ifndef MINDIRSEL
32 +                                /* minimum directory seek length */
33 + #define MINDIRSEL       (4*BUFSIZ/sizeof(BEAMI))
34 + #endif
35  
36   #ifndef BSD
37   #define write   writebuf        /* safe i/o routines */
38   #define read    readbuf
39   #endif
40  
41 < #define FRAGBLK         256     /* number of fragments to allocate at a time */
41 > #define FRAGBLK         512     /* number of fragments to allocate at a time */
42  
43 + int     hdfragflags = FF_DEFAULT;       /* tells when to free fragments */
44   unsigned        hdcachesize = CACHESIZE*1024*1024;      /* target cache size */
45 < unsigned long   hdclock;        /* clock value */
45 > unsigned long   hdclock;                /* clock value */
46  
47   HOLO    *hdlist[HDMAX+1];       /* holodeck pointers (NULL term.) */
48  
# Line 204 | Line 213 | int    i;
213                  }
214                  copystruct(hp->dirseg+j, hp->dirseg+(j-1));
215          }
216 <        mindist = nbeams(hp);                   /* find closest neighbors */
217 <        for (j = hp->dirty; --j; )
218 <                if (hp->dirseg[j].s - (hp->dirseg[j-1].s + hp->dirseg[j-1].n)
219 <                                < mindist) {
220 <                        mindist = hp->dirseg[j].s -
216 >        do {                            /* check neighbors */
217 >                mindist = nbeams(hp);           /* find closest */
218 >                for (j = hp->dirty; --j; )
219 >                        if (hp->dirseg[j].s -
220 >                                        (hp->dirseg[j-1].s + hp->dirseg[j-1].n)
221 >                                        < mindist) {
222 >                                minpos = j;
223 >                                mindist = hp->dirseg[j].s -
224                                          (hp->dirseg[j-1].s + hp->dirseg[j-1].n);
225 <                        minpos = j;
226 <                }
227 <        if (hp->dirty > MAXDIRSE || mindist <= BUFSIZ/sizeof(BEAMI)) {
225 >                        }
226 >                                                /* good enough? */
227 >                if (hp->dirty <= MAXDIRSE && mindist > MINDIRSEL)
228 >                        break;
229                  j = minpos - 1;                 /* coalesce neighbors */
230                  if (hp->dirseg[j].s + hp->dirseg[j].n <
231 <                                hp->dirseg[minpos].s + hp->dirseg[minpos].n) {
231 >                                hp->dirseg[minpos].s + hp->dirseg[minpos].n)
232                          hp->dirseg[j].n = hp->dirseg[minpos].s +
233                                          hp->dirseg[minpos].n - hp->dirseg[j].s;
234 <                }
235 <                hp->dirty--;                    /* close the gap */
223 <                for (j = minpos; j < hp->dirty; j++)
234 >                hp->dirty--;
235 >                while (++j < hp->dirty)         /* close the gap */
236                          copystruct(hp->dirseg+j, hp->dirseg+(j+1));
237 <        }
237 >        } while (mindist <= MINDIRSEL);
238   }
239  
240  
# Line 342 | Line 354 | int    nr;                     /* number of new rays desired */
354          RAYVAL  *p;
355          int     n;
356  
357 <        if (nr <= 0)
358 <                return(NULL);
359 <        if (i < 1 | i > nbeams(hp))
348 <                error(CONSISTENCY, "bad beam index given to hdnewrays");
357 >        if (nr <= 0) return(NULL);
358 >        CHECK(i < 1 | i > nbeams(hp),
359 >                        CONSISTENCY, "bad beam index given to hdnewrays");
360          if (hp->bl[i] != NULL)
361                  hp->bl[i]->tick = hdclock;      /* preempt swap */
362          if (hdcachesize > 0 && hdmemuse(0) >= hdcachesize)
# Line 354 | Line 365 | int    nr;                     /* number of new rays desired */
365                  n = hp->bi[i].nrd + nr;
366                  hp->bl[i] = (BEAM *)hdrealloc(NULL, hdbsiz(n), "hdnewrays");
367                  blglob(hp)->nrm += n;
368 <                if (n = hp->bl[i]->nrm = hp->bi[i].nrd) {
368 >                if ((n = hp->bl[i]->nrm = hp->bi[i].nrd)) {
369                          errno = 0;
370                          if (lseek(hp->fd, hp->bi[i].fo, 0) < 0)
371                                  error(SYSTEM, "seek error on holodeck file");
# Line 368 | Line 379 | int    nr;                     /* number of new rays desired */
379                                  hdbsiz(hp->bl[i]->nrm + nr), "hdnewrays");
380                  blglob(hp)->nrm += nr;
381          }
382 +        if (hdfragflags&FF_ALLOC && hp->bi[i].nrd)
383 +                hdfreefrag(hp, i);              /* relinquish old fragment */
384          p = hdbray(hp->bl[i]) + hp->bl[i]->nrm;
385          hp->bl[i]->nrm += nr;                   /* update in-core structure */
386          bzero((char *)p, nr*sizeof(RAYVAL));
# Line 383 | Line 396 | register int   i;
396   {
397          register int    n;
398  
399 <        if (i < 1 | i > nbeams(hp))
400 <                error(CONSISTENCY, "bad beam index given to hdgetbeam");
399 >        CHECK(i < 1 | i > nbeams(hp),
400 >                        CONSISTENCY, "bad beam index given to hdgetbeam");
401          if (hp->bl[i] == NULL) {                /* load from disk */
402                  if (!(n = hp->bi[i].nrd))
403                          return(NULL);
# Line 398 | Line 411 | register int   i;
411                  n *= sizeof(RAYVAL);
412                  if (read(hp->fd, (char *)hdbray(hp->bl[i]), n) != n)
413                          error(SYSTEM, "error reading beam from holodeck file");
414 +                if (hdfragflags&FF_READ)
415 +                        hdfreefrag(hp, i);      /* relinquish old frag. */
416          }
417          blglob(hp)->tick = hp->bl[i]->tick = hdclock++; /* update LRU clock */
418          return(hp->bl[i]);
# Line 435 | Line 450 | int    (*bf)();                /* callback function (optional) */
450                                          /* precheck consistency */
451          if (n <= 0) return;
452          for (i = n; i--; )
453 <                if (hb[i].h == NULL || hb[i].b < 1 | hb[i].b > nbeams(hb[i].h))
453 >                if (hb[i].h==NULL || hb[i].b<1 | hb[i].b>nbeams(hb[i].h))
454                          error(CONSISTENCY, "bad beam in hdloadbeams");
455                                          /* sort list for optimal access */
456          qsort((char *)hb, n, sizeof(HDBEAMI), hdfilord);
# Line 468 | Line 483 | int    (*bf)();                /* callback function (optional) */
483   }
484  
485  
486 < hdfreefrag(fd, bi)                      /* free a file fragment */
487 < int     fd;
488 < register BEAMI  *bi;
486 > int
487 > hdfreefrag(hp, i)                       /* free a file fragment */
488 > HOLO    *hp;
489 > int     i;
490   {
491 +        register BEAMI  *bi = &hp->bi[i];
492          register struct fraglist        *f;
493          register int    j, k;
494  
495 <        if (bi->nrd == 0)
496 <                return;
497 < #ifdef DEBUG
498 <        if (fd < 0 | fd >= nhdfragls || !hdfragl[fd].nlinks)
499 <                error(CONSISTENCY, "bad file descriptor in hdfreefrag");
483 < #endif
484 <        f = &hdfragl[fd];
495 >        if (bi->nrd <= 0)
496 >                return(0);
497 >        DCHECK(hp->fd < 0 | hp->fd >= nhdfragls || !hdfragl[hp->fd].nlinks,
498 >                        CONSISTENCY, "bad file descriptor in hdfreefrag");
499 >        f = &hdfragl[hp->fd];
500          if (f->nfrags % FRAGBLK == 0) { /* delete empty remnants */
501                  for (j = k = 0; k < f->nfrags; j++, k++) {
502                          while (f->fi[k].nrd == 0)
# Line 494 | Line 509 | register BEAMI *bi;
509                  f->nfrags = j;
510          }
511          j = f->nfrags++;                /* allocate a slot in free list */
512 < #if MAXFRAG
513 <        if (j >= MAXFRAG-1)
514 <                f->nfrags--;
512 > #if MAXFRAGB
513 >        if (j >= MAXFRAGB*FRAGBLK) {
514 >                f->nfrags = j--;        /* stop list growth */
515 >                if (bi->nrd <= f->fi[j].nrd)
516 >                        return(0);      /* new one no better than discard */
517 >        }
518   #endif
519 <        if (j % FRAGBLK == 0) {         /* more free list space */
519 >        if (j % FRAGBLK == 0) {         /* more (or less) free list space */
520                  register BEAMI  *newp;
521                  if (f->fi == NULL)
522                          newp = (BEAMI *)malloc((j+FRAGBLK)*sizeof(BEAMI));
# Line 507 | Line 525 | register BEAMI *bi;
525                                          (j+FRAGBLK)*sizeof(BEAMI));
526                  if (newp == NULL) {
527                          f->nfrags--;    /* graceful failure */
528 <                        return;
528 >                        return(0);
529                  }
530                  f->fi = newp;
531          }
# Line 534 | Line 552 | register BEAMI *bi;
552                          }
553                          break;
554                  }
555 +        biglob(hp)->nrd -= bi->nrd;             /* tell fragment it's free */
556 +        bi->nrd = 0;
557 +        bi->fo = 0;
558 +        return(1);
559   }
560  
561  
562 + int
563 + hdfragOK(fd, listlen, listsiz)  /* get fragment list status for file */
564 + int     fd;
565 + int     *listlen;
566 + register int4   *listsiz;
567 + {
568 +        register struct fraglist        *f;
569 +        register int    i;
570 +
571 +        if (fd < 0 | fd >= nhdfragls || !(f = &hdfragl[fd])->nlinks)
572 +                return(0);              /* listless */
573 +        if (listlen != NULL)
574 +                *listlen = f->nfrags;
575 +        if (listsiz != NULL)
576 +                for (i = f->nfrags, *listsiz = 0; i--; )
577 +                        *listsiz += f->fi[i].nrd;
578 + #if MAXFRAGB
579 +        return(f->nfrags < MAXFRAGB*FRAGBLK);
580 + #else
581 +        return(1);                      /* list never fills up */
582 + #endif
583 + }
584 +
585 +
586   long
587   hdallocfrag(fd, nrays)          /* allocate a file fragment */
588   int     fd;
589   unsigned int4   nrays;
590   {
591          register struct fraglist        *f;
592 <        register int    j, k;
592 >        register int    j;
593          long    nfo;
594  
595          if (nrays == 0)
596                  return(-1L);
597 < #ifdef DEBUG
598 <        if (fd < 0 | fd >= nhdfragls || !hdfragl[fd].nlinks)
553 <                error(CONSISTENCY, "bad file descriptor in hdallocfrag");
554 < #endif
597 >        DCHECK(fd < 0 | fd >= nhdfragls || !hdfragl[fd].nlinks,
598 >                        CONSISTENCY, "bad file descriptor in hdallocfrag");
599          f = &hdfragl[fd];
600 <        k = -1;                         /* find closest-sized fragment */
601 <        for (j = f->nfrags; j-- > 0; )
602 <                if (f->fi[j].nrd >= nrays &&
603 <                                (k < 0 || f->fi[j].nrd < f->fi[k].nrd))
560 <                        if (f->fi[k=j].nrd == nrays)
561 <                                break;
562 <        if (k < 0) {                    /* no fragment -- extend file */
600 >        for (j = f->nfrags; j-- > 0; )  /* first fit algorithm */
601 >                if (f->fi[j].nrd >= nrays)
602 >                        break;
603 >        if (j < 0) {                    /* no fragment -- extend file */
604                  nfo = f->flen;
605                  f->flen += nrays*sizeof(RAYVAL);
606          } else {                        /* else use fragment */
607 <                nfo = f->fi[k].fo;
608 <                f->fi[k].fo += nrays*sizeof(RAYVAL);
609 <                f->fi[k].nrd -= nrays;
607 >                nfo = f->fi[j].fo;
608 >                f->fi[j].fo += nrays*sizeof(RAYVAL);
609 >                f->fi[j].nrd -= nrays;
610          }
611          return(nfo);
612   }
# Line 582 | Line 623 | register int   i;
623                                          /* check file status */
624          if (hdfragl[hp->fd].writerr)
625                  return(-1);
626 < #ifdef DEBUG
627 <        if (i < 1 | i > nbeams(hp))
587 <                error(CONSISTENCY, "bad beam index in hdsyncbeam");
588 < #endif
626 >        DCHECK(i < 1 | i > nbeams(hp),
627 >                        CONSISTENCY, "bad beam index in hdsyncbeam");
628                                          /* is current fragment OK? */
629          if (hp->bl[i] == NULL || (nrays = hp->bl[i]->nrm) == hp->bi[i].nrd)
630                  return(0);
631 <        if (hp->bi[i].nrd)              /* relinquish old fragment */
632 <                hdfreefrag(hp->fd, &hp->bi[i]);
631 >        if (hdfragflags&FF_WRITE && hp->bi[i].nrd)
632 >                hdfreefrag(hp, i);      /* relinquish old fragment */
633          if (nrays) {                    /* get and write new fragment */
634                  nfo = hdallocfrag(hp->fd, nrays);
635                  errno = 0;
# Line 605 | Line 644 | register int   i;
644                  hp->bi[i].fo = nfo;
645          } else
646                  hp->bi[i].fo = 0L;
647 <        biglob(hp)->nrd += nrays - hp->bi[i].nrd;
609 <        hp->bi[i].nrd = nrays;
647 >        biglob(hp)->nrd += hp->bi[i].nrd = nrays;
648          markdirty(hp, i);               /* section directory now out of date */
649          return(1);
650   }
# Line 634 | Line 672 | register int   i;
672                                  nchanged += hdfreebeam(hp, i);
673                  return(nchanged);
674          }
675 < #ifdef DEBUG
676 <        if (i < 1 | i > nbeams(hp))
639 <                error(CONSISTENCY, "bad beam index to hdfreebeam");
640 < #endif
675 >        DCHECK(i < 1 | i > nbeams(hp),
676 >                        CONSISTENCY, "bad beam index to hdfreebeam");
677          if (hp->bl[i] == NULL)
678                  return(0);
679                                          /* check for additions */
# Line 656 | Line 692 | hdkillbeam(hp, i)              /* delete beam from holodeck */
692   register HOLO   *hp;
693   register int    i;
694   {
659 –        static BEAM     emptybeam;
695          int     nchanged;
696  
697          if (hp == NULL) {               /* clobber all holodecks */
# Line 670 | Line 705 | register int   i;
705                  for (i = nbeams(hp); i > 0; i--)
706                          if (hp->bi[i].nrd > 0 || hp->bl[i] != NULL)
707                                  nchanged += hdkillbeam(hp, i);
708 < #ifdef DEBUG
709 <                if (biglob(hp)->nrd != 0 | blglob(hp)->nrm != 0)
675 <                        error(CONSISTENCY, "bad beam count in hdkillbeam");
676 < #endif
708 >                DCHECK(biglob(hp)->nrd != 0 | blglob(hp)->nrm != 0,
709 >                                CONSISTENCY, "bad beam count in hdkillbeam");
710                  return(nchanged);
711          }
712 < #ifdef DEBUG
713 <        if (i < 1 | i > nbeams(hp))
681 <                error(CONSISTENCY, "bad beam index to hdkillbeam");
682 < #endif
712 >        DCHECK(i < 1 | i > nbeams(hp),
713 >                        CONSISTENCY, "bad beam index to hdkillbeam");
714          if (hp->bl[i] != NULL) {        /* free memory */
715                  blglob(hp)->nrm -= nchanged = hp->bl[i]->nrm;
716                  free((char *)hp->bl[i]);
717          } else
718                  nchanged = hp->bi[i].nrd;
719 <        if (hp->bi[i].nrd) {            /* free file fragment */
720 <                hp->bl[i] = &emptybeam;
721 <                hdsyncbeam(hp, i);
719 >        if (hp->bi[i].nrd) {
720 >                if (hdfragflags&FF_KILL)
721 >                        hdfreefrag(hp, i);
722 >                hp->bi[i].nrd = 0;      /* make sure it's gone */
723 >                hp->bi[i].fo = 0L;
724          }
725          hp->bl[i] = NULL;
726          return(nchanged);

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines