Christopher Rackauckas

Generate intermediate tuples instead of arrays, i.e. ```julia δ₁, δ₂, δ₃, δ₄ = (0.003/100, 0.004/100, (0.015+0.030+0.064+0.213+0.718)/(5*100), (2.384+8.466+12.497+1.117)/(4*100)) ``` This code: ```julia δ = vcat(repeat([δ₁],1),repeat([δ₂],1),repeat([δ₃],1),repeat([δ₄],4-1-1-1)) ``` is a similar issue: get...

@ajkeller34 an MWE for what I am talking about goes like this. We have united arrays `A` and `B` which exactly match types. We want to define scalars `a` and...

In the example above `a` and `b` should actually be unitless, so it's easier in this case. But I'm still not sure how to represent the Array properly. Use union...

That tutorial was written in 2016 (the time stamps are just newer because of the auto-updating infrastructure). In 202(1), I think the way that Unitful.jl expresses units is a really...

xref https://github.com/luchr/ODEInterface.jl/issues/27 It seems like you don't have a working Fortran compiler? This is why we don't include ODEInterface by default: it requires specific compilers so you need to do...

https://github.com/SciML/ODEInterfaceDiffEq.jl#installation

So the workaround would be to set an initial dt like `dt=0.01`. It's a bit odd that this broadcast doesn't make a LabelledArray though.

The 2020 way is to use DiffEqFlux because it's much more flexible. https://diffeqflux.sciml.ai/dev/examples/optimization_ode/ This tutorial solves with a ton of initial conditions and parallelizes over them: https://diffeqflux.sciml.ai/dev/examples/optimization_sde/ Just remake `u0`...

If you do `sol=nothing` and then call `GC.gc()` does it clear the memory? It should, and if that doesn't there's a weird error to investigate. Another thing that's helpful is...

My guess is that it doesn't clear the memory for `sol` until after the computation is done, since in theory it could be modified in the function?