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

update Singular spkg to 4.3.2+ #35676

Closed
1 task done
dimpase opened this issue May 23, 2023 · 24 comments · Fixed by #35934
Closed
1 task done

update Singular spkg to 4.3.2+ #35676

dimpase opened this issue May 23, 2023 · 24 comments · Fixed by #35934

Comments

@dimpase
Copy link
Member

dimpase commented May 23, 2023

Is there an existing issue for this?

  • I have searched the existing issues for a bug report that matches the one I want to file, without success.

Problem Description

currently we ship 4.3.1 with patches, the latest is 4.3.2p1

Experiments indicate that 4.3.2 will need work, though.

Proposed Solution

do the work

Alternatives Considered

n/a

Additional Information

No response

@antonio-rojas
Copy link
Contributor

I wasn't aware of tht 4.3.2p1 tag, as there is no corresponging tarballs in their server. Do you know if they're planning to release from github from now on, or is it just going to be random?

@antonio-rojas
Copy link
Contributor

Updated downstream, there are no new issues that need fixing in sagelib compared to the current 4.3.1p3 version.

The only problem is that they don't provide a source tarball, so it needs to be downloaded from github and docs need to be built from source. Depending on the texinfo version used at build time, this can cause issues when Sage tries to parse the info file.

I've opened #35678 to handle this.

@dimpase
Copy link
Member Author

dimpase commented May 24, 2023

I've asked whether there will be an official source tarball for 4.3.2p1: Singular/Singular#1175

@dimpase
Copy link
Member Author

dimpase commented Jun 5, 2023

@antonio-rojas
Copy link
Contributor

4.3.2p2 is broken

Singular/Singular#1176

@tornaria
Copy link
Contributor

After fixing 4.3.2p2 with Singular/Singular@8d54773, I get:

sage -t --random-seed=129530826087015995108246630021130802663 src/sage/rings/asymptotic/asymptotics_multivariate_generating_functions.py
**********************************************************************
File "src/sage/rings/asymptotic/asymptotics_multivariate_generating_functions.py", line 1254, in sage.rings.asymptotic.asymptotics_multivariate_generating_functions.FractionWithFactoredDenominator.leinartas_decomposition
Failed example:
    decomp
Expected:
    (0, []) +
    (-(x*y^2*sin(x) + x^2*y + x*y + y*sin(x) + x)*y, [(y, 1)]) +
    ((x*y^2*sin(x) + x^2*y + x*y + y*sin(x) + x)*x*y, [(x*y + 1, 1)]) +
    (x*y^2*sin(x) + x^2*y + x*y + y*sin(x) + x, [(y, 1), (x, 1)])
Got:
    (0, []) + (-(x*y^2*sin(x) + x^2*y + x*y + y*sin(x) + x)*x, [(x, 1)]) + ((x*y^2*sin(x) + x^2*y + x*y + y*sin(x) + x)*x*y, [(x*y + 1, 1)]) + (x*y^2*sin(x) + x^2*y + x*y + y*sin(x) + x, [(y, 1), (x, 1)])
**********************************************************************
File "src/sage/rings/asymptotic/asymptotics_multivariate_generating_functions.py", line 1613, in sage.rings.asymptotic.asymptotics_multivariate_generating_functions.FractionWithFactoredDenominator.?
Failed example:
    decomp = F.asymptotic_decomposition(alpha); decomp
Expected:
    (0, []) +
    (16*r*(3/x - 2/z) + 16/x - 16/z,
     [(x + 2*y + z - 4, 1), (2*x + y + z - 4, 1)])
Got:
    (0, []) + (16*r*(2/x - 1/y) + 32/3/x - 16/3/y, [(x + 2*y + z - 4, 1), (2*x + y + z - 4, 1)])
**********************************************************************
2 items had failures:
   1 of  96 in sage.rings.asymptotic.asymptotics_multivariate_generating_functions.FractionWithFactoredDenominator.?
   1 of  32 in sage.rings.asymptotic.asymptotics_multivariate_generating_functions.FractionWithFactoredDenominator.leinartas_decomposition
    [793 tests, 2 failures, 4.94 s]
