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/mana/starcontour.c

    r7917 r20936  
    2222  N = 0;
    2323  Npts = 100;
    24   REALLOCATE (vecx[0].elements, float, MAX (Npts, 1));
    25   REALLOCATE (vecy[0].elements, float, MAX (Npts, 1));
     24  ResetVector (vecx, OPIHI_FLT, MAX (Npts, 1));
     25  ResetVector (vecy, OPIHI_FLT, MAX (Npts, 1));
    2626
    2727  Nx = buf[0].matrix.Naxis[0];
    … …  
    4141        min = TRUE;
    4242        xmin = x + (zo - zt)/(v[x + 1 + y*Nx] - zt);
    43         vecx[0].elements[N] = xmin;
    44         vecy[0].elements[N] = y;
     43        vecx[0].elements.Flt[N] = xmin;
     44        vecy[0].elements.Flt[N] = y;
    4545        N ++;
    4646        if (N >= Npts) {
    4747          Npts += 100;
    48           REALLOCATE (vecx[0].elements, float, MAX (Npts, 1));
    49           REALLOCATE (vecy[0].elements, float, MAX (Npts, 1));
     48          REALLOCATE (vecx[0].elements.Flt, opihi_flt, MAX (Npts, 1));
     49          REALLOCATE (vecy[0].elements.Flt, opihi_flt, MAX (Npts, 1));
    5050        }
    5151      }
    … …  
    5959        max = TRUE;
    6060        xmax = x - (zo - zt)/(v[x - 1 + y*Nx] - zt);
    61         vecx[0].elements[N] = xmax;
    62         vecy[0].elements[N] = y;
     61        vecx[0].elements.Flt[N] = xmax;
     62        vecy[0].elements.Flt[N] = y;
    6363        N ++;
    6464        if (N >= Npts) {
    6565          Npts += 100;
    66           REALLOCATE (vecx[0].elements, float, MAX (Npts, 1));
    67           REALLOCATE (vecy[0].elements, float, MAX (Npts, 1));
     66          REALLOCATE (vecx[0].elements.Flt, opihi_flt, MAX (Npts, 1));
     67          REALLOCATE (vecy[0].elements.Flt, opihi_flt, MAX (Npts, 1));
    6868        }
    6969      }
    … …  
    8686        min = TRUE;
    8787        xmin = x + (zo - zt)/(v[x + 1 + y*Nx] - zt);
    88         vecx[0].elements[N] = xmin;
    89         vecy[0].elements[N] = y;
     88        vecx[0].elements.Flt[N] = xmin;
     89        vecy[0].elements.Flt[N] = y;
    9090        N ++;
    9191        if (N >= Npts) {
    9292          Npts += 100;
    93           REALLOCATE (vecx[0].elements, float, MAX (Npts, 1));
    94           REALLOCATE (vecy[0].elements, float, MAX (Npts, 1));
     93          REALLOCATE (vecx[0].elements.Flt, opihi_flt, MAX (Npts, 1));
     94          REALLOCATE (vecy[0].elements.Flt, opihi_flt, MAX (Npts, 1));
    9595        }
    9696      }
    … …  
    104104        max = TRUE;
    105105        xmax = x - (zo - zt)/(v[x - 1 + y*Nx] - zt);
    106         vecx[0].elements[N] = xmax;
    107         vecy[0].elements[N] = y;
     106        vecx[0].elements.Flt[N] = xmax;
     107        vecy[0].elements.Flt[N] = y;
    108108        N ++;
    109109        if (N >= Npts) {
    110110          Npts += 100;
    111           REALLOCATE (vecx[0].elements, float, MAX (Npts, 1));
    112           REALLOCATE (vecy[0].elements, float, MAX (Npts, 1));
     111          REALLOCATE (vecx[0].elements.Flt, opihi_flt, MAX (Npts, 1));
     112          REALLOCATE (vecy[0].elements.Flt, opihi_flt, MAX (Npts, 1));
    113113        }
    114114      }
Note: See TracChangeset for help on using the changeset viewer.