ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/rv2.c
Revision: 2.62
Committed: Wed Oct 5 17:20:55 2011 UTC (12 years, 6 months ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad4R1
Changes since 2.61: +11 -19 lines
Log Message:
Windows compatibility changes from Bill Hoffman

File Contents

# User Rev Content
1 greg 1.1 #ifndef lint
2 greg 2.62 static const char RCSid[] = "$Id: rv2.c,v 2.61 2010/10/20 16:36:10 greg Exp $";
3 greg 1.1 #endif
4     /*
5     * rv2.c - command routines used in tracing a view.
6     *
7 greg 2.38 * External symbols declared in rpaint.h
8     */
9    
10 schorsch 2.41 #include "copyright.h"
11    
12 schorsch 2.40 #include <ctype.h>
13 schorsch 2.41 #include <string.h>
14 greg 1.1
15 schorsch 2.40 #include "platform.h"
16 greg 2.62 #include "rtprocess.h" /* win_popen() */
17     #include "paths.h"
18 greg 1.1 #include "ray.h"
19 greg 2.46 #include "source.h"
20 schorsch 2.43 #include "ambient.h"
21 greg 1.1 #include "otypes.h"
22     #include "rpaint.h"
23    
24 greg 2.38 extern int psample; /* pixel sample size */
25     extern double maxdiff; /* max. sample difference */
26 greg 1.1
27 greg 2.4 #define CTRL(c) ((c)-'@')
28 greg 1.1
29 greg 2.6 #ifdef SMLFLT
30     #define sscanvec(s,v) (sscanf(s,"%f %f %f",v,v+1,v+2)==3)
31     #else
32     #define sscanvec(s,v) (sscanf(s,"%lf %lf %lf",v,v+1,v+2)==3)
33     #endif
34    
35 greg 2.20 extern char rifname[128]; /* rad input file name */
36 greg 2.18
37 greg 1.1 extern char *progname;
38 greg 1.29 extern char *octname;
39 greg 1.1
40    
41 greg 2.55 void
42 schorsch 2.43 getframe( /* get a new frame */
43     char *s
44     )
45 greg 1.1 {
46 greg 1.12 if (getrect(s, &pframe) < 0)
47 greg 1.1 return;
48     pdepth = 0;
49     }
50    
51    
52 greg 2.55 void
53 schorsch 2.43 getrepaint( /* get area and repaint */
54     char *s
55     )
56 greg 1.12 {
57     RECT box;
58    
59     if (getrect(s, &box) < 0)
60     return;
61 greg 2.55 paintrect(&ptrunk, &box);
62 greg 1.12 }
63    
64    
65 greg 2.55 void
66 greg 2.62 getview( /* get/show/save view parameters */
67 schorsch 2.43 char *s
68     )
69 greg 1.1 {
70     FILE *fp;
71     char buf[128];
72 greg 1.10 char *fname;
73 greg 1.3 int change = 0;
74 greg 2.44 VIEW nv = ourview;
75 greg 1.1
76 greg 2.9 while (isspace(*s))
77     s++;
78     if (*s == '-') { /* command line parameters */
79     if (sscanview(&nv, s))
80     newview(&nv);
81     else
82     error(COMMAND, "bad view option(s)");
83     return;
84     }
85 greg 2.62 if (nextword(buf, sizeof(buf), s) != NULL) { /* write to a file */
86 greg 1.10 if ((fname = getpath(buf, NULL, 0)) == NULL ||
87 greg 1.1 (fp = fopen(fname, "a")) == NULL) {
88     sprintf(errmsg, "cannot open \"%s\"", buf);
89     error(COMMAND, errmsg);
90     return;
91     }
92 greg 1.2 fputs(progname, fp);
93 greg 1.1 fprintview(&ourview, fp);
94 greg 1.19 fputs(sskip(s), fp);
95 greg 2.18 putc('\n', fp);
96 greg 1.1 fclose(fp);
97     return;
98     }
99     sprintf(buf, "view type (%c): ", ourview.type);
100     (*dev->comout)(buf);
101 greg 1.18 (*dev->comin)(buf, NULL);
102 greg 2.4 if (buf[0] == CTRL('C')) return;
103 greg 1.3 if (buf[0] && buf[0] != ourview.type) {
104 greg 1.1 nv.type = buf[0];
105 greg 1.3 change++;
106 greg 2.44 }
107 greg 1.1 sprintf(buf, "view point (%.6g %.6g %.6g): ",
108     ourview.vp[0], ourview.vp[1], ourview.vp[2]);
109     (*dev->comout)(buf);
110 greg 1.18 (*dev->comin)(buf, NULL);
111 greg 2.4 if (buf[0] == CTRL('C')) return;
112 greg 2.6 if (sscanvec(buf, nv.vp))
113 greg 1.3 change++;
114 greg 1.1 sprintf(buf, "view direction (%.6g %.6g %.6g): ",
115 greg 2.45 ourview.vdir[0]*ourview.vdist,
116     ourview.vdir[1]*ourview.vdist,
117     ourview.vdir[2]*ourview.vdist);
118 greg 1.1 (*dev->comout)(buf);
119 greg 1.18 (*dev->comin)(buf, NULL);
120 greg 2.4 if (buf[0] == CTRL('C')) return;
121 greg 2.45 if (sscanvec(buf, nv.vdir)) {
122     nv.vdist = 1.;
123 greg 1.3 change++;
124 greg 2.45 }
125 greg 1.1 sprintf(buf, "view up (%.6g %.6g %.6g): ",
126     ourview.vup[0], ourview.vup[1], ourview.vup[2]);
127     (*dev->comout)(buf);
128 greg 1.18 (*dev->comin)(buf, NULL);
129 greg 2.4 if (buf[0] == CTRL('C')) return;
130 greg 2.6 if (sscanvec(buf, nv.vup))
131 greg 1.3 change++;
132 greg 1.1 sprintf(buf, "view horiz and vert size (%.6g %.6g): ",
133     ourview.horiz, ourview.vert);
134     (*dev->comout)(buf);
135 greg 1.18 (*dev->comin)(buf, NULL);
136 greg 2.4 if (buf[0] == CTRL('C')) return;
137 greg 1.3 if (sscanf(buf, "%lf %lf", &nv.horiz, &nv.vert) == 2)
138     change++;
139 greg 2.26 sprintf(buf, "fore and aft clipping plane (%.6g %.6g): ",
140     ourview.vfore, ourview.vaft);
141     (*dev->comout)(buf);
142     (*dev->comin)(buf, NULL);
143     if (buf[0] == CTRL('C')) return;
144     if (sscanf(buf, "%lf %lf", &nv.vfore, &nv.vaft) == 2)
145     change++;
146 greg 1.14 sprintf(buf, "view shift and lift (%.6g %.6g): ",
147     ourview.hoff, ourview.voff);
148 greg 1.1 (*dev->comout)(buf);
149 greg 1.18 (*dev->comin)(buf, NULL);
150 greg 2.4 if (buf[0] == CTRL('C')) return;
151 greg 1.14 if (sscanf(buf, "%lf %lf", &nv.hoff, &nv.voff) == 2)
152 greg 1.3 change++;
153     if (change)
154     newview(&nv);
155 greg 1.1 }
156    
157    
158 greg 2.55 void
159 schorsch 2.43 lastview( /* return to a previous view */
160     char *s
161     )
162 greg 1.1 {
163     char buf[128];
164 greg 1.10 char *fname;
165 greg 1.1 int success;
166     VIEW nv;
167 greg 2.62 /* get parameters from a file */
168     if (nextword(buf, sizeof(buf), s) != NULL) {
169 schorsch 2.42 nv = stdview;
170 greg 2.19 if ((fname = getpath(buf, "", R_OK)) == NULL ||
171 greg 1.33 (success = viewfile(fname, &nv, NULL)) == -1) {
172 greg 1.1 sprintf(errmsg, "cannot open \"%s\"", buf);
173     error(COMMAND, errmsg);
174     return;
175     }
176     if (!success)
177     error(COMMAND, "wrong file format");
178     else
179     newview(&nv);
180     return;
181     }
182 greg 1.20 if (oldview.type == 0) { /* no old view! */
183 greg 1.1 error(COMMAND, "no previous view");
184     return;
185     }
186 schorsch 2.42 nv = ourview;
187     ourview = oldview;
188     oldview = nv;
189 greg 2.55 newimage(NULL);
190 greg 2.18 }
191    
192    
193 greg 2.55 void
194 schorsch 2.43 saveview( /* save view to rad file */
195     char *s
196     )
197 greg 2.18 {
198     char view[64];
199     char *fname;
200     FILE *fp;
201    
202 greg 2.19 if (*atos(view, sizeof(view), s)) {
203     if (isint(view)) {
204     error(COMMAND, "cannot write view by number");
205     return;
206     }
207 greg 2.18 s = sskip(s);
208 greg 2.19 }
209 greg 2.62 if (nextword(rifname, sizeof(rifname), s) == NULL && !rifname[0]) {
210 greg 2.18 error(COMMAND, "no previous rad file");
211     return;
212     }
213     if ((fname = getpath(rifname, NULL, 0)) == NULL ||
214     (fp = fopen(fname, "a")) == NULL) {
215     sprintf(errmsg, "cannot open \"%s\"", rifname);
216     error(COMMAND, errmsg);
217     return;
218     }
219     fputs("view= ", fp);
220     fputs(view, fp);
221     fprintview(&ourview, fp);
222     putc('\n', fp);
223     fclose(fp);
224     }
225    
226    
227 greg 2.55 void
228 schorsch 2.43 loadview( /* load view from rad file */
229     char *s
230     )
231 greg 2.18 {
232     char buf[512];
233     char *fname;
234     FILE *fp;
235     VIEW nv;
236    
237     strcpy(buf, "rad -n -s -V -v ");
238     if (sscanf(s, "%s", buf+strlen(buf)) == 1)
239     s = sskip(s);
240     else
241     strcat(buf, "1");
242 greg 2.62 if (nextword(rifname, sizeof(rifname), s) == NULL && !rifname[0]) {
243 greg 2.18 error(COMMAND, "no previous rad file");
244     return;
245     }
246 greg 2.19 if ((fname = getpath(rifname, "", R_OK)) == NULL) {
247 greg 2.18 sprintf(errmsg, "cannot access \"%s\"", rifname);
248     error(COMMAND, errmsg);
249     return;
250     }
251     sprintf(buf+strlen(buf), " %s", fname);
252     if ((fp = popen(buf, "r")) == NULL) {
253     error(COMMAND, "cannot run rad");
254     return;
255     }
256     buf[0] = '\0';
257     fgets(buf, sizeof(buf), fp);
258     pclose(fp);
259 schorsch 2.42 nv = stdview;
260 greg 2.18 if (!sscanview(&nv, buf)) {
261     error(COMMAND, "rad error -- no such view?");
262     return;
263     }
264     newview(&nv);
265 greg 1.1 }
266    
267    
268 greg 2.55 void
269 schorsch 2.43 getaim( /* aim camera */
270     char *s
271     )
272 greg 1.1 {
273 greg 2.44 VIEW nv = ourview;
274 greg 1.1 double zfact;
275    
276     if (getinterest(s, 1, nv.vdir, &zfact) < 0)
277     return;
278 greg 1.21 zoomview(&nv, zfact);
279 greg 1.1 newview(&nv);
280     }
281    
282    
283 greg 2.55 void
284 greg 2.45 getfocus( /* set focus distance */
285     char *s
286     )
287     {
288 greg 2.53 char buf[64];
289 greg 2.45 double dist;
290    
291     if (sscanf(s, "%lf", &dist) < 1) {
292     int x, y;
293     RAY thisray;
294     if (dev->getcur == NULL)
295     return;
296     (*dev->comout)("Pick focus point\n");
297     if ((*dev->getcur)(&x, &y) == ABORT)
298     return;
299     if ((thisray.rmax = viewray(thisray.rorg, thisray.rdir,
300     &ourview, (x+.5)/hresolu, (y+.5)/vresolu)) < -FTINY) {
301     error(COMMAND, "not on image");
302     return;
303     }
304 greg 2.49 rayorigin(&thisray, PRIMARY, NULL, NULL);
305 greg 2.45 if (!localhit(&thisray, &thescene)) {
306     error(COMMAND, "not a local object");
307     return;
308     }
309     dist = thisray.rot;
310     } else if (dist <= .0) {
311     error(COMMAND, "focus distance must be positive");
312     return;
313     }
314     ourview.vdist = dist;
315 greg 2.53 sprintf(buf, "Focus distance set to %f\n", dist);
316     (*dev->comout)(buf);
317 greg 2.45 }
318    
319    
320 greg 2.55 void
321 schorsch 2.43 getmove( /* move camera */
322     char *s
323     )
324 greg 1.1 {
325     FVECT vc;
326     double mag;
327    
328     if (getinterest(s, 0, vc, &mag) < 0)
329     return;
330 greg 1.4 moveview(0.0, 0.0, mag, vc);
331 greg 1.1 }
332    
333    
334 greg 2.55 void
335 schorsch 2.43 getrotate( /* rotate camera */
336     char *s
337     )
338 greg 1.1 {
339 greg 2.44 VIEW nv = ourview;
340 greg 1.1 FVECT v1;
341 greg 1.4 double angle, elev, zfact;
342 greg 1.1
343 greg 1.4 elev = 0.0; zfact = 1.0;
344     if (sscanf(s, "%lf %lf %lf", &angle, &elev, &zfact) < 1) {
345 greg 1.1 error(COMMAND, "missing angle");
346     return;
347     }
348     spinvector(nv.vdir, ourview.vdir, ourview.vup, angle*(PI/180.));
349     if (elev != 0.0) {
350     fcross(v1, nv.vdir, ourview.vup);
351     normalize(v1);
352     spinvector(nv.vdir, nv.vdir, v1, elev*(PI/180.));
353     }
354 greg 1.21 zoomview(&nv, zfact);
355 greg 1.1 newview(&nv);
356     }
357    
358    
359 greg 2.55 void
360 schorsch 2.43 getpivot( /* pivot viewpoint */
361 greg 2.55 char *s
362 schorsch 2.43 )
363 greg 1.1 {
364     FVECT vc;
365 greg 1.4 double angle, elev, mag;
366 greg 1.1
367 greg 1.4 elev = 0.0;
368     if (sscanf(s, "%lf %lf", &angle, &elev) < 1) {
369 greg 1.1 error(COMMAND, "missing angle");
370     return;
371     }
372 greg 2.33 if (getinterest(sskip2(s,2), 0, vc, &mag) < 0)
373 greg 1.1 return;
374 greg 1.4 moveview(angle, elev, mag, vc);
375 greg 1.1 }
376    
377    
378 greg 2.55 void
379 schorsch 2.43 getexposure( /* get new exposure */
380     char *s
381     )
382 greg 1.1 {
383     char buf[128];
384 greg 2.55 char *cp;
385 greg 1.1 int x, y;
386 greg 2.55 PNODE *p = &ptrunk;
387 greg 2.16 int adapt = 0;
388     double e = 1.0;
389 greg 2.25
390 greg 1.1 for (cp = s; isspace(*cp); cp++)
391     ;
392 greg 2.16 if (*cp == '@') {
393     adapt++;
394 greg 2.25 while (isspace(*++cp))
395     ;
396 greg 2.16 }
397 greg 1.1 if (*cp == '\0') { /* normalize to point */
398     if (dev->getcur == NULL)
399     return;
400     (*dev->comout)("Pick point for exposure\n");
401     if ((*dev->getcur)(&x, &y) == ABORT)
402     return;
403 greg 2.55 p = findrect(x, y, &ptrunk, -1);
404 greg 1.1 } else {
405     if (*cp == '=') { /* absolute setting */
406     p = NULL;
407     e = 1.0/exposure;
408     for (cp++; isspace(*cp); cp++)
409     ;
410     if (*cp == '\0') { /* interactive */
411 greg 2.15 sprintf(buf, "exposure (%f): ", exposure);
412 greg 1.1 (*dev->comout)(buf);
413 greg 1.18 (*dev->comin)(buf, NULL);
414 greg 1.1 for (cp = buf; isspace(*cp); cp++)
415     ;
416     if (*cp == '\0')
417     return;
418     }
419     }
420     if (*cp == '+' || *cp == '-') /* f-stops */
421     e *= pow(2.0, atof(cp));
422     else /* multiplier */
423     e *= atof(cp);
424     }
425     if (p != NULL) { /* relative setting */
426 greg 2.58 compavg(p);
427 greg 2.14 if (bright(p->v) < 1e-15) {
428 greg 1.1 error(COMMAND, "cannot normalize to zero");
429     return;
430     }
431 greg 2.16 if (adapt)
432     e *= 106./pow(1.219+pow(luminance(p->v)/exposure,.4),2.5)/exposure;
433     else
434     e *= 0.5 / bright(p->v);
435 greg 1.1 }
436     if (e <= FTINY || fabs(1.0 - e) <= FTINY)
437     return;
438     scalepict(&ptrunk, e);
439     exposure *= e;
440     redraw();
441     }
442    
443 greg 2.38 typedef union {int i; double d; COLOR C;} *MyUptr;
444 greg 2.57 #define FEQ(x,y) (fabs((x)-(y)) <= FTINY)
445 greg 1.1
446 greg 2.55 int
447 schorsch 2.43 getparam( /* get variable from user */
448     char *str,
449     char *dsc,
450     int typ,
451     void *p
452     )
453 greg 1.24 {
454 greg 2.55 MyUptr ptr = (MyUptr)p;
455 greg 1.24 int i0;
456     double d0, d1, d2;
457     char buf[48];
458    
459     switch (typ) {
460     case 'i': /* integer */
461     if (sscanf(str, "%d", &i0) != 1) {
462     (*dev->comout)(dsc);
463     sprintf(buf, " (%d): ", ptr->i);
464     (*dev->comout)(buf);
465     (*dev->comin)(buf, NULL);
466     if (sscanf(buf, "%d", &i0) != 1)
467 greg 1.31 return(0);
468 greg 1.24 }
469 greg 2.57 if (ptr->i == i0)
470     return(0);
471 greg 1.24 ptr->i = i0;
472 greg 2.55 break;
473 greg 1.24 case 'r': /* real */
474     if (sscanf(str, "%lf", &d0) != 1) {
475     (*dev->comout)(dsc);
476     sprintf(buf, " (%.6g): ", ptr->d);
477     (*dev->comout)(buf);
478     (*dev->comin)(buf, NULL);
479     if (sscanf(buf, "%lf", &d0) != 1)
480 greg 1.31 return(0);
481 greg 1.24 }
482 greg 2.57 if (FEQ(ptr->d, d0))
483     return(0);
484 greg 1.24 ptr->d = d0;
485 greg 2.55 break;
486 greg 1.24 case 'b': /* boolean */
487     if (sscanf(str, "%1s", buf) != 1) {
488     (*dev->comout)(dsc);
489 greg 2.2 sprintf(buf, "? (%c): ", ptr->i ? 'y' : 'n');
490 greg 1.24 (*dev->comout)(buf);
491     (*dev->comin)(buf, NULL);
492 greg 2.57 if (buf[0] == '\0')
493 greg 1.31 return(0);
494 greg 1.24 }
495 greg 2.57 if (strchr("yY+1tTnN-0fF", buf[0]) == NULL)
496     return(0);
497     i0 = strchr("yY+1tT", buf[0]) != NULL;
498     if (ptr->i == i0)
499     return(0);
500     ptr->i = i0;
501 greg 2.55 break;
502 greg 1.24 case 'C': /* color */
503     if (sscanf(str, "%lf %lf %lf", &d0, &d1, &d2) != 3) {
504     (*dev->comout)(dsc);
505     sprintf(buf, " (%.6g %.6g %.6g): ",
506     colval(ptr->C,RED),
507     colval(ptr->C,GRN),
508     colval(ptr->C,BLU));
509     (*dev->comout)(buf);
510     (*dev->comin)(buf, NULL);
511     if (sscanf(buf, "%lf %lf %lf", &d0, &d1, &d2) != 3)
512 greg 1.31 return(0);
513 greg 1.24 }
514 greg 2.57 if (FEQ(colval(ptr->C,RED), d0) &&
515     FEQ(colval(ptr->C,GRN), d1) &&
516     FEQ(colval(ptr->C,BLU), d2))
517     return(0);
518 greg 1.24 setcolor(ptr->C, d0, d1, d2);
519 greg 2.55 break;
520     default:
521     return(0); /* shouldn't happen */
522 greg 1.24 }
523 greg 2.55 newparam++;
524     return(1);
525 greg 1.24 }
526    
527    
528 greg 2.55 void
529 schorsch 2.43 setparam( /* get/set program parameter */
530 greg 2.55 char *s
531 schorsch 2.43 )
532 greg 1.1 {
533 greg 2.59 int prev_newp = newparam;
534 greg 1.24 char buf[128];
535 greg 1.1
536     if (s[0] == '\0') {
537 greg 1.23 (*dev->comout)(
538 greg 2.60 "aa ab ad ar as av aw b bv dc dv dj ds dt i lr lw me ma mg ms ps pt ss st u: ");
539 greg 1.18 (*dev->comin)(buf, NULL);
540 greg 1.1 s = buf;
541     }
542     switch (s[0]) {
543 greg 2.51 case 'u': /* uncorrelated sampling */
544     getparam(s+1, "uncorrelated sampling", 'b',
545 greg 2.50 (void *)&rand_samp);
546     break;
547 greg 1.1 case 'l': /* limit */
548     switch (s[1]) {
549     case 'w': /* weight */
550 greg 2.38 getparam(s+2, "limit weight", 'r',
551     (void *)&minweight);
552 greg 1.1 break;
553     case 'r': /* reflection */
554 greg 2.38 getparam(s+2, "limit reflection", 'i',
555     (void *)&maxdepth);
556 greg 1.1 break;
557     default:
558     goto badparam;
559     }
560     break;
561 greg 1.5 case 'd': /* direct */
562     switch (s[1]) {
563     case 'j': /* jitter */
564 greg 2.38 getparam(s+2, "direct jitter", 'r',
565     (void *)&dstrsrc);
566 greg 1.9 break;
567     case 'c': /* certainty */
568 greg 2.38 getparam(s+2, "direct certainty", 'r',
569     (void *)&shadcert);
570 greg 1.5 break;
571 greg 1.8 case 't': /* threshold */
572 greg 2.38 getparam(s+2, "direct threshold", 'r',
573     (void *)&shadthresh);
574 greg 1.5 break;
575 greg 2.13 case 'v': /* visibility */
576 greg 2.38 getparam(s+2, "direct visibility", 'b',
577     (void *)&directvis);
578 greg 1.30 break;
579     case 's': /* sampling */
580 greg 2.38 getparam(s+2, "direct sampling", 'r',
581     (void *)&srcsizerat);
582 greg 1.23 break;
583 greg 1.5 default:
584 greg 1.1 goto badparam;
585     }
586 greg 1.23 break;
587 greg 2.27 case 'b': /* back faces or black and white */
588     switch (s[1]) {
589     case 'v': /* back face visibility */
590 greg 2.38 getparam(s+2, "back face visibility", 'b',
591     (void *)&backvis);
592 greg 2.27 break;
593     case '\0': /* black and white */
594 greg 2.29 case ' ':
595 greg 2.27 case 'y': case 'Y': case 't': case 'T': case '1': case '+':
596     case 'n': case 'N': case 'f': case 'F': case '0': case '-':
597 greg 2.38 getparam(s+1, "black and white", 'b',
598     (void *)&greyscale);
599 greg 2.59 newparam = prev_newp;
600 greg 2.27 break;
601     default:
602     goto badparam;
603     }
604 greg 1.23 break;
605     case 'i': /* irradiance */
606 greg 2.38 getparam(s+1, "irradiance", 'b',
607     (void *)&do_irrad);
608 greg 1.1 break;
609     case 'a': /* ambient */
610     switch (s[1]) {
611     case 'v': /* value */
612 greg 2.38 getparam(s+2, "ambient value", 'C',
613     (void *)ambval);
614 greg 2.34 break;
615     case 'w': /* weight */
616 greg 2.38 getparam(s+2, "ambient value weight", 'i',
617     (void *)&ambvwt);
618 greg 1.1 break;
619     case 'a': /* accuracy */
620 greg 2.38 if (getparam(s+2, "ambient accuracy", 'r',
621     (void *)&ambacc))
622 greg 2.17 setambacc(ambacc);
623 greg 1.1 break;
624     case 'd': /* divisions */
625 greg 2.38 getparam(s+2, "ambient divisions", 'i',
626     (void *)&ambdiv);
627 greg 1.1 break;
628     case 's': /* samples */
629 greg 2.38 getparam(s+2, "ambient super-samples", 'i',
630     (void *)&ambssamp);
631 greg 1.1 break;
632     case 'b': /* bounces */
633 greg 2.38 getparam(s+2, "ambient bounces", 'i',
634     (void *)&ambounce);
635 greg 1.1 break;
636     case 'r':
637 greg 2.38 if (getparam(s+2, "ambient resolution", 'i',
638     (void *)&ambres))
639 greg 2.8 setambres(ambres);
640 greg 2.32 break;
641     default:
642     goto badparam;
643     }
644     break;
645     case 'm': /* medium */
646     switch (s[1]) {
647     case 'e': /* extinction coefficient */
648     getparam(s+2, "extinction coefficient", 'C',
649 greg 2.38 (void *)cextinction);
650 greg 2.32 break;
651     case 'a': /* scattering albedo */
652 greg 2.35 getparam(s+2, "scattering albedo", 'C',
653 greg 2.38 (void *)salbedo);
654 greg 2.32 break;
655     case 'g': /* scattering eccentricity */
656 greg 2.38 getparam(s+2, "scattering eccentricity", 'r',
657     (void *)&seccg);
658 greg 2.32 break;
659     case 's': /* sampling distance */
660     getparam(s+2, "mist sampling distance", 'r',
661 greg 2.38 (void *)&ssampdist);
662 greg 1.1 break;
663     default:
664     goto badparam;
665     }
666     break;
667 greg 2.5 case 'p': /* pixel */
668 greg 1.1 switch (s[1]) {
669 greg 2.5 case 's': /* sample */
670 greg 2.38 if (getparam(s+2, "pixel sample", 'i',
671     (void *)&psample))
672 greg 1.31 pdepth = 0;
673 greg 1.1 break;
674 greg 1.8 case 't': /* threshold */
675 greg 2.38 if (getparam(s+2, "pixel threshold", 'r',
676     (void *)&maxdiff))
677 greg 1.31 pdepth = 0;
678 greg 2.5 break;
679     default:
680     goto badparam;
681     }
682 greg 2.59 newparam = prev_newp;
683 greg 2.5 break;
684     case 's': /* specular */
685     switch (s[1]) {
686 greg 2.60 case 's': /* sampling */
687     getparam(s+2, "specular sampling", 'r',
688 greg 2.38 (void *)&specjitter);
689 greg 2.5 break;
690     case 't': /* threshold */
691 greg 2.38 getparam(s+2, "specular threshold", 'r',
692     (void *)&specthresh);
693 greg 1.1 break;
694     default:
695     goto badparam;
696     }
697     break;
698     case '\0': /* nothing */
699     break;
700     default:;
701     badparam:
702 greg 1.26 *sskip(s) = '\0';
703 greg 1.1 sprintf(errmsg, "%s: unknown variable", s);
704     error(COMMAND, errmsg);
705     break;
706     }
707     }
708    
709    
710 greg 2.55 void
711     traceray( /* trace a single ray */
712     char *s
713     )
714 greg 1.1 {
715 greg 2.61 RAY thisray;
716     char buf[512];
717 greg 1.1
718 greg 2.26 thisray.rmax = 0.0;
719    
720 greg 2.6 if (!sscanvec(s, thisray.rorg) ||
721 greg 2.33 !sscanvec(sskip2(s,3), thisray.rdir)) {
722 greg 2.61 int x, y;
723 greg 1.1
724     if (dev->getcur == NULL)
725     return;
726     (*dev->comout)("Pick ray\n");
727     if ((*dev->getcur)(&x, &y) == ABORT)
728     return;
729    
730 greg 2.26 if ((thisray.rmax = viewray(thisray.rorg, thisray.rdir,
731     &ourview, (x+.5)/hresolu, (y+.5)/vresolu)) < -FTINY) {
732 greg 1.21 error(COMMAND, "not on image");
733     return;
734     }
735    
736 greg 1.1 } else if (normalize(thisray.rdir) == 0.0) {
737     error(COMMAND, "zero ray direction");
738     return;
739     }
740    
741 greg 2.56 ray_trace(&thisray);
742 greg 1.1
743     if (thisray.ro == NULL)
744     (*dev->comout)("ray hit nothing");
745     else {
746 greg 2.47 OBJREC *mat = NULL;
747 greg 2.52 OBJREC *mod = NULL;
748 greg 2.47 char matspec[256];
749 greg 2.61 OBJREC *ino;
750    
751 greg 2.47 matspec[0] = '\0';
752     if (thisray.ro->omod != OVOID) {
753     mod = objptr(thisray.ro->omod);
754     mat = findmaterial(mod);
755     }
756     if (thisray.rod < 0.0)
757     strcpy(matspec, "back of ");
758     if (mod != NULL) {
759     strcat(matspec, mod->oname);
760 greg 2.48 if (mat != mod && mat != NULL)
761     sprintf(matspec+strlen(matspec),
762     " (%s)", mat->oname);
763 greg 2.47 } else
764     strcat(matspec, VOIDID);
765     sprintf(buf, "ray hit %s %s \"%s\"", matspec,
766 greg 1.1 ofun[thisray.ro->otype].funame,
767     thisray.ro->oname);
768 greg 2.37 if ((ino = objptr(thisray.robj)) != thisray.ro)
769     sprintf(buf+strlen(buf), " in %s \"%s\"",
770     ofun[ino->otype].funame, ino->oname);
771 greg 1.1 (*dev->comout)(buf);
772 greg 1.18 (*dev->comin)(buf, NULL);
773 greg 1.1 if (thisray.rot >= FHUGE)
774     (*dev->comout)("at infinity");
775     else {
776 greg 2.30 sprintf(buf, "at (%.6g %.6g %.6g) (%.6g)",
777     thisray.rop[0], thisray.rop[1],
778     thisray.rop[2], thisray.rt);
779 greg 1.1 (*dev->comout)(buf);
780     }
781 greg 1.18 (*dev->comin)(buf, NULL);
782 greg 2.38 sprintf(buf, "value (%.5g %.5g %.5g) (%.3gL)",
783 greg 1.1 colval(thisray.rcol,RED),
784     colval(thisray.rcol,GRN),
785 greg 1.32 colval(thisray.rcol,BLU),
786     luminance(thisray.rcol));
787 greg 1.1 (*dev->comout)(buf);
788     }
789 greg 1.18 (*dev->comin)(buf, NULL);
790 greg 1.1 }
791    
792    
793 greg 2.55 void
794     writepict( /* write the picture to a file */
795     char *s
796     )
797 greg 1.1 {
798     static char buf[128];
799 greg 1.10 char *fname;
800 greg 1.1 FILE *fp;
801     COLR *scanline;
802     int y;
803 greg 2.62 /* XXX relies on words.c 2.11 behavior */
804     if (nextword(buf, sizeof(buf), s) == NULL && !buf[0]) {
805 greg 1.1 error(COMMAND, "no file");
806     return;
807     }
808 greg 1.10 if ((fname = getpath(buf, NULL, 0)) == NULL ||
809 greg 1.1 (fp = fopen(fname, "w")) == NULL) {
810     sprintf(errmsg, "cannot open \"%s\"", buf);
811     error(COMMAND, errmsg);
812     return;
813     }
814 schorsch 2.40 SET_FILE_BINARY(fp);
815 greg 1.1 (*dev->comout)("writing \"");
816     (*dev->comout)(fname);
817     (*dev->comout)("\"...\n");
818     /* write header */
819 greg 2.24 newheader("RADIANCE", fp);
820 greg 1.1 fputs(progname, fp);
821     fprintview(&ourview, fp);
822 greg 1.29 if (octname != NULL)
823     fprintf(fp, " %s\n", octname);
824     else
825     putc('\n', fp);
826 greg 1.28 fprintf(fp, "SOFTWARE= %s\n", VersionID);
827 greg 2.38 fputnow(fp);
828 greg 1.1 if (exposure != 1.0)
829 greg 1.13 fputexpos(exposure, fp);
830 greg 1.15 if (dev->pixaspect != 1.0)
831     fputaspect(dev->pixaspect, fp);
832 greg 1.22 fputformat(COLRFMT, fp);
833 greg 1.11 putc('\n', fp);
834 greg 1.33 fprtresolu(hresolu, vresolu, fp);
835 greg 1.1
836 greg 1.14 scanline = (COLR *)malloc(hresolu*sizeof(COLR));
837 greg 1.27 if (scanline == NULL) {
838     error(COMMAND, "not enough memory!");
839     fclose(fp);
840     unlink(fname);
841     return;
842     }
843 greg 1.14 for (y = vresolu-1; y >= 0; y--) {
844     getpictcolrs(y, scanline, &ptrunk, hresolu, vresolu);
845     if (fwritecolrs(scanline, hresolu, fp) < 0)
846 greg 1.1 break;
847     }
848 greg 2.38 free((void *)scanline);
849 greg 1.1 if (fclose(fp) < 0)
850     error(COMMAND, "write error");
851     }