Yet another adaptative maxsteps

PerezHz · Sep 3, 2023 · 326d79c · 326d79c
1 parent 4a4e323
commit 326d79c
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/orbit_determination/gauss_method.jl b/src/orbit_determination/gauss_method.jl
@@ -400,7 +400,11 @@ function adaptative_maxsteps(radec::Vector{RadecMPC{T}}) where {T <: AbstractFlo
     # Time difference [days]
     Δ_day = Δ_ms / 86_400_000
     # Adaptative maxsteps
-    return 35 + ceil(Int, 5*Δ_day/160)
+    if Δ_day <= 30
+        return ceil(Int, -7*Δ_day/3) + 100
+    else
+        return ceil(Int, (Δ_day + 360)/13)
+    end
 end
 
 for i in 1:2