The original code used the + operator to join the initial value and the integration resuslts. This led to incorrect array shapes and numpy treating the data inside as dtype=object. Thus scipy interp1d did not work directly with the output of ddeint.
Tested with:
1D DDE (provided examples and own models)
The original code used the
+
operator to join the initial value and the integration resuslts. This led to incorrect array shapes and numpy treating the data inside as dtype=object. Thus scipy interp1d did not work directly with the output of ddeint.Tested with: 1D DDE (provided examples and own models)
Not tested with: Systems of DDEs