Changeset 358 in Sophya for trunk/Poubelle/archTOI.old/toiiter.cc

Timestamp:

Aug 6, 1999, 3:46:37 PM (26 years ago)

Author:

ansari

Message:

auxilliary GPS

File:

• trunk/Poubelle/archTOI.old/toiiter.cc (modified) (9 diffs)

trunk/Poubelle/archTOI.old/toiiter.cc

@@ -38 +38 @@
  lastSample = -1;
  maxLookAhead = 10000;
+ 
+ auxGPS = NULL;
 
 }
@@ -72 +74 @@
   lastSample = x.lastSample;
   maxLookAhead = x.maxLookAhead;
+  
+  auxGPS = x.auxGPS; 
+  if (auxGPS) auxGPS = auxGPS->clone();
 }
 
@@ -78 +83 @@
   delete rawIter;
   delete[] interp;
+  delete auxGPS;
 }
 
@@ -311 +317 @@
   while (1) {
     if (!NextSample()) return false; // end of files
-    double t = tBlock0+(file->blockNum() * file->nEchBlock() + imes) * perEch/86400.;
+    double t = getMJD();
     if (t < tStart) continue;
     if (t > tEnd) return false;
@@ -378 +384 @@
 bool TOIIter::canGetValue(int column) {
    if (column < 0 || column >= infos.size()) return false;
+   TOIKind kind = infos[column].kind;
+   if (auxGPS && 
+       (kind == longitude || kind == latitude || kind == altitude)) {
+         double dummy;
+         return auxGPS->getLocation(getMJD(), dummy, dummy, dummy) == 0;
+   }
    if (rawIter) {
        return interp[column].canGet(lastSample);
    }
-   TOIKind kind = infos[column].kind;
    int index = infos[column].index;
    switch (kind) {
@@ -418 +429 @@
 double TOIIter::getValue(int column) {
    if (column < 0 || column >= infos.size()) return -1;
+   TOIKind kind = infos[column].kind;
+   if (auxGPS && 
+       (kind == longitude || kind == latitude || kind == altitude)) {
+         double lat,lon,alt;
+         if (auxGPS->getLocation(getMJD(), lat, lon, alt)) return -99999;
+         if (kind == longitude) return lon;
+         if (kind == latitude)  return lat;
+         if (kind == altitude)  return alt;         
+   }
    if (rawIter) {
      if (infos[column].interpolated)
@@ -424 +444 @@
        return interp[column].getEValue(lastSample);
    }
-   TOIKind kind = infos[column].kind;
    int index = infos[column].index;
    switch (kind) {
      case sampleNum:
-       return file->blockNum() * file->nEchBlock() + imes;
+       return getSampleIndex();
      case internalTime:
-       //return (file->blockNum() * file->nEchBlock() + imes) * file->perEchant();
-       return (file->blockNum() * file->nEchBlock() + imes) * perEch;
+       return getSampleIndex() * perEch;
      case mjd:
-      /* printf("mjd: %d %d %g %g %g\n",file->blockNum(),
-         (file->blockNum() * file->nEchBlock() + imes),
-         file->perEchant(),
-         (file->blockNum() * file->nEchBlock() + imes) * file->perEchant()/86400.,
-         tBlock0+(file->blockNum() * file->nEchBlock() + imes) * file->perEchant()/86400.); */
-       //return tBlock0+(file->blockNum() * file->nEchBlock() + imes) * file->perEchant()/86400.;
-       return tBlock0+(file->blockNum() * file->nEchBlock() + imes) * perEch/86400.;
+       return getMJD();
      case boloTens:
        return file->getMuVBolo(index, imes);
@@ -474 +486 @@
 bool   TOIIter::newValue(int column) {
    if (column < 0 || column >= infos.size()) return false;
+   TOIKind kind = infos[column].kind;
+   if (auxGPS && 
+       (kind == longitude || kind == latitude || kind == altitude)) {
+      return true;
+   }
    if (rawIter) {
        return interp[column].isNewValue(lastSample);
    }
-   TOIKind kind = infos[column].kind;
    switch (kind) {
      case sampleNum:
@@ -567 +583 @@
    }
 }
+
+double TOIIter::getMJD() {
+  int sample = getSampleIndex();
+  return tBlock0 + sample*perEch/86400.;
+}
  
 bool TOIIter::fetchAhead() { // Seulement si delegation