----------------------------------------------------------------------
sage -t --random-seed=129530826087015995108246630021130802663 src/sage/rings/asymptotic/asymptotics_multivariate_generating_functions.py  # 2 doctests failed
----------------------------------------------------------------------

@kiwifb
Copy link
Member

kiwifb commented Jun 13, 2023

I can confirm 4.3.2_p1 does work and do not break any tests here. I probably will not get to play with _p2 if I wait for it in Gentoo main tree. If install is broken we may wait for _p3 or 4.3.3 whichever comes first.

@tornaria
Copy link
Contributor

I tracked the first doctest failure down to

sage: R.<x,y> = PolynomialRing(QQ)
sage: J = R.ideal([x, y, x*y + 1])
sage: R.one().lift(J)  # with singular 4.3.2p1
[-y, 0, 1]
sage: R.one().lift()  # with singular 4.3.2p2
[0, -x, 1]

The second doctest failure seems to arise from something similar since the output goes through a Leinartas decomposition at some point.

Easy fix:

--- a/src/sage/rings/asymptotic/asymptotics_multivariate_generating_functions.py
+++ b/src/sage/rings/asymptotic/asymptotics_multivariate_generating_functions.py
@@ -1251,7 +1251,7 @@ class FractionWithFactoredDenominator(RingElement):
             sage: H = R(f.denominator())
             sage: ff = FFPD(G, H.factor())
             sage: decomp = ff.leinartas_decomposition()
-            sage: decomp
+            sage: decomp  # random - non canonical depends on singular version
             (0, []) +
             (-(x*y^2*sin(x) + x^2*y + x*y + y*sin(x) + x)*y, [(y, 1)]) +
             ((x*y^2*sin(x) + x^2*y + x*y + y*sin(x) + x)*x*y, [(x*y + 1, 1)]) +
@@ -1611,9 +1611,7 @@ class FractionWithFactoredDenominator(RingElement):
             (-16, [(x + 2*y + z - 4, 1), (2*x + y + z - 4, 2)])
             sage: alpha = [3, 3, 2]
             sage: decomp = F.asymptotic_decomposition(alpha); decomp
-            (0, []) +
-            (16*r*(3/x - 2/z) + 16/x - 16/z,
-             [(x + 2*y + z - 4, 1), (2*x + y + z - 4, 1)])
+            (0, []) + (..., [(x + 2*y + z - 4, 1), (2*x + y + z - 4, 1)])
             sage: F1 = decomp[1]
             sage: p = {x: 1, y: 1, z: 1}
             sage: asy = F1.asymptotics(p, alpha, 2, verbose=True) # long time

Note that the first test then goes on to check the computed decomposition actually satisfies the definition. And the second test is for asymptotics which is still tested.

@kiwifb
Copy link
Member

kiwifb commented Jun 20, 2023

Really, in the first case since there is a check that the result is correct, we do not need to display and check that output. The fact that it is consistent should be enough. And I do not really think we should display decomp in the second case either, it is just an intermediate result.

@jhpalmieri
Copy link
Member

On OS X, homebrew's Singular 4.3.2p2 causes several doctest failures — see #35273. Has anyone tested out p1 or building your own p2 on OS X to see if those failures are present?

@mkoeppe
Copy link
Contributor

mkoeppe commented Jul 11, 2023

Has anyone tested out p1 or building your own p2 on OS X to see if those failures are present?

Upgrade in #35934 passes these tests on macOS

@antonio-rojas
Copy link
Contributor

4.3.2p3 is out, with one test regression

File "/usr/lib/python3.11/site-packages/sage/rings/polynomial/hilbert.pyx", line 584, in sage.rings.polynomial.hilbert.hilbert_poincare_series
Failed example:
    J.hilbert_numerator(algorithm='singular')
