-
-
Notifications
You must be signed in to change notification settings - Fork 206
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
Unknowns in the initialization_eqs
should be added to the variable list
#2787
Comments
y is missing a guess. The error message needs to be handled but that is trying to throw an error about missing a guess on y which would be required to solve the nonlinear system. |
does |
This is my doubt too. (and I wonder if is possible to avoid that) |
As @ChrisRackauckas mentioned, in this case the initialization equations are not satisfiable (x=>1, y=>1, x^2 + y^2 ~ 1). |
Yes, the equations are not satisfiable. Note that it's not But also, |
Ok. That makes sense. Thanks |
You also invented a new feature, which is easier to just support. #2789 |
Describe the bug 🐞
The unknowns in the
initialization_eqs
aren't automatically added to the variable map/list. ATM, it should additionally be added to theguesses
. See the example (from docs) below.Expected behavior
Unknowns that appear in the lhs should be just added to the list.
Minimal Reproducible Example 👇
(The example is modified from here: https://docs.sciml.ai/ModelingToolkit/stable/tutorials/initialization/
prob in the docs:
prob = ODEProblem(pend, [x => 1], (0.0, 1.5), [g => 1], guesses = [λ => 0, y => 1], initialization_eqs = [y ~ 1])
)Error & Stacktrace⚠️
Environment (please complete the following information):
using Pkg; Pkg.status()
versioninfo()
Additional context
One might have to account for the
Differential
operator inintializtion_eqs
ofD(w) ~ 0
sort.The text was updated successfully, but these errors were encountered: