I’m working on model using Differential Equation that has two state variables, current and temperature. So when using
DifferentialEquations.jl I define something like
f(du, u, p, t) = ...
u is a vector with two elements, each of which I’d like to have different units and therefore different types. And du has yet different units from
u since we add
1/s when we take the time derivative.
Any advice or packages on how to deal with this in a type stable way? For now I strip all the units off and pass floats to
f, but I’d rather propagate the units all the way through if possible.