Fix iod

PerezHz · Aug 28, 2024 · fb4aa41 · fb4aa41
1 parent 8f166bc
commit fb4aa41
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/orbitdetermination/orbitdetermination.jl b/src/orbitdetermination/orbitdetermination.jl
@@ -361,15 +361,15 @@ function iod(radec::Vector{RadecMPC{T}}, params::NEOParameters{T};
     end
     # Reduce tracklets by polynomial regression
     tracklets = reduce_tracklets(radec)
+    # Set jet transport variables
+    varorder = max(params.tsaorder, params.gaussorder, params.jtlsorder)
+    scaled_variables("dx", ones(T, 6); order = varorder)
     # Gauss method
     if gauss
         sol = min(sol, gaussinitcond(radec, tracklets, params; dynamics))
         # Termination condition
         nrms(sol) <= params.gaussQmax && length(sol.res) == length(radec) && return sol
     end
-    # Set jet transport variables
-    varorder = max(params.tsaorder, params.gaussorder, params.jtlsorder)
-    scaled_variables("dx", ones(T, 6); order = varorder)
     # Iterate tracklets
     for i in eachindex(tracklets)
         # Admissible region