Expected:
    120*t^33 - 3465*t^32 + 48180*t^31 - 429374*t^30 + 2753520*t^29 - 13522410*t^28 + 52832780*t^27 - 168384150*t^26 + 445188744*t^25 - 987193350*t^24 + 1847488500*t^23 + 1372406746*t^22 - 403422496*t^21 - 8403314*t^20 - 471656596*t^19 + 1806623746*t^18 + 752776200*t^17 + 752776200*t^16 - 1580830020*t^15 + 1673936550*t^14 - 1294246800*t^13 + 786893250*t^12 - 382391100*t^11 + 146679390*t^10 - 42299400*t^9 + 7837830*t^8 - 172260*t^7 - 468930*t^6 + 183744*t^5 - 39270*t^4 + 5060*t^3 - 330*t^2 + 1  
Got:
    overflow at t^22
    overflow at t^21
    overflow at t^20
    overflow at t^19
    overflow at t^18
    120*t^33 - 3465*t^32 + 48180*t^31 - 429374*t^30 + 2753520*t^29 - 13522410*t^28 + 52832780*t^27 - 168384150*t^26 + 445188744*t^25 - 987193350*t^24 + 1847488500*t^23 + 1372406746*t^22 - 403422496*t^21 - 8403314*t^20 - 471656596*t^19 + 1806623746*t^18 + 752776200*t^17 + 752776200*t^16 - 1580830020*t^15 + 1673936550*t^14 - 1294246800*t^13 + 786893250*t^12 - 382391100*t^11 + 146679390*t^10 - 42299400*t^9 + 7837830*t^8 - 172260*t^7 - 468930*t^6 + 183744*t^5 - 39270*t^4 + 5060*t^3 - 330*t^2 + 1

I don't understand what this test is doing here anyway - it's testing a method from another class which is already tested in its class. Should we just drop it?

@tornaria
Copy link
Contributor

4.3.2p3 is out, with one test regression

File "/usr/lib/python3.11/site-packages/sage/rings/polynomial/hilbert.pyx", line 584, in sage.rings.polynomial.hilbert.hilbert_poincare_series
Failed example:
    J.hilbert_numerator(algorithm='singular')
Expected:
    120*t^33 - 3465*t^32 + 48180*t^31 - 429374*t^30 + 2753520*t^29 - 13522410*t^28 + 52832780*t^27 - 168384150*t^26 + 445188744*t^25 - 987193350*t^24 + 1847488500*t^23 + 1372406746*t^22 - 403422496*t^21 - 8403314*t^20 - 471656596*t^19 + 1806623746*t^18 + 752776200*t^17 + 752776200*t^16 - 1580830020*t^15 + 1673936550*t^14 - 1294246800*t^13 + 786893250*t^12 - 382391100*t^11 + 146679390*t^10 - 42299400*t^9 + 7837830*t^8 - 172260*t^7 - 468930*t^6 + 183744*t^5 - 39270*t^4 + 5060*t^3 - 330*t^2 + 1  
Got:
    overflow at t^22
    overflow at t^21
    overflow at t^20
    overflow at t^19
    overflow at t^18
    120*t^33 - 3465*t^32 + 48180*t^31 - 429374*t^30 + 2753520*t^29 - 13522410*t^28 + 52832780*t^27 - 168384150*t^26 + 445188744*t^25 - 987193350*t^24 + 1847488500*t^23 + 1372406746*t^22 - 403422496*t^21 - 8403314*t^20 - 471656596*t^19 + 1806623746*t^18 + 752776200*t^17 + 752776200*t^16 - 1580830020*t^15 + 1673936550*t^14 - 1294246800*t^13 + 786893250*t^12 - 382391100*t^11 + 146679390*t^10 - 42299400*t^9 + 7837830*t^8 - 172260*t^7 - 468930*t^6 + 183744*t^5 - 39270*t^4 + 5060*t^3 - 330*t^2 + 1

I don't understand what this test is doing here anyway - it's testing a method from another class which is already tested in its class. Should we just drop it?

