IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Dec 7, 2008, 3:31:01 PM (18 years ago)
Author:
eugene
Message:

big update from eam_branch_20081124 with updates to Opihi math

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Ohana/src/opihi/cmd.data/vgrid.c

    r9275 r20936  
    55  int i, Nx, Ny, Xb, Yb;
    66  float Xmin, Xmax, dX, Ymin, Ymax, dY;
    7   float *buf, *val, *x, *y, *z;
     7  float *buf;
     8  opihi_flt *val, *x, *y, *z;
    89  int *Nval;
    910  Buffer *bf;
     
    1213  if (argc != 11) {
    1314    gprint (GP_ERR, "USAGE: vgrid x y z buffer Xmin Xmax dX Ymin Ymax dY\n");
     15    gprint (GP_ERR, "  re-grid values from a triplet of vectors (x,y,z) into an image\n");
     16    gprint (GP_ERR, "  the vectors must be floating-point type\n");
    1417    return (FALSE);
    1518  }
     
    2225  if (vx[0].Nelements != vy[0].Nelements) return (FALSE);
    2326  if (vx[0].Nelements != vz[0].Nelements) return (FALSE);
     27
     28  REQUIRE_VECTOR_FLT (vx, FALSE);
     29  REQUIRE_VECTOR_FLT (vy, FALSE);
     30  REQUIRE_VECTOR_FLT (vz, FALSE);
    2431
    2532  Xmin = atof (argv[5]);
     
    3946  strcpy (bf[0].file, "(empty)");
    4047
    41   ALLOCATE (val, float, Nx*Ny);
    42   bzero (val, Nx*Ny*sizeof(float));
     48  ALLOCATE (val, opihi_flt, Nx*Ny);
     49  bzero (val, Nx*Ny*sizeof(opihi_flt));
    4350  ALLOCATE (Nval, int, Nx*Ny);
    4451  bzero (Nval, Nx*Ny*sizeof(int));
    4552
    46   x = vx[0].elements;
    47   y = vy[0].elements;
    48   z = vz[0].elements;
     53  x = vx[0].elements.Flt;
     54  y = vy[0].elements.Flt;
     55  z = vz[0].elements.Flt;
    4956  for (i = 0; i < vx[0].Nelements; i++, x++, y++, z++) {
    5057    Xb = (*x - Xmin) / dX;
Note: See TracChangeset for help on using the changeset viewer.