ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/util/ranimate.c
Revision: 2.27
Committed: Wed Jun 3 16:37:51 1998 UTC (25 years, 10 months ago) by gwlarson
Content type: text/plain
Branch: MAIN
Changes since 2.26: +41 -41 lines
Log Message:
made variable ordering alphabetical for easier parsing

File Contents

# Content
1 /* Copyright (c) 1996 Regents of the University of California */
2
3 #ifndef lint
4 static char SCCSid[] = "$SunId$ LBL";
5 #endif
6
7 /*
8 * Radiance animation control program
9 */
10
11 #include "standard.h"
12 #include <ctype.h>
13 #include <sys/types.h>
14 #include <sys/stat.h>
15 #include "view.h"
16 #include "vars.h"
17 #include "netproc.h"
18 /* default remote shell */
19 #ifdef _AUX_SOURCE
20 #define REMSH "remsh"
21 #else
22 #define REMSH "rsh"
23 #endif
24 /* input variables (alphabetical by name) */
25 #define ANIMATE 0 /* animation command */
26 #define ARCHIVE 1 /* archiving command */
27 #define BASENAME 2 /* output image base name */
28 #define DIRECTORY 3 /* working (sub)directory */
29 #define DISKSPACE 4 /* how much disk space to use */
30 #define END 5 /* ending frame number */
31 #define EXPOSURE 6 /* how to compute exposure */
32 #define HOST 7 /* rendering host machine */
33 #define INTERP 8 /* # frames to interpolate */
34 #define MBLUR 9 /* samples for motion blur */
35 #define NEXTANIM 10 /* next animation file */
36 #define OCTREE 11 /* octree file name */
37 #define OVERSAMP 12 /* # times to oversample image */
38 #define PFILT 13 /* pfilt options */
39 #define PINTERP 14 /* pinterp options */
40 #define RENDER 15 /* rendering options */
41 #define RESOLUTION 16 /* desired final resolution */
42 #define RIF 17 /* rad input file */
43 #define RSH 18 /* remote shell script or program */
44 #define RTRACE 19 /* use rtrace with pinterp? */
45 #define START 20 /* starting frame number */
46 #define TRANSFER 21 /* frame transfer command */
47 #define VIEWFILE 22 /* animation frame views */
48
49 int NVARS = 23; /* total number of variables */
50
51 VARIABLE vv[] = { /* variable-value pairs */
52 {"ANIMATE", 2, 0, NULL, onevalue},
53 {"ARCHIVE", 2, 0, NULL, onevalue},
54 {"BASENAME", 3, 0, NULL, onevalue},
55 {"DIRECTORY", 3, 0, NULL, onevalue},
56 {"DISKSPACE", 3, 0, NULL, fltvalue},
57 {"END", 2, 0, NULL, intvalue},
58 {"EXPOSURE", 3, 0, NULL, onevalue},
59 {"host", 4, 0, NULL, NULL},
60 {"INTERPOLATE", 3, 0, NULL, intvalue},
61 {"MBLUR", 2, 0, NULL, onevalue},
62 {"NEXTANIM", 3, 0, NULL, onevalue},
63 {"OCTREE", 3, 0, NULL, onevalue},
64 {"OVERSAMPLE", 2, 0, NULL, fltvalue},
65 {"pfilt", 2, 0, NULL, catvalues},
66 {"pinterp", 2, 0, NULL, catvalues},
67 {"render", 3, 0, NULL, catvalues},
68 {"RESOLUTION", 3, 0, NULL, onevalue},
69 {"RIF", 3, 0, NULL, onevalue},
70 {"RSH", 3, 0, NULL, onevalue},
71 {"RTRACE", 2, 0, NULL, boolvalue},
72 {"START", 2, 0, NULL, intvalue},
73 {"TRANSFER", 2, 0, NULL, onevalue},
74 {"VIEWFILE", 2, 0, NULL, onevalue},
75 };
76
77 #define SFNAME "STATUS" /* status file name */
78
79 struct {
80 char host[64]; /* control host name */
81 int pid; /* control process id */
82 char cfname[128]; /* control file name */
83 int rnext; /* next frame to render */
84 int fnext; /* next frame to filter */
85 int tnext; /* next frame to transfer */
86 } astat; /* animation status */
87
88 char *progname; /* our program name */
89 char *cfname; /* our control file name */
90
91 int nowarn = 0; /* turn warnings off? */
92 int silent = 0; /* silent mode? */
93 int noaction = 0; /* take no action? */
94
95 char *remsh; /* remote shell program/script */
96 char rendopt[2048]; /* rendering options */
97 char rresopt[32]; /* rendering resolution options */
98 char fresopt[32]; /* filter resolution options */
99 int pfiltalways; /* always use pfilt? */
100
101 char arcargs[10240]; /* files to archive */
102 char *arcfirst, *arcnext; /* pointers to first and next argument */
103
104 struct pslot {
105 int pid; /* process ID (0 if empty) */
106 int fout; /* output frame number */
107 int (*rcvf)(); /* recover function */
108 } *pslot; /* process slots */
109 int npslots; /* number of process slots */
110
111 #define phostname(ps) ((ps)->hostname[0] ? (ps)->hostname : astat.host)
112
113 struct pslot *findpslot();
114
115 PSERVER *lastpserver; /* last process server with error */
116
117 VIEW *getview();
118 char *getexp(), *dirfile();
119
120 extern time_t fdate(), time();
121
122
123 main(argc, argv)
124 int argc;
125 char *argv[];
126 {
127 int explicate = 0;
128 int i;
129
130 progname = argv[0]; /* get arguments */
131 for (i = 1; i < argc && argv[i][0] == '-'; i++)
132 switch (argv[i][1]) {
133 case 'e': /* print variables */
134 explicate++;
135 break;
136 case 'w': /* turn off warnings */
137 nowarn++;
138 break;
139 case 's': /* silent mode */
140 silent++;
141 break;
142 case 'n': /* take no action */
143 noaction++;
144 break;
145 default:
146 goto userr;
147 }
148 if (i != argc-1)
149 goto userr;
150 cfname = argv[i];
151 /* load variables */
152 loadvars(cfname);
153 /* check variables */
154 checkvalues();
155 /* did we get DIRECTORY? */
156 checkdir();
157 /* check status */
158 if (getastat() < 0) {
159 fprintf(stderr, "%s: exiting\n", progname);
160 quit(1);
161 }
162 /* pfilt always if options given */
163 pfiltalways = vdef(PFILT);
164 /* load RIF if any */
165 if (vdef(RIF))
166 getradfile(vval(RIF));
167 /* set defaults */
168 setdefaults();
169 /* print variables */
170 if (explicate)
171 printvars(stdout);
172 /* set up process servers */
173 sethosts();
174 /* run animation */
175 animate();
176 /* all done */
177 if (vdef(NEXTANIM)) {
178 argv[i] = vval(NEXTANIM); /* just change input file */
179 if (!silent)
180 printargs(argc, argv, stdout);
181 if ((argv[0] = getpath(progname,getenv("PATH"),X_OK)) == NULL)
182 fprintf(stderr, "%s: command not found\n", progname);
183 else
184 execv(progname, argv);
185 quit(1);
186 }
187 quit(0);
188 userr:
189 fprintf(stderr, "Usage: %s [-s][-n][-w][-e] anim_file\n", progname);
190 quit(1);
191 }
192
193
194 getastat() /* check/set animation status */
195 {
196 char sfname[256];
197 FILE *fp;
198
199 sprintf(sfname, "%s/%s", vval(DIRECTORY), SFNAME);
200 if ((fp = fopen(sfname, "r")) == NULL) {
201 if (errno != ENOENT) {
202 perror(sfname);
203 return(-1);
204 }
205 astat.rnext = astat.fnext = astat.tnext = 0;
206 goto setours;
207 }
208 if (fscanf(fp, "Control host: %s\n", astat.host) != 1)
209 goto fmterr;
210 if (fscanf(fp, "Control PID: %d\n", &astat.pid) != 1)
211 goto fmterr;
212 if (fscanf(fp, "Control file: %s\n", astat.cfname) != 1)
213 goto fmterr;
214 if (fscanf(fp, "Next render: %d\n", &astat.rnext) != 1)
215 goto fmterr;
216 if (fscanf(fp, "Next filter: %d\n", &astat.fnext) != 1)
217 goto fmterr;
218 if (fscanf(fp, "Next transfer: %d\n", &astat.tnext) != 1)
219 goto fmterr;
220 fclose(fp);
221 if (astat.pid != 0) { /* thinks it's still running */
222 if (strcmp(myhostname(), astat.host)) {
223 fprintf(stderr,
224 "%s: process %d may still be running on host %s\n",
225 progname, astat.pid, astat.host);
226 return(-1);
227 }
228 if (kill(astat.pid, 0) != -1 || errno != ESRCH) {
229 fprintf(stderr, "%s: process %d is still running\n",
230 progname, astat.pid);
231 return(-1);
232 }
233 /* assume it is dead */
234 }
235 if (strcmp(cfname, astat.cfname) && astat.pid != 0) { /* other's */
236 fprintf(stderr, "%s: unfinished job \"%s\"\n",
237 progname, astat.cfname);
238 return(-1);
239 }
240 /* check control file mods. */
241 if (!nowarn && fdate(cfname) > fdate(sfname))
242 fprintf(stderr,
243 "%s: warning - control file modified since last run\n",
244 progname);
245 setours: /* set our values */
246 strcpy(astat.host, myhostname());
247 astat.pid = getpid();
248 strcpy(astat.cfname, cfname);
249 return(0);
250 fmterr:
251 fprintf(stderr, "%s: format error in status file \"%s\"\n",
252 progname, sfname);
253 fclose(fp);
254 return(-1);
255 }
256
257
258 putastat() /* put out current status */
259 {
260 char buf[256];
261 FILE *fp;
262
263 if (noaction)
264 return;
265 sprintf(buf, "%s/%s", vval(DIRECTORY), SFNAME);
266 if ((fp = fopen(buf, "w")) == NULL) {
267 perror(buf);
268 quit(1);
269 }
270 fprintf(fp, "Control host: %s\n", astat.host);
271 fprintf(fp, "Control PID: %d\n", astat.pid);
272 fprintf(fp, "Control file: %s\n", astat.cfname);
273 fprintf(fp, "Next render: %d\n", astat.rnext);
274 fprintf(fp, "Next filter: %d\n", astat.fnext);
275 fprintf(fp, "Next transfer: %d\n", astat.tnext);
276 fclose(fp);
277 }
278
279
280 checkdir() /* make sure we have our directory */
281 {
282 struct stat stb;
283
284 if (!vdef(DIRECTORY)) {
285 fprintf(stderr, "%s: %s undefined\n",
286 progname, vnam(DIRECTORY));
287 quit(1);
288 }
289 if (stat(vval(DIRECTORY), &stb) == -1) {
290 if (errno == ENOENT && mkdir(vval(DIRECTORY), 0777) == 0)
291 return;
292 perror(vval(DIRECTORY));
293 quit(1);
294 }
295 if (!(stb.st_mode & S_IFDIR)) {
296 fprintf(stderr, "%s: not a directory\n", vval(DIRECTORY));
297 quit(1);
298 }
299 }
300
301
302 setdefaults() /* set default values */
303 {
304 extern char *atos();
305 int decades;
306 char buf[256];
307
308 if (vdef(ANIMATE)) {
309 vval(OCTREE) = NULL;
310 vdef(OCTREE) = 0;
311 } else if (!vdef(OCTREE)) {
312 fprintf(stderr, "%s: either %s or %s must be defined\n",
313 progname, vnam(OCTREE), vnam(ANIMATE));
314 quit(1);
315 }
316 if (!vdef(VIEWFILE)) {
317 fprintf(stderr, "%s: %s undefined\n", progname, vnam(VIEWFILE));
318 quit(1);
319 }
320 if (!vdef(HOST)) {
321 vval(HOST) = LHOSTNAME;
322 vdef(HOST)++;
323 }
324 if (!vdef(START)) {
325 vval(START) = "1";
326 vdef(START)++;
327 }
328 if (!vdef(END)) {
329 sprintf(buf, "%d", countviews()+vint(START)-1);
330 vval(END) = savqstr(buf);
331 vdef(END)++;
332 }
333 if (vint(END) < vint(START)) {
334 fprintf(stderr, "%s: ending frame less than starting frame\n",
335 progname);
336 quit(1);
337 }
338 if (!vdef(BASENAME)) {
339 decades = (int)log10((double)vint(END)) + 1;
340 if (decades < 3) decades = 3;
341 sprintf(buf, "%s/frame%%0%dd", vval(DIRECTORY), decades);
342 vval(BASENAME) = savqstr(buf);
343 vdef(BASENAME)++;
344 }
345 if (!vdef(RESOLUTION)) {
346 vval(RESOLUTION) = "640";
347 vdef(RESOLUTION)++;
348 }
349 if (!vdef(OVERSAMP)) {
350 vval(OVERSAMP) = "2";
351 vdef(OVERSAMP)++;
352 }
353 if (!vdef(INTERP)) {
354 vval(INTERP) = "0";
355 vdef(INTERP)++;
356 }
357 if (!vdef(MBLUR)) {
358 vval(MBLUR) = "0";
359 vdef(MBLUR)++;
360 }
361 if (!vdef(RTRACE)) {
362 vval(RTRACE) = "F";
363 vdef(RTRACE)++;
364 }
365 if (!vdef(DISKSPACE)) {
366 if (!nowarn)
367 fprintf(stderr,
368 "%s: warning - no %s setting, assuming 100 Mbytes available\n",
369 progname, vnam(DISKSPACE));
370 vval(DISKSPACE) = "100";
371 vdef(DISKSPACE)++;
372 }
373 if (!vdef(RSH)) {
374 vval(RSH) = REMSH;
375 vdef(RSH)++;
376 }
377 /* locate remote shell program */
378 atos(buf, sizeof(buf), vval(RSH));
379 if ((remsh = getpath(buf, getenv("PATH"), X_OK)) != NULL)
380 remsh = savqstr(remsh);
381 else
382 remsh = vval(RSH); /* will generate error if used */
383
384 /* append rendering options */
385 if (vdef(RENDER))
386 sprintf(rendopt+strlen(rendopt), " %s", vval(RENDER));
387 }
388
389
390 sethosts() /* set up process servers */
391 {
392 extern char *iskip();
393 char buf[256], *dir, *uname;
394 int np;
395 register char *cp;
396 int i;
397
398 npslots = 0;
399 if (noaction)
400 return;
401 for (i = 0; i < vdef(HOST); i++) { /* add each host */
402 dir = uname = NULL;
403 np = 1;
404 strcpy(cp=buf, nvalue(HOST, i)); /* copy to buffer */
405 cp = sskip(cp); /* skip host name */
406 while (isspace(*cp))
407 *cp++ = '\0';
408 if (*cp) { /* has # processes? */
409 np = atoi(cp);
410 if ((cp = iskip(cp)) == NULL || (*cp && !isspace(*cp)))
411 badvalue(HOST);
412 while (isspace(*cp))
413 cp++;
414 if (*cp) { /* has directory? */
415 dir = cp;
416 cp = sskip(cp); /* skip dir. */
417 while (isspace(*cp))
418 *cp++ = '\0';
419 if (*cp) { /* has user? */
420 uname = cp;
421 if (*sskip(cp))
422 badvalue(HOST);
423 }
424 }
425 }
426 if (addpserver(buf, dir, uname, np) == NULL) {
427 if (!nowarn)
428 fprintf(stderr,
429 "%s: cannot execute on host \"%s\"\n",
430 progname, buf);
431 } else
432 npslots += np;
433 }
434 if (npslots == 0) {
435 fprintf(stderr, "%s: no working process servers\n", progname);
436 quit(1);
437 }
438 pslot = (struct pslot *)calloc(npslots, sizeof(struct pslot));
439 if (pslot == NULL) {
440 perror("malloc");
441 quit(1);
442 }
443 }
444
445
446 getradfile(rfargs) /* run rad and get needed variables */
447 char *rfargs;
448 {
449 static short mvar[] = {OCTREE,PFILT,RESOLUTION,EXPOSURE,-1};
450 char combuf[256];
451 register int i;
452 register char *cp;
453 char *pippt;
454 /* create rad command */
455 sprintf(rendopt, " @%s/render.opt", vval(DIRECTORY));
456 sprintf(combuf,
457 "rad -v 0 -s -e -w %s OPTFILE=%s | egrep '^[ \t]*(NOMATCH",
458 rfargs, rendopt+2);
459 cp = combuf;
460 while (*cp) {
461 if (*cp == '|') pippt = cp;
462 cp++;
463 } /* match unset variables */
464 for (i = 0; mvar[i] >= 0; i++)
465 if (!vdef(mvar[i])) {
466 *cp++ = '|';
467 strcpy(cp, vnam(mvar[i]));
468 while (*cp) cp++;
469 pippt = NULL;
470 }
471 if (pippt != NULL)
472 strcpy(pippt, "> /dev/null"); /* nothing to match */
473 else {
474 sprintf(cp, ")[ \t]*=' > %s/radset.var", vval(DIRECTORY));
475 cp += 11; /* point to file name */
476 }
477 if (system(combuf)) {
478 fprintf(stderr, "%s: error executing rad command:\n\t%s\n",
479 progname, combuf);
480 quit(1);
481 }
482 if (pippt == NULL) { /* load variables and remove file */
483 loadvars(cp);
484 unlink(cp);
485 }
486 }
487
488
489 animate() /* run animation */
490 {
491 int xres, yres;
492 float pa, mult;
493 int frames_batch;
494 register int i;
495 double d1, d2;
496 /* compute rpict resolution */
497 i = sscanf(vval(RESOLUTION), "%d %d %f", &xres, &yres, &pa);
498 mult = vflt(OVERSAMP);
499 if (i == 3) {
500 sprintf(rresopt, "-x %d -y %d -pa %.3f", (int)(mult*xres),
501 (int)(mult*yres), pa);
502 sprintf(fresopt, "-x %d -y %d -pa %.3f", xres, yres, pa);
503 } else if (i) {
504 if (i == 1) yres = xres;
505 sprintf(rresopt, "-x %d -y %d", (int)(mult*xres),
506 (int)(mult*yres));
507 sprintf(fresopt, "-x %d -y %d -pa 1", xres, yres);
508 } else
509 badvalue(RESOLUTION);
510 /* consistency checks */
511 if (vdef(ANIMATE)) {
512 if (vint(INTERP)) {
513 if (!nowarn)
514 fprintf(stderr,
515 "%s: resetting %s=0 for animation\n",
516 progname, vnam(INTERP));
517 vval(INTERP) = "0";
518 }
519 if (atoi(vval(MBLUR))) { /* can't handle this yet */
520 if (!nowarn)
521 fprintf(stderr,
522 "%s: resetting %s=0 for animation\n",
523 progname, vnam(MBLUR));
524 vval(MBLUR) = "0";
525 }
526 }
527 /* figure # frames per batch */
528 d1 = mult*xres*mult*yres*4; /* space for orig. picture */
529 if ((i=vint(INTERP)) || atoi(vval(MBLUR)))
530 d1 += mult*xres*mult*yres*sizeof(float); /* Z-buffer */
531 d2 = xres*yres*4; /* space for final picture */
532 frames_batch = (i+1)*(vflt(DISKSPACE)*1048576.-d1)/(d1+i*d2);
533 if (frames_batch < i+2) {
534 fprintf(stderr, "%s: insufficient disk space allocated\n",
535 progname);
536 quit(1);
537 }
538 /* initialize archive argument list */
539 i = vdef(ARCHIVE) ? strlen(vval(ARCHIVE))+132 : 132;
540 arcnext = arcfirst = arcargs + i;
541 /* initialize status file */
542 if (astat.rnext == 0)
543 astat.rnext = astat.fnext = astat.tnext = vint(START);
544 putastat();
545 /* render in batches */
546 while (astat.tnext <= vint(END)) {
547 renderframes(frames_batch);
548 filterframes();
549 transferframes();
550 }
551 /* mark status as finished */
552 astat.pid = 0;
553 putastat();
554 /* close open files */
555 getview(0);
556 getexp(0);
557 }
558
559
560 renderframes(nframes) /* render next nframes frames */
561 int nframes;
562 {
563 static char vendbuf[16];
564 VIEW *vp;
565 FILE *fp = NULL;
566 char vfname[128];
567 int lastframe;
568 register int i;
569
570 if (astat.tnext < astat.rnext) /* other work to do first */
571 return;
572 /* create batch view file */
573 if (!vdef(ANIMATE)) {
574 sprintf(vfname, "%s/anim.vf", vval(DIRECTORY));
575 if ((fp = fopen(vfname, "w")) == NULL) {
576 perror(vfname);
577 quit(1);
578 }
579 }
580 /* bound batch properly */
581 lastframe = astat.rnext + nframes - 1;
582 if ((lastframe-1) % (vint(INTERP)+1)) /* need even interval */
583 lastframe += vint(INTERP)+1 - ((lastframe-1)%(vint(INTERP)+1));
584 if (lastframe > vint(END)) /* check for end */
585 lastframe = vint(END);
586 /* render each view */
587 for (i = astat.rnext; i <= lastframe; i++) {
588 if ((vp = getview(i)) == NULL) {
589 if (!nowarn)
590 fprintf(stderr,
591 "%s: ran out of views before last frame\n",
592 progname);
593 sprintf(vval(END)=vendbuf, "%d", i-1);
594 lastframe = i - 1;
595 break;
596 }
597 if (vdef(ANIMATE)) /* animate frame */
598 animrend(i, vp);
599 else { /* else record it */
600 fputs(VIEWSTR, fp);
601 fprintview(vp, fp);
602 putc('\n', fp);
603 }
604 }
605 if (vdef(ANIMATE)) /* wait for renderings to finish */
606 bwait(0);
607 else { /* else if walk-through */
608 fclose(fp); /* close view file */
609 walkwait(astat.rnext, lastframe, vfname); /* walk it */
610 unlink(vfname); /* remove view file */
611 }
612 astat.rnext = i; /* update status */
613 putastat();
614 }
615
616
617 filterframes() /* catch up with filtering */
618 {
619 VIEW *vp;
620 register int i;
621
622 if (astat.tnext < astat.fnext) /* other work to do first */
623 return;
624 /* filter each view */
625 for (i = astat.fnext; i < astat.rnext; i++) {
626 if ((vp = getview(i)) == NULL) { /* get view i */
627 fprintf(stderr,
628 "%s: unexpected error reading view for frame %d\n",
629 progname, i);
630 quit(1);
631 }
632 dofilt(i, vp, getexp(i), 0); /* filter frame */
633 }
634 bwait(0); /* wait for filter processes */
635 archive(); /* archive originals */
636 astat.fnext = i; /* update status */
637 putastat();
638 }
639
640
641 transferframes() /* catch up with picture transfers */
642 {
643 char combuf[10240], *fbase;
644 register char *cp;
645 register int i;
646
647 if (astat.tnext >= astat.fnext) /* nothing to do, yet */
648 return;
649 if (!vdef(TRANSFER)) { /* no transfer function -- leave 'em */
650 astat.tnext = astat.fnext;
651 putastat(); /* update status */
652 return;
653 }
654 strcpy(combuf, "cd "); /* start transfer command */
655 fbase = dirfile(cp = combuf+3, vval(BASENAME));
656 if (*cp) {
657 while (*++cp) ;
658 *cp++ = ';'; *cp++ = ' ';
659 } else
660 cp = combuf;
661 strcpy(cp, vval(TRANSFER));
662 while (*cp) cp++;
663 /* make argument list */
664 for (i = astat.tnext; i < astat.fnext; i++) {
665 *cp++ = ' ';
666 sprintf(cp, fbase, i);
667 while (*cp) cp++;
668 strcpy(cp, ".pic");
669 cp += 4;
670 }
671 if (runcom(combuf)) { /* transfer frames */
672 fprintf(stderr, "%s: error running transfer command\n",
673 progname);
674 quit(1);
675 }
676 astat.tnext = i; /* update status */
677 putastat();
678 }
679
680
681 animrend(frame, vp) /* start animation frame */
682 int frame;
683 VIEW *vp;
684 {
685 extern int recover();
686 char combuf[2048];
687 char fname[128];
688
689 sprintf(fname, vval(BASENAME), frame);
690 strcat(fname, ".unf");
691 if (access(fname, F_OK) == 0)
692 return;
693 sprintf(combuf, "%s %d | rpict%s%s -w0 %s > %s", vval(ANIMATE), frame,
694 rendopt, viewopt(vp), rresopt, fname);
695 bruncom(combuf, frame, recover); /* run in background */
696 }
697
698
699 walkwait(first, last, vfn) /* walk-through frames */
700 int first, last;
701 char *vfn;
702 {
703 char combuf[2048];
704 char *inspoint;
705 register int i;
706
707 if (!noaction && vint(INTERP)) /* create dummy frames */
708 for (i = first; i <= last; i++)
709 if (i < vint(END) && (i-1) % (vint(INTERP)+1)) {
710 sprintf(combuf, vval(BASENAME), i);
711 strcat(combuf, ".unf");
712 close(open(combuf, O_RDONLY|O_CREAT, 0666));
713 }
714 /* create command */
715 sprintf(combuf, "rpict%s%s -w0", rendopt,
716 viewopt(getview(first>1 ? first-1 : 1)));
717 if (vint(INTERP) || atoi(vval(MBLUR)))
718 sprintf(combuf+strlen(combuf), " -z %s.zbf", vval(BASENAME));
719 sprintf(combuf+strlen(combuf), " -o %s.unf %s -S %d",
720 vval(BASENAME), rresopt, first);
721 inspoint = combuf + strlen(combuf);
722 sprintf(inspoint, " %s < %s", vval(OCTREE), vfn);
723 /* run in parallel */
724 i = (last-first+1)/(vint(INTERP)+1);
725 if (i < 1) i = 1;
726 if (pruncom(combuf, inspoint, i)) {
727 fprintf(stderr, "%s: error rendering frames %d through %d\n",
728 progname, first, last);
729 quit(1);
730 }
731 if (!noaction && vint(INTERP)) /* remove dummy frames */
732 for (i = first; i <= last; i++)
733 if (i < vint(END) && (i-1) % (vint(INTERP)+1)) {
734 sprintf(combuf, vval(BASENAME), i);
735 strcat(combuf, ".unf");
736 unlink(combuf);
737 }
738 }
739
740
741 int
742 recover(frame) /* recover the specified frame */
743 int frame;
744 {
745 static int *rfrm; /* list of recovered frames */
746 static int nrfrms = 0;
747 char combuf[2048];
748 char fname[128];
749 register char *cp;
750 register int i;
751 /* check to see if recovered already */
752 for (i = nrfrms; i--; )
753 if (rfrm[i] == frame)
754 return(0);
755 /* build command */
756 sprintf(fname, vval(BASENAME), frame);
757 if (vdef(ANIMATE))
758 sprintf(combuf, "%s %d | rpict%s -w0",
759 vval(ANIMATE), frame, rendopt);
760 else
761 sprintf(combuf, "rpict%s -w0", rendopt);
762 cp = combuf + strlen(combuf);
763 if (vint(INTERP) || atoi(vval(MBLUR))) {
764 sprintf(cp, " -z %s.zbf", fname);
765 while (*cp) cp++;
766 }
767 sprintf(cp, " -ro %s.unf", fname);
768 while (*cp) cp++;
769 if (!vdef(ANIMATE)) {
770 *cp++ = ' ';
771 strcpy(cp, vval(OCTREE));
772 }
773 if (runcom(combuf)) /* run command */
774 return(1);
775 /* add frame to recovered list */
776 if (nrfrms)
777 rfrm = (int *)realloc((char *)rfrm, (nrfrms+1)*sizeof(int));
778 else
779 rfrm = (int *)malloc(sizeof(int));
780 if (rfrm == NULL) {
781 perror("malloc");
782 quit(1);
783 }
784 rfrm[nrfrms++] = frame;
785 return(0);
786 }
787
788
789 int
790 frecover(frame) /* recover filtered frame */
791 int frame;
792 {
793 VIEW *vp;
794 char *ex;
795
796 vp = getview(frame);
797 ex = getexp(frame);
798 if (dofilt(frame, vp, ex, 2) && dofilt(frame, vp, ex, 1))
799 return(1);
800 return(0);
801 }
802
803
804 archive() /* archive and remove renderings */
805 {
806 #define RMCOML (sizeof(rmcom)-1)
807 static char rmcom[] = "rm -f";
808 char basedir[128];
809 int dlen, alen;
810 register int j;
811
812 if (arcnext == arcfirst)
813 return; /* nothing to do */
814 dirfile(basedir, vval(BASENAME));
815 dlen = strlen(basedir);
816 if (vdef(ARCHIVE)) { /* run archive command */
817 alen = strlen(vval(ARCHIVE));
818 if (dlen) {
819 j = alen + dlen + 5;
820 strncpy(arcfirst-j, "cd ", 3);
821 strncpy(arcfirst-j+3, basedir, dlen);
822 (arcfirst-j)[dlen+3] = ';'; (arcfirst-j)[dlen+4] = ' ';
823 } else
824 j = alen;
825 strncpy(arcfirst-alen, vval(ARCHIVE), alen);
826 if (runcom(arcfirst-j)) {
827 fprintf(stderr, "%s: error running archive command\n",
828 progname);
829 quit(1);
830 }
831 }
832 if (dlen) {
833 j = RMCOML + dlen + 5;
834 strncpy(arcfirst-j, "cd ", 3);
835 strncpy(arcfirst-j+3, basedir, dlen);
836 (arcfirst-j)[dlen+3] = ';'; (arcfirst-j)[dlen+4] = ' ';
837 } else
838 j = RMCOML;
839 /* run remove command */
840 strncpy(arcfirst-RMCOML, rmcom, RMCOML);
841 runcom(arcfirst-j);
842 arcnext = arcfirst; /* reset argument list */
843 #undef RMCOML
844 }
845
846
847 int
848 dofilt(frame, vp, ep, rvr) /* filter frame */
849 int frame;
850 VIEW *vp;
851 char *ep;
852 int rvr;
853 {
854 extern int frecover();
855 static int iter = 0;
856 char fnbefore[128], fnafter[128], *fbase;
857 char combuf[1024], fname0[128], fname1[128];
858 int usepinterp, usepfilt, nora_rgbe;
859 int frseq[2];
860 /* check what is needed */
861 usepinterp = atoi(vval(MBLUR));
862 usepfilt = pfiltalways | ep==NULL;
863 if (ep != NULL && !strcmp(ep, "1"))
864 ep = "+0";
865 nora_rgbe = strcmp(vval(OVERSAMP),"1") || ep==NULL ||
866 *ep != '+' || *ep != '-' || !isint(ep);
867 /* compute rendered views */
868 frseq[0] = frame - ((frame-1) % (vint(INTERP)+1));
869 frseq[1] = frseq[0] + vint(INTERP) + 1;
870 fbase = dirfile(NULL, vval(BASENAME));
871 if (frseq[1] > vint(END))
872 frseq[1] = vint(END);
873 if (frseq[1] == frame) { /* pfilt only */
874 frseq[0] = frseq[1];
875 usepinterp = 0; /* update what's needed */
876 usepfilt |= nora_rgbe;
877 } else if (frseq[0] == frame) { /* no interpolation needed */
878 if (!rvr && frame > 1+vint(INTERP)) { /* archive previous */
879 *arcnext++ = ' ';
880 sprintf(arcnext, fbase, frame-vint(INTERP)-1);
881 while (*arcnext) arcnext++;
882 strcpy(arcnext, ".unf");
883 arcnext += 4;
884 if (usepinterp || vint(INTERP)) { /* and Z-buf */
885 *arcnext++ = ' ';
886 sprintf(arcnext, fbase, frame-vint(INTERP)-1);
887 while (*arcnext) arcnext++;
888 strcpy(arcnext, ".zbf");
889 arcnext += 4;
890 }
891 }
892 if (!usepinterp) /* update what's needed */
893 usepfilt |= nora_rgbe;
894 } else /* interpolation needed */
895 usepinterp++;
896 if (frseq[1] >= astat.rnext) /* next batch unavailable */
897 frseq[1] = frseq[0];
898 sprintf(fnbefore, vval(BASENAME), frseq[0]);
899 sprintf(fnafter, vval(BASENAME), frseq[1]);
900 if (rvr == 1 && recover(frseq[0])) /* recover before frame? */
901 return(1);
902 /* generate command */
903 if (usepinterp) { /* using pinterp */
904 if (rvr == 2 && recover(frseq[1])) /* recover after? */
905 return(1);
906 if (atoi(vval(MBLUR))) {
907 sprintf(fname0, "%s/vw0%c", vval(DIRECTORY),
908 'a'+(iter%26));
909 sprintf(fname1, "%s/vw1%c", vval(DIRECTORY),
910 'a'+(iter%26));
911 if (!noaction) {
912 FILE *fp; /* motion blurring */
913 if ((fp = fopen(fname0, "w")) == NULL) {
914 perror(fname0); quit(1);
915 }
916 fputs(VIEWSTR, fp);
917 fprintview(vp, fp);
918 putc('\n', fp); fclose(fp);
919 if ((vp = getview(frame+1)) == NULL) {
920 fprintf(stderr,
921 "%s: unexpected error reading view for frame %d\n",
922 progname, frame+1);
923 quit(1);
924 }
925 if ((fp = fopen(fname1, "w")) == NULL) {
926 perror(fname1); quit(1);
927 }
928 fputs(VIEWSTR, fp);
929 fprintview(vp, fp);
930 putc('\n', fp); fclose(fp);
931 }
932 sprintf(combuf,
933 "(pmblur %s %d %s %s; rm -f %s %s) | pinterp -B",
934 *sskip(vval(MBLUR)) ? sskip2(vval(MBLUR),1) : "1",
935 atoi(vval(MBLUR)),
936 fname0, fname1, fname0, fname1);
937 iter++;
938 } else /* no blurring */
939 strcpy(combuf, "pinterp");
940 strcat(combuf, viewopt(vp));
941 if (vbool(RTRACE))
942 sprintf(combuf+strlen(combuf), " -ff -fr '%s -w0 %s'",
943 rendopt+1, vval(OCTREE));
944 if (vdef(PINTERP))
945 sprintf(combuf+strlen(combuf), " %s", vval(PINTERP));
946 if (usepfilt)
947 sprintf(combuf+strlen(combuf), " %s", rresopt);
948 else
949 sprintf(combuf+strlen(combuf), " %s -e %s",
950 fresopt, ep);
951 sprintf(combuf+strlen(combuf), " %s.unf %s.zbf",
952 fnbefore, fnbefore);
953 if (frseq[1] != frseq[0])
954 sprintf(combuf+strlen(combuf), " %s.unf %s.zbf",
955 fnafter, fnafter);
956 if (usepfilt) { /* also pfilt */
957 if (vdef(PFILT))
958 sprintf(combuf+strlen(combuf), " | pfilt %s",
959 vval(PFILT));
960 else
961 strcat(combuf, " | pfilt");
962 if (ep != NULL)
963 sprintf(combuf+strlen(combuf), " -1 -e %s %s",
964 ep, fresopt);
965 else
966 sprintf(combuf+strlen(combuf), " %s", fresopt);
967 }
968 } else if (usepfilt) { /* pfilt only */
969 if (rvr == 2)
970 return(1);
971 if (vdef(PFILT))
972 sprintf(combuf, "pfilt %s", vval(PFILT));
973 else
974 strcpy(combuf, "pfilt");
975 if (ep != NULL)
976 sprintf(combuf+strlen(combuf), " -1 -e %s %s %s.unf",
977 ep, fresopt, fnbefore);
978 else
979 sprintf(combuf+strlen(combuf), " %s %s.unf",
980 fresopt, fnbefore);
981 } else { /* else just check it */
982 if (rvr == 2)
983 return(1);
984 sprintf(combuf, "ra_rgbe -e %s -r %s.unf", ep, fnbefore);
985 }
986 /* output file name */
987 sprintf(fname0, vval(BASENAME), frame);
988 sprintf(combuf+strlen(combuf), " > %s.pic", fname0);
989 if (rvr) /* in recovery */
990 return(runcom(combuf));
991 bruncom(combuf, frame, frecover); /* else run in background */
992 return(0);
993 }
994
995
996 VIEW *
997 getview(n) /* get view number n */
998 int n;
999 {
1000 static FILE *viewfp = NULL; /* view file pointer */
1001 static int viewnum = 0; /* current view number */
1002 static VIEW curview = STDVIEW; /* current view */
1003 char linebuf[256];
1004
1005 if (n == 0) { /* signal to close file and clean up */
1006 if (viewfp != NULL) {
1007 fclose(viewfp);
1008 viewfp = NULL;
1009 viewnum = 0;
1010 copystruct(&curview, &stdview);
1011 }
1012 return(NULL);
1013 }
1014 if (viewfp == NULL) { /* open file */
1015 if ((viewfp = fopen(vval(VIEWFILE), "r")) == NULL) {
1016 perror(vval(VIEWFILE));
1017 quit(1);
1018 }
1019 } else if (n > 0 && n < viewnum) { /* rewind file */
1020 if (viewnum == 1 && feof(viewfp))
1021 return(&curview); /* just one view */
1022 if (fseek(viewfp, 0L, 0) == EOF) {
1023 perror(vval(VIEWFILE));
1024 quit(1);
1025 }
1026 copystruct(&curview, &stdview);
1027 viewnum = 0;
1028 }
1029 if (n < 0) { /* get next view */
1030 register int c = getc(viewfp);
1031 if (c == EOF)
1032 return((VIEW *)NULL); /* that's it */
1033 ungetc(c, viewfp);
1034 n = viewnum + 1;
1035 }
1036 while (n > viewnum) { /* scan to desired view */
1037 if (fgets(linebuf, sizeof(linebuf), viewfp) == NULL)
1038 return(viewnum==1 ? &curview : (VIEW *)NULL);
1039 if (isview(linebuf) && sscanview(&curview, linebuf) > 0)
1040 viewnum++;
1041 }
1042 return(&curview); /* return it */
1043 }
1044
1045
1046 int
1047 countviews() /* count views in view file */
1048 {
1049 int n;
1050
1051 if (getview(n=1) == NULL)
1052 return(0);
1053 while (getview(-1) != NULL)
1054 n++;
1055 return(n);
1056 }
1057
1058
1059 char *
1060 getexp(n) /* get exposure for nth frame */
1061 int n;
1062 {
1063 extern char *fskip();
1064 static char expval[32];
1065 static FILE *expfp = NULL;
1066 static long *exppos;
1067 static int curfrm;
1068 register char *cp;
1069
1070 if (n == 0) { /* signal to close file */
1071 if (expfp != NULL) {
1072 fclose(expfp);
1073 free((char *)exppos);
1074 expfp = NULL;
1075 }
1076 return(NULL);
1077 } else if (n > vint(END)) /* request past end (error?) */
1078 return(NULL);
1079 if (!vdef(EXPOSURE)) /* no setting (auto) */
1080 return(NULL);
1081 if (isflt(vval(EXPOSURE))) /* always the same */
1082 return(vval(EXPOSURE));
1083 if (expfp == NULL) { /* open exposure file */
1084 if ((expfp = fopen(vval(EXPOSURE), "r")) == NULL) {
1085 fprintf(stderr,
1086 "%s: cannot open exposure file \"%s\"\n",
1087 progname, vval(EXPOSURE));
1088 quit(1);
1089 }
1090 curfrm = vint(END) + 1; /* init lookup tab. */
1091 exppos = (long *)malloc(curfrm*sizeof(long *));
1092 if (exppos == NULL) {
1093 perror(progname);
1094 quit(1);
1095 }
1096 while (curfrm--)
1097 exppos[curfrm] = -1L;
1098 curfrm = 0;
1099 }
1100 /* find position in file */
1101 if (n-1 != curfrm && n != curfrm && exppos[n-1] >= 0 &&
1102 fseek(expfp, exppos[curfrm=n-1], 0) == EOF) {
1103 fprintf(stderr, "%s: seek error on exposure file\n", progname);
1104 quit(1);
1105 }
1106 while (n > curfrm) { /* read exposure */
1107 if (exppos[curfrm] < 0)
1108 exppos[curfrm] = ftell(expfp);
1109 if (fgets(expval, sizeof(expval), expfp) == NULL) {
1110 fprintf(stderr, "%s: too few exposures\n",
1111 vval(EXPOSURE));
1112 quit(1);
1113 }
1114 curfrm++;
1115 cp = fskip(expval); /* check format */
1116 if (cp != NULL)
1117 while (isspace(*cp))
1118 *cp++ = '\0';
1119 if (cp == NULL || *cp) {
1120 fprintf(stderr,
1121 "%s: exposure format error on line %d\n",
1122 vval(EXPOSURE), curfrm);
1123 quit(1);
1124 }
1125 }
1126 return(expval); /* return value */
1127 }
1128
1129
1130 struct pslot *
1131 findpslot(pid) /* find or allocate a process slot */
1132 int pid;
1133 {
1134 register struct pslot *psempty = NULL;
1135 register int i;
1136
1137 for (i = 0; i < npslots; i++) { /* look for match */
1138 if (pslot[i].pid == pid)
1139 return(pslot+i);
1140 if (psempty == NULL && pslot[i].pid == 0)
1141 psempty = pslot+i;
1142 }
1143 return(psempty); /* return emtpy slot (error if NULL) */
1144 }
1145
1146
1147 int
1148 donecom(ps, pn, status) /* clean up after finished process */
1149 PSERVER *ps;
1150 int pn;
1151 int status;
1152 {
1153 register PROC *pp;
1154 register struct pslot *psl;
1155
1156 pp = ps->proc + pn;
1157 if (pp->elen) { /* pass errors */
1158 if (ps->hostname[0])
1159 fprintf(stderr, "%s: ", ps->hostname);
1160 fprintf(stderr, "Error output from: %s\n", pp->com);
1161 fputs(pp->errs, stderr);
1162 fflush(stderr);
1163 if (ps->hostname[0])
1164 status = 1; /* because rsh doesn't return status */
1165 }
1166 lastpserver = NULL;
1167 psl = findpslot(pp->pid); /* check for bruncom() slot */
1168 if (psl->pid) {
1169 if (status) {
1170 if (psl->rcvf != NULL) /* attempt recovery */
1171 status = (*psl->rcvf)(psl->fout);
1172 if (status) {
1173 fprintf(stderr,
1174 "%s: error rendering frame %d\n",
1175 progname, psl->fout);
1176 quit(1);
1177 }
1178 lastpserver = ps;
1179 }
1180 psl->pid = 0; /* free process slot */
1181 } else if (status)
1182 lastpserver = ps;
1183 freestr(pp->com); /* free command string */
1184 return(status);
1185 }
1186
1187
1188 int
1189 serverdown() /* check status of last process server */
1190 {
1191 if (lastpserver == NULL || !lastpserver->hostname[0])
1192 return(0);
1193 if (pserverOK(lastpserver)) /* server still up? */
1194 return(0);
1195 delpserver(lastpserver); /* else delete it */
1196 if (pslist == NULL) {
1197 fprintf(stderr, "%s: all process servers are down\n",
1198 progname);
1199 quit(1);
1200 }
1201 return(1);
1202 }
1203
1204
1205 int
1206 bruncom(com, fout, rf) /* run a command in the background */
1207 char *com;
1208 int fout;
1209 int (*rf)();
1210 {
1211 int pid;
1212 register struct pslot *psl;
1213
1214 if (noaction) {
1215 if (!silent)
1216 printf("\t%s\n", com); /* echo command */
1217 return(0);
1218 }
1219 com = savestr(com); /* else start it when we can */
1220 while ((pid = startjob(NULL, com, donecom)) == -1)
1221 bwait(1);
1222 if (!silent) { /* echo command */
1223 PSERVER *ps;
1224 int psn = pid;
1225 ps = findjob(&psn);
1226 printf("\t%s\n", com);
1227 printf("\tProcess started on %s\n", phostname(ps));
1228 fflush(stdout);
1229 }
1230 psl = findpslot(pid); /* record info. in appropriate slot */
1231 psl->pid = pid;
1232 psl->fout = fout;
1233 psl->rcvf = rf;
1234 return(pid);
1235 }
1236
1237
1238 bwait(ncoms) /* wait for batch job(s) to finish */
1239 int ncoms;
1240 {
1241 int status;
1242
1243 if (noaction)
1244 return;
1245 while ((status = wait4job(NULL, -1)) != -1) {
1246 serverdown(); /* update server status */
1247 if (--ncoms == 0)
1248 break; /* done enough */
1249 }
1250 }
1251
1252
1253 int
1254 pruncom(com, ppins, maxcopies) /* run a command in parallel over network */
1255 char *com, *ppins;
1256 int maxcopies;
1257 {
1258 int retstatus = 0;
1259 int hostcopies;
1260 char buf[10240], *com1, *s;
1261 int status;
1262 int pfd;
1263 register int n;
1264 register PSERVER *ps;
1265
1266 if (!silent)
1267 printf("\t%s\n", com); /* echo command */
1268 if (noaction)
1269 return(0);
1270 fflush(stdout);
1271 /* start jobs on each server */
1272 for (ps = pslist; ps != NULL; ps = ps->next) {
1273 hostcopies = 0;
1274 if (maxcopies > 1 && ps->nprocs > 1 && ppins != NULL) {
1275 strcpy(com1=buf, com); /* build -PP command */
1276 sprintf(com1+(ppins-com), " -PP %s/%s.persist",
1277 vval(DIRECTORY), phostname(ps));
1278 strcat(com1, ppins);
1279 } else
1280 com1 = com;
1281 while (maxcopies > 0) {
1282 s = savestr(com1);
1283 if (startjob(ps, s, donecom) != -1) {
1284 sleep(20);
1285 hostcopies++;
1286 maxcopies--;
1287 } else {
1288 freestr(s);
1289 break;
1290 }
1291 }
1292 if (!silent && hostcopies) {
1293 if (hostcopies > 1)
1294 printf("\t%d duplicate processes", hostcopies);
1295 else
1296 printf("\tProcess");
1297 printf(" started on %s\n", phostname(ps));
1298 fflush(stdout);
1299 }
1300 }
1301 /* wait for jobs to finish */
1302 while ((status = wait4job(NULL, -1)) != -1)
1303 retstatus += status && !serverdown();
1304 /* terminate parallel rpict's */
1305 for (ps = pslist; ps != NULL; ps = ps->next) {
1306 sprintf(buf, "%s/%s.persist", vval(DIRECTORY), phostname(ps));
1307 if ((pfd = open(buf, O_RDONLY)) >= 0) {
1308 n = read(pfd, buf, sizeof(buf)-1); /* get PID */
1309 buf[n] = '\0';
1310 close(pfd);
1311 for (n = 0; buf[n] && !isspace(buf[n]); n++)
1312 ;
1313 /* terminate */
1314 sprintf(buf, "kill -ALRM %d", atoi(buf+n));
1315 wait4job(ps, startjob(ps, buf, NULL));
1316 }
1317 }
1318 return(retstatus);
1319 }
1320
1321
1322 runcom(cs) /* run a command locally and wait for it */
1323 char *cs;
1324 {
1325 if (!silent) /* echo it */
1326 printf("\t%s\n", cs);
1327 if (noaction)
1328 return(0);
1329 fflush(stdout); /* flush output and pass to shell */
1330 return(system(cs));
1331 }
1332
1333
1334 rmfile(fn) /* remove a file */
1335 char *fn;
1336 {
1337 if (!silent)
1338 #ifdef MSDOS
1339 printf("\tdel %s\n", fn);
1340 #else
1341 printf("\trm -f %s\n", fn);
1342 #endif
1343 if (noaction)
1344 return(0);
1345 return(unlink(fn));
1346 }
1347
1348
1349 badvalue(vc) /* report bad variable value and exit */
1350 int vc;
1351 {
1352 fprintf(stderr, "%s: bad value for variable '%s'\n",
1353 progname, vnam(vc));
1354 quit(1);
1355 }
1356
1357
1358 char *
1359 dirfile(df, path) /* separate path into directory and file */
1360 char *df;
1361 register char *path;
1362 {
1363 register int i;
1364 int psep;
1365
1366 for (i = 0, psep = -1; path[i]; i++)
1367 if (path[i] == '/')
1368 psep = i;
1369 if (df != NULL)
1370 if (psep == 0) {
1371 df[0] = '/';
1372 df[1] = '\0';
1373 } else if (psep > 0) {
1374 strncpy(df, path, psep);
1375 df[psep] = '\0';
1376 } else
1377 df[0] = '\0';
1378 return(path+psep+1);
1379 }