It doesn't belong here, but it's catching some error. Indeed the coefficients for t^22 to t^18 are different between J.hilbert_numerator(algorithm='singular') and J.hilbert_numerator().

Based on

sage: factor(J.hilbert_numerator())
(t - 1)^30 * (120*t^3 + 135*t^2 + 30*t + 1)

vs. J.hilbert_numerator(algorithm='singular') being irreducible, I'm incline to believe the singular result is a bug. And it has been a bug since singular 4.2.1p2 (see 720d10e).

At least 4.3.2p2 is warning that there was overflow in those coefficients. Indeed the five incorrect coefficients differ by 2^32 from the correct ones:

sage: J.hilbert_numerator(algorithm='singular') - J.hilbert_numerator()
overflow at t^22
overflow at t^21
overflow at t^20
overflow at t^19
overflow at t^18
4294967296*t^22 - 4294967296*t^21 + 4294967296*t^20 - 4294967296*t^19 + 4294967296*t^18

@dimpase
Copy link
Member Author

dimpase commented Jul 17, 2023

it's insane that coefficients aren't arbitrary precision internally.
Or it could be a fallout from Sage using undocumented libsingular stuff - in such a case no promises on stability of the API...

@antonio-rojas
Copy link
Contributor

it's insane that coefficients aren't arbitrary precision internally. Or it could be a fallout from Sage using undocumented libsingular stuff - in such a case no promises on stability of the API...

No, the issue is in Singular:

