Opened 2 years ago

Last modified 8 weeks ago

#29671 new defect

Minimal associated primes for polynomial ideals can hang unexpectedly

Reported by: kedlaya Owned by:
Priority: major Milestone: sage-9.7
Component: commutative algebra Keywords: minimal associated primes, singular
Cc: Merged in:
Authors: Reviewers:
Report Upstream: Reported upstream. No feedback yet. Work issues:
Branch: Commit:
Dependencies: Stopgaps:

Status badges

Description

This might be an upstream issue. I don't know if finding a more minimal example would help to isolate this, but note that this is only in two variables over a quadratic field, and the degrees are not absurdly large.

sage: L.<zeta3> = CyclotomicField(3)                                                        
sage: Q.<x5,x6> = L[] 
sage: I = Q.ideal((x5*x6^16 - 3*x5*x6^15 + 3*x5*x6^14 + x5*x6^13 - 6*x5*x6^12 + 7*x5*x6^11 - 3*x5*x6^10 - 3*x5*x6^9 + 7*x5*x6^8 - 6*x5*x6^7 + x5*x6^6 + 3*x5*x6^5 - 3*x5*x6^4 + x5*x6^3 + (-zeta3)*x6^13 + (3*zeta3)*x6^12 + (-3*zeta3)*x6^11 + (-zeta3)*x6^10 + (6*zeta3)*x6^9 + (-7*zeta3)*x6^8 + (3*zeta3)*x6^7 + (3*zeta3)*x6^6 + (-7*zeta3)*x6^5 + (6*zeta3)*x6^4 + (-zeta3)*x6^3 + (-3*zeta3)*x6^2 + (3*zeta3)*x6 + (-zeta3), x5^2*x6^5 - x5^2*x6^3 + (3*zeta3)*x5*x6^15 + (-7*zeta3)*x5*x6^14 + (4*zeta3)*x5*x6^13 + (7*zeta3)*x5*x6^12 + (-14*zeta3)*x5*x6^11 + (10*zeta3)*x5*x6^10 + (-10*zeta3)*x5*x6^8 + (14*zeta3)*x5*x6^7 + (-8*zeta3)*x5*x6^6 + (-4*zeta3)*x5*x6^5 + (8*zeta3)*x5*x6^4 + (-3*zeta3)*x5*x6^3 + (-zeta3)*x5*x6^2 + (zeta3)*x5 + (3*zeta3 + 3)*x6^12 + (-7*zeta3 - 7)*x6^11 + (4*zeta3 + 4)*x6^10 + (7*zeta3 + 7)*x6^9 + (-14*zeta3 - 14)*x6^8 + (10*zeta3 + 10)*x6^7 + (-10*zeta3 - 10)*x6^5 + (14*zeta3 + 14)*x6^4 + (-8*zeta3 - 8)*x6^3 + (-4*zeta3 - 4)*x6^2 + (8*zeta3 + 8)*x6 + (-3*zeta3 - 3), x5^4*x6^4 - x5^4*x6^3 + (2*zeta3)*x5^3*x6^4 + (-2*zeta3)*x5^3*x6^3 + (-zeta3)*x5^3*x6 + (zeta3)*x5^3 + (-2*zeta3 - 2)*x5^2*x6^4 + (2*zeta3 + 2)*x5^2*x6^3 + (2*zeta3 + 2)*x5^2*x6 + (-2*zeta3 - 2)*x5^2 + 7*x5*x6^15 - 18*x5*x6^14 + 11*x5*x6^13 + 18*x5*x6^12 - 36*x5*x6^11 + 25*x5*x6^10 - 25*x5*x6^8 + 35*x5*x6^7 - 19*x5*x6^6 - 11*x5*x6^5 + 20*x5*x6^4 - 7*x5*x6^3 - 2*x5*x6 + 2*x5 + (-7*zeta3)*x6^12 + (18*zeta3)*x6^11 + (-11*zeta3)*x6^10 + (-18*zeta3)*x6^9 + (36*zeta3)*x6^8 + (-25*zeta3)*x6^7 + (25*zeta3)*x6^5 + (-35*zeta3)*x6^4 + (19*zeta3)*x6^3 + (11*zeta3)*x6^2 + (-20*zeta3)*x6 + (7*zeta3)))

Now lots of things behave normally and run instantly:

sage: gb = I.groebner_basis()
sage: l = I.associated_primes()

But try:

sage: l = I.minimal_associated_primes()

and Sage hangs, apparently within sage.libs.singular.function_factory.ff.primdec__lib.minAssGTZ. (Sage reports that the Singular version is 4.1.1p2.p0.)

Change History (8)

comment:1 Changed 2 years ago by kedlaya

Looks like this is an upstream issue:

sage: singular(I).minAssGTZ()

also hangs, whereas similar commands like

sage: singular(I).primdecGTZ()

work just fine.

comment:2 Changed 2 years ago by kedlaya

  • Keywords singular added

comment:3 Changed 2 years ago by kedlaya

  • Report Upstream changed from N/A to Reported upstream. No feedback yet.

comment:4 Changed 20 months ago by mkoeppe

  • Milestone changed from sage-9.2 to sage-9.3

comment:5 Changed 14 months ago by mkoeppe

  • Milestone changed from sage-9.3 to sage-9.4

Moving to 9.4, as 9.3 has been released.

comment:6 Changed 10 months ago by mkoeppe

  • Milestone changed from sage-9.4 to sage-9.5

comment:7 Changed 6 months ago by mkoeppe

  • Milestone changed from sage-9.5 to sage-9.6

comment:8 Changed 8 weeks ago by mkoeppe

  • Milestone changed from sage-9.6 to sage-9.7
Note: See TracTickets for help on using tickets.