QuantumOptics.jl-documentation
QuantumOptics.jl-documentation copied to clipboard
qojulia
Missing details about default solver options
I'm wondering about the choice of default solver options;
DP5() for the algorithm,
https://github.com/qojulia/QuantumOptics.jl/blob/8012aa8cb224df75dcc0e1301b24c740c616d330/src/timeevolution_base.jl#L14-L18
and for the tolerance,
https://github.com/qojulia/QuantumOptics.jl/blob/8012aa8cb224df75dcc0e1301b24c740c616d330/src/timeevolution_base.jl#L59-L66
I can relate to the tolerance choice, as the default in DiffEqBase is quite low in my opinion (reltol=1e-3 & abstol=1e-6).
But not sure about the algorithm, why DP5 ? and not Tsit5 (as recommended for non-stiff problems)? or letting DiffEq choose automatically (with alg=nothing)?
The Time-evolution doesn't offer much insight. I would love for the reasoning for this to be documented 🙏