> ring R=0,(x0,x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35,x36,x37,x38,x39,x40,x41,x42,x43),dp;
> ideal i=-x1*x11+x0*x12,-x2*x11+x0*x13,-x3*x11+x0*x14,-x4*x11+x0*x15,-x5*x11+x0*x16,-x6*x11+x0*x17,-x7*x11+x0*x18,-x8*x11+x0*x19,-x9*x11+x0*x20,-x10*x11+x0*x21,-x2*x12+x1*x13,-x3*x12+x1*x14,-x4*x12+x1*x15,-x5*x12+x1*x16,-x6*x12+x1*x17,-x7*x12+x1*x18,-x8*x12+x1*x19,-x9*x12+x1*x20,-x10*x12+x1*x21,-x3*x13+x2*x14,-x4*x13+x2*x15,-x5*x13+x2*x16,-x6*x13+x2*x17,-x7*x13+x2*x18,-x8*x13+x2*x19,-x9*x13+x2*x20,-x10*x13+x2*x21,-x4*x14+x3*x15,-x5*x14+x3*x16,-x6*x14+x3*x17,-x7*x14+x3*x18,-x8*x14+x3*x19,-x9*x14+x3*x20,-x10*x14+x3*x21,-x5*x15+x4*x16,-x6*x15+x4*x17,-x7*x15+x4*x18,-x8*x15+x4*x19,-x9*x15+x4*x20,-x10*x15+x4*x21,-x6*x16+x5*x17,-x7*x16+x5*x18,-x8*x16+x5*x19,-x9*x16+x5*x20,-x10*x16+x5*x21,-x7*x17+x6*x18,-x8*x17+x6*x19,-x9*x17+x6*x20,-x10*x17+x6*x21,-x8*x18+x7*x19,-x9*x18+x7*x20,-x10*x18+x7*x21,-x9*x19+x8*x20,-x10*x19+x8*x21,-x10*x20+x9*x21,-x1*x22+x0*x23,-x2*x22+x0*x24,-x3*x22+x0*x25,-x4*x22+x0*x26,-x5*x22+x0*x27,-x6*x22+x0*x28,-x7*x22+x0*x29,-x8*x22+x0*x30,-x9*x22+x0*x31,-x10*x22+x0*x32,-x2*x23+x1*x24,-x3*x23+x1*x25,-x4*x23+x1*x26,-x5*x23+x1*x27,-x6*x23+x1*x28,-x7*x23+x1*x29,-x8*x23+x1*x30,-x9*x23+x1*x31,-x10*x23+x1*x32,-x3*x24+x2*x25,-x4*x24+x2*x26,-x5*x24+x2*x27,-x6*x24+x2*x28,-x7*x24+x2*x29,-x8*x24+x2*x30,-x9*x24+x2*x31,-x10*x24+x2*x32,-x4*x25+x3*x26,-x5*x25+x3*x27,-x6*x25+x3*x28,-x7*x25+x3*x29,-x8*x25+x3*x30,-x9*x25+x3*x31,-x10*x25+x3*x32,-x5*x26+x4*x27,-x6*x26+x4*x28,-x7*x26+x4*x29,-x8*x26+x4*x30,-x9*x26+x4*x31,-x10*x26+x4*x32,-x6*x27+x5*x28,-x7*x27+x5*x29,-x8*x27+x5*x30,-x9*x27+x5*x31,-x10*x27+x5*x32,-x7*x28+x6*x29,-x8*x28+x6*x30,-x9*x28+x6*x31,-x10*x28+x6*x32,-x8*x29+x7*x30,-x9*x29+x7*x31,-x10*x29+x7*x32,-x9*x30+x8*x31,-x10*x30+x8*x32,-x10*x31+x9*x32,-x1*x33+x0*x34,-x2*x33+x0*x35,-x3*x33+x0*x36,-x4*x33+x0*x37,-x5*x33+x0*x38,-x6*x33+x0*x39,-x7*x33+x0*x40,-x8*x33+x0*x41,-x9*x33+x0*x42,-x10*x33+x0*x43,-x2*x34+x1*x35,-x3*x34+x1*x36,-x4*x34+x1*x37,-x5*x34+x1*x38,-x6*x34+x1*x39,-x7*x34+x1*x40,-x8*x34+x1*x41,-x9*x34+x1*x42,-x10*x34+x1*x43,-x3*x35+x2*x36,-x4*x35+x2*x37,-x5*x35+x2*x38,-x6*x35+x2*x39,-x7*x35+x2*x40,-x8*x35+x2*x41,-x9*x35+x2*x42,-x10*x35+x2*x43,-x4*x36+x3*x37,-x5*x36+x3*x38,-x6*x36+x3*x39,-x7*x36+x3*x40,-x8*x36+x3*x41,-x9*x36+x3*x42,-x10*x36+x3*x43,-x5*x37+x4*x38,-x6*x37+x4*x39,-x7*x37+x4*x40,-x8*x37+x4*x41,-x9*x37+x4*x42,-x10*x37+x4*x43,-x6*x38+x5*x39,-x7*x38+x5*x40,-x8*x38+x5*x41,-x9*x38+x5*x42,-x10*x38+x5*x43,-x7*x39+x6*x40,-x8*x39+x6*x41,-x9*x39+x6*x42,-x10*x39+x6*x43,-x8*x40+x7*x41,-x9*x40+x7*x42,-x10*x40+x7*x43,-x9*x41+x8*x42,-x10*x41+x8*x43,-x10*x42+x9*x43,-x12*x22+x11*x23,-x13*x22+x11*x24,-x14*x22+x11*x25,-x15*x22+x11*x26,-x16*x22+x11*x27,-x17*x22+x11*x28,-x18*x22+x11*x29,-x19*x22+x11*x30,-x20*x22+x11*x31,-x21*x22+x11*x32,-x13*x23+x12*x24,-x14*x23+x12*x25,-x15*x23+x12*x26,-x16*x23+x12*x27,-x17*x23+x12*x28,-x18*x23+x12*x29,-x19*x23+x12*x30,-x20*x23+x12*x31,-x21*x23+x12*x32,-x14*x24+x13*x25,-x15*x24+x13*x26,-x16*x24+x13*x27,-x17*x24+x13*x28,-x18*x24+x13*x29,-x19*x24+x13*x30,-x20*x24+x13*x31,-x21*x24+x13*x32,-x15*x25+x14*x26,-x16*x25+x14*x27,-x17*x25+x14*x28,-x18*x25+x14*x29,-x19*x25+x14*x30,-x20*x25+x14*x31,-x21*x25+x14*x32,-x16*x26+x15*x27,-x17*x26+x15*x28,-x18*x26+x15*x29,-x19*x26+x15*x30,-x20*x26+x15*x31,-x21*x26+x15*x32,-x17*x27+x16*x28,-x18*x27+x16*x29,-x19*x27+x16*x30,-x20*x27+x16*x31,-x21*x27+x16*x32,-x18*x28+x17*x29,-x19*x28+x17*x30,-x20*x28+x17*x31,-x21*x28+x17*x32,-x19*x29+x18*x30,-x20*x29+x18*x31,-x21*x29+x18*x32,-x20*x30+x19*x31,-x21*x30+x19*x32,-x21*x31+x20*x32,-x12*x33+x11*x34,-x13*x33+x11*x35,-x14*x33+x11*x36,-x15*x33+x11*x37,-x16*x33+x11*x38,-x17*x33+x11*x39,-x18*x33+x11*x40,-x19*x33+x11*x41,-x20*x33+x11*x42,-x21*x33+x11*x43,-x13*x34+x12*x35,-x14*x34+x12*x36,-x15*x34+x12*x37,-x16*x34+x12*x38,-x17*x34+x12*x39,-x18*x34+x12*x40,-x19*x34+x12*x41,-x20*x34+x12*x42,-x21*x34+x12*x43,-x14*x35+x13*x36,-x15*x35+x13*x37,-x16*x35+x13*x38,-x17*x35+x13*x39,-x18*x35+x13*x40,-x19*x35+x13*x41,-x20*x35+x13*x42,-x21*x35+x13*x43,-x15*x36+x14*x37,-x16*x36+x14*x38,-x17*x36+x14*x39,-x18*x36+x14*x40,-x19*x36+x14*x41,-x20*x36+x14*x42,-x21*x36+x14*x43,-x16*x37+x15*x38,-x17*x37+x15*x39,-x18*x37+x15*x40,-x19*x37+x15*x41,-x20*x37+x15*x42,-x21*x37+x15*x43,-x17*x38+x16*x39,-x18*x38+x16*x40,-x19*x38+x16*x41,-x20*x38+x16*x42,-x21*x38+x16*x43,-x18*x39+x17*x40,-x19*x39+x17*x41,-x20*x39+x17*x42,-x21*x39+x17*x43,-x19*x40+x18*x41,-x20*x40+x18*x42,-x21*x40+x18*x43,-x20*x41+x19*x42,-x21*x41+x19*x43,-x21*x42+x20*x43,-x23*x33+x22*x34,-x24*x33+x22*x35,-x25*x33+x22*x36,-x26*x33+x22*x37,-x27*x33+x22*x38,-x28*x33+x22*x39,-x29*x33+x22*x40,-x30*x33+x22*x41,-x31*x33+x22*x42,-x32*x33+x22*x43,-x24*x34+x23*x35,-x25*x34+x23*x36,-x26*x34+x23*x37,-x27*x34+x23*x38,-x28*x34+x23*x39,-x29*x34+x23*x40,-x30*x34+x23*x41,-x31*x34+x23*x42,-x32*x34+x23*x43,-x25*x35+x24*x36,-x26*x35+x24*x37,-x27*x35+x24*x38,-x28*x35+x24*x39,-x29*x35+x24*x40,-x30*x35+x24*x41,-x31*x35+x24*x42,-x32*x35+x24*x43,-x26*x36+x25*x37,-x27*x36+x25*x38,-x28*x36+x25*x39,-x29*x36+x25*x40,-x30*x36+x25*x41,-x31*x36+x25*x42,-x32*x36+x25*x43,-x27*x37+x26*x38,-x28*x37+x26*x39,-x29*x37+x26*x40,-x30*x37+x26*x41,-x31*x37+x26*x42,-x32*x37+x26*x43,-x28*x38+x27*x39,-x29*x38+x27*x40,-x30*x38+x27*x41,-x31*x38+x27*x42,-x32*x38+x27*x43,-x29*x39+x28*x40,-x30*x39+x28*x41,-x31*x39+x28*x42,-x32*x39+x28*x43,-x30*x40+x29*x41,-x31*x40+x29*x42,-x32*x40+x29*x43,-x31*x41+x30*x42,-x32*x41+x30*x43,-x32*x42+x31*x43;
> hilb(i);
// ** i is no standard basis
overflow at t^22
overflow at t^21
overflow at t^20
overflow at t^19
overflow at t^18
[...]
//  1806623746 t^18
//  -471656596 t^19
//  -8403314 t^20
//  -403422496 t^21
[...]

