IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Nov 25, 2008, 4:51:39 PM (18 years ago)
Author:
eugene
Message:

convert opihi vectors to double type, add option for int type as well

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/eam_branch_20081124/Ohana/src/opihi/dvo/imdata.c

    r20820 r20839  
    125125  /* create output vector */
    126126  NPTS = 1000;
    127   REALLOCATE (vec[0].elements, float, NPTS);
     127  REALLOCATE (vec[0].elements.Flt, opihi_flt, NPTS);
    128128  vec[0].Nelements = N = 0;
    129129
     
    153153          if ((catalog.measure[i].t < start) || (catalog.measure[i].t > stop)) continue;
    154154          n = catalog.measure[i].averef;
    155           vec[0].elements[N] = catalog.average[n].R - catalog.measure[i].dR / 3600.0;
    156           N++;
    157           CHECK_REALLOCATE (vec[0].elements, float, NPTS, N, 1000);
     155          vec[0].elements.Flt[N] = catalog.average[n].R - catalog.measure[i].dR / 3600.0;
     156          N++;
     157          CHECK_REALLOCATE (vec[0].elements.Flt, opihi_flt, NPTS, N, 1000);
    158158        }
    159159        break;
     
    162162          if ((catalog.measure[i].t < start) || (catalog.measure[i].t > stop)) continue;
    163163          n = catalog.measure[i].averef;
    164           vec[0].elements[N] = catalog.average[n].D - catalog.measure[i].dD / 3600.0;
    165           N++;
    166           CHECK_REALLOCATE (vec[0].elements, float, NPTS, N, 1000);
     164          vec[0].elements.Flt[N] = catalog.average[n].D - catalog.measure[i].dD / 3600.0;
     165          N++;
     166          CHECK_REALLOCATE (vec[0].elements.Flt, opihi_flt, NPTS, N, 1000);
    167167        }
    168168        break;
     
    170170        for (i = 0; i < catalog.Nmeasure; i++) {
    171171          if ((catalog.measure[i].t < start) || (catalog.measure[i].t > stop)) continue;
    172           vec[0].elements[N] = catalog.measure[i].M;
    173           N++;
    174           CHECK_REALLOCATE (vec[0].elements, float, NPTS, N, 1000);
     172          vec[0].elements.Flt[N] = catalog.measure[i].M;
     173          N++;
     174          CHECK_REALLOCATE (vec[0].elements.Flt, opihi_flt, NPTS, N, 1000);
    175175        }
    176176        break;
     
    178178        for (i = 0; i < catalog.Nmeasure; i++) {
    179179          if ((catalog.measure[i].t < start) || (catalog.measure[i].t > stop)) continue;
    180           vec[0].elements[N] = catalog.measure[i].dM;
    181           N++;
    182           CHECK_REALLOCATE (vec[0].elements, float, NPTS, N, 1000);
     180          vec[0].elements.Flt[N] = catalog.measure[i].dM;
     181          N++;
     182          CHECK_REALLOCATE (vec[0].elements.Flt, opihi_flt, NPTS, N, 1000);
    183183        }
    184184        break;
     
    186186        for (i = 0; i < catalog.Nmeasure; i++) {
    187187          if ((catalog.measure[i].t < start) || (catalog.measure[i].t > stop)) continue;
    188           vec[0].elements[N] = catalog.measure[i].Mcal;
    189           N++;
    190           CHECK_REALLOCATE (vec[0].elements, float, NPTS, N, 1000);
     188          vec[0].elements.Flt[N] = catalog.measure[i].Mcal;
     189          N++;
     190          CHECK_REALLOCATE (vec[0].elements.Flt, opihi_flt, NPTS, N, 1000);
    191191        }
    192192        break;
     
    195195          if ((catalog.measure[i].t < start) || (catalog.measure[i].t > stop)) continue;
    196196          n = catalog.measure[i].averef;
    197           // vec[0].elements[N] = catalog.average[n].M;
     197          // vec[0].elements.Flt[N] = catalog.average[n].M;
    198198          N++;
    199199        }
     
    202202        for (i = 0; i < catalog.Nmeasure; i++) {
    203203          if ((catalog.measure[i].t < start) || (catalog.measure[i].t > stop)) continue;
    204           vec[0].elements[N] = catalog.measure[i].photcode;
    205           N++;
    206           CHECK_REALLOCATE (vec[0].elements, float, NPTS, N, 1000);
     204          vec[0].elements.Flt[N] = catalog.measure[i].photcode;
     205          N++;
     206          CHECK_REALLOCATE (vec[0].elements.Flt, opihi_flt, NPTS, N, 1000);
    207207        }
    208208        break;
     
    210210        for (i = 0; i < catalog.Nmeasure; i++) {
    211211          if ((catalog.measure[i].t < start) || (catalog.measure[i].t > stop)) continue;
    212           vec[0].elements[N] = TimeValue (catalog.measure[i].t, TimeReference, TimeFormat);
    213           N++;
    214           CHECK_REALLOCATE (vec[0].elements, float, NPTS, N, 1000);
     212          vec[0].elements.Flt[N] = TimeValue (catalog.measure[i].t, TimeReference, TimeFormat);
     213          N++;
     214          CHECK_REALLOCATE (vec[0].elements.Flt, opihi_flt, NPTS, N, 1000);
    215215        }
    216216        break;
     
    220220 
    221221  vec[0].Nelements = N;
    222   REALLOCATE (vec[0].elements, float, MAX(1,N));
     222  REALLOCATE (vec[0].elements.Flt, opihi_flt, MAX(1,N));
    223223  return (TRUE);
    224224}
Note: See TracChangeset for help on using the changeset viewer.