ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/ambio.c
Revision: 2.6
Committed: Fri Apr 11 20:31:37 2014 UTC (10 years ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.5: +101 -14 lines
Log Message:
Partial implementation of Hessian gradient calculation (-DNEWAMB)

File Contents

# Content
1 #ifndef lint
2 static const char RCSid[] = "$Id: ambio.c,v 2.5 2003/02/25 02:47:22 greg Exp $";
3 #endif
4 /*
5 * Read and write portable ambient values
6 *
7 * Declarations of external symbols in ambient.h
8 */
9
10 #include "copyright.h"
11
12 #include "ray.h"
13 #include "ambient.h"
14
15
16 #define badflt(x) ((x) < -FHUGE || (x) > FHUGE)
17
18 #define badvec(v) (badflt((v)[0]) | badflt((v)[1]) | badflt((v)[2]))
19
20
21 void
22 putambmagic(fp) /* write out ambient value magic number */
23 FILE *fp;
24 {
25 putint(AMBMAGIC, 2, fp);
26 }
27
28
29 int
30 hasambmagic(fp) /* read in and check validity of magic # */
31 FILE *fp;
32 {
33 int magic;
34
35 magic = getint(2, fp);
36 if (feof(fp))
37 return(0);
38 return(magic == AMBMAGIC);
39 }
40
41
42 #ifdef NEWAMB
43
44 #define putpos(v,fp) putflt((v)[0],fp);putflt((v)[1],fp);putflt((v)[2],fp)
45
46 #define getpos(v,fp) (v)[0]=getflt(fp);(v)[1]=getflt(fp);(v)[2]=getflt(fp)
47
48 #define putv2(v2,fp) putflt((v2)[0],fp);putflt((v2)[1],fp)
49
50 #define getv2(v2,fp) (v2)[0]=getflt(fp);(v2)[1]=getflt(fp)
51
52 int
53 writambval( /* write ambient value to stream */
54 AMBVAL *av,
55 FILE *fp
56 )
57 {
58 COLR col;
59
60 putint(av->lvl, 1, fp);
61 putflt(av->weight, fp);
62 putpos(av->pos, fp);
63 putint(av->ndir, sizeof(av->ndir), fp);
64 putint(av->udir, sizeof(av->udir), fp);
65 setcolr(col, colval(av->val,RED),
66 colval(av->val,GRN), colval(av->val,BLU));
67 fwrite((char *)col, sizeof(col), 1, fp);
68 putv2(av->rad, fp);
69 putv2(av->gpos, fp);
70 putv2(av->gdir, fp);
71 return(ferror(fp) ? -1 : 0);
72 }
73
74
75 int
76 ambvalOK( /* check consistency of ambient value */
77 AMBVAL *av
78 )
79 {
80 double d;
81
82 if (badvec(av->pos)) return(0);
83 if (!av->ndir | !av->udir) return(0);
84 if (av->lvl < 0 || av->lvl > 100) return(0);
85 if (av->weight <= 0. || av->weight > 1.) return(0);
86 if ((av->rad[0] <= 0.) | (av->rad[0] >= FHUGE)) return(0);
87 if ((av->rad[1] <= 0.) | (av->rad[1] >= FHUGE)) return(0);
88 if (colval(av->val,RED) < 0. ||
89 colval(av->val,RED) > FHUGE ||
90 colval(av->val,GRN) < 0. ||
91 colval(av->val,GRN) > FHUGE ||
92 colval(av->val,BLU) < 0. ||
93 colval(av->val,BLU) > FHUGE) return(0);
94 if (badflt(av->gpos[0]) || badflt(av->gpos[1])) return(0);
95 if (badflt(av->gdir[0]) || badflt(av->gdir[1])) return(0);
96 return(1);
97 }
98
99
100 int
101 readambval( /* read ambient value from stream */
102 AMBVAL *av,
103 FILE *fp
104 )
105 {
106 COLR col;
107
108 av->lvl = getint(1, fp);
109 if (feof(fp))
110 return(0);
111 av->weight = getflt(fp);
112 getpos(av->pos, fp);
113 av->ndir = getint(sizeof(av->ndir), fp);
114 av->udir = getint(sizeof(av->udir), fp);
115 if (fread((char *)col, sizeof(col), 1, fp) != 1)
116 return(0);
117 colr_color(av->val, col);
118 getv2(av->rad, fp);
119 getv2(av->gpos, fp);
120 getv2(av->gdir, fp);
121 return(feof(fp) ? 0 : ambvalOK(av));
122 }
123
124
125 #else /* ! NEWAMB */
126
127
128 #define putvec(v,fp) putflt((v)[0],fp);putflt((v)[1],fp);putflt((v)[2],fp)
129
130 #define getvec(v,fp) (v)[0]=getflt(fp);(v)[1]=getflt(fp);(v)[2]=getflt(fp)
131
132
133 int
134 writambval(av, fp) /* write ambient value to stream */
135 AMBVAL *av;
136 FILE *fp;
137 {
138 COLR col;
139
140 putint(av->lvl, 1, fp);
141 putflt(av->weight, fp);
142 putvec(av->pos, fp);
143 putvec(av->dir, fp);
144 setcolr(col, colval(av->val,RED),
145 colval(av->val,GRN), colval(av->val,BLU));
146 fwrite((char *)col, sizeof(col), 1, fp);
147 putflt(av->rad, fp);
148 putvec(av->gpos, fp);
149 putvec(av->gdir, fp);
150 return(ferror(fp) ? -1 : 0);
151 }
152
153
154 int
155 ambvalOK(av) /* check consistency of ambient value */
156 AMBVAL *av;
157 {
158 double d;
159
160 if (badvec(av->pos)) return(0);
161 if (badvec(av->dir)) return(0);
162 d = DOT(av->dir,av->dir);
163 if (d < 0.9999 || d > 1.0001) return(0);
164 if (av->lvl < 0 || av->lvl > 100) return(0);
165 if (av->weight <= 0. || av->weight > 1.) return(0);
166 if (av->rad <= 0. || av->rad >= FHUGE) return(0);
167 if (colval(av->val,RED) < 0. ||
168 colval(av->val,RED) > FHUGE ||
169 colval(av->val,GRN) < 0. ||
170 colval(av->val,GRN) > FHUGE ||
171 colval(av->val,BLU) < 0. ||
172 colval(av->val,BLU) > FHUGE) return(0);
173 if (badvec(av->gpos)) return(0);
174 if (badvec(av->gdir)) return(0);
175 return(1);
176 }
177
178
179 int
180 readambval(av, fp) /* read ambient value from stream */
181 AMBVAL *av;
182 FILE *fp;
183 {
184 COLR col;
185
186 av->lvl = getint(1, fp);
187 if (feof(fp))
188 return(0);
189 av->weight = getflt(fp);
190 getvec(av->pos, fp);
191 getvec(av->dir, fp);
192 if (fread((char *)col, sizeof(col), 1, fp) != 1)
193 return(0);
194 colr_color(av->val, col);
195 av->rad = getflt(fp);
196 getvec(av->gpos, fp);
197 getvec(av->gdir, fp);
198 return(feof(fp) ? 0 : ambvalOK(av));
199 }
200
201 #endif /* ! NEWAMB */