dimpase referenced this issue in Singular/Singular Jul 18, 2023
@dimpase
Copy link
Member Author

dimpase commented Jul 18, 2023

Is it reported upstream? In p3 there were some tweaks to hilb(), but apparently not enough to fix regression from p2.

@antonio-rojas
Copy link
Contributor

Is it reported upstream? In p3 there were some tweaks to hilb(), but apparently not enough to fix regression from p2.

Note that this has been the case since 4.2.1.p2 - Before, it aborted with an overflow error. From 4.2.1.p2 to 4.3.2.p2 it gave a wrong answer due to overflow. And now it still gives the same wrong answer, but warns about the overflow. I'm not sure this would be considered a bug upstream, given they explicitly changed the code to warn about the overflow (but still return a possibly wrong answer).

@dimpase
Copy link
Member Author

dimpase commented Jul 18, 2023

IMHO it should at least error out, not just warn. I expressed my opinion here: Singular/Singular@d1ba061
and would like to open a proper Issue - cause it's just silly of them...

@dimpase
Copy link
Member Author

dimpase commented Jul 18, 2023

See Singular/Singular#1135 for the issue tracking this example

@antonio-rojas
Copy link
Contributor

Regardless of the upstream report outcome, we need to make it clear that this is giving a wrong answer (and fix the test for singular 4.3.2.p3), done in #35977

