From 7576e5bbd69d42b7b35a7eea803fd704f72e2916 Mon Sep 17 00:00:00 2001 From: Herman Sletmoen Date: Wed, 13 Aug 2025 11:40:03 +0200 Subject: [PATCH 1/2] Update default perturbation solver to Rodas5P --- src/solve.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/solve.jl b/src/solve.jl index 442cc967..364a4819 100644 --- a/src/solve.jl +++ b/src/solve.jl @@ -263,7 +263,7 @@ end function solve( prob::CosmologyProblem, ks::Union{Nothing, AbstractArray} = nothing; bgopts = (alg = Rodas4P(), reltol = 1e-9, abstol = 1e-9), - ptopts = (alg = KenCarp4(), reltol = 1e-8, abstol = 1e-8), + ptopts = (alg = Rodas5P(), reltol = 1e-7, abstol = 1e-7), shootopts = (alg = NewtonRaphson(), abstol = 1e-5), thread = true, verbose = false, kwargs... ) @@ -276,7 +276,7 @@ If `threads`, integration over independent perturbation modes are parallellized. function solve( prob::CosmologyProblem, ks::Union{Nothing, AbstractArray} = nothing; bgopts = (alg = Rodas4P(), reltol = 1e-9, abstol = 1e-9), - ptopts = (alg = KenCarp4(), reltol = 1e-8, abstol = 1e-8), + ptopts = (alg = Rodas5P(), reltol = 1e-7, abstol = 1e-7), shootopts = (alg = NewtonRaphson(), abstol = 1e-5), thread = true, verbose = false, kwargs... ) From 0b1b06e5a0a83ceba61d39a4d91760c364dcfea6 Mon Sep 17 00:00:00 2001 From: Herman Sletmoen Date: Wed, 13 Aug 2025 17:23:17 +0200 Subject: [PATCH 2/2] Update solve.jl --- src/solve.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/solve.jl b/src/solve.jl index 364a4819..564f2c36 100644 --- a/src/solve.jl +++ b/src/solve.jl @@ -263,7 +263,7 @@ end function solve( prob::CosmologyProblem, ks::Union{Nothing, AbstractArray} = nothing; bgopts = (alg = Rodas4P(), reltol = 1e-9, abstol = 1e-9), - ptopts = (alg = Rodas5P(), reltol = 1e-7, abstol = 1e-7), + ptopts = (alg = Rodas5P(), reltol = 1e-5, abstol = 1e-5), shootopts = (alg = NewtonRaphson(), abstol = 1e-5), thread = true, verbose = false, kwargs... ) @@ -276,7 +276,7 @@ If `threads`, integration over independent perturbation modes are parallellized. function solve( prob::CosmologyProblem, ks::Union{Nothing, AbstractArray} = nothing; bgopts = (alg = Rodas4P(), reltol = 1e-9, abstol = 1e-9), - ptopts = (alg = Rodas5P(), reltol = 1e-7, abstol = 1e-7), + ptopts = (alg = Rodas5P(), reltol = 1e-5, abstol = 1e-5), shootopts = (alg = NewtonRaphson(), abstol = 1e-5), thread = true, verbose = false, kwargs... )