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.42 by greg, Sat Feb 22 02:07:24 2003 UTC vs.
Revision 3.52 by schorsch, Sun Jul 27 22:12:02 2003 UTC

# Line 7 | Line 7 | static const char      RCSid[] = "$Id$";
7   *      9/30/97 GWLarson
8   */
9  
10 + #include "copyright.h"
11 +
12 + #include <string.h>
13 +
14   #include "holo.h"
15  
16   #ifndef CACHESIZE
17 < #ifdef BIGMEM
14 < #define CACHESIZE       17      /* default cache size (Mbytes, 0==inf) */
15 < #else
17 > #ifdef SMLMEM
18   #define CACHESIZE       5
19 + #else
20 + #define CACHESIZE       17      /* default cache size (Mbytes, 0==inf) */
21   #endif
22   #endif
23   #ifndef FREEBEAMS
# Line 52 | Line 56 | static struct fraglist {
56          short   writable;       /* 0 read-only, <0 write error encountered */
57          int     nfrags;         /* number of known fragments */
58          BEAMI   *fi;            /* fragments, descending file position */
59 <        long    flen;           /* last known file length */
59 >        off_t   flen;           /* last known file length */
60   } *hdfragl;             /* fragment lists, indexed by file descriptor */
61  
62   static int      nhdfragls;      /* size of hdfragl array */
# Line 66 | Line 70 | HDGRID *hproto;
70          register HOLO   *hp;
71          int     n;
72                                  /* copy grid to temporary header */
73 <        bcopy((char *)hproto, (char *)&hdhead, sizeof(HDGRID));
73 >        memcpy((void *)&hdhead, (void *)hproto, sizeof(HDGRID));
74                                  /* compute grid vectors and sizes */
75          hdcompgrid(&hdhead);
76                                  /* allocate header with directory */
# Line 74 | Line 78 | HDGRID *hproto;
78          if ((hp = (HOLO *)malloc(n)) == NULL)
79                  return(NULL);
80                                  /* copy header information */
81 <        copystruct(hp, &hdhead);
81 >        *hp = hdhead;
82                                  /* allocate and clear beam list */
83          hp->bl = (BEAM **)malloc((nbeams(hp)+1)*sizeof(BEAM *)+sizeof(BEAM));
84          if (hp->bl == NULL) {
85                  free((void *)hp);
86                  return(NULL);
87          }
88 <        bzero((char *)hp->bl, (nbeams(hp)+1)*sizeof(BEAM *)+sizeof(BEAM));
88 >        memset((void *)hp->bl, '\0', (nbeams(hp)+1)*sizeof(BEAM *)+sizeof(BEAM));
89          hp->bl[0] = (BEAM *)(hp->bl+nbeams(hp)+1);      /* set blglob(hp) */
90          hp->fd = -1;
91          hp->dirty = 0;
92          hp->priv = NULL;
93                                  /* clear beam directory */
94 <        bzero((char *)hp->bi, (nbeams(hp)+1)*sizeof(BEAMI));
94 >        memset((void *)hp->bi, '\0', (nbeams(hp)+1)*sizeof(BEAMI));
95          return(hp);             /* all is well */
96   }
97  
# Line 101 | Line 105 | char   *rout;
105          register char   *newp;
106                                          /* call malloc/realloc */
107          if (ptr == NULL) newp = (char *)malloc(siz);
108 <        else newp = (char *)realloc(ptr, siz);
108 >        else newp = (char *)realloc((void *)ptr, siz);
109                                          /* check success */
110          if (newp == NULL && rout != NULL) {
111                  hdfreecache(25, NULL);  /* free some memory */
# Line 123 | Line 127 | int    wr;
127          if (fd >= nhdfragls) {
128                  hdfragl = (struct fraglist *)hdrealloc((char *)hdfragl,
129                                  (fd+1)*sizeof(struct fraglist), "hdattach");
130 <                bzero((char *)(hdfragl+nhdfragls),
131 <                                (fd+1-nhdfragls)*sizeof(struct fraglist));
130 >                memset((void *)(hdfragl+nhdfragls),
131 >                                '\0', (fd+1-nhdfragls)*sizeof(struct fraglist));
132                  nhdfragls = fd+1;
133          }
134          hdfragl[fd].nlinks++;
135          hdfragl[fd].writable = wr;              /* set writable flag */
136 <        hdfragl[fd].flen = lseek(fd, (off_t)0L, 2);     /* get file length */
136 >        hdfragl[fd].flen = lseek(fd, (off_t)0, 2);      /* get file length */
137   }
138  
139  
# Line 139 | Line 143 | int    wr;
143   hdrelease(fd)           /* stop tracking file fragments for some section */
144   register int    fd;
145   {
146 <        if (fd < 0 | fd >= nhdfragls || !hdfragl[fd].nlinks)
146 >        if ((fd < 0) | (fd >= nhdfragls) || !hdfragl[fd].nlinks)
147                  return;
148          if (!--hdfragl[fd].nlinks && hdfragl[fd].nfrags) {
149                  free((void *)hdfragl[fd].fi);
# Line 154 | Line 158 | hdinit(fd, hproto)     /* initialize a holodeck section in
158   int     fd;                     /* corresponding file descriptor */
159   HDGRID  *hproto;                /* holodeck section grid */
160   {
161 <        long    rtrunc;
162 <        long    fpos;
161 >        off_t   rtrunc;
162 >        off_t   fpos;
163          int     writable;
164          register HOLO   *hp;
165          register int    n;
166                                          /* prepare for system errors */
167          errno = 0;
168 <        if ((fpos = lseek(fd, (off_t)0L, 1)) < 0)
168 >        if ((fpos = lseek(fd, (off_t)0, 1)) < 0)
169                  error(SYSTEM, "cannot determine holodeck file position");
170          if (hproto == NULL) {           /* assume we're loading it */
171                  HDGRID  hpr;
# Line 186 | Line 190 | HDGRID *hproto;                /* holodeck section grid */
190                  if (fd < nhdfragls && hdfragl[fd].nlinks)
191                          writable = hdfragl[fd].writable;
192                  else
193 <                        writable = lseek(fd, (off_t)fpos, 0) == fpos &&
193 >                        writable = lseek(fd, fpos, 0) == fpos &&
194                                  write(fd, (char *)hp, sizeof(HDGRID)) ==
195                                                          sizeof(HDGRID);
196          } else {                        /* else assume we're creating it */
# Line 209 | Line 213 | HDGRID *hproto;                /* holodeck section grid */
213          fpos = hdfilen(fd);
214          biglob(hp)->nrd = rtrunc = 0;
215          for (n = hproto == NULL ? nbeams(hp) : 0; n > 0; n--)
216 <                if (hp->bi[n].nrd)
216 >                if (hp->bi[n].nrd) {
217                          if (hp->bi[n].fo+hp->bi[n].nrd*sizeof(RAYVAL) > fpos) {
218                                  rtrunc += hp->bi[n].nrd;
219                                  hp->bi[n].nrd = 0;
220                          } else
221                                  biglob(hp)->nrd += hp->bi[n].nrd;
222 +                }
223          if (rtrunc) {
224                  sprintf(errmsg, "truncated section, %ld rays lost (%.1f%%)",
225                                  rtrunc, 100.*rtrunc/(rtrunc+biglob(hp)->nrd));
# Line 242 | Line 247 | int    i;
247          register int    j;
248  
249          if (!hp->dirty++) {                     /* write smudge first time */
250 <                if (lseek(hp->fd, (off_t)(biglob(hp)->fo+(i-1)*sizeof(BEAMI)), 0) < 0
250 >                if (lseek(hp->fd, biglob(hp)->fo+(i-1)*sizeof(BEAMI), 0) < 0
251                                  || write(hp->fd, (char *)&smudge,
252                                          sizeof(BEAMI)) != sizeof(BEAMI))
253                          error(SYSTEM, "seek/write error in hdmarkdirty");
# Line 257 | Line 262 | int    i;
262                          hp->dirseg[j].n = 1;
263                          break;
264                  }
265 <                copystruct(hp->dirseg+j, hp->dirseg+(j-1));
265 >                *(hp->dirseg+j) = *(hp->dirseg+(j-1));
266          }
267          do {                            /* check neighbors */
268                  mindist = nbeams(hp);           /* find closest */
# Line 279 | Line 284 | int    i;
284                                          hp->dirseg[minpos].n - hp->dirseg[j].s;
285                  hp->dirty--;
286                  while (++j < hp->dirty)         /* close the gap */
287 <                        copystruct(hp->dirseg+j, hp->dirseg+(j+1));
287 >                        *(hp->dirseg+j) = *(hp->dirseg+(j+1));
288          } while (mindist <= MINDIRSEL);
289   }
290  
# Line 305 | Line 310 | int    all;
310                  return(0);
311          errno = 0;                      /* write dirty segments */
312          for (j = 0; j < hp->dirty; j++) {
313 <                if (lseek(hp->fd, (off_t)(biglob(hp)->fo +
314 <                                (hp->dirseg[j].s-1)*sizeof(BEAMI)), 0) < 0)
313 >                if (lseek(hp->fd, biglob(hp)->fo +
314 >                                (hp->dirseg[j].s-1)*sizeof(BEAMI), 0) < 0)
315                          error(SYSTEM, "cannot seek on holodeck file");
316                  n = hp->dirseg[j].n * sizeof(BEAMI);
317                  if (write(hp->fd, (char *)(hp->bi+hp->dirseg[j].s), n) != n)
# Line 346 | Line 351 | int    all;                    /* include overhead (painful) */
351   }
352  
353  
354 < long
354 > off_t
355   hdfilen(fd)             /* return file length for fd */
356   int     fd;
357   {
358 <        long    fpos, flen;
358 >        off_t   fpos, flen;
359  
360          if (fd < 0)
361                  return(-1);
362          if (fd >= nhdfragls || !hdfragl[fd].nlinks) {
363 <                if ((fpos = lseek(fd, (off_t)0L, 1)) < 0)
363 >                if ((fpos = lseek(fd, (off_t)0, 1)) < 0)
364                          return(-1);
365 <                flen = lseek(fd, (off_t)0L, 2);
366 <                lseek(fd, (off_t)fpos, 0);
365 >                flen = lseek(fd, (off_t)0, 2);
366 >                lseek(fd, fpos, 0);
367                  return(flen);
368          }
369          return(hdfragl[fd].flen);
370   }
371  
372  
373 < long
373 > off_t
374   hdfiluse(fd, all)       /* compute file usage (in bytes) */
375   int     fd;                     /* open file descriptor to check */
376   int     all;                    /* include overhead and unflushed data */
377   {
378 <        long    total = 0;
378 >        off_t   total = 0;
379          register int    i, j;
380  
381          for (j = 0; hdlist[j] != NULL; j++) {
# Line 401 | Line 406 | int    nr;                     /* number of new rays desired */
406          int     n;
407  
408          if (nr <= 0) return(NULL);
409 <        CHECK(i < 1 | i > nbeams(hp),
409 >        CHECK((i < 1) | (i > nbeams(hp)),
410                          CONSISTENCY, "bad beam index given to hdnewrays");
411          if (hp->bl[i] != NULL)
412                  hp->bl[i]->tick = hdclock;      /* preempt swap */
# Line 413 | Line 418 | int    nr;                     /* number of new rays desired */
418                  blglob(hp)->nrm += n;
419                  if ((n = hp->bl[i]->nrm = hp->bi[i].nrd)) {
420                          errno = 0;
421 <                        if (lseek(hp->fd, (off_t)hp->bi[i].fo, 0) < 0)
421 >                        if (lseek(hp->fd, hp->bi[i].fo, 0) < 0)
422                                  error(SYSTEM, "seek error on holodeck file");
423                          n *= sizeof(RAYVAL);
424                          if (read(hp->fd, (char *)hdbray(hp->bl[i]), n) != n)
# Line 429 | Line 434 | int    nr;                     /* number of new rays desired */
434                  hdfreefrag(hp, i);              /* relinquish old fragment */
435          p = hdbray(hp->bl[i]) + hp->bl[i]->nrm;
436          hp->bl[i]->nrm += nr;                   /* update in-core structure */
437 <        bzero((char *)p, nr*sizeof(RAYVAL));
437 >        memset((void *)p, '\0', nr*sizeof(RAYVAL));
438          blglob(hp)->tick = hp->bl[i]->tick = hdclock++; /* update LRU clock */
439          return(p);                              /* point to new rays */
440   }
# Line 442 | Line 447 | register int   i;
447   {
448          register int    n;
449  
450 <        CHECK(i < 1 | i > nbeams(hp),
450 >        CHECK((i < 1) | (i > nbeams(hp)),
451                          CONSISTENCY, "bad beam index given to hdgetbeam");
452          if (hp->bl[i] == NULL) {                /* load from disk */
453                  if (!(n = hp->bi[i].nrd))
# Line 452 | Line 457 | register int   i;
457                  hp->bl[i] = (BEAM *)hdrealloc(NULL, hdbsiz(n), "hdgetbeam");
458                  blglob(hp)->nrm += hp->bl[i]->nrm = n;
459                  errno = 0;
460 <                if (lseek(hp->fd, (off_t)hp->bi[i].fo, 0) < 0)
460 >                if (lseek(hp->fd, hp->bi[i].fo, 0) < 0)
461                          error(SYSTEM, "seek error on holodeck file");
462                  n *= sizeof(RAYVAL);
463                  if (read(hp->fd, (char *)hdbray(hp->bl[i]), n) != n)
# Line 469 | Line 474 | int
474   hdfilord(hb1, hb2)      /* order beams for quick loading */
475   register HDBEAMI        *hb1, *hb2;
476   {
477 <        register long   c;
477 >        register off_t  c;
478                                  /* residents go first */
479          if (hb2->h->bl[hb2->b] != NULL)
480                  return(hb1->h->bl[hb1->b] == NULL);
481          if (hb1->h->bl[hb1->b] != NULL)
482                  return(-1);
483                                  /* otherwise sort by file descriptor */
484 <        if ((c = hb1->h->fd - hb2->h->fd))
485 <                return(c);
484 >        if (hb1->h->fd != hb2->h->fd)
485 >                return(hb1->h->fd - hb2->h->fd);
486                                  /* then by position in file */
487          c = hb1->h->bi[hb1->b].fo - hb2->h->bi[hb2->b].fo;
488          return(c > 0 ? 1 : c < 0 ? -1 : 0);
# Line 487 | Line 492 | register HDBEAMI       *hb1, *hb2;
492   hdloadbeams(hb, n, bf)  /* load a list of beams in optimal order */
493   register HDBEAMI        *hb;    /* list gets sorted by hdfilord() */
494   int     n;                      /* list length */
495 < int     (*bf)();                /* callback function (optional) */
495 > void    (*bf)();                /* callback function (optional) */
496   {
497          unsigned        origcachesize, memuse;
498          int     bytesloaded, needbytes, bytes2free;
# Line 496 | Line 501 | int    (*bf)();                /* callback function (optional) */
501                                          /* precheck consistency */
502          if (n <= 0) return;
503          for (i = n; i--; )
504 <                if (hb[i].h==NULL || hb[i].b<1 | hb[i].b>nbeams(hb[i].h))
504 >                if (hb[i].h==NULL || (hb[i].b<1) | (hb[i].b>nbeams(hb[i].h)))
505                          error(CONSISTENCY, "bad beam in hdloadbeams");
506                                          /* sort list for optimal access */
507 <        qsort((char *)hb, n, sizeof(HDBEAMI), hdfilord);
507 >        qsort((void *)hb, n, sizeof(HDBEAMI), hdfilord);
508          bytesloaded = 0;                /* run through loaded beams */
509          for ( ; n && (bp = hb->h->bl[hb->b]) != NULL; n--, hb++) {
510                  bp->tick = hdclock;     /* preempt swap */
# Line 551 | Line 556 | int    i;
556                                  if (++k >= f->nfrags)
557                                          goto endloop;
558                          if (k > j)
559 <                                copystruct(f->fi+j, f->fi+k);
559 >                                *(f->fi+j) = *(f->fi+k);
560                  }
561          endloop:
562                  f->nfrags = j;
# Line 569 | Line 574 | int    i;
574                  if (f->fi == NULL)
575                          newp = (BEAMI *)malloc((j+FRAGBLK)*sizeof(BEAMI));
576                  else
577 <                        newp = (BEAMI *)realloc((char *)f->fi,
577 >                        newp = (BEAMI *)realloc((void *)f->fi,
578                                          (j+FRAGBLK)*sizeof(BEAMI));
579                  if (newp == NULL) {
580                          f->nfrags--;    /* graceful failure */
# Line 583 | Line 588 | int    i;
588                          f->fi[j].nrd = bi->nrd;
589                          break;
590                  }
591 <                copystruct(f->fi+j, f->fi+(j-1));
591 >                *(f->fi+j) = *(f->fi+(j-1));
592          }
593                                          /* coalesce adjacent fragments */
594                                                  /* successors never empty */
# Line 602 | Line 607 | int    i;
607                  }
608          biglob(hp)->nrd -= bi->nrd;             /* tell fragment it's free */
609          bi->nrd = 0;
610 <        bi->fo = 0L;
610 >        bi->fo = 0;
611          hdmarkdirty(hp, i);                     /* assume we'll reallocate */
612          return(1);
613   }
# Line 612 | Line 617 | int
617   hdfragOK(fd, listlen, listsiz)  /* get fragment list status for file */
618   int     fd;
619   int     *listlen;
620 < register int4   *listsiz;
620 > register int32  *listsiz;
621   {
622          register struct fraglist        *f;
623          register int    i;
624  
625 <        if (fd < 0 | fd >= nhdfragls || !(f = &hdfragl[fd])->nlinks)
625 >        if ((fd < 0) | (fd >= nhdfragls) || !(f = &hdfragl[fd])->nlinks)
626                  return(0);              /* listless */
627          if (listlen != NULL)
628                  *listlen = f->nfrags;
# Line 632 | Line 637 | register int4  *listsiz;
637   }
638  
639  
640 < long
640 > off_t
641   hdallocfrag(fd, nrays)          /* allocate a file fragment */
642   int     fd;
643 < unsigned int4   nrays;
643 > uint32  nrays;
644   {
645          register struct fraglist        *f;
646          register int    j;
647 <        long    nfo;
647 >        off_t   nfo;
648  
649          if (nrays == 0)
650                  return(-1L);
# Line 667 | Line 672 | register HOLO  *hp;
672   register int    i;
673   {
674          int     fragfreed;
675 <        unsigned int4   nrays;
675 >        uint32  nrays;
676          unsigned int    n;
677 <        long    nfo;
677 >        off_t   nfo;
678                                          /* check file status */
679          if (hdfragl[hp->fd].writable <= 0)
680                  return(hdfragl[hp->fd].writable);
# Line 683 | Line 688 | register int   i;
688          if (nrays) {                    /* get and write new fragment */
689                  nfo = hdallocfrag(hp->fd, nrays);
690                  errno = 0;
691 <                if (lseek(hp->fd, (off_t)nfo, 0) < 0)
691 >                if (lseek(hp->fd, nfo, 0) < 0)
692                          error(SYSTEM, "cannot seek on holodeck file");
693                  n = hp->bl[i]->nrm * sizeof(RAYVAL);
694                  if (write(hp->fd, (char *)hdbray(hp->bl[i]), n) != n) {
# Line 693 | Line 698 | register int   i;
698                  }
699                  hp->bi[i].fo = nfo;
700          } else
701 <                hp->bi[i].fo = 0L;
701 >                hp->bi[i].fo = 0;
702          biglob(hp)->nrd += nrays - hp->bi[i].nrd;
703          hp->bi[i].nrd = nrays;
704          if (!fragfreed)
# Line 757 | Line 762 | register int   i;
762                  return(nchanged);
763          }
764          if (i == 0) {                   /* clobber entire holodeck */
765 <                if (biglob(hp)->nrd == 0 & blglob(hp)->nrm == 0)
765 >                if ((biglob(hp)->nrd == 0) & (blglob(hp)->nrm == 0))
766                          return(0);              /* already empty */
767                  nchanged = 0;
768                  nchanged = 0;
# Line 781 | Line 786 | register int   i;
786          if (hp->bi[i].nrd && !(hdfragflags&FF_KILL && hdfreefrag(hp,i))) {
787                  biglob(hp)->nrd -= hp->bi[i].nrd;       /* free failed */
788                  hp->bi[i].nrd = 0;
789 <                hp->bi[i].fo = 0L;
789 >                hp->bi[i].fo = 0;
790                  hdmarkdirty(hp, i);
791          }
792          return(nchanged);
# Line 813 | Line 818 | register HOLO  *hp;                    /* section we're adding from */
818                                  hb[j].b = i;
819                                  break;
820                          }
821 <                        copystruct(hb+j, hb+(j-1));
821 >                        *(hb+j) = *(hb+(j-1));
822                  }
823          }
824          return(nents);                  /* return new list length */

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines