Update default parameters

src/orbitdetermination/tooshortarc.jl

@@ -1,18 +1,28 @@
 @doc raw"""
     adam(radec::Vector{RadecMPC{T}}, A::AdmissibleRegion{T}, ρ::T, v_ρ::T,
-         params::NEOParameters{T}; scale::Symbol = :linear, η::T = 25.0,
-         μ::T = 0.75, ν::T = 0.9, ϵ::T = 1e-8, Qtol::T = 0.001, adamorder::Int = 2,
-         dynamics::D = newtonian!) where {T <: Real, D}
+         params::NEOParameters{T}; kwargs...) where {T <: Real, D}
 
 Adaptative moment estimation (ADAM) minimizer of normalized mean square
-residual over the manifold of variations of `A`.
+residual of `radec` over the manifold of variations of `A`, starting
+from `(ρ, v_ρ)`.
+
+## Keyword arguments
+
+- `scale::Symbol`: horizontal scale, either `:log` (default) or `:linear`.
+- `η::T`: learning rate (default: `25.0`).
+- `μ::T`: first moment (default: `0.75`).
+- `ν::T`: second moment (default: `0.9`).
+- `ϵ::T`: numerical stability constant (default: `1e-8`).
+- `Qtol::T`: target function relative tolerance (default: `1e-5`).
+- `adamorder::Int`: jet transport order (default: `2`).
+- `dynamics::D`: dynamical model (default: `newtonian!`).
 
 !!! reference
     See Algorithm 1 of https://doi.org/10.48550/arXiv.1412.6980.
 """
 function adam(radec::Vector{RadecMPC{T}}, A::AdmissibleRegion{T}, ρ::T, v_ρ::T,
               params::NEOParameters{T}; scale::Symbol = :linear, η::T = 25.0,
-              μ::T = 0.75, ν::T = 0.9, ϵ::T = 1e-8, Qtol::T = 0.001, adamorder::Int = 2,
+              μ::T = 0.75, ν::T = 0.9, ϵ::T = 1e-8, Qtol::T = 1e-5, adamorder::Int = 2,
               dynamics::D = newtonian!) where {T <: Real, D}
     # Initial time of integration [julian days]
     jd0 = datetime2julian(A.date)

src/propagation/parameters.jl

@@ -89,7 +89,7 @@ Parameters for all orbit determination functions.
 
 - `H_max::T`: maximum absolute magnitude (default: `34.5`).
 - `a_max::T`: maximum semimajor axis (default: `100.0`).
-- `adamiter::Int`: maximum number of iterations for `ADAM` optimizer (default: `200`).
+- `adamiter::Int`: maximum number of iterations for `ADAM` optimizer (default: `500`).
 - `tsaorder::Int`: order of the jet transport perturbation (default: `6`).
 - `tsaQmax::T`: nrms threshold (default: `1.5`).
 
@@ -108,7 +108,7 @@ function NEOParameters(;
     abstol::T = 1e-20, parse_eqs::Bool = true, bwdoffset::T = 0.5, fwdoffset::T = 0.5,
     coeffstol::T = 10.0, debias_table::String = "2018", max_triplets::Int = 10,
     gaussorder::Int = 5, adamhelp::Bool = false, gaussQmax::T = 5.0, H_max::T = 34.5,
-    a_max::T = 100.0, adamiter::Int = 200, tsaorder::Int = 6, tsaQmax::T = 1.5,
+    a_max::T = 100.0, adamiter::Int = 500, tsaorder::Int = 6, tsaQmax::T = 1.5,
     jtlsiter::Int = 5, newtoniter::Int = 5, jtlsorder::Int = 5, max_per::T = 18.0
     ) where {T <: Real}