-
-
Notifications
You must be signed in to change notification settings - Fork 4.3k
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
Incorrect behaviour for solve_poly_system #26682
Comments
As a workaround you can use |
Somehow the Groebner basis comes out as
That shouldn't happen because if 1 is in the Groebner basis then it should be the whole basis. |
The problem is here: sympy/sympy/polys/groebnertools.py Lines 103 to 107 in 73ae5e9
We get a small h that should be zero but is not quite zero because of floating point rounding error:
|
Possibly the algorithms used for |
Using solve_poly_system(eqns,s,c,domain=QQ) seems slow because of all the rational computations, and I need to solve many of these systems. Any work around for that? |
The algorithm needs to use exact calculations to work properly. Please open a separate issue about slowness and give an example of a case that is slow. |
This code raises NotImplementedError: only zero-dimensional systems supported (finite number of solutions), but the system is indeed zero-dimensional. Furthermore, if the order of the generators is changed to c,s instead, the code works as expected.
I've been generally running into a lot of trouble with this function claiming the system is not zero-dimensional or returning no solutions when there are solutions for systems generated by intersecting a circle with a random affine line like the one above.
The text was updated successfully, but these errors were encountered: