ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/px/pinterp.c
(Generate patch)

Comparing ray/src/px/pinterp.c (file contents):
Revision 2.27 by greg, Mon Oct 16 11:40:16 1995 UTC vs.
Revision 2.37 by schorsch, Sun Jul 27 22:12:03 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1995 Regents of the University of California */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   * Interpolate and extrapolate pictures with different view parameters.
6   *
7   *      Greg Ward       09Dec89
8   */
9  
10 < #include "standard.h"
10 > #include "copyright.h"
11  
12   #include <ctype.h>
13 + #include <string.h>
14  
15 + #include "standard.h"
16 + #include "rtprocess.h" /* Windows: must come before color.h */
17   #include "view.h"
18 –
18   #include "color.h"
19  
21 – #include "resolu.h"
22 –
20   #define LOG2            0.69314718055994530942
21  
22   #define pscan(y)        (ourpict+(y)*hresolu)
# Line 30 | Line 27 | static char SCCSid[] = "$SunId$ LBL";
27   #define averaging       (ourweigh != NULL)
28   #define blurring        (ourbpict != NULL)
29   #define usematrix       (hasmatrix & !averaging)
30 < #define zisnorm         (!usematrix | ourview.type != VT_PER)
30 > #define zisnorm         ((!usematrix) | (ourview.type != VT_PER))
31  
32   #define MAXWT           1000.           /* maximum pixel weight (averaging) */
33  
# Line 39 | Line 36 | static char SCCSid[] = "$SunId$ LBL";
36  
37   #define PACKSIZ         256             /* max. calculation packet size */
38  
39 < #define RTCOM           "rtrace -h- -ovl -fff "
39 > #define RTCOM           "rtrace -h- -ovl -fff -ld- -i- -I- "
40  
41   #define ABS(x)          ((x)>0?(x):-(x))
42  
# Line 89 | Line 86 | double theirexp;                       /* input picture exposure */
86   MAT4    theirs2ours;                    /* transformation matrix */
87   int     hasmatrix = 0;                  /* has transformation matrix */
88  
89 < int     PDesc[3] = {-1,-1,-1};          /* rtrace process descriptor */
93 < #define childpid        (PDesc[2])
89 > static SUBPROC PDesc = SP_INACTIVE; /* rtrace process descriptor */
90   unsigned short  queue[PACKSIZ][2];      /* pending pixels */
91   int     packsiz;                        /* actual packet size */
92   int     queuesiz = 0;                   /* number of pixels pending */
# Line 213 | Line 209 | char   *argv[];
209                          if (argv[an][2] != 'f')
210                                  goto badopt;
211                          check(3,"s");
212 <                        gotvfile = viewfile(argv[++an], &ourview, 0, 0);
212 >                        gotvfile = viewfile(argv[++an], &ourview, NULL);
213                          if (gotvfile < 0)
214                                  syserror(argv[an]);
215                          else if (gotvfile == 0) {
# Line 236 | Line 232 | char   *argv[];
232                  fillo &= ~F_BACK;
233          if (doavg < 0)
234                  doavg = (argc-an) > 2;
235 <        if (expcomp != NULL)
236 <                if (expcomp[0] == '+' | expcomp[0] == '-') {
235 >        if (expcomp != NULL) {
236 >                if ((expcomp[0] == '+') | (expcomp[0] == '-')) {
237                          expadj = atof(expcomp) + (expcomp[0]=='+' ? .5 : -.5);
238                          if (doavg | doblur)
239                                  rexpadj = pow(2.0, atof(expcomp));
# Line 251 | Line 247 | char   *argv[];
247                          if (!(doavg | doblur))
248                                  rexpadj = pow(2.0, (double)expadj);
249                  }
250 +        }
251                                                  /* set view */
252 <        if (nextview(doblur ? stdin : NULL) == EOF) {
252 >        if (nextview(doblur ? stdin : (FILE *)NULL) == EOF) {
253                  fprintf(stderr, "%s: no view on standard input!\n",
254                                  progname);
255                  exit(1);
# Line 262 | Line 259 | char   *argv[];
259          if (doavg) {
260                  ourspict = (COLOR *)bmalloc(hresolu*vresolu*sizeof(COLOR));
261                  ourweigh = (float *)bmalloc(hresolu*vresolu*sizeof(float));
262 <                if (ourspict == NULL | ourweigh == NULL)
262 >                if ((ourspict == NULL) | (ourweigh == NULL))
263                          syserror(progname);
264          } else {
265                  ourpict = (COLR *)bmalloc(hresolu*vresolu*sizeof(COLR));
# Line 279 | Line 276 | char   *argv[];
276                  syserror(progname);
277                                                          /* new header */
278          newheader("RADIANCE", stdout);
279 +        fputnow(stdout);
280                                                          /* run pictures */
281          do {
282 <                bzero((char *)ourzbuf, hresolu*vresolu*sizeof(float));
282 >                memset((char *)ourzbuf, '\0', hresolu*vresolu*sizeof(float));
283                  for (i = an; i < argc; i += 2)
284                          addpicture(argv[i], argv[i+1]);
285                  if (fillo&F_BACK)                       /* fill in spaces */
# Line 301 | Line 299 | char   *argv[];
299                  fprintview(&avgview, stdout);
300                  putc('\n', stdout);
301          }
302 <        if (pixaspect < .99 | pixaspect > 1.01)
302 >        if ((pixaspect < .99) | (pixaspect > 1.01))
303                  fputaspect(pixaspect, stdout);
304          if (ourexp > 0)
305                  ourexp *= rexpadj;
306          else
307                  ourexp = rexpadj;
308 <        if (ourexp < .995 | ourexp > 1.005)
308 >        if ((ourexp < .995) | (ourexp > 1.005))
309                  fputexpos(ourexp, stdout);
310          if (strcmp(ourfmt, PICFMT))             /* print format if known */
311                  fputformat(ourfmt, stdout);
# Line 328 | Line 326 | userr:
326   }
327  
328  
329 + int
330   headline(s)                             /* process header string */
331   char    *s;
332   {
333          char    fmt[32];
334  
335          if (isheadid(s))
336 <                return;
336 >                return(0);
337          if (formatval(fmt, s)) {
338                  if (globmatch(ourfmt, fmt)) {
339                          wrongformat = 0;
340                          strcpy(ourfmt, fmt);
341                  } else
342                          wrongformat = 1;
343 <                return;
343 >                return(0);
344          }
345          if (nvavg < 2) {
346                  putc('\t', stdout);
# Line 349 | Line 348 | char   *s;
348          }
349          if (isexpos(s)) {
350                  theirexp *= exposval(s);
351 <                return;
351 >                return(0);
352          }
353          if (isview(s) && sscanview(&theirview, s) > 0)
354                  gotview++;
355 +        return(0);
356   }
357  
358  
# Line 375 | Line 375 | FILE   *fp;
375                  exit(1);
376          }
377          if (!nvavg) {                   /* first view */
378 <                copystruct(&avgview, &ourview);
378 >                avgview = ourview;
379                  return(nvavg++);
380          }
381                                          /* add to average view */
# Line 414 | Line 414 | compavgview()                          /* compute average view */
414          avgview.vfore *= f;
415          avgview.vaft *= f;
416          if (setview(&avgview) != NULL)          /* in case of emergency... */
417 <                copystruct(&avgview, &ourview);
417 >                avgview = ourview;
418          pixaspect = viewaspect(&avgview) * hresolu / vresolu;
419   }
420  
# Line 435 | Line 435 | char   *pfile, *zspec;
435                  syserror(pfile);
436                                          /* get header with exposure and view */
437          theirexp = 1.0;
438 <        copystruct(&theirview, &stdview);
438 >        theirview = stdview;
439          gotview = 0;
440          if (nvavg < 2)
441                  printf("%s:\n", pfile);
# Line 448 | Line 448 | char   *pfile, *zspec;
448                  ourexp = theirexp;
449          else if (ABS(theirexp-ourexp) > .01*ourexp)
450                  fprintf(stderr, "%s: different exposure (warning)\n", pfile);
451 <        if (err = setview(&theirview)) {
451 >        if ( (err = setview(&theirview)) ) {
452                  fprintf(stderr, "%s: %s\n", pfile, err);
453                  exit(1);
454          }
# Line 471 | Line 471 | char   *pfile, *zspec;
471          if (xlim == NULL)
472                  syserror(progname);
473          if (!getperim(xlim, &ylim, zin, zfd)) { /* overlapping area? */
474 <                free((char *)zin);
475 <                free((char *)xlim);
474 >                free((void *)zin);
475 >                free((void *)xlim);
476                  if (zfd != -1)
477                          close(zfd);
478                  fclose(pfp);
# Line 482 | Line 482 | char   *pfile, *zspec;
482          scanin = (COLR *)malloc(scanlen(&tresolu)*sizeof(COLR));
483          plast = (struct position *)calloc(scanlen(&tresolu),
484                          sizeof(struct position));
485 <        if (scanin == NULL | plast == NULL)
485 >        if ((scanin == NULL) | (plast == NULL))
486                  syserror(progname);
487                                          /* skip to starting point */
488          for (y = 0; y < ylim.min; y++)
# Line 491 | Line 491 | char   *pfile, *zspec;
491                          exit(1);
492                  }
493          if (zfd != -1 && lseek(zfd,
494 <                        (long)ylim.min*scanlen(&tresolu)*sizeof(float), 0) < 0)
494 >                        (off_t)ylim.min*scanlen(&tresolu)*sizeof(float), 0) < 0)
495                  syserror(zspec);
496                                          /* load image */
497          for (y = ylim.min; y <= ylim.max; y++) {
# Line 506 | Line 506 | char   *pfile, *zspec;
506                  addscanline(xlim+y, y, scanin, zin, plast);
507          }
508                                          /* clean up */
509 <        free((char *)xlim);
510 <        free((char *)scanin);
511 <        free((char *)zin);
512 <        free((char *)plast);
509 >        free((void *)xlim);
510 >        free((void *)scanin);
511 >        free((void *)zin);
512 >        free((void *)plast);
513          fclose(pfp);
514          if (zfd != -1)
515                  close(zfd);
# Line 522 | Line 522 | register VIEW  *vw1, *vw2;
522   {
523          double  m4t[4][4];
524  
525 <        if (vw1->type != VT_PER & vw1->type != VT_PAR)
525 >        if ((vw1->type != VT_PER) & (vw1->type != VT_PAR))
526                  return(0);
527 <        if (vw2->type != VT_PER & vw2->type != VT_PAR)
527 >        if ((vw2->type != VT_PER) & (vw2->type != VT_PAR))
528                  return(0);
529          setident4(xfmat);
530          xfmat[0][0] = vw1->hvec[0];
# Line 608 | Line 608 | double z;
608                  s1x = p1->x - p0->x;
609                  s1y = p1->y - p0->y;
610                  l1 = ABS(s1x);
611 <                if (p1isy = (ABS(s1y) > l1))
611 >                if ( (p1isy = (ABS(s1y) > l1)) )
612                          l1 = ABS(s1y);
613                  else if (l1 < 1)
614                          l1 = 1;
# Line 641 | Line 641 | double z;
641                  y1 = p0->y + c1*s1y/l1;
642                  for (c2 = l2; c2-- > 0; ) {
643                          x = x1 + c2*s2x/l2;
644 <                        if (x < 0 | x >= hresolu)
644 >                        if ((x < 0) | (x >= hresolu))
645                                  continue;
646                          y = y1 + c2*s2y/l2;
647 <                        if (y < 0 | y >= vresolu)
647 >                        if ((y < 0) | (y >= vresolu))
648                                  continue;
649                          if (averaging) {
650                                  if (zscan(y)[x] <= 0 || zscan(y)[x]-z
# Line 672 | Line 672 | register FVECT pos;
672   {
673          FVECT   pt, tdir, odir;
674          double  d;
675 –        register int    i;
675          
676          if (pos[2] <= 0)                /* empty pixel */
677                  return(0);
678          if (usematrix) {
679                  pos[0] += theirview.hoff - .5;
680                  pos[1] += theirview.voff - .5;
681 <                if (normdist & theirview.type == VT_PER)
681 >                if (normdist & (theirview.type == VT_PER))
682                          d = sqrt(1. + pos[0]*pos[0]*theirview.hn2
683                                          + pos[1]*pos[1]*theirview.vn2);
684                  else
# Line 703 | Line 702 | register FVECT pos;
702          } else {
703                  if (viewray(pt, tdir, &theirview, pos[0], pos[1]) < -FTINY)
704                          return(0);
705 <                if (!normdist & theirview.type == VT_PER)       /* adjust */
705 >                if ((!normdist) & (theirview.type == VT_PER))   /* adjust */
706                          pos[2] *= sqrt(1. + pos[0]*pos[0]*theirview.hn2
707                                          + pos[1]*pos[1]*theirview.vn2);
708                  pt[0] += tdir[0]*pos[2];
# Line 713 | Line 712 | register FVECT pos;
712                  if (pos[2] <= 0)
713                          return(0);
714          }
715 <        if (pos[0] < 0 | pos[0] >= 1-FTINY | pos[1] < 0 | pos[1] >= 1-FTINY)
715 >        if ((pos[0] < 0) | (pos[0] >= 1-FTINY) | (pos[1] < 0) | (pos[1] >= 1-FTINY))
716                  return(0);
717          if (!averaging)
718                  return(1);
719 <        if (ourview.type == VT_PAR)             /* compute our direction */
720 <                VCOPY(odir, ourview.vdir);
721 <        else
722 <                for (i = 0; i < 3; i++)
723 <                        odir[i] = (pt[i] - ourview.vp[i])/pos[2];
724 <        d = DOT(odir,tdir);                     /* compute pixel weight */
725 <        if (d >= 1.-1./MAXWT/MAXWT)
719 >                                                /* compute pixel weight */
720 >        if (ourview.type == VT_PAR) {
721 >                d = DOT(ourview.vdir,tdir);
722 >                d = 1. - d*d;
723 >        } else {
724 >                VSUB(odir, pt, ourview.vp);
725 >                d = DOT(odir,tdir);
726 >                d = 1. - d*d/DOT(odir,odir);
727 >        }
728 >        if (d <= 1./MAXWT/MAXWT)
729                  return(MAXWT);          /* clip to maximum weight */
730 <        return(1./sqrt(1.-d));
730 >        return(1./sqrt(d));
731   }
732  
733  
# Line 756 | Line 758 | int    zfd;
758          yl->min = 32000; yl->max = 0;           /* search for points on image */
759          for (y = step - 1; y < numscans(&tresolu); y += step) {
760                  if (zfd != -1) {
761 <                        if (lseek(zfd, (long)y*scanlen(&tresolu)*sizeof(float),
761 >                        if (lseek(zfd, (off_t)y*scanlen(&tresolu)*sizeof(float),
762                                          0) < 0)
763                                  syserror("lseek");
764                          if (read(zfd, (char *)zline,
# Line 802 | Line 804 | int    zfd;
804                                  xl[y-1].max = xl[y-step].max;
805                  }
806                  for (x = 2; x < step; x++)
807 <                        copystruct(xl+y-x, xl+y-1);
807 >                        *(xl+y-x) = *(xl+y-1);
808          }
809          if (yl->max >= numscans(&tresolu))
810                  yl->max = numscans(&tresolu) - 1;
811          y -= step;
812          for (x = numscans(&tresolu) - 1; x > y; x--)    /* fill bottom rows */
813 <                copystruct(xl+x, xl+y);
813 >                *(xl+x) = *(xl+y);
814          return(yl->max >= yl->min);
815   }
816  
# Line 917 | Line 919 | int    samp;
919                                  xback = -2;
920                          }
921          }
922 <        free((char *)yback);
922 >        free((void *)yback);
923   }
924  
925  
# Line 938 | Line 940 | int    (*fill)();
940   clipaft()                       /* perform aft clipping as indicated */
941   {
942          register int    x, y;
943 <        int     adjtest = ourview.type == VT_PER & zisnorm;
943 >        int     adjtest = (ourview.type == VT_PER) & zisnorm;
944          double  tstdist;
945          double  yzn2, vx;
946  
# Line 961 | Line 963 | clipaft()                      /* perform aft clipping as indicated */
963                                                  continue;
964                                  }
965                                  if (averaging)
966 <                                        bzero(sscan(y)[x], sizeof(COLOR));
966 >                                        memset(sscan(y)[x], '\0', sizeof(COLOR));
967                                  else
968 <                                        bzero(pscan(y)[x], sizeof(COLR));
968 >                                        memset(pscan(y)[x], '\0', sizeof(COLR));
969                                  zscan(y)[x] = 0.0;
970                          }
971          }
# Line 1046 | Line 1048 | writedistance(fname)                   /* write out z file (alters buf
1048   char    *fname;
1049   {
1050          int     donorm = normdist & !zisnorm ? 1 :
1051 <                        ourview.type == VT_PER & !normdist & zisnorm ? -1 : 0;
1051 >                        (ourview.type == VT_PER) & !normdist & zisnorm ? -1 : 0;
1052          int     fd;
1053          int     y;
1054  
# Line 1095 | Line 1097 | char   *prog, *args;
1097          int     rval;
1098          register char   **wp, *cp;
1099  
1100 <        if (childpid != -1) {
1100 >        if (PDesc.running) {
1101                  fprintf(stderr, "%s: too many calculations\n", progname);
1102                  exit(1);
1103          }
# Line 1114 | Line 1116 | char   *prog, *args;
1116          }
1117          *wp = NULL;
1118                                                  /* start process */
1119 <        if ((rval = open_process(PDesc, argv)) < 0)
1119 >        if ((rval = open_process(&PDesc, argv)) < 0)
1120                  syserror(progname);
1121          if (rval == 0) {
1122                  fprintf(stderr, "%s: command not found\n", argv[0]);
# Line 1129 | Line 1131 | char   *prog, *args;
1131  
1132   caldone()                               /* done with calculation */
1133   {
1134 <        if (childpid == -1)
1134 >        if (!PDesc.running)
1135                  return;
1136          clearqueue();
1137 <        close_process(PDesc);
1136 <        childpid = -1;
1137 >        close_process(&PDesc);
1138   }
1139  
1140  
# Line 1168 | Line 1169 | clearqueue()                           /* process queue */
1169                  *fbp++ = dir[0]; *fbp++ = dir[1]; *fbp++ = dir[2];
1170          }
1171                                          /* mark end and get results */
1172 <        bzero((char *)fbp, 6*sizeof(float));
1173 <        if (process(PDesc, fbuf, fbuf, 4*sizeof(float)*queuesiz,
1172 >        memset((char *)fbp, '\0', 6*sizeof(float));
1173 >        if (process(&PDesc, (char *)fbuf, (char *)fbuf,
1174 >                        4*sizeof(float)*(queuesiz+1),
1175                          6*sizeof(float)*(queuesiz+1)) !=
1176 <                        4*sizeof(float)*queuesiz) {
1176 >                        4*sizeof(float)*(queuesiz+1)) {
1177                  fprintf(stderr, "%s: error reading from rtrace process\n",
1178                                  progname);
1179                  exit(1);

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines