Improve inference

Improve OrbitalData inference
CliMA · Sep 10, 2024 · 4b284c7 · 4b284c7
1 parent cc49082
commit 4b284c7
Showing 1 changed file with 12 additions and 5 deletions.
diff --git a/src/Insolation.jl b/src/Insolation.jl
@@ -28,14 +28,21 @@ end
 
 function OrbitalData()
     datapath = joinpath(artifact"laskar2004", "INSOL.LA2004.BTL.csv")
-    x, _ = readdlm(datapath, ',', Float64, header = true)
-    t_range = (x[1, 1] * 1e3):1e3:(x[end, 1] * 1e3) # array of every 1 kyr to range of years
+    Tx = Tuple{Matrix{Float64}, Matrix{AbstractString}}
+    x, _ = readdlm(datapath, ',', Float64, header = true)::Tx
+    T = StepRangeLen{
+        Float64,
+        Base.TwicePrecision{Float64},
+        Base.TwicePrecision{Float64},
+        Int,
+    }
+    t_range = ((x[1, 1] * 1e3):1e3:(x[end, 1] * 1e3))::T # array of every 1 kyr to range of years
     e_spline_etp =
-        cubic_spline_interpolation(t_range, x[:, 2], extrapolation_bc = NaN)
+        cubic_spline_interpolation(t_range, x[:, 2]; extrapolation_bc = NaN)
     γ_spline_etp =
-        cubic_spline_interpolation(t_range, x[:, 3], extrapolation_bc = NaN)
+        cubic_spline_interpolation(t_range, x[:, 3]; extrapolation_bc = NaN)
     ϖ_spline_etp =
-        cubic_spline_interpolation(t_range, x[:, 4], extrapolation_bc = NaN)
+        cubic_spline_interpolation(t_range, x[:, 4]; extrapolation_bc = NaN)
 
     E = typeof(e_spline_etp)
     G = typeof(γ_spline_etp)