Yingbo Ma

Add benchmarks for special matrices

1

Aside from the differences of BLAS and LAPACK in Julia and Matlab, it maybe meaningful to compare the algorithms between Julia and Matlab. For instance, `sqrtm` function requires to do...

#1683

#1675

structural_simplify for SDESystem

1

```julia @variables t u[1:3](t) params = @parameters ß=0.2 µ=0.1 D = Differential(t) eqs = [D(u[1]) ~ -ß * u[1] * u[2], D(u[2]) ~ ß * u[1] * u[2] - µ...

Implicit DAE observed indexing is broken

2

https://github.com/SciML/ModelingToolkitStandardLibrary.jl/pull/56 @sharanry could you take a look?

State selection

4

Without state selection, `ODAEProblem` w/ `Tsit5` gives terrible numerical drift on the pendulum model. ```julia using ModelingToolkit, LinearAlgebra, OrdinaryDiffEq function Pendulum(;name, g=9.81, r=1.) gval, rval = g, r @parameters g...

Not connected inputs are "time-dependent parameters", and connected inputs should be reduced. Outputs can always be reduced. Derivative of inputs should be reduced as well. E.g. ```julia using ModelingToolkit @variables...