ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/rv2.c
Revision: 2.65
Committed: Fri Apr 11 20:27:23 2014 UTC (9 years, 11 months ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad4R2P2, rad5R0, rad5R1, rad4R2, rad4R2P1
Changes since 2.64: +2 -6 lines
Log Message:
Made reading of vectors more consistent with FVFORMAT macro

File Contents

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