I am attempting to fit several SIR models at once to model a Bayesian Hierarchical Model.
The current signature for the ODE integrators is
real ODEFunc(real t, real y, real params, real x_r, int x_i)
Is there any way to solve an ODE which looks something like this?
real[,] ODEFunc(real t, real[,] y, real[,] params, real x_r, int x_i)
Essentially, if we have a differential equation that looks like this:
dydt = params * y
in the 1d case, dydt, params, and y are all reals. Can the ODE solver work in the case where dydt, params, and y are real[n], where n is the number of different fits I want?
Sorry if this was asked somewhere else or if I’m ignorant of something. I’m not really well versed in how the ODE solvers work , and I couldn’t find a related topic in the forum that answered this question.