Return to contour.c CVS log

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

Annotation of researchv10no/cmd/view2d/contour.c, revision 1.1

1.1     ! root        1: /* undocumented feature:
        !             2:    -clx,x,x,x  sets contour levels at x, x, x, and x.   up to 8 levels allowed.
        !             3: */
        !             4: #include <stdio.h>
        !             5: #include <math.h>
        !             6: #include "view2d.h"
        !             7: 
        !             8: extern Rd2d rd;
        !             9: char *progname;
        !            10: short timewarp;
        !            11: double ts, te;
        !            12: 
        !            13: main(argc,argv)
        !            14:   int argc;
        !            15:   char **argv;
        !            16: {
        !            17:   float *p, *pict;
        !            18:   float u2, v2, pow2();
        !            19:   short *v, *q;
        !            20:   char *map;
        !            21:   int prange;
        !            22:   int ncontour = 6;
        !            23:   int i, j, r;
        !            24:   short nx, ny, fru, frv;
        !            25:   int nf;
        !            26:   double frtime;
        !            27:   int zero;
        !            28:   int g2nx, g2ny;
        !            29:   float g2s[6], g2d[4], g2c[300];
        !            30:   float NaN;
        !            31:   int explcont = 0;
        !            32:   timewarp = 0;
        !            33:   progname=argv[0];
        !            34: 
        !            35:   for(argc--, argv++; *argv && (**argv == '-' ); argv++){
        !            36:     switch(argv[0][1]){
        !            37:      case 'c':
        !            38:        if(argv[0][2]=='l'){
        !            39:          ncontour = sscanf(&(argv[0][3]),"%f,%f,%f,%f,%f,%f,%f,%f",
        !            40:             g2c+0,g2c+1,g2c+2,g2c+3,g2c+4,g2c+5,g2c+6,g2c+7,g2c+8);
        !            41:          explcont = 1;
        !            42:          break;
        !            43:        }
        !            44:        ncontour = atoi(&argv[0][2]);
        !            45:        break;
        !            46:      case 't':
        !            47:        i = sscanf(&argv[0][2], "%hd, %hd", &ts, &te);
        !            48:        if(i==0) error("bad TS,TE");
        !            49:        timewarp = i;
        !            50:        break;
        !            51:      case 'T':
        !            52:        break;
        !            53:      default:
        !            54:      err:
        !            55:        fprintf(stderr,"unrecognized option: %s\n",argv[0]);
        !            56:        exit(1);
        !            57:      }
        !            58:    }
        !            59: 
        !            60:   rd2dh(0,&nx,&ny);
        !            61:   v = (short *)Malloc(nx*ny*sizeof(short));
        !            62:   pict = (float *)Malloc(nx*ny*sizeof(float));
        !            63: 
        !            64:   g2open_();
        !            65:   g2d[0] = 0;
        !            66:   g2d[1] = 0;
        !            67:   g2d[2] = nx-1;
        !            68:   g2d[3] = ny-1;
        !            69:   i = 4;
        !            70:   g2sce_(g2s,g2d,&i);
        !            71:   g2nx = nx;
        !            72:   g2ny = ny;
        !            73:   fprintf(stderr,"contours at:\n");
        !            74:   for( i=1; i<=ncontour; i++){
        !            75:     if(!explcont)
        !            76:       g2c[i-1] = rd.fmin + i*(rd.fmax-rd.fmin)/(ncontour+1);
        !            77:     fprintf(stderr," %g\n",g2c[i-1]);
        !            78:   }
        !            79: 
        !            80:   /* for a not-a-number flag, choose a value below all the data */
        !            81:   NaN = rd.fmin;
        !            82:   if(NaN > -1) NaN = -1;
        !            83:   NaN *= 2;
        !            84: 
        !            85:   u2 = rd.u;
        !            86:   v2 = pow2(rd.v);
        !            87:   while(rd2di(&frtime,v)){
        !            88:     p = pict;
        !            89:     q = v;
        !            90:     for(i=0; i<nx*ny; i++){
        !            91:       r = *q;
        !            92:       if( r < -BIG ){ *p = NaN; }
        !            93:        else{
        !            94:         *p = (r-u2)*v2;
        !            95:        }
        !            96:       p++; q++;
        !            97:     }
        !            98:     g2ff_();
        !            99:     g2box_(g2s,g2d);
        !           100:     g2co_(g2s,g2d,&g2nx,&g2ny,pict,&ncontour,g2c,&NaN);
        !           101:   }
        !           102:   g2clos_();
        !           103: }

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