DAE Problems

Missing docstring.

Missing docstring for SciMLBase.DAEProblem. Check Documenter's build log for details.

Missing docstring.

Missing docstring for SciMLBase.DAEFunction. Check Documenter's build log for details.

Solution Type

Missing docstring.

Missing docstring for SciMLBase.DAESolution. Check Documenter's build log for details.

Alias Specifier

SciMLBase.DAEAliasSpecifier — Type
DAEAliasSpecifier(;
    alias_p = nothing, alias_f = nothing, alias_u0 = nothing,
    alias_du0 = nothing, alias_tstops = nothing, alias = nothing
)

Control which DAEProblem inputs and solver option arrays may be aliased.

alias_u0 controls the initial state, alias_du0 controls the initial derivative array, alias_p controls the parameter object, alias_f controls the DAE function object, and alias_tstops controls the tstops vector. A value of nothing delegates to the solver default. Set alias = true or alias = false to apply the same policy to all fields.

Keywords

  • alias_p::Union{Bool, Nothing}: alias the parameter object.
  • alias_f::Union{Bool, Nothing}: alias the DAE function object.
  • alias_u0::Union{Bool, Nothing}: alias the u0 array.
  • alias_du0::Union{Bool, Nothing}: alias the du0 array.
  • alias_tstops::Union{Bool, Nothing}: alias the tstops array.
  • alias::Union{Bool, Nothing}: set every field of the DAEAliasSpecifier.

Example Problems

Examples problems can be found in DiffEqProblemLibrary.jl.

To use a sample problem, such as prob_dae_resrob, you can do something like:

import DiffEqProblemLibrary.DAEProblemLibrary, Sundials
prob = DAEProblemLibrary.prob_dae_resrob
sol = solve(prob, Sundials.IDA())
DAEProblemLibrary.prob_dae_resrob — Constant

The Robertson biochemical reactions in DAE form

\[\begin{align*} \frac{dy₁}{dt} &= -k₁y₁+k₃y₂y₃ \\ \frac{dy₂}{dt} &= k₁y₁-k₂y₂^2-k₃y₂y₃ \\ 1 &= y₁ + y₂ + y₃ \end{align*}\]

where $k₁=0.04$, $k₂=3×10^7$, $k₃=10^4$. For details, see: Hairer Norsett Wanner Solving Ordinary Differential Equations I - Nonstiff Problems Page 129 Usually solved on $[0,10^{11}]$