@tornaria
Copy link
Contributor

I found two more regressions in singular 4.3.2p3:

sage -t --random-seed=46572117214005167953262287166108255484 src/sage/rings/polynomial/multi_polynomial_ideal.py
**********************************************************************
File "src/sage/rings/polynomial/multi_polynomial_ideal.py", line 2448, in sage.rings.polynomial.multi_polynomial_ideal.?.saturation
Failed example:
    I.saturation(J)
Expected:
    (Ideal (y, x^5) of Multivariate Polynomial Ring in x, y, z over Rational Field, 4)
Got:
    (Ideal (y, x^5) of Multivariate Polynomial Ring in x, y, z over Rational Field,
     0)
**********************************************************************
File "src/sage/rings/polynomial/multi_polynomial_ideal.py", line 2459, in sage.rings.polynomial.multi_polynomial_ideal.?.saturation
Failed example:
    I.saturation(other = J)
Expected:
    (Ideal (y, x^5) of Multivariate Polynomial Ring in x, y, z over Algebraic Field, 4)
Got:
    (Ideal (y, x^5) of Multivariate Polynomial Ring in x, y, z over Algebraic Field,
     0)
**********************************************************************

@antonio-rojas
Copy link
Contributor

I found two more regressions in singular 4.3.2p3:

Reported at Singular/Singular#1181

@antonio-rojas
Copy link
Contributor

I found two more regressions in singular 4.3.2p3:

Reported at Singular/Singular#1181

Fixed, unfortunately by changing the API. The corresponding Sage change is now #35980

@antonio-rojas
Copy link
Contributor

And opened #36018 to fix test failures in 4.3.2.p4 due to changes in hilb output.

@vbraun vbraun closed this as completed in bfb6b3c Sep 1, 2023
@mkoeppe mkoeppe added this to the sage-10.2 milestone Sep 1, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment