Sunode pymc, need for dummy parameter

Questions v5
1

On the github readme for sunode ( GitHub - pymc-devs/sunode: Solve ODEs fast, with support for PyMC) they have an ‘extra’ parameter that is set to np.zeros(1). (This is in the call to sunode.wrappers.as_pytensor.solve_ivp)

This doesn’t seem to be optional, and if you do not include it the sampling fails ValueError: When changing to a smaller dtype... I don’t understand why this is needed? Is this working as intended?

2

@aseyboldt

3

Hello,

That isn’t really intended behavior, but a known bug that I didn’t get around to fixing yet.
I might have a good excuse to work a bit more on sunode soon, so hopefully I’ll find the time then…

2 Likes