- Timestamp:
- Nov 30, 2008, 7:45:21 AM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branch_20081124/Ohana/src/opihi/dvo/paverage.c
r20839 r20857 10 10 double Radius, Rmin, Rmax, R, D; 11 11 unsigned IDclip, IDchoice, LimExclude; 12 opihi_flt *Xvec, *Yvec, *Zvec;12 float *Xvec, *Yvec, *Zvec; 13 13 void *Signal; 14 14 … … 38 38 photcode = GetPhotcodebyName (argv[Narg]); 39 39 if (!photcode) { 40 fprintf (stderr, "unknown photcode %s\n", argv[Narg]);41 return (FALSE);40 fprintf (stderr, "unknown photcode %s\n", argv[Narg]); 41 return (FALSE); 42 42 } 43 43 remove_argument (Narg, &argc, argv); 44 44 Nsec = GetPhotcodeNsec (photcode[0].code); 45 45 if (Nsec == -1) { 46 fprintf (stderr, "photcode %s is not an AVERAGE photcode\n", argv[Narg]);47 return (FALSE);46 fprintf (stderr, "photcode %s is not an AVERAGE photcode\n", argv[Narg]); 47 return (FALSE); 48 48 } 49 49 } … … 96 96 Npts = 0; 97 97 NPTS = 1000; 98 ALLOCATE (Xvec, opihi_flt, NPTS);99 ALLOCATE (Yvec, opihi_flt, NPTS);100 ALLOCATE (Zvec, opihi_flt, NPTS);98 ALLOCATE (Xvec, float, NPTS); 99 ALLOCATE (Yvec, float, NPTS); 100 ALLOCATE (Zvec, float, NPTS); 101 101 102 102 // prepare to handle interrupt signals … … 133 133 R = average[i].R; 134 134 D = average[i].D; 135 status = RD_to_XY (&Xvec[Npts], &Yvec[Npts], R, D, &graphmode.coords);135 status = fRD_to_XY (&Xvec[Npts], &Yvec[Npts], R, D, &graphmode.coords); 136 136 if (!status) continue; 137 137 Npts ++; 138 138 139 139 if (Npts == NPTS - 1) { 140 NPTS += 1000;141 REALLOCATE (Xvec, opihi_flt, NPTS);142 REALLOCATE (Yvec, opihi_flt, NPTS);143 REALLOCATE (Zvec, opihi_flt, NPTS);140 NPTS += 1000; 141 REALLOCATE (Xvec, float, NPTS); 142 REALLOCATE (Yvec, float, NPTS); 143 REALLOCATE (Zvec, float, NPTS); 144 144 } 145 145 if ((Npts > NCHUNK) || (Nloaded >= 25)) { 146 PlotVectorTriplet (kapa, Npts, Xvec, Yvec, Zvec, &graphmode); 147 Npts = 0; 148 Nloaded = 0; 146 KapaPrepPlot (kapa, Npts, graphmode); 147 KapaPlotVector (kapa, Npts, Xvec, "x"); 148 KapaPlotVector (kapa, Npts, Yvec, "y"); 149 KapaPlotVector (kapa, Npts, Zvec, "z"); 150 Npts = 0; 151 Nloaded = 0; 149 152 } 150 153 } … … 152 155 dvo_catalog_free (&catalog); 153 156 } 154 if (Npts > 0) PlotVectorTriplet (kapa, Npts, Xvec, Yvec, Zvec, &graphmode); 157 if (Npts > 0) { 158 KapaPrepPlot (kapa, Npts, graphmode); 159 KapaPlotVector (kapa, Npts, Xvec, "x"); 160 KapaPlotVector (kapa, Npts, Yvec, "y"); 161 KapaPlotVector (kapa, Npts, Zvec, "z"); 162 } 155 163 156 164 free (Xvec);
Note:
See TracChangeset
for help on using the changeset viewer.
