Return to moviefil.c CVS log

Up to [Research Unix] / researchv10no / cmd / view2d

Annotation of researchv10no/cmd/view2d/moviefil.c, revision 1.1.1.1

1.1       root        1: /* From alice!doyle Fri Feb  1 17:19 EST 1985 */
                      2: /*  modified slightly by ehg */
                      3: /*  based on mvefil.f by rksmith */
                      4: #include <stdio.h>
                      5: #include <math.h>
                      6: #include "view2d.h"
                      7: 
                      8: #define max(a,b)               (a<b ? b : a)
                      9: #define min(a,b)               (a>b ? b : a)
                     10: #define abs(x)                 (x>=0 ? x : -(x))
                     11: #define inf 1.e30
                     12: 
                     13: moviefil(fd, nx, ny, time, outside, fa)
                     14: int fd, nx, ny;
                     15: float time, outside, fa[];
                     16: {
                     17:   int i;
                     18:   int npts;
                     19:   float f, v2;
                     20:   short u, v;
                     21:   float fmin, fmax, grange, gmax, gmin;
                     22:   double tim;
                     23:   short *data, *d;
                     24:   char *malloc();
                     25: 
                     26:     tim = time;
                     27:     npts = nx*ny;
                     28:     data = (short *)malloc((unsigned)npts*sizeof(short));
                     29:     if(data==NULL){fprintf(stderr,"out of memory!\n");exit(9);}
                     30: 
                     31:     fmin = inf;
                     32:     fmax = -inf;
                     33:     for (i=0;i<npts;i++) {
                     34:       if(fa[i]>outside){
                     35:        fmin = min(fa[i], fmin);
                     36:        fmax = max(fa[i], fmax);
                     37:       }
                     38:     }
                     39:        
                     40:     gmax = (fmax<0)?0:fmax;
                     41:     gmin = (fmin>0)?0:fmin;
                     42:     grange = gmax - gmin;
                     43:     if(grange<=0){ v = 0; u=0; }
                     44:     else{
                     45:       f = 2*BIG-2;
                     46:       v = log(grange)/log(2.)-14;
                     47:       v2 = pow(2.,(double)v);
                     48:       while( 2*grange < f*v2 ){ v--; v2 /= 2; }
                     49:       f = -(gmin+gmax)/(2*v2);
                     50:       u = f;
                     51:       if( (f<-BIG)||(f>BIG) ){
                     52:         fprintf(stderr,"fmin=%g fmax=%g u=%d v=%d\n",fmin,fmax,u,v);
                     53:         fprintf(stderr,"u out of bounds f=%g",f);
                     54:        exit(1);
                     55:       }
                     56:     }
                     57: 
                     58:     for(d = data, i = 0; i < npts; i++, d++){
                     59:       if(fa[i]<=outside){
                     60:         *d = -BIG-1;  /* out of bounds */
                     61:       }else{
                     62:         *d = u + fa[i]/v2;
                     63:       }
                     64:     }
                     65:     view2d(fd,nx,ny,tim,u,v,0,0,0,data);
                     66:     free((char *)data);
                     67: }
                     68: 
                     69: 
                     70: fmovie_(nx,ny,time,outside,f)
                     71:   int *nx, *ny;
                     72:   float *time, *outside, *f;
                     73: {
                     74:   moviefil(1,*nx,*ny,*time,*outside,f);
                     75: }

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.