Use faster Julian date calculations for offsets

JuliaAstro · Aug 27, 2018 · d6e476f · d6e476f
1 parent b9ac603
commit d6e476f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/offsets.jl b/src/offsets.jl
@@ -44,13 +44,13 @@ the quantity TDB-TT can differ by as much as about 4 microseconds.
 
 """
 function tai_offset(::BarycentricDynamicalTime, ep)
-    dt = get(days(ep - J2000_EPOCH))
+    dt = j2000(TT, ep)
     g = 357.53 + 0.9856003dt
     tai_offset(TT, ep) + 0.001658sind(g) + 0.000014sind(2g)
 end
 
 function tai_offset(::BarycentricDynamicalTime, ep, ut, elong, u, v)
-    t = get(centuries(ep - J2000_EPOCH)) / 10.0
+    t = get(centuries(j2000(TT, ep) * days)) / 10.0
     # Convert UT to local solar time in radians.
     tsol = mod(ut, 1.0) * 2π  + elong