ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/modobject.c
Revision: 2.15
Committed: Sun Dec 8 19:25:36 2013 UTC (10 years, 5 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.14: +28 -28 lines
Log Message:
Made eqobjects() a little more general and exposed eqreal() to library

File Contents

# User Rev Content
1 greg 1.1 #ifndef lint
2 greg 2.15 static const char RCSid[] = "$Id: modobject.c,v 2.14 2013/12/08 18:59:53 greg Exp $";
3 greg 1.1 #endif
4     /*
5     * Routines for tracking object modifiers
6 greg 2.7 *
7     * External symbols declared in object.h
8     */
9    
10 greg 2.8 #include "copyright.h"
11 greg 1.1
12     #include "standard.h"
13    
14     #include "object.h"
15    
16 greg 1.2 #include "otypes.h"
17    
18 greg 1.1
19     static struct ohtab {
20     int hsiz; /* current table size */
21     OBJECT *htab; /* table, if allocated */
22     } modtab = {100, NULL}, objtab = {1000, NULL}; /* modifiers and objects */
23    
24 schorsch 2.10 static int otndx(char *, struct ohtab *);
25 greg 1.3
26 greg 1.1
27 greg 2.9 OBJECT
28 greg 2.13 objndx( /* get object number from pointer */
29     OBJREC *op
30     )
31 greg 1.1 {
32 greg 2.13 int i, j;
33 greg 1.1
34 gwlarson 2.6 for (i = nobjects>>OBJBLKSHFT; i >= 0; i--) {
35 greg 2.3 j = op - objblock[i];
36 greg 2.12 if ((j >= 0) & (j < OBJBLKSIZ))
37 gwlarson 2.6 return((i<<OBJBLKSHFT) + j);
38 greg 2.3 }
39     return(OVOID);
40 greg 1.1 }
41    
42    
43 greg 2.9 OBJECT
44 greg 2.13 lastmod( /* find modifier definition before obj */
45     OBJECT obj,
46     char *mname
47     )
48 greg 2.3 {
49 greg 2.13 OBJREC *op;
50     int i;
51 greg 2.3
52 greg 2.4 i = modifier(mname); /* try hash table first */
53 greg 2.12 if ((obj == OVOID) | (i < obj))
54 greg 2.4 return(i);
55     for (i = obj; i-- > 0; ) { /* need to search */
56 greg 2.3 op = objptr(i);
57 greg 2.13 if (ismodifier(op->otype) && op->oname[0] == mname[0] &&
58     !strcmp(op->oname, mname))
59 greg 2.3 return(i);
60     }
61     return(OVOID);
62     }
63    
64    
65 greg 2.9 OBJECT
66 greg 2.13 modifier( /* get a modifier number from its name */
67     char *mname
68     )
69 greg 1.1 {
70 greg 2.13 int ndx;
71 greg 1.1
72     ndx = otndx(mname, &modtab);
73     return(modtab.htab[ndx]);
74     }
75    
76    
77 greg 2.3 #ifdef GETOBJ
78 greg 2.9 OBJECT
79 greg 2.13 object( /* get an object number from its name */
80     char *oname
81     )
82 greg 2.3 {
83 greg 2.13 int ndx;
84 greg 2.3
85     ndx = otndx(oname, &objtab);
86     return(objtab.htab[ndx]);
87     }
88     #endif
89    
90    
91 greg 2.15 int
92 greg 2.14 eqreal( /* are two real values close enough to equal? */
93     double d1,
94     double d2
95     )
96     {
97     if (d2 != 0.0)
98     d1 = d1/d2 - 1.0;
99     return((-FTINY <= d1) & (d1 <= FTINY));
100     }
101    
102    
103     int
104     eqobjects( /* check if two objects are equal */
105     OBJECT obj1,
106     OBJECT obj2
107     )
108     {
109     OBJREC *op1, *op2;
110     int i;
111    
112 greg 2.15 while (obj1 != obj2) {
113     if (obj1 == OVOID)
114     return(0);
115     if (obj2 == OVOID)
116     return(0);
117     op1 = objptr(obj1);
118     op2 = objptr(obj2);
119     if (op1->otype != op2->otype)
120     return(0);
121     if (op1->oargs.nsargs != op2->oargs.nsargs)
122     return(0);
123     if (op1->oargs.nfargs != op2->oargs.nfargs)
124     return(0);
125 greg 2.14 #ifdef IARGS
126 greg 2.15 if (op1->oargs.niargs != op2->oargs.niargs)
127 greg 2.14 return(0);
128 greg 2.15 for (i = op1->oargs.niargs; i-- > 0; )
129     if (op1->oargs.iarg[i] != op2->oargs.iarg[i])
130     return(0);
131 greg 2.14 #endif
132 greg 2.15 for (i = op1->oargs.nfargs; i-- > 0; )
133     if (!eqreal(op1->oargs.farg[i], op2->oargs.farg[i]))
134     return(0);
135     for (i = op1->oargs.nsargs; i-- > 0; )
136     if (strcmp(op1->oargs.sarg[i], op2->oargs.sarg[i]))
137     return(0);
138     obj1 = op1->omod;
139     obj2 = op2->omod;
140     }
141 greg 2.14 return(1);
142     }
143    
144    
145 greg 2.7 void
146 greg 2.13 insertobject( /* insert new object into our list */
147     OBJECT obj
148     )
149 greg 1.1 {
150 greg 2.13 int i;
151 greg 1.1
152     if (ismodifier(objptr(obj)->otype)) {
153     i = otndx(objptr(obj)->oname, &modtab);
154 greg 2.14 if (eqobjects(obj, modtab.htab[i]))
155     return;
156 greg 1.1 modtab.htab[i] = obj;
157     }
158 greg 2.3 #ifdef GETOBJ
159     else {
160     i = otndx(objptr(obj)->oname, &objtab);
161     objtab.htab[i] = obj;
162     }
163     #endif
164 greg 1.1 for (i = 0; addobjnotify[i] != NULL; i++)
165     (*addobjnotify[i])(obj);
166     }
167    
168    
169 greg 2.7 void
170 greg 2.13 clearobjndx(void) /* clear object hash tables */
171 greg 2.7 {
172     if (modtab.htab != NULL) {
173     free((void *)modtab.htab);
174     modtab.htab = NULL;
175     modtab.hsiz = 100;
176     }
177     if (objtab.htab != NULL) {
178     free((void *)objtab.htab);
179     objtab.htab = NULL;
180     objtab.hsiz = 100;
181     }
182     }
183    
184    
185 greg 1.1 static int
186 greg 2.13 nexthsiz( /* return next hash table size */
187     int oldsiz
188     )
189 greg 1.1 {
190     static int hsiztab[] = {
191     251, 509, 1021, 2039, 4093, 8191, 16381, 0
192     };
193 greg 2.13 int *hsp;
194 greg 1.1
195     for (hsp = hsiztab; *hsp; hsp++)
196     if (*hsp > oldsiz)
197     return(*hsp);
198     return(oldsiz*2 + 1); /* not always prime */
199     }
200    
201    
202     static int
203 greg 2.13 otndx( /* get object table index for name */
204     char *name,
205     struct ohtab *tab
206     )
207 greg 1.1 {
208     OBJECT *oldhtab;
209     int hval, i;
210 greg 2.13 int ndx;
211 greg 1.1
212     if (tab->htab == NULL) { /* new table */
213     tab->hsiz = nexthsiz(tab->hsiz);
214     tab->htab = (OBJECT *)malloc(tab->hsiz*sizeof(OBJECT));
215     if (tab->htab == NULL)
216     error(SYSTEM, "out of memory in otndx");
217     ndx = tab->hsiz;
218     while (ndx--) /* empty it */
219     tab->htab[ndx] = OVOID;
220     }
221     /* look up object */
222     hval = shash(name);
223     tryagain:
224     for (i = 0; i < tab->hsiz; i++) {
225 greg 2.11 ndx = (hval + (unsigned long)i*i) % tab->hsiz;
226 greg 1.1 if (tab->htab[ndx] == OVOID ||
227     !strcmp(objptr(tab->htab[ndx])->oname, name))
228     return(ndx);
229     }
230     /* table is full, reallocate */
231     oldhtab = tab->htab;
232     ndx = tab->hsiz;
233     tab->htab = NULL;
234     while (ndx--)
235     if (oldhtab[ndx] != OVOID) {
236     i = otndx(objptr(oldhtab[ndx])->oname, tab);
237     tab->htab[i] = oldhtab[ndx];
238     }
239 greg 2.7 free((void *)oldhtab);
240 greg 1.1 goto tryagain; /* should happen only once! */
241     }