def testEllab(self):
"""tpm.ellab => apply elliptic aberration."""
# pytpm/tests/c_tests/ellab_test.c
v6 = tpm.V6S(r=1e9, alpha=tpm.h2r(20), delta=tpm.d2r(40.0))
v6 = v6.s2c()
v6 = tpm.ellab(tpm.J2000, v6, -1)
v6 = v6.c2s()
self.assertAlmostEqual(v6.r, 1e9, 5)
self.assertAlmostEqual(tpm.r2h(tpm.r2r(v6.alpha)), 20.000007838, 8)
self.assertAlmostEqual(tpm.r2d(tpm.r2r(v6.delta)), 39.999987574, 8)
def testEllab(self):
"""tpm.ellab => apply elliptic aberration."""
# pytpm/tests/c_tests/ellab_test.c
v6 = tpm.V6S(r=1e9, alpha=tpm.h2r(20), delta=tpm.d2r(40.0))
v6 = v6.s2c()
v6 = tpm.ellab(tpm.J2000, v6, -1)
v6 = v6.c2s()
self.assertAlmostEqual(v6.r, 1e9,5)
self.assertAlmostEqual(tpm.r2h(tpm.r2r(v6.alpha)), 20.000007838,8)
self.assertAlmostEqual(tpm.r2d(tpm.r2r(v6.delta)), 39.999987574,8)
def testR2H(self):
"""tpm.r2h() => convert radians into hours."""
import math
self.assertAlmostEqual(tpm.r2h(math.pi), 12.0)
self.assertAlmostEqual(tpm.r2h(2*math.pi), 24.0)
self.assertAlmostEqual(tpm.r2h(-(12.5/12.0)*math.pi), -12.5)
def testR2H(self):
"""tpm.r2h() => convert radians into hours."""
import math
self.assertAlmostEqual(tpm.r2h(math.pi), 12.0)
self.assertAlmostEqual(tpm.r2h(2 * math.pi), 24.0)
self.assertAlmostEqual(tpm.r2h(-(12.5 / 12.0) * math.pi), -12.5)
