ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/rv2.c
Revision: 2.58
Committed: Fri Sep 5 19:45:41 2008 UTC (15 years, 7 months ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.57: +2 -1 lines
Log Message:
Fixed bug that broke "exposure" command

File Contents

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