Skip to content

Commit

Permalink
Minor fixes
Browse files Browse the repository at this point in the history
  • Loading branch information
LuEdRaMo committed May 29, 2024
1 parent e3af3bf commit f12c881
Show file tree
Hide file tree
Showing 2 changed files with 2 additions and 2 deletions.
2 changes: 1 addition & 1 deletion src/orbitdetermination/gaussinitcond.jl
Original file line number Diff line number Diff line change
Expand Up @@ -467,7 +467,7 @@ function gaussinitcond(radec::Vector{RadecMPC{T}}, tracklets::Vector{Tracklet{T}
# Gauss method solution
sol = gauss_method(observatories, dates, α .+ dq[1:3], δ .+ dq[4:6], params)
# Filter non-physical (negative) rho solutions
filter!(x -> cte(cte(getindex(getproperty(x, ), 2))) > 0, sol)
filter!(x -> all(cte.(x.ρ) .> 0), sol)
# Filter Gauss solutions by heliocentric energy
filter!(x -> cte(cte(heliocentric_energy(x.statevect))) <= 0, sol)
# Iterate over Gauss solutions
Expand Down
2 changes: 1 addition & 1 deletion src/orbitdetermination/orbitdetermination.jl
Original file line number Diff line number Diff line change
Expand Up @@ -59,7 +59,7 @@ function propres!(res::Vector{OpticalResidual{T, U}}, radec::Vector{RadecMPC{T}}
fwd = _propagate(dynamics, jd0, nyears_fwd, q0, buffer, params)
if !issuccessfulprop(bwd, t0 - _jd0_; tol = params.coeffstol) ||
!issuccessfulprop(fwd, tf - _jd0_; tol = params.coeffstol)
res = Vector{OpticalResidual{T, U}}(undef, 0)
empty!(res)
return bwd, fwd
end
# O-C residuals
Expand Down

0 comments on commit f12c881

Please sign in to comment.