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

Comparing ray/src/hd/rhoptimize.c (file contents):
Revision 3.5 by gwlarson, Thu Nov 12 15:07:49 1998 UTC vs.
Revision 3.11 by greg, Fri Jun 13 15:27:04 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1998 Silicon Graphics, Inc. */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ SGI";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   * Optimize holodeck for quick access.
6   *
# Line 20 | Line 17 | static char SCCSid[] = "$SunId$ SGI";
17  
18   char    *progname;
19   char    tempfile[128];
20 + int     dupchecking = 0;
21  
22   extern char     *rindex();
25 – extern int      quit();
23   extern long     rhinitcopy();
24  
25  
# Line 35 | Line 32 | char   *argv[];
32          long    nextipos, lastopos, thisopos;
33  
34          progname = argv[0];
35 <        if (argc < 2 | argc > 3) {
36 <                fprintf(stderr, "Usage: %s input.hdk [output.hdk]\n", progname);
35 >        argv++; argc--;                 /* duplicate checking flag? */
36 >        if (argc > 1 && !strcmp(argv[0], "-u")) {
37 >                dupchecking++;
38 >                argv++; argc--;
39 >        }
40 >        if (argc < 1 | argc > 2) {
41 >                fprintf(stderr, "Usage: %s [-u] input.hdk [output.hdk]\n",
42 >                                progname);
43                  exit(1);
44          }
45 <        inpname = argv[1];
46 <        if (argc == 3)                  /* use given output file */
47 <                outname = argv[2];
45 >        inpname = argv[0];              /* get input file */
46 >        argv++; argc--;
47 >        if (argc == 1)                  /* use given output file */
48 >                outname = argv[0];
49          else {                          /* else use temporary file */
50                  if (access(inpname, R_OK|W_OK) < 0) {   /* check permissions */
51                          sprintf(errmsg, "cannot access \"%s\"", inpname);
# Line 60 | Line 64 | char   *argv[];
64          lastopos = 0L;                  /* copy sections one by one */
65          while (nextipos != 0L) {
66                                          /* set input position; get next */
67 <                lseek(hdfd[0], nextipos, 0);
67 >                lseek(hdfd[0], (off_t)nextipos, 0);
68                  read(hdfd[0], (char *)&nextipos, sizeof(nextipos));
69                                          /* get output position; set last */
70 <                thisopos = lseek(hdfd[1], 0L, 2);
70 >                thisopos = lseek(hdfd[1], (off_t)0, 2);
71                  if (lastopos > 0L) {
72 <                        lseek(hdfd[1], lastopos, 0);
72 >                        lseek(hdfd[1], (off_t)lastopos, 0);
73                          write(hdfd[1], (char *)&thisopos, sizeof(thisopos));
74 <                        lseek(hdfd[1], 0L, 2);
74 >                        lseek(hdfd[1], (off_t)0, 2);
75                  }
76                  lastopos = thisopos;
77                  thisopos = 0L;          /* write place holder */
# Line 99 | Line 103 | char   *infn, *outfn;
103                  sprintf(errmsg, "cannot open \"%s\" for reading", infn);
104                  error(SYSTEM, errmsg);
105          }
106 +        if (access(outfn, F_OK) == 0) {
107 +                sprintf(errmsg, "output file \"%s\" already exists!", outfn);
108 +                error(USER, errmsg);
109 +        }
110          if ((outfp = fopen(outfn, "w+")) == NULL) {
111                  sprintf(errmsg, "cannot open \"%s\" for writing", outfn);
112                  error(SYSTEM, errmsg);
# Line 133 | Line 141 | char   *infn, *outfn;
141   }
142  
143  
144 + int
145 + nuniq(rva, n)                   /* sort unique rays to front of beam list */
146 + register RAYVAL *rva;
147 + int     n;
148 + {
149 +        register int    i, j;
150 +        RAYVAL  rtmp;
151 +
152 +        for (j = 0; j < n; j++)
153 +                for (i = j+1; i < n; i++)
154 +                        if ( rva[i].d == rva[j].d &&
155 +                                        rva[i].r[0][0]==rva[j].r[0][0] &&
156 +                                        rva[i].r[0][1]==rva[j].r[0][1] &&
157 +                                        rva[i].r[1][0]==rva[j].r[1][0] &&
158 +                                        rva[i].r[1][1]==rva[j].r[1][1] ) {
159 +                                n--;            /* swap duplicate with end */
160 +                                copystruct(&rtmp, rva+n);
161 +                                copystruct(rva+n, rva+i);
162 +                                copystruct(rva+i, &rtmp);
163 +                                i--;            /* recheck one we swapped */
164 +                        }
165 +        return(n);
166 + }
167 +
168 +
169   static BEAMI    *beamdir;
170  
171   static int
172   bpcmp(b1p, b2p)                 /* compare beam positions on disk */
173   int     *b1p, *b2p;
174   {
175 <        register long   pdif = beamdir[*b1p].fo - beamdir[*b2p].fo;
175 >        register off_t  pdif = beamdir[*b1p].fo - beamdir[*b2p].fo;
176  
144 –        if (pdif > 0L) return(1);
177          if (pdif < 0L) return(-1);
178 <        return(0);
178 >        return(pdif > 0L);
179   }
180  
181   static HOLO     *hout;
# Line 155 | Line 187 | int    *bq, nb;
187   {
188          register int    i;
189          register BEAM   *bp;
190 +        int     n;
191  
192          beamdir = hp->bi;               /* sort based on file position */
193 <        qsort((char *)bq, nb, sizeof(*bq), bpcmp);
193 >        qsort((void *)bq, nb, sizeof(*bq), bpcmp);
194                                          /* transfer and free each beam */
195          for (i = 0; i < nb; i++) {
196                  bp = hdgetbeam(hp, bq[i]);
197 <                bcopy((char *)hdbray(bp), (char *)hdnewrays(hout,bq[i],bp->nrm),
198 <                                bp->nrm*sizeof(RAYVAL));
197 >                DCHECK(bp==NULL, CONSISTENCY, "empty beam in xferclump");
198 >                n = dupchecking ? nuniq(hdbray(bp),bp->nrm) : bp->nrm;
199 >                bcopy((void *)hdbray(bp), (void *)hdnewrays(hout,bq[i],n),
200 >                                n*sizeof(RAYVAL));
201                  hdfreebeam(hp, bq[i]);
202          }
203 <        hdflush(hout);                  /* write & free clump */
203 >        hdfreebeam(hout, 0);            /* write & free clump */
204          return(0);
205   }
206  
# Line 185 | Line 220 | int    ifd, ofd;
220   }
221  
222  
223 + void
224   eputs(s)                        /* put error message to stderr */
225   register char  *s;
226   {
# Line 204 | Line 240 | register char  *s;
240   }
241  
242  
243 + void
244   quit(code)                      /* exit the program gracefully */
245   int     code;
246   {

Diff Legend

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