Skip to content
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

SageMath 10.3 Qepcad crashed -- automatically restarting #38302

Closed
2 tasks done
anirjoshi opened this issue Jun 28, 2024 · 7 comments
Closed
2 tasks done

SageMath 10.3 Qepcad crashed -- automatically restarting #38302

anirjoshi opened this issue Jun 28, 2024 · 7 comments
Labels

Comments

@anirjoshi
Copy link

anirjoshi commented Jun 28, 2024

Steps To Reproduce

  1. I installed Sage-Math 10.3 for my Mac M2 from https://github.com/3-manifolds/Sage_macOS/releases/
  2. I tried running a program containing QEPCAD
  3. Sage outputs Qepcad crashed -- automatically restarting
  4. Previously asked about this, but the response was that this is resolved in 10.3 previous issue
  5. The input program is the following: Note that even without memcells='10000000 +L5000' it crashes.
var('x_5_0')
var('x_5_1')
var('d')
var('x_0_0')
var('x_0_1')
var('x_1_0')
var('x_1_1')
var('x_2_0')
var('x_2_1')
var('x_3_0')
var('x_3_1')
var('x_4_0')
var('x_4_1')
var('tmp')
qf = qepcad_formula
constr0 = (tmp - tmp +(Rational(-1/1)* (d)) + (Rational(1/1)* (x_0_0^2)) + (Rational(-3/4))< 0)
constr1 = (tmp - tmp +(Rational(-1/2)) + (Rational(-1/1)* (d))< 0)
constr2 = (tmp - tmp +(Rational(1/1)) + (Rational(-1/1)* (d))< 0)
constr3 = (tmp - tmp +(Rational(-1/2)) + (Rational(-1/1)* (d))< 0)
constr4 = (tmp - tmp +(Rational(1/4)) + (Rational(-1/1)* (d))< 0)
constr5 = (tmp - tmp +(Rational(-1/1)* (d)) + (Rational(-1/1)) + (Rational(1/1)* (x_5_1^2)) + (Rational(1/1)* (x_5_0^2))< 0)
constr6 = (tmp - tmp +(Rational(-1/1)* (x_0_0)) + (Rational(-1/1)* (d)) + (Rational(-1/1)* (x_0_0^2)) + (Rational(-1/4))< 0)
constr7 = (tmp - tmp +(Rational(2/1)* (x_0_0)) + (Rational(-1/1)* (d)) + (Rational(-1/1)* (x_0_0^2)) + (Rational(-9/4))< 0)
constr8 = (tmp - tmp +(Rational(1/1)* (x_0_0)) + (Rational(-1/1)* (d)) + (Rational(-1/1)* (x_0_0^2)) + (Rational(-1/4))< 0)
constr9 = (tmp - tmp +(Rational(-1/1)* (x_0_0)) + (Rational(-1/1)* (d)) + (Rational(-1/1)* (x_0_0^2)) + (Rational(-3/2))< 0)
constr10 = (tmp - tmp +(Rational(1/1)* (x_5_1)) + (Rational(-1/1)* (x_5_0^2)) + (Rational(-1/1)* (d)) + (Rational(3/4)) + (Rational(2/1)* (x_0_0)* (x_5_0)) + (Rational(-1/1)* (x_5_1^2)) + (Rational(-1/1)* (x_0_0^2))< 0)
constr11 = (tmp - tmp +(Rational(-3/2)) + (Rational(-1/1)* (d))< 0)
constr12 = (tmp - tmp +(Rational(0/1)) + (Rational(-1/1)* (d))< 0)
constr13 = (tmp - tmp +(Rational(3/4)) + (Rational(-1/1)* (d))< 0)
constr14 = (tmp - tmp +(Rational(-1/1)* (x_5_1)) + (Rational(-1/1)* (x_5_1^2)) + (Rational(1/2)) + (Rational(-1/1)* (x_5_0)) + (Rational(-1/1)* (x_5_0^2)) + (Rational(-1/1)* (d))< 0)
constr15 = (tmp - tmp +(Rational(1/2)) + (Rational(-1/1)* (d))< 0)
constr16 = (tmp - tmp +(Rational(-5/4)) + (Rational(-1/1)* (d))< 0)
constr17 = (tmp - tmp +(Rational(-2/1)* (x_5_1)) + (Rational(-1/1)* (x_5_1^2)) + (Rational(-1/1)) + (Rational(2/1)* (x_5_0)) + (Rational(-1/1)* (x_5_0^2)) + (Rational(-1/1)* (d))< 0)
constr18 = (tmp - tmp +(Rational(-1/4)) + (Rational(-1/1)* (d))< 0)
constr19 = (tmp - tmp +(Rational(-1/1)* (x_5_1)) + (Rational(-1/1)* (x_5_1^2)) + (Rational(1/2)) + (Rational(1/1)* (x_5_0)) + (Rational(-1/1)* (x_5_0^2)) + (Rational(-1/1)* (d))< 0)
constr20 = (tmp - tmp +(Rational(-2/1)* (x_5_1)) + (Rational(-1/1)* (x_5_1^2)) + (Rational(-1/4)) + (Rational(-1/1)* (x_5_0)) + (Rational(-1/1)* (x_5_0^2)) + (Rational(-1/1)* (d))< 0)
constr21 = (tmp - tmp +(Rational(-1/1)* (d))< 0)
F = qf.and_(constr0, constr1, constr2, constr3, constr4, constr5, constr6, constr7, constr8, constr9, constr10, constr11, constr12, constr13, constr14, constr15, constr16, constr17, constr18, constr19, constr20, constr21)
print(qepcad(qf.exists(x_0_0,F) ,memcells='10000000 +L5000'))

