Yingbo Ma

`ndims(x)` isn't the chunk size though. Do you mean `axes(x, dims)`? (so that it doesn't even assume 1-based indexing.)

The following code produces this error, too. ```julia using StochasticDiffEq, BSON function lorenz(du,u,p,t) du[1] = 10.0(u[2]-u[1]) du[2] = u[1]*(28.0-u[3]) - u[2] du[3] = u[1]*u[2] - (8/3)*u[3] end function σ_lorenz(du,u,p,t) du[1]...

Maybe force to print the deprecation? Julia now hides depwarns by default.

Sometimes unit normalization is necessary to get a stable model. I wonder if we could support that as well. For instance, convert all farads to microfarads or second to micro...

Unitful has `preferunits`, so we can use ```julia julia> using Unitful julia> Unitful.preferunits(u"ms") julia> upreferred(1e-12u"F") 1.0 A² ms⁴ kg⁻¹ m⁻² ``` to scale everything to a consistent time scale.

MTK doesn't really support array variables right now, and I think we should be careful on designing it.

You should remove `.vscode` and `.DS_Store`, then the PR is good to go after formatting.

Format failure.

The current MTK code example of scalarized code is ```julia using ModelingToolkit l = 1 n = 5 h = l / n @variables t u[1:n](t) D = Differential(t) #...

Are you sure you are using JuliaFormater 1.0?