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/dvo/procks.c

    r19823 r20936  
    9595
    9696  /* data has been loaded, get ready to plot it */
    97   Yvec.Nelements = Xvec.Nelements = 3*Nrocks;
    98   ALLOCATE (Xvec.elements, float, Xvec.Nelements);
    99   ALLOCATE (Yvec.elements, float, Yvec.Nelements);
     97  SetVector (&Xvec, OPIHI_FLT, 3*Nrocks);
     98  SetVector (&Yvec, OPIHI_FLT, 3*Nrocks);
    10099 
    101100  /* project stars to screen display coords */
     
    106105      while (rocks[i].ra[j] < Rmin) rocks[i].ra[j] += 360.0;
    107106      while (rocks[i].ra[j] > Rmax) rocks[i].ra[j] -= 360.0;
    108       fRD_to_XY (&Xvec.elements[N], &Yvec.elements[N], rocks[i].ra[j], rocks[i].dec[j], &graphmode.coords);
     107      RD_to_XY (&Xvec.elements.Flt[N], &Yvec.elements.Flt[N], rocks[i].ra[j], rocks[i].dec[j], &graphmode.coords);
    109108      N ++;
    110109    }
     
    114113  graphmode.style = 2; /* set style to points */
    115114  graphmode.etype = 0; /* no errorbars */
    116   PlotVectorPair (kapa, N, Xvec.elements, Yvec.elements, &graphmode);
     115  PlotVectorPair (kapa, &Xvec, &Yvec, &graphmode);
    117116
    118117  /* now plot vectors between two extrema */
    119118  Yvec.Nelements = Xvec.Nelements = 2*Nrocks;
    120   REALLOCATE (Xvec.elements, float, Xvec.Nelements);
    121   REALLOCATE (Yvec.elements, float, Yvec.Nelements);
     119  REALLOCATE (Xvec.elements.Flt, opihi_flt, Xvec.Nelements);
     120  REALLOCATE (Yvec.elements.Flt, opihi_flt, Yvec.Nelements);
    122121 
    123122  /* project stars to screen display coords */
     
    136135    while (rocks[i].ra[N1] < Rmin) rocks[i].ra[N1] += 360.0;
    137136    while (rocks[i].ra[N1] > Rmax) rocks[i].ra[N1] -= 360.0;
    138     fRD_to_XY (&Xvec.elements[N], &Yvec.elements[N], rocks[i].ra[N0], rocks[i].dec[N0], &graphmode.coords);
     137    RD_to_XY (&Xvec.elements.Flt[N], &Yvec.elements.Flt[N], rocks[i].ra[N0], rocks[i].dec[N0], &graphmode.coords);
    139138    N ++;
    140     fRD_to_XY (&Xvec.elements[N], &Yvec.elements[N], rocks[i].ra[N1], rocks[i].dec[N1], &graphmode.coords);
     139    RD_to_XY (&Xvec.elements.Flt[N], &Yvec.elements.Flt[N], rocks[i].ra[N1], rocks[i].dec[N1], &graphmode.coords);
    141140    N ++;
    142141  }
     
    147146  graphmode.etype = 0; /* no errorbars */
    148147
    149   PlotVectorPair (kapa, N, Xvec.elements, Yvec.elements, &graphmode);
     148  PlotVectorPair (kapa, &Xvec, &Yvec, &graphmode);
    150149
    151   free (Xvec.elements);
    152   free (Yvec.elements);
     150  free (Xvec.elements.Flt);
     151  free (Yvec.elements.Flt);
    153152  return (TRUE);
    154153
Note: See TracChangeset for help on using the changeset viewer.