Expected Behavior

It should output the formula the we get after quantifier elimination.

Actual Behavior

Sage crashes and outputs Qepcad crashed -- automatically restarting

Additional Information

Similar issues were raised previously in a lot of places including:
link 1
link 2
link 3
link 4

Environment

- **OS**: Mac M2 Ventura 13.4
- **Sage Version**: 10.3

Checklist

  • I have searched the existing issues for a bug report that matches the one I want to file, without success.
  • I have read the documentation and troubleshoot guide
@maxale
Copy link
Contributor

maxale commented Jun 29, 2024

This can be marked as a duplicate of #37365
A close look reveals that QEPCAD in this case complains Reason for the failure: Too few cells reclaimed. and the issue is resolved by adding an extra 0 to the memecells parameter, that is, memcells='100000000 +L5000'

@anirjoshi
Copy link
Author

anirjoshi commented Jun 29, 2024

@maxale Thank you for the pointer. But it crashes even after adding an extra 0 in the memcells parameter.

@maxale
Copy link
Contributor

maxale commented Jun 29, 2024

It works fine at Sagecell with memcells='100000000 +L5000':
https://sagecell.sagemath.org/?q=uivsfv

@anirjoshi
Copy link
Author

@maxale Thank you! Unfortunately, it still seems to crash in my system. However, the output d - 1 > 0 /\ x51^2 + x50^2 - d - 1 < 0 of Sagecell doesn't seem to make sense because there is no variable with the name x51 and x50.

@maxale
Copy link
Contributor

maxale commented Jun 29, 2024

Whether or not it makes sense is a different story. The point is that it works fine at Sagecell. Try to debug your QEPCAD with .go() commands.

@anirjoshi
Copy link
Author

The output has a variable that is not present in the input. Probably, I should file it in another issue? Anyways @maxale Thank you for the quick response.

@DaveWitteMorris
Copy link
Member

As explained above, this can be closed as a duplicate of #37365

@mkoeppe mkoeppe closed this as not planned Won't fix, can't repro, duplicate, stale Jul 8, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

4 participants