@@ -62,7 +62,7 @@ prob = ODEProblem{true, SciMLBase.FullSpecialize}(sys, [], (0, 1.0), [])
6262
6363model_true =  structural_simplify (lotka_true ())
6464prob_true =  ODEProblem {true, SciMLBase.FullSpecialize} (model_true, [], (0 , 1.0 ), [])
65- sol_ref =  solve (prob_true, Rodas5P (), abstol =  1e-8  , reltol =  1e-8 )
65+ sol_ref =  solve (prob_true, Rodas5P (), abstol =  1e-10  , reltol =  1e-8 )
6666
6767x0 =  default_values (sys)[nn. p]
6868
@@ -74,7 +74,7 @@ function loss(x, (prob, sol_ref, get_vars, get_refs, set_x))
7474    new_p =  set_x (prob, x)
7575    new_prob =  remake (prob, p =  new_p, u0 =  eltype (x).(prob. u0))
7676    ts =  sol_ref. t
77-     new_sol =  solve (new_prob, Rodas5P (), abstol =  1e-8  , reltol =  1e-8 , saveat =  ts)
77+     new_sol =  solve (new_prob, Rodas5P (), abstol =  1e-10  , reltol =  1e-8 , saveat =  ts)
7878
7979    loss =  zero (eltype (x))
8080
@@ -122,7 +122,7 @@ op = OptimizationProblem(of, x0, ps)
122122#      false
123123#  end
124124
125- res =  solve (op, Adam (), maxiters =  5000 )# , callback = plot_cb)
125+ res =  solve (op, Adam (), maxiters =  10000 )# , callback = plot_cb)
126126
127127@test  res. objective <  1 
128128
0 commit comments