Changeset 6268 for trunk/psLib/src/astro/psTime.c

Timestamp:

Jan 31, 2006, 1:24:21 PM (20 years ago)

Author:

drobbin

Message:

Fixed Time conversions and ToMJD & ToJD fxns to use TAI type. Updated transforms tests.

File:

• trunk/psLib/src/astro/psTime.c (modified) (9 diffs)

trunk/psLib/src/astro/psTime.c

@@ -10 +10 @@
  *  @author Ross Harman, MHPCC
  *
- *  @version $Revision: 1.79 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-01-26 21:10:22 $
+ *  @version $Revision: 1.80 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-01-31 23:24:21 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -93 +93 @@
 static psTime* convertTimeTTTAI(psTime* time);
 static psTime* convertTimeUTCUT1(psTime* time);
+
 
 /** Removes leading and trailing whitespace and # characters from a string. The cleaned string is a new null
@@ -581 +582 @@
     time->type = PS_TIME_TAI;
 
+    //XXX: Set leapseconds to TRUE
+    //    time->leapsecond = true;
     return time;
 }
@@ -676 +679 @@
     if(time->type == PS_TIME_UT1) {
         psError(PS_ERR_BAD_PARAMETER_VALUE,true,"Cannot convert from UT1 time type");
+        //        return NULL;
         return time;
     }
@@ -1161 +1165 @@
 
     // Determine Julian and modified Julian dates used in table lookup and time delta calculation
-    jd = psTimeToJD(time);
-    mjd = psTimeToMJD(time);
+    if(time->sec < 0) {
+        // psTime earlier than epoch
+        jd = time->sec / SEC_PER_DAY - time->nsec / NSEC_PER_DAY + JD_EPOCH_OFFSET;
+        mjd = time->sec / SEC_PER_DAY - time->nsec / NSEC_PER_DAY + MJD_EPOCH_OFFSET;
+    } else {
+        // psTime greater than epoch
+        jd = time->sec / SEC_PER_DAY + time->nsec / NSEC_PER_DAY + JD_EPOCH_OFFSET;
+        mjd = time->sec / SEC_PER_DAY + time->nsec / NSEC_PER_DAY + MJD_EPOCH_OFFSET;
+    }
 
     // Set ceiling of the julian date to the last entry in the lookup table
@@ -1205 +1216 @@
     PS_ASSERT_INT_WITHIN_RANGE(time2->nsec,0,(psU32)((1e9)-1),0);
     diff = abs((psS64)p_psTimeGetTAIDelta((psTime*)time1)-(psS64)p_psTimeGetTAIDelta((psTime*)time2));
-
     return diff;
 }
@@ -1248 +1258 @@
     PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,(psU32)((1e9)-1),NAN);
 
+    psTime *time2 = p_psTimeCopy(time);
+    //XXX: ADD says that this formula works only for PS_TIME_TAI, so adding the following:
+    if (time->type == PS_TIME_UTC || time->type == PS_TIME_TT) {
+        time2 = psTimeConvert(time2, PS_TIME_TAI);
+    }
+
     // Julian date conversion
-    if(time->sec < 0) {
+    if(time2->sec < 0) {
         // psTime earlier than epoch
-        jd = time->sec / SEC_PER_DAY - time->nsec / NSEC_PER_DAY + JD_EPOCH_OFFSET;
+        jd = time2->sec / SEC_PER_DAY - time2->nsec / NSEC_PER_DAY + JD_EPOCH_OFFSET;
     } else {
         // psTime greater than epoch
-        jd = time->sec / SEC_PER_DAY + time->nsec / NSEC_PER_DAY + JD_EPOCH_OFFSET;
-    }
+        jd = time2->sec / SEC_PER_DAY + time2->nsec / NSEC_PER_DAY + JD_EPOCH_OFFSET;
+    }
+
+    psFree(time2);
     return jd;
 }
@@ -1266 +1284 @@
     PS_ASSERT_PTR_NON_NULL(time,NAN);
     PS_ASSERT_INT_WITHIN_RANGE(time->nsec,0,(psU32)((1e9)-1),NAN);
+    psTime *time2 = p_psTimeCopy(time);
+    //XXX: ADD says that this formula works only for PS_TIME_TAI, so adding the following:
+    if (time->type == PS_TIME_UTC || time->type == PS_TIME_TT) {
+        time2 = psTimeConvert(time2, PS_TIME_TAI);
+    }
 
     // Modified Julian date conversion
-    if(time->sec < 0) {
+    if(time2->sec < 0) {
         // psTime earlier than epoch
-        mjd = time->sec / SEC_PER_DAY - time->nsec / NSEC_PER_DAY + MJD_EPOCH_OFFSET;
+        mjd = time2->sec / SEC_PER_DAY - time2->nsec / NSEC_PER_DAY + MJD_EPOCH_OFFSET;
     } else {
         // psTime greater than epoch
-        mjd = time->sec / SEC_PER_DAY + time->nsec / NSEC_PER_DAY + MJD_EPOCH_OFFSET;
-    }
+        mjd = time2->sec / SEC_PER_DAY + time2->nsec / NSEC_PER_DAY + MJD_EPOCH_OFFSET;
+    }
+    psFree(time2);
     return mjd;
 }
@@ -1810 +1834 @@
     PS_ASSERT_PTR_NON_NULL(inTime, NULL);
     psTime *outTime = psTimeAlloc(inTime->type);
-    *outTime = *inTime;
+    //    *outTime = *inTime;
+    outTime->sec = inTime->sec;
+    outTime->nsec = inTime->nsec;
+    outTime->leapsecond = inTime->leapsecond;
     return outTime;
 }
 
+