Arno Strouwen
Arno Strouwen
```jl using ModelingToolkit @variables t x1(t) x2(t) D = Differential(t) @named fol = ODESystem([ D(x1) ~ -x2,D(x2) ~ -x1]) prob = ODEProblem(fol, [x1 => 1.0,x2 => 2.0], (0.0,10.0)) sparsepatternprob =...
https://mtk.sciml.ai/dev/tutorials/acausal_components/#Simplifying-and-Solving-this-System Docs say 2 states for ODEProblem, 1 state for ODAEProblem, but it looks to me like both are reduced to 1 state?
```jl using OrdinaryDiffEq using DiffEqSensitivity using ForwardDiff function f(du, u, p, t) du[1] = dx = p[1] * u[1] - p[2] * u[1] * u[2] du[2] = dy = -p[3]...
Is it worthwhile to set up sensitivity problems with `DiffEqArrayOperator` ? I made a forward prototype that does not exploit the structure of A, i.e. it is a dense matrix....
Issues with JacVecOperator and LinSolveGMRES() when performing sensitivity analysis. Also a more general ReverseDiff issue? ```jl using OrdinaryDiffEq, DiffEqSensitivity, DiffEqOperators, Zygote, ReverseDiff, ForwardDiff, ModelingToolkit function fiip(du,u,p,t) du[1] = dx =...
```jl using DiffEqSensitivity using OrdinaryDiffEq using Zygote using ForwardDiff using BenchmarkTools function fiip(du,u,p,t) du[1] = dx = p[1]*u[1] - p[2]*u[1]*u[2] du[2] = dy = -p[3]*u[2] + p[4]*u[1]*u[2] end u0 =...
https://symbolics.juliasymbolics.org/dev/manual/build_function/ seems to be written back when there was only scalar. How does build function now work for functions that contain vector/matrix operations? Sentences like "If ex is an array...
https://symbolics.juliasymbolics.org/dev/tutorials/converting_to_C/ The generated function is named diffeqf by default. Maybe this should be generalized since the scope of this project has increased?
Currently, if you go to the documentation of a Julia package e.g. https://juliadocs.github.io/Documenter.jl/ the `` element of the page is `Home · Documenter.jl`, where the "Home" comes from the first...
Should the pages keyword argument not be in here? https://juliadocs.github.io/Documenter.jl/dev/lib/public/#Documenter.makedocs