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

Comparing ray/src/common/objset.c (file contents):
Revision 1.1 by greg, Thu Feb 2 10:34:37 1989 UTC vs.
Revision 2.9 by greg, Tue Feb 25 02:47:21 2003 UTC

# Line 1 | Line 1
1 /* Copyright (c) 1986 Regents of the University of California */
2
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6
4   /*
5   *  objset.c - routines for maintaining object sets.
6   *
7 < *      7/28/85
7 > *  External symbols declared in object.h
8   */
9  
10 + #include "copyright.h"
11 +
12   #include  "standard.h"
13  
14   #include  "octree.h"
15  
16   #include  "object.h"
17  
18 < #define  OSTSIZ         3037            /* object table size (a prime!) */
18 > #ifndef  OSTSIZ
19 > #ifdef  BIGMEM
20 > #define  OSTSIZ         262139          /* object table size (a prime!) */
21 > #else
22 > #define  OSTSIZ         32749           /* object table size (a prime!) */
23 > #endif
24 > #endif
25  
26   static OBJECT  *ostable[OSTSIZ];        /* the object set table */
27  
28  
29 + void
30   insertelem(os, obj)             /* insert obj into os, no questions */
31   register OBJECT  *os;
32   OBJECT  obj;
# Line 36 | Line 42 | OBJECT  obj;
42   }
43  
44  
45 + void
46   deletelem(os, obj)              /* delete obj from os, no questions */
47   register OBJECT  *os;
48   OBJECT  obj;
49   {
50          register int  i;
51  
52 <        for (i = (*os++)--; i > 0 && *os < obj; i--, os++)
53 <                ;
52 >        i = (*os)--;
53 >        os++;
54 >        while (i > 0 && *os < obj) {
55 >                i--;
56 >                os++;
57 >        }
58          while (--i > 0) {
59                  os[0] = os[1];
60                  os++;
# Line 51 | Line 62 | OBJECT  obj;
62   }
63  
64  
65 + int
66   inset(os, obj)                  /* determine if object is in set */
67   register OBJECT  *os;
68   OBJECT  obj;
# Line 58 | Line 70 | OBJECT  obj;
70          int  upper, lower;
71          register int  cm, i;
72  
73 +        if ((i = os[0]) <= 6) {         /* linear search algorithm */
74 +                cm = obj;
75 +                while (i-- > 0)
76 +                        if (*++os == cm)
77 +                                return(1);
78 +                return(0);
79 +        }
80          lower = 1;
81 <        upper = cm = os[0] + 1;
82 <
81 >        upper = cm = i + 1;
82 >                                        /* binary search algorithm */
83          while ((i = (lower + upper) >> 1) != cm) {
84                  cm = obj - os[i];
85                  if (cm > 0)
# Line 75 | Line 94 | OBJECT  obj;
94   }
95  
96  
97 + int
98   setequal(os1, os2)              /* determine if two sets are equal */
99   register OBJECT  *os1, *os2;
100   {
# Line 87 | Line 107 | register OBJECT  *os1, *os2;
107   }
108  
109  
110 + void
111   setcopy(os1, os2)               /* copy object set os2 into os1 */
112   register OBJECT  *os1, *os2;
113   {
# Line 97 | Line 118 | register OBJECT  *os1, *os2;
118   }
119  
120  
121 + OBJECT *
122 + setsave(os)                     /* allocate space and save set */
123 + register OBJECT  *os;
124 + {
125 +        OBJECT  *osnew;
126 +        register OBJECT  *oset;
127 +        register int  i;
128 +
129 +        if ((osnew = oset = (OBJECT *)malloc((*os+1)*sizeof(OBJECT))) == NULL)
130 +                error(SYSTEM, "out of memory in setsave\n");
131 +        for (i = *os; i-- >= 0; )       /* inline setcopy */
132 +                *oset++ = *os++;
133 +        return(osnew);
134 + }
135 +
136 +
137 + void
138 + setunion(osr, os1, os2)         /* osr = os1 Union os2 */
139 + register OBJECT  *osr, *os1, *os2;
140 + {
141 +        register int    i1, i2;
142 +
143 +        osr[0] = 0;
144 +        for (i1 = i2 = 1; i1 <= os1[0] || i2 <= os2[0]; ) {
145 +                while (i1 <= os1[0] && (i2 > os2[0] || os1[i1] <= os2[i2])) {
146 +                        osr[++osr[0]] = os1[i1];
147 +                        if (i2 <= os2[0] && os2[i2] == os1[i1])
148 +                                i2++;
149 +                        i1++;
150 +                }
151 +                while (i2 <= os2[0] && (i1 > os1[0] || os2[i2] < os1[i1]))
152 +                        osr[++osr[0]] = os2[i2++];
153 +        }
154 + }
155 +
156 +
157 + void
158 + setintersect(osr, os1, os2)     /* osr = os1 Intersect os2 */
159 + register OBJECT  *osr, *os1, *os2;
160 + {
161 +        register int    i1, i2;
162 +
163 +        osr[0] = 0;
164 +        if (os1[0] <= 0)
165 +                return;
166 +        for (i1 = i2 = 1; i2 <= os2[0]; ) {
167 +                while (os1[i1] < os2[i2])
168 +                        if (++i1 > os1[0])
169 +                                return;
170 +                while (os2[i2] < os1[i1])
171 +                        if (++i2 > os2[0])
172 +                                return;
173 +                if (os1[i1] == os2[i2])
174 +                        osr[++osr[0]] = os2[i2++];
175 +        }
176 + }
177 +
178 +
179   OCTREE
180   fullnode(oset)                  /* return octree for object set */
181   OBJECT  *oset;
# Line 153 | Line 232 | memerr:
232   }
233  
234  
235 + void
236   objset(oset, ot)                /* get object set for full node */
237   register OBJECT  *oset;
238   OCTREE  ot;
# Line 162 | Line 242 | OCTREE  ot;
242  
243          if (!isfull(ot))
244                  goto noderr;
245 <        i = oseti(ot);
246 <        if ((os = ostable[i%OSTSIZ]) == NULL)
245 >        ot = oseti(ot);
246 >        if ((os = ostable[ot%OSTSIZ]) == NULL)
247                  goto noderr;
248 <        for (i /= OSTSIZ; i--; os += *os + 1)
248 >        for (i = ot/OSTSIZ; i--; os += *os + 1)
249                  if (*os <= 0)
250                          goto noderr;
251          for (i = *os; i-- >= 0; )               /* copy set here */
# Line 173 | Line 253 | OCTREE  ot;
253          return;
254   noderr:
255          error(CONSISTENCY, "bad node in objset");
256 + }
257 +
258 +
259 + int
260 + dosets(f)                               /* loop through all sets */
261 + int     (*f)();
262 + {
263 +        int  res = 0;
264 +        int  n;
265 +        register OBJECT  *os;
266 +
267 +        for (n = 0; n < OSTSIZ; n++) {
268 +                if ((os = ostable[n]) == NULL)
269 +                        continue;
270 +                while (*os > 0) {
271 +                        res += (*f)(os);
272 +                        os += *os + 1;
273 +                }
274 +        }
275 +        return(res);
276 + }
277 +
278 +
279 + void
280 + donesets()                      /* free ALL SETS in our table */
281 + {
282 +        register int  n;
283 +
284 +        for (n = 0; n < OSTSIZ; n++)
285 +                if (ostable[n] != NULL) {
286 +                        free((void *)ostable[n]);
287 +                        ostable[n] = NULL;
288 +                }
289   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines