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.32 by gwlarson, Fri Jul 16 16:11:31 1999 UTC vs.
Revision 2.36 by schorsch, Mon Jul 21 22:30:18 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1999 Silicon Graphics, Inc. */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ SGI";
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 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)
235 >        if (expcomp != NULL) {
236                  if (expcomp[0] == '+' | expcomp[0] == '-') {
237                          expadj = atof(expcomp) + (expcomp[0]=='+' ? .5 : -.5);
238                          if (doavg | doblur)
# 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 : (FILE *)NULL) == EOF) {
253                  fprintf(stderr, "%s: no view on standard input!\n",
# 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 377 | 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 416 | 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 437 | 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 473 | 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 493 | 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 508 | 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 760 | 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 806 | 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 921 | Line 919 | int    samp;
919                                  xback = -2;
920                          }
921          }
922 <        free((char *)yback);
922 >        free((void *)yback);
923   }
924  
925  
# Line 965 | 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 1099 | 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 1118 | 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 1133 | 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);
1140 <        childpid = -1;
1137 >        close_process(&PDesc);
1138   }
1139  
1140  
# Line 1172 | 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+1),
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+1)) {
1177                  fprintf(stderr, "%s: error reading from rtrace process\n",

Diff Legend

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