-
Notifications
You must be signed in to change notification settings - Fork 48
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Integrate the entire dynamics for higher-order ODEs #1139
Conversation
Bencher Report
Click to view all benchmark results
|
fig.savefig("/tmp/test_integrate_odes_higher_order.png", dpi=300) | ||
|
||
# verify | ||
np.testing.assert_allclose(x[-1], 0.10737970490959549) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
How do we know these values are correct (are they from just running the NESTML model once, or from an independent implementation of the ODEs)?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
These values are from repeated simulations of the NESTML-generated code, equivalent to providing integrate_odes(x',x)
in the NESTML model. I have changed the test now to include calls to both integrate_odes
and compare their output.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
With the latest changes, integrate_odes(x, x')
is not supported and hence I changed the test again to compare it against numerical values. Do you have any other ideas to check this?
Thanks for the PR! Could you add/update the documentation for the |
integrate_odes(x) | ||
# integrate_odes(x', x) # equivalent to integrate_odes(x) | ||
|
||
Here, ``integrate_odes(x)`` integrates the entire dynamics of ``x(t)`` and is equivalent to ``integrate_odes(x', x)``. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Should it be allowed to write variables with derivatives inside integrate_odes()
at all? As in, x
is allowed, but x'
is not. Otherwise, what is the semantics of saving integrate_odes(x')
?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I have now updated one of the existing cocos to check this condition.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Great, thank you!
No description provided.