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

Comparing ray/src/util/vwrays.c (file contents):
Revision 3.7 by schorsch, Sun Jul 27 22:12:04 2003 UTC vs.
Revision 3.14 by greg, Mon Jun 22 00:12:04 2009 UTC

# Line 5 | Line 5 | static const char      RCSid[] = "$Id$";
5   * Compute rays corresponding to a given picture or view.
6   */
7  
8 –
9 – #include "standard.h"
10 –
8   #include "platform.h"
9 + #include "standard.h"
10 + #include "random.h"
11   #include "view.h"
12  
13 < extern int      putf(), putd(), puta();
13 > typedef void putfunc(FVECT ro, FVECT rd);
14 > static putfunc puta;
15 > static putfunc putf;
16 > static putfunc putd;
17 > static void pix2rays(FILE *fp);
18 > static void putrays(void);
19  
20 < int     (*putr)() = puta;
20 > static putfunc *putr = puta;
21  
22   VIEW    vw = STDVIEW;
23  
# Line 21 | Line 25 | RESOLU rs = {PIXSTANDARD, 512, 512};
25  
26   double  pa = 1.;
27  
28 + double  pj = 0.;
29 +
30   int     zfd = -1;
31  
32   int     fromstdin = 0;
# Line 28 | Line 34 | int    fromstdin = 0;
34   char    *progname;
35  
36  
37 < main(argc, argv)
38 < int     argc;
39 < char    *argv[];
37 > int
38 > main(
39 >        int     argc,
40 >        char    *argv[]
41 > )
42   {
43          char    *err;
44          int     rval, getdim = 0;
# Line 48 | Line 56 | char   *argv[];
56                                  break;
57                          case 'f':                       /* float */
58                                  putr = putf;
59 +                                SET_FILE_BINARY(stdout);
60                                  break;
61                          case 'd':                       /* double */
62                                  putr = putd;
63 +                                SET_FILE_BINARY(stdout);
64                                  break;
65                          default:
66                                  goto userr;
# Line 91 | Line 101 | char   *argv[];
101                                  exit(1);
102                          }
103                          break;
104 <                case 'p':                       /* pixel aspect ratio */
105 <                        pa = atof(argv[++i]);
104 >                case 'p':                       /* pixel aspect or jitter */
105 >                        if (argv[i][2] == 'a')
106 >                                pa = atof(argv[++i]);
107 >                        else if (argv[i][2] == 'j')
108 >                                pj= atof(argv[++i]);
109 >                        else
110 >                                goto userr;
111                          break;
112                  case 'i':                       /* get pixels from stdin */
113                          fromstdin = 1;
# Line 142 | Line 157 | userr:
157   }
158  
159  
160 < pix2rays(FILE *fp)
160 > static void
161 > jitterloc(
162 >        RREAL   loc[2]
163 > )
164   {
165 +        if (pj > FTINY) {
166 +                loc[0] += pj*(.5 - frandom())/rs.xr;
167 +                loc[1] += pj*(.5 - frandom())/rs.yr;
168 +        }
169 + }
170 +
171 +
172 + static void
173 + pix2rays(
174 +        FILE *fp
175 + )
176 + {
177          static FVECT    rorg, rdir;
178          float   zval;
179          double  px, py;
180 +        RREAL   loc[2];
181          int     pp[2];
182          double  d;
183          register int    i;
184  
185          while (fscanf(fp, "%lf %lf", &px, &py) == 2) {
186 <                if (px < 0 || px >= rs.xr ||
187 <                                py < 0 || py >= rs.yr) {
157 <                        fprintf(stderr,
158 <                                "%s: (x,y) pair (%.0f,%.0f) out of range\n",
159 <                                        progname, px, py);
160 <                        exit(1);
161 <                }
186 >                px += .5; py += .5;
187 >                loc[0] = px/rs.xr; loc[1] = py/rs.yr;
188                  if (zfd >= 0) {
189 <                        loc2pix(pp, &rs, px/rs.xr, py/rs.yr);
189 >                        loc2pix(pp, &rs, loc[0], loc[1]);
190                          if (lseek(zfd,
191 <                                (pp[1]*scanlen(&rs)+pp[0])*sizeof(float), 0)
192 <                                        < 0 ||
191 >                                (pp[1]*scanlen(&rs)+pp[0])*sizeof(float),
192 >                                                SEEK_SET) < 0 ||
193                                          read(zfd, &zval, sizeof(float))
194 <                                        < sizeof(float)) {
194 >                                                < sizeof(float)) {
195                                  fprintf(stderr, "%s: depth buffer read error\n",
196                                                  progname);
197                                  exit(1);
198                          }
199                  }
200 <                d = viewray(rorg, rdir, &vw, px/rs.xr, py/rs.yr);
200 >                jitterloc(loc);
201 >                d = viewray(rorg, rdir, &vw, loc[0], loc[1]);
202                  if (d < -FTINY)
203                          rorg[0] = rorg[1] = rorg[2] =
204                          rdir[0] = rdir[1] = rdir[2] = 0.;
# Line 192 | Line 219 | pix2rays(FILE *fp)
219   }
220  
221  
222 < putrays()
222 > static void
223 > putrays(void)
224   {
225 <        static RREAL    loc[2];
226 <        static FVECT    rorg, rdir;
227 <        float   *zbuf;
225 >        RREAL   loc[2];
226 >        FVECT   rorg, rdir;
227 >        float   *zbuf = NULL;
228          int     sc;
229          double  d;
230          register int    si, i;
# Line 219 | Line 247 | putrays()
247                  }
248                  for (si = 0; si < scanlen(&rs); si++) {
249                          pix2loc(loc, &rs, si, sc);
250 +                        jitterloc(loc);
251                          d = viewray(rorg, rdir, &vw, loc[0], loc[1]);
252                          if (d < -FTINY)
253                                  rorg[0] = rorg[1] = rorg[2] =
# Line 239 | Line 268 | putrays()
268   }
269  
270  
271 < puta(ro, rd)            /* put out ray in ASCII format */
272 < FVECT   ro, rd;
271 > static void
272 > puta(           /* put out ray in ASCII format */
273 >        FVECT   ro,
274 >        FVECT   rd
275 > )
276   {
277          printf("%.5e %.5e %.5e %.5e %.5e %.5e\n",
278                          ro[0], ro[1], ro[2],
# Line 248 | Line 280 | FVECT  ro, rd;
280   }
281  
282  
283 < putf(ro, rd)            /* put out ray in float format */
284 < FVECT   ro, rd;
283 > static void
284 > putf(           /* put out ray in float format */
285 >        FVECT   ro,
286 >        FVECT   rd
287 > )
288   {
289          float v[6];
290  
# Line 259 | Line 294 | FVECT  ro, rd;
294   }
295  
296  
297 < putd(ro, rd)            /* put out ray in double format */
298 < FVECT   ro, rd;
297 > static void
298 > putd(           /* put out ray in double format */
299 >        FVECT   ro,
300 >        FVECT   rd
301 > )
302   {
303          double v[6];
304  

Diff Legend

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