ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/src/rt/rv2.c
Revision: 2.12
Committed: Thu Oct 15 21:13:14 1992 UTC (33 years ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.11: +3 -0 lines
Log Message: 
added setmode(fd,O_BINARY) where needed for MSDOS

File Contents

# User Rev Content
1 greg 2.9 /* Copyright (c) 1992 Regents of the University of California */
2 greg 1.1
3     #ifndef lint
4     static char SCCSid[] = "$SunId$ LBL";
5     #endif
6    
7     /*
8     * rv2.c - command routines used in tracing a view.
9     *
10     * 3/24/87
11     */
12    
13     #include "ray.h"
14    
15     #include "octree.h"
16    
17     #include "otypes.h"
18    
19     #include "rpaint.h"
20    
21 greg 1.33 #include "resolu.h"
22    
23 greg 1.1 #include <ctype.h>
24    
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    
34 greg 1.29 extern char VersionID[];
35 greg 1.1 extern char *progname;
36 greg 1.29 extern char *octname;
37 greg 1.1
38    
39     getframe(s) /* get a new frame */
40     char *s;
41     {
42 greg 1.12 if (getrect(s, &pframe) < 0)
43 greg 1.1 return;
44     pdepth = 0;
45     }
46    
47    
48 greg 1.12 getrepaint(s) /* get area and repaint */
49     char *s;
50     {
51     RECT box;
52    
53     if (getrect(s, &box) < 0)
54     return;
55 greg 1.14 paintrect(&ptrunk, 0, 0, hresolu, vresolu, &box);
56 greg 1.12 }
57    
58    
59 greg 1.1 getview(s) /* get/show view parameters */
60     char *s;
61     {
62     FILE *fp;
63     char buf[128];
64 greg 1.10 char *fname;
65 greg 1.3 int change = 0;
66 greg 1.1 VIEW nv;
67    
68 greg 2.9 while (isspace(*s))
69     s++;
70     if (*s == '-') { /* command line parameters */
71     copystruct(&nv, &ourview);
72     if (sscanview(&nv, s))
73     newview(&nv);
74     else
75     error(COMMAND, "bad view option(s)");
76     return;
77     }
78 greg 1.1 if (sscanf(s, "%s", buf) == 1) { /* write parameters to a file */
79 greg 1.10 if ((fname = getpath(buf, NULL, 0)) == NULL ||
80 greg 1.1 (fp = fopen(fname, "a")) == NULL) {
81     sprintf(errmsg, "cannot open \"%s\"", buf);
82     error(COMMAND, errmsg);
83     return;
84     }
85 greg 1.2 fputs(progname, fp);
86 greg 1.1 fprintview(&ourview, fp);
87 greg 1.19 fputs(sskip(s), fp);
88 greg 1.1 fputs("\n", fp);
89     fclose(fp);
90     return;
91     }
92     sprintf(buf, "view type (%c): ", ourview.type);
93     (*dev->comout)(buf);
94 greg 1.18 (*dev->comin)(buf, NULL);
95 greg 2.4 if (buf[0] == CTRL('C')) return;
96 greg 1.3 if (buf[0] && buf[0] != ourview.type) {
97 greg 1.1 nv.type = buf[0];
98 greg 1.3 change++;
99     } else
100 greg 1.1 nv.type = ourview.type;
101     sprintf(buf, "view point (%.6g %.6g %.6g): ",
102     ourview.vp[0], ourview.vp[1], ourview.vp[2]);
103     (*dev->comout)(buf);
104 greg 1.18 (*dev->comin)(buf, NULL);
105 greg 2.4 if (buf[0] == CTRL('C')) return;
106 greg 2.6 if (sscanvec(buf, nv.vp))
107 greg 1.3 change++;
108     else
109 greg 1.1 VCOPY(nv.vp, ourview.vp);
110     sprintf(buf, "view direction (%.6g %.6g %.6g): ",
111     ourview.vdir[0], ourview.vdir[1], ourview.vdir[2]);
112     (*dev->comout)(buf);
113 greg 1.18 (*dev->comin)(buf, NULL);
114 greg 2.4 if (buf[0] == CTRL('C')) return;
115 greg 2.6 if (sscanvec(buf, nv.vdir))
116 greg 1.3 change++;
117     else
118 greg 1.1 VCOPY(nv.vdir, ourview.vdir);
119     sprintf(buf, "view up (%.6g %.6g %.6g): ",
120     ourview.vup[0], ourview.vup[1], ourview.vup[2]);
121     (*dev->comout)(buf);
122 greg 1.18 (*dev->comin)(buf, NULL);
123 greg 2.4 if (buf[0] == CTRL('C')) return;
124 greg 2.6 if (sscanvec(buf, nv.vup))
125 greg 1.3 change++;
126     else
127 greg 1.1 VCOPY(nv.vup, ourview.vup);
128     sprintf(buf, "view horiz and vert size (%.6g %.6g): ",
129     ourview.horiz, ourview.vert);
130     (*dev->comout)(buf);
131 greg 1.18 (*dev->comin)(buf, NULL);
132 greg 2.4 if (buf[0] == CTRL('C')) return;
133 greg 1.3 if (sscanf(buf, "%lf %lf", &nv.horiz, &nv.vert) == 2)
134     change++;
135     else {
136 greg 1.1 nv.horiz = ourview.horiz; nv.vert = ourview.vert;
137     }
138 greg 1.14 sprintf(buf, "view shift and lift (%.6g %.6g): ",
139     ourview.hoff, ourview.voff);
140 greg 1.1 (*dev->comout)(buf);
141 greg 1.18 (*dev->comin)(buf, NULL);
142 greg 2.4 if (buf[0] == CTRL('C')) return;
143 greg 1.14 if (sscanf(buf, "%lf %lf", &nv.hoff, &nv.voff) == 2)
144 greg 1.3 change++;
145     else {
146 greg 1.14 nv.hoff = ourview.hoff; nv.voff = ourview.voff;
147 greg 1.1 }
148 greg 1.3 if (change)
149     newview(&nv);
150 greg 1.1 }
151    
152    
153     lastview(s) /* return to a previous view */
154     char *s;
155     {
156     char buf[128];
157 greg 1.10 char *fname;
158 greg 1.1 int success;
159     VIEW nv;
160    
161     if (sscanf(s, "%s", buf) == 1) { /* get parameters from a file */
162 greg 1.17 copystruct(&nv, &stdview);
163 greg 1.10 if ((fname = getpath(buf, NULL, 0)) == NULL ||
164 greg 1.33 (success = viewfile(fname, &nv, NULL)) == -1) {
165 greg 1.1 sprintf(errmsg, "cannot open \"%s\"", buf);
166     error(COMMAND, errmsg);
167     return;
168     }
169     if (!success)
170     error(COMMAND, "wrong file format");
171     else
172     newview(&nv);
173     return;
174     }
175 greg 1.20 if (oldview.type == 0) { /* no old view! */
176 greg 1.1 error(COMMAND, "no previous view");
177     return;
178     }
179 greg 1.17 copystruct(&nv, &ourview);
180     copystruct(&ourview, &oldview);
181     copystruct(&oldview, &nv);
182 greg 1.1 newimage();
183     }
184    
185    
186     getaim(s) /* aim camera */
187     char *s;
188     {
189     double zfact;
190     VIEW nv;
191    
192     if (getinterest(s, 1, nv.vdir, &zfact) < 0)
193     return;
194 greg 1.21 nv.type = ourview.type;
195 greg 1.1 VCOPY(nv.vp, ourview.vp);
196     VCOPY(nv.vup, ourview.vup);
197 greg 1.14 nv.hoff = ourview.hoff; nv.voff = ourview.voff;
198 greg 1.21 nv.horiz = ourview.horiz; nv.vert = ourview.vert;
199     zoomview(&nv, zfact);
200 greg 1.1 newview(&nv);
201     }
202    
203    
204     getmove(s) /* move camera */
205     char *s;
206     {
207     FVECT vc;
208     double mag;
209    
210     if (getinterest(s, 0, vc, &mag) < 0)
211     return;
212 greg 1.4 moveview(0.0, 0.0, mag, vc);
213 greg 1.1 }
214    
215    
216     getrotate(s) /* rotate camera */
217     char *s;
218     {
219     VIEW nv;
220     FVECT v1;
221 greg 1.4 double angle, elev, zfact;
222 greg 1.1
223 greg 1.4 elev = 0.0; zfact = 1.0;
224     if (sscanf(s, "%lf %lf %lf", &angle, &elev, &zfact) < 1) {
225 greg 1.1 error(COMMAND, "missing angle");
226     return;
227     }
228 greg 1.21 nv.type = ourview.type;
229 greg 1.1 VCOPY(nv.vp, ourview.vp);
230     VCOPY(nv.vup, ourview.vup);
231 greg 1.14 nv.hoff = ourview.hoff; nv.voff = ourview.voff;
232 greg 1.1 spinvector(nv.vdir, ourview.vdir, ourview.vup, angle*(PI/180.));
233     if (elev != 0.0) {
234     fcross(v1, nv.vdir, ourview.vup);
235     normalize(v1);
236     spinvector(nv.vdir, nv.vdir, v1, elev*(PI/180.));
237     }
238 greg 1.21 nv.horiz = ourview.horiz; nv.vert = ourview.vert;
239     zoomview(&nv, zfact);
240 greg 1.1 newview(&nv);
241     }
242    
243    
244     getpivot(s) /* pivot viewpoint */
245     register char *s;
246     {
247     FVECT vc;
248 greg 1.4 double angle, elev, mag;
249 greg 1.1
250 greg 1.4 elev = 0.0;
251     if (sscanf(s, "%lf %lf", &angle, &elev) < 1) {
252 greg 1.1 error(COMMAND, "missing angle");
253     return;
254     }
255 greg 1.4 if (getinterest(sskip(sskip(s)), 0, vc, &mag) < 0)
256 greg 1.1 return;
257 greg 1.4 moveview(angle, elev, mag, vc);
258 greg 1.1 }
259    
260    
261     getexposure(s) /* get new exposure */
262     char *s;
263     {
264     char buf[128];
265     register char *cp;
266     register PNODE *p;
267     RECT r;
268     int x, y;
269     double e;
270    
271     for (cp = s; isspace(*cp); cp++)
272     ;
273     if (*cp == '\0') { /* normalize to point */
274     if (dev->getcur == NULL)
275     return;
276     (*dev->comout)("Pick point for exposure\n");
277     if ((*dev->getcur)(&x, &y) == ABORT)
278     return;
279     r.l = r.d = 0;
280 greg 1.14 r.r = hresolu; r.u = vresolu;
281 greg 1.1 p = findrect(x, y, &ptrunk, &r, -1);
282     e = 1.0;
283     } else {
284     if (*cp == '=') { /* absolute setting */
285     p = NULL;
286     e = 1.0/exposure;
287     for (cp++; isspace(*cp); cp++)
288     ;
289     if (*cp == '\0') { /* interactive */
290     sprintf(buf, "exposure (%lf): ", exposure);
291     (*dev->comout)(buf);
292 greg 1.18 (*dev->comin)(buf, NULL);
293 greg 1.1 for (cp = buf; isspace(*cp); cp++)
294     ;
295     if (*cp == '\0')
296     return;
297     }
298     } else { /* normalize to average */
299     p = &ptrunk;
300     e = 1.0;
301     }
302     if (*cp == '+' || *cp == '-') /* f-stops */
303     e *= pow(2.0, atof(cp));
304     else /* multiplier */
305     e *= atof(cp);
306     }
307     if (p != NULL) { /* relative setting */
308 greg 2.7 if (bright(p->v) < 1e-25) {
309 greg 1.1 error(COMMAND, "cannot normalize to zero");
310     return;
311     }
312 greg 1.6 e *= 0.5 / bright(p->v);
313 greg 1.1 }
314     if (e <= FTINY || fabs(1.0 - e) <= FTINY)
315     return;
316     scalepict(&ptrunk, e);
317     exposure *= e;
318     redraw();
319     }
320    
321    
322 greg 1.24 getparam(str, dsc, typ, ptr) /* get variable from user */
323     char *str, *dsc;
324     int typ;
325     register union {int i; double d; COLOR C;} *ptr;
326     {
327 greg 1.25 extern char *index();
328 greg 1.24 int i0;
329     double d0, d1, d2;
330     char buf[48];
331    
332     switch (typ) {
333     case 'i': /* integer */
334     if (sscanf(str, "%d", &i0) != 1) {
335     (*dev->comout)(dsc);
336     sprintf(buf, " (%d): ", ptr->i);
337     (*dev->comout)(buf);
338     (*dev->comin)(buf, NULL);
339     if (sscanf(buf, "%d", &i0) != 1)
340 greg 1.31 return(0);
341 greg 1.24 }
342     ptr->i = i0;
343 greg 1.31 return(1);
344 greg 1.24 case 'r': /* real */
345     if (sscanf(str, "%lf", &d0) != 1) {
346     (*dev->comout)(dsc);
347     sprintf(buf, " (%.6g): ", ptr->d);
348     (*dev->comout)(buf);
349     (*dev->comin)(buf, NULL);
350     if (sscanf(buf, "%lf", &d0) != 1)
351 greg 1.31 return(0);
352 greg 1.24 }
353     ptr->d = d0;
354 greg 1.31 return(1);
355 greg 1.24 case 'b': /* boolean */
356     if (sscanf(str, "%1s", buf) != 1) {
357     (*dev->comout)(dsc);
358 greg 2.2 sprintf(buf, "? (%c): ", ptr->i ? 'y' : 'n');
359 greg 1.24 (*dev->comout)(buf);
360     (*dev->comin)(buf, NULL);
361 greg 1.25 if (buf[0] == '\0' ||
362     index("yY+1tTnN-0fF", buf[0]) == NULL)
363 greg 1.31 return(0);
364 greg 1.24 }
365 greg 1.25 ptr->i = index("yY+1tT", buf[0]) != NULL;
366 greg 1.31 return(1);
367 greg 1.24 case 'C': /* color */
368     if (sscanf(str, "%lf %lf %lf", &d0, &d1, &d2) != 3) {
369     (*dev->comout)(dsc);
370     sprintf(buf, " (%.6g %.6g %.6g): ",
371     colval(ptr->C,RED),
372     colval(ptr->C,GRN),
373     colval(ptr->C,BLU));
374     (*dev->comout)(buf);
375     (*dev->comin)(buf, NULL);
376     if (sscanf(buf, "%lf %lf %lf", &d0, &d1, &d2) != 3)
377 greg 1.31 return(0);
378 greg 1.24 }
379     setcolor(ptr->C, d0, d1, d2);
380 greg 1.31 return(1);
381 greg 1.24 }
382     }
383    
384    
385 greg 1.1 setparam(s) /* get/set program parameter */
386     register char *s;
387     {
388     extern int psample;
389     extern double maxdiff;
390     extern double minweight;
391     extern int maxdepth;
392     extern double dstrsrc;
393 greg 1.5 extern double shadthresh;
394 greg 1.9 extern double shadcert;
395 greg 1.1 extern COLOR ambval;
396     extern double ambacc;
397     extern int ambres;
398     extern int ambdiv;
399     extern int ambssamp;
400     extern int ambounce;
401 greg 1.23 extern int directinvis;
402 greg 1.30 extern double srcsizerat;
403 greg 1.23 extern int do_irrad;
404 greg 2.5 extern double specjitter;
405     extern double specthresh;
406 greg 1.24 char buf[128];
407 greg 1.1
408     if (s[0] == '\0') {
409 greg 1.23 (*dev->comout)(
410 greg 2.5 "aa ab ad ar as av b dc di dj ds dt i lr lw ps pt sj st: ");
411 greg 1.18 (*dev->comin)(buf, NULL);
412 greg 1.1 s = buf;
413     }
414     switch (s[0]) {
415     case 'l': /* limit */
416     switch (s[1]) {
417     case 'w': /* weight */
418 greg 1.24 getparam(s+2, "limit weight", 'r', &minweight);
419 greg 1.1 break;
420     case 'r': /* reflection */
421 greg 1.24 getparam(s+2, "limit reflection", 'i', &maxdepth);
422 greg 1.1 break;
423     default:
424     goto badparam;
425     }
426     break;
427 greg 1.5 case 'd': /* direct */
428     switch (s[1]) {
429     case 'j': /* jitter */
430 greg 1.24 getparam(s+2, "direct jitter", 'r', &dstrsrc);
431 greg 1.9 break;
432     case 'c': /* certainty */
433 greg 1.24 getparam(s+2, "direct certainty", 'r', &shadcert);
434 greg 1.5 break;
435 greg 1.8 case 't': /* threshold */
436 greg 1.24 getparam(s+2, "direct threshold", 'r', &shadthresh);
437 greg 1.5 break;
438 greg 1.23 case 'i': /* invisibility */
439 greg 1.24 getparam(s+2, "direct invisibility",
440     'b', &directinvis);
441 greg 1.30 break;
442     case 's': /* sampling */
443     getparam(s+2, "direct sampling", 'r', &srcsizerat);
444 greg 1.23 break;
445 greg 1.5 default:
446 greg 1.1 goto badparam;
447     }
448 greg 1.23 break;
449     case 'b': /* black and white */
450 greg 1.24 getparam(s+1, "black and white", 'b', &greyscale);
451 greg 1.23 break;
452     case 'i': /* irradiance */
453 greg 1.24 getparam(s+1, "irradiance", 'b', &do_irrad);
454 greg 1.1 break;
455     case 'a': /* ambient */
456     switch (s[1]) {
457     case 'v': /* value */
458 greg 2.11 getparam(s+2, "ambient value", 'C', (COLOR *)ambval);
459 greg 1.1 break;
460     case 'a': /* accuracy */
461 greg 1.24 getparam(s+2, "ambient accuracy", 'r', &ambacc);
462 greg 1.1 break;
463     case 'd': /* divisions */
464 greg 1.24 getparam(s+2, "ambient divisions", 'i', &ambdiv);
465 greg 1.1 break;
466     case 's': /* samples */
467 greg 1.24 getparam(s+2, "ambient super-samples", 'i', &ambssamp);
468 greg 1.1 break;
469     case 'b': /* bounces */
470 greg 1.24 getparam(s+2, "ambient bounces", 'i', &ambounce);
471 greg 1.1 break;
472     case 'r':
473 greg 1.31 if (getparam(s+2, "ambient resolution", 'i', &ambres))
474 greg 2.8 setambres(ambres);
475 greg 1.1 break;
476     default:
477     goto badparam;
478     }
479     break;
480 greg 2.5 case 'p': /* pixel */
481 greg 1.1 switch (s[1]) {
482 greg 2.5 case 's': /* sample */
483     if (getparam(s+2, "pixel sample", 'i', &psample))
484 greg 1.31 pdepth = 0;
485 greg 1.1 break;
486 greg 1.8 case 't': /* threshold */
487 greg 2.5 if (getparam(s+2, "pixel threshold", 'r', &maxdiff))
488 greg 1.31 pdepth = 0;
489 greg 2.5 break;
490     default:
491     goto badparam;
492     }
493     break;
494     case 's': /* specular */
495     switch (s[1]) {
496     case 'j': /* jitter */
497     getparam(s+2, "specular jitter", 'r', &specjitter);
498     break;
499     case 't': /* threshold */
500     getparam(s+2, "specular threshold", 'r', &specthresh);
501 greg 1.1 break;
502     default:
503     goto badparam;
504     }
505     break;
506     case '\0': /* nothing */
507     break;
508     default:;
509     badparam:
510 greg 1.26 *sskip(s) = '\0';
511 greg 1.1 sprintf(errmsg, "%s: unknown variable", s);
512     error(COMMAND, errmsg);
513     break;
514     }
515     }
516    
517    
518     traceray(s) /* trace a single ray */
519     char *s;
520     {
521     char buf[128];
522     int x, y;
523     RAY thisray;
524    
525 greg 2.6 if (!sscanvec(s, thisray.rorg) ||
526     !sscanvec(sskip(sskip(sskip(s))), thisray.rdir)) {
527 greg 1.1
528     if (dev->getcur == NULL)
529     return;
530     (*dev->comout)("Pick ray\n");
531     if ((*dev->getcur)(&x, &y) == ABORT)
532     return;
533    
534 greg 1.21 if (viewray(thisray.rorg, thisray.rdir, &ourview,
535     (x+.5)/hresolu, (y+.5)/vresolu) < 0) {
536     error(COMMAND, "not on image");
537     return;
538     }
539    
540 greg 1.1 } else if (normalize(thisray.rdir) == 0.0) {
541     error(COMMAND, "zero ray direction");
542     return;
543     }
544    
545     rayorigin(&thisray, NULL, PRIMARY, 1.0);
546    
547     rayvalue(&thisray);
548    
549     if (thisray.ro == NULL)
550     (*dev->comout)("ray hit nothing");
551     else {
552     sprintf(buf, "ray hit %s %s \"%s\"",
553     objptr(thisray.ro->omod)->oname,
554     ofun[thisray.ro->otype].funame,
555     thisray.ro->oname);
556     (*dev->comout)(buf);
557 greg 1.18 (*dev->comin)(buf, NULL);
558 greg 1.1 if (thisray.rot >= FHUGE)
559     (*dev->comout)("at infinity");
560     else {
561     sprintf(buf, "at (%.6g %.6g %.6g)", thisray.rop[0],
562     thisray.rop[1], thisray.rop[2]);
563     (*dev->comout)(buf);
564     }
565 greg 1.18 (*dev->comin)(buf, NULL);
566 greg 1.32 sprintf(buf, "value (%.5g %.5g %.5g) (%.1fL)",
567 greg 1.1 colval(thisray.rcol,RED),
568     colval(thisray.rcol,GRN),
569 greg 1.32 colval(thisray.rcol,BLU),
570     luminance(thisray.rcol));
571 greg 1.1 (*dev->comout)(buf);
572     }
573 greg 1.18 (*dev->comin)(buf, NULL);
574 greg 1.1 }
575    
576    
577     writepict(s) /* write the picture to a file */
578     char *s;
579     {
580     static char buf[128];
581 greg 1.10 char *fname;
582 greg 1.1 FILE *fp;
583     COLR *scanline;
584     int y;
585    
586     if (sscanf(s, "%s", buf) != 1 && buf[0] == '\0') {
587     error(COMMAND, "no file");
588     return;
589     }
590 greg 1.10 if ((fname = getpath(buf, NULL, 0)) == NULL ||
591 greg 1.1 (fp = fopen(fname, "w")) == NULL) {
592     sprintf(errmsg, "cannot open \"%s\"", buf);
593     error(COMMAND, errmsg);
594     return;
595     }
596 greg 2.12 #ifdef MSDOS
597     setmode(fileno(fp), O_BINARY);
598     #endif
599 greg 1.1 (*dev->comout)("writing \"");
600     (*dev->comout)(fname);
601     (*dev->comout)("\"...\n");
602     /* write header */
603     fputs(progname, fp);
604     fprintview(&ourview, fp);
605 greg 1.29 if (octname != NULL)
606     fprintf(fp, " %s\n", octname);
607     else
608     putc('\n', fp);
609 greg 1.28 fprintf(fp, "SOFTWARE= %s\n", VersionID);
610 greg 1.1 if (exposure != 1.0)
611 greg 1.13 fputexpos(exposure, fp);
612 greg 1.15 if (dev->pixaspect != 1.0)
613     fputaspect(dev->pixaspect, fp);
614 greg 1.22 fputformat(COLRFMT, fp);
615 greg 1.11 putc('\n', fp);
616 greg 1.33 fprtresolu(hresolu, vresolu, fp);
617 greg 1.1
618 greg 1.14 scanline = (COLR *)malloc(hresolu*sizeof(COLR));
619 greg 1.27 if (scanline == NULL) {
620     error(COMMAND, "not enough memory!");
621     fclose(fp);
622     unlink(fname);
623     return;
624     }
625 greg 1.14 for (y = vresolu-1; y >= 0; y--) {
626     getpictcolrs(y, scanline, &ptrunk, hresolu, vresolu);
627     if (fwritecolrs(scanline, hresolu, fp) < 0)
628 greg 1.1 break;
629     }
630 greg 1.27 free((char *)scanline);
631 greg 1.1 if (fclose(fp) < 0)
632     error(COMMAND, "write error");
633     }