def testJulianEpochFromTAI(self):
"""Test julianEpochFromTAI and inverse
"""
for tai in (4232.89, 20000.32, 56350.03, 74222.9):
self.assertAlmostEqual(coordConv.julianEpochFromTAI(tai), julianEpochFromTAI(tai), places=7)
self.assertAlmostEqual(coordConv.taiFromJulianEpoch(coordConv.julianEpochFromTAI(tai)), tai, places=5)
def testJulianEpochFromTAI(self):
"""Test julianEpochFromTAI and inverse
"""
for tai in (4232.89, 20000.32, 56350.03, 74222.9):
self.assertAlmostEqual(coordConv.julianEpochFromTAI(tai),
julianEpochFromTAI(tai),
places=7)
self.assertAlmostEqual(coordConv.taiFromJulianEpoch(
coordConv.julianEpochFromTAI(tai)),
tai,
places=5)
def getCoordSys(coordSysCode, date, tai):
if coordSysCode < -1:
# the old TCC use LAST computed from TAI; the new TCC uses TAI
date = tai
elif coordSysCode == -1 and date == 0:
date = coordConv.julianEpochFromTAI(tai)
return CSysDict[coordSysCode](date)
def getCoordSys(coordSysCode, date, tai):
if coordSysCode < -1:
# the old TCC use LAST computed from TAI; the new TCC uses TAI
date = tai
elif coordSysCode == -1 and date == 0:
date = coordConv.julianEpochFromTAI(tai)
return CSysDict[coordSysCode](date)
def testDateFromTAI(self):
"""Test dateFromTAI
"""
for cls in JulianSysList:
coordSys = cls(2000)
for tai in (4232.89, 20000.32, 56350.03, 74222.9):
predDate = coordConv.julianEpochFromTAI(tai)
self.assertAlmostEqual(predDate, coordSys.dateFromTAI(tai))
coordSys = coordConv.FK4CoordSys(2000)
for tai in (4232.89, 20000.32, 56350.03, 74222.9):
ttDays = (tai + coordConv.TT_TAI) / coordConv.SecPerDay
predDate = 1900.0 + ((ttDays - 15019.81352 ) / 365.242198781)
self.assertAlmostEqual(predDate, coordSys.dateFromTAI(tai), places=5)
for cls in AzAltSysList:
for tai in (4232.89, 20000.32, 56350.03, 74222.9):
coordSys = cls(2000)
self.assertAlmostEqual(tai, coordSys.dateFromTAI(tai))
def testDateFromTAI(self):
"""Test dateFromTAI
"""
for cls in JulianSysList:
coordSys = cls(2000)
for tai in (4232.89, 20000.32, 56350.03, 74222.9):
predDate = coordConv.julianEpochFromTAI(tai)
self.assertAlmostEqual(predDate, coordSys.dateFromTAI(tai))
coordSys = coordConv.FK4CoordSys(2000)
for tai in (4232.89, 20000.32, 56350.03, 74222.9):
ttDays = (tai + coordConv.TT_TAI) / coordConv.SecPerDay
predDate = 1900.0 + ((ttDays - 15019.81352) / 365.242198781)
self.assertAlmostEqual(predDate,
coordSys.dateFromTAI(tai),
places=5)
for cls in AzAltSysList:
for tai in (4232.89, 20000.32, 56350.03, 74222.9):
coordSys = cls(2000)
self.assertAlmostEqual(tai, coordSys.dateFromTAI(tai))
