DifferentialEquations.jl
DifferentialEquations.jl copied to clipboard
SciML
Multi-language suite for high-performance solvers of differential equations and scientific machine learning (SciML) components. Ordinary differential equations (ODEs), stochastic differential equation...
The first run of a large (>1000) system of ODEs, using an algorithm that has `autodiff=true ` as the default (i.e `AutoTsit5(Rosenbrock23())`) results in a very long compilation time. Setting...
MWE: ``` using DifferentialEquations function usode!(du,u,p,t) C1 = 1.63e-9 RS = 790.0 C2 = 3.3e-9 C0 = 0.3e-9 L1 = 1e-3 L2 = 82e-3 CM = 4e-12 LM = 8.0...
Hi, It seems that there is an issue solving a `SplitODEProblem` with initial conditions of type `u0::Array{Float64,2}`. I get the error : `BoundsError: attempt to access 200-element view(::Array{Float64,2}, :, 1)...
Trying to be lazy and doing `SavedValues(Float64, Array{any})` leads to an Internal Error. ``` _ _ _ _(_)_ | Documentation: https://docs.julialang.org (_) | (_) (_) | _ _ _| |_...
The function evaluation count in the `.destats` substruct is useful as part of gauging computational effort taken over a solution. My use case is simple non-stiff ODEs. I have not...
Per Slack discussion with @ChrisRackauckas, opening an issue about the need to have an interface for users to get at the `W` operator components. In my case I have a...
Specifying the solver arg `linsolve=LinSolveGMRES` when solving a in-place function ```julia A = [1. 0 0 -5 4 -2 4 -3 -4 0 0 1 5 -2 2 3] u0...
Hi, I couldn't figure out from the docs how to use derivative-based methods for complex problems. Some things are mentioned in the news blog, others I still don't know. It...
DifferentialEquations supports storing additional values (besides u, du) at a pre-defined time-grid. However, **it seems not possible to store additional values also at event points**. Note, this is important for...
The callback `save_func(u, t, integrator)` can be provided to compute additional outputs `y` based on the current time t and the current state u (and the model parameters integrator.p): `y=f(u,t,p)`....