Hello, I encountered the following error with Sage 7.2 (running Pari 2.8.0 - development git 61b65cc) on 64-bit Linux:
sage:
factor(10168294586031875795858081963666896068309279739028149015822994808470409*x^40
+ 54637623902063579051764164063785862273611
227214908262321295448827189084*x^39 +
70277533614869515097962450874425962549220558005689745585991279744832844*x^38
- 8376262189094092
345214251474588345374416636909696030473181959426070676*x^37 +
242436471429285303967664599557556624957448729852485841450904112819960*x
^36 -
16467308196223131257423201642300320128598078516427018512063667037294*x^35 +
994670440764331398827492095544600212856465140498468
589958417872072*x^34 +
1669447355892151599814976161752696223640623068834337285730125664*x^33 +
47878211805201490953500075842624444622
1319502304181118993796382*x^32 -
11027053324867589949389559654243101141839981063096910359753562*x^31 -
650300612550426855523171806453
147668269344451836561774976635*x^30 +
1570684208380806232824734619938834236093091252040919342718*x^29 +
34145278060467227277351650969
0661978537006252836870160018*x^28 -
966249275643573722337167379239879222114833277088010990*x^27 -
11930600081984079242665842322914049
3263927594569597794*x^26 +
1713994524513493848614719052055059346207593390535634*x^25 +
4235662454113565448888950245312909262295898883
3235*x^24 - 996720692002969642895416146776609782878203402522*x^23 -
13192515526064968724761611985177806405225265171*x^22 + 3432294679
31700304107325597744650482081790788*x^21 +
2664889270854122998123116283199878096309356*x^20 -
983892014425892473122607648686510399940
82*x^19 - 150060150030105789105698925850838409429*x^18 +
25211033715397630291436024039238871516*x^17 - 870382582902471525325447717394
57079*x^16 - 5390615133300676018637015393468804*x^15 +
29306521731548203900579062036148*x^14 + 898735965174867625096564403276*x^13
-
4155515638403008599438631634*x^12 - 111947828289343132957097300*x^11 +
236978606843233982987143*x^10 + 9962170815201528010504*x^9 + 1
2452485804904321512*x^8 - 546338689137860232*x^7 - 2961407384736921*x^6 +
7430129362962*x^5 + 136659878406*x^4 + 634377186*x^3 + 1511
809*x^2 + 1898*x + 1)
(100837962028354558549359124161681597*x^20 +
270917930127842441064021001407727686*x^19 -
15465362385090532027681012126697308*x^18 + 1
6986375609134831570968842775350*x^17 -
29476942979162692929636792384016*x^16 + 147550926306681751374468608357*x^15
+ 1333976962730521
5432298110948*x^14 + 33478195196624326088400980*x^13 -
3231442293652186023750739*x^12 + 24170655490731720935025*x^11 +
86595567527548
7498612*x^10 - 20060803083331090107*x^9 - 243338441403803107*x^8 +
5265846272236422*x^7 + 34263040880952*x^6 - 648382533165*x^5 - 415
1669958*x^4 + 27170397*x^3 + 305604*x^2 + 949*x + 1)^2
sage: factor(256*x^40 + 1024*x^35 + 1792*x^30 + 1792*x^25 + 1120*x^20 +
448*x^15 + 112*x^10 + 16*x^5 + 1)
---------------------------------------------------------------------------
PariError Traceback (most recent call last)
<ipython-input-2-74fc882135d8> in <module>()
----> 1 factor(Integer(256)*x**Integer(40) + Integer(1024)*x**Integer(35) +
Integer(1792)*x**Integer(30) + Integer(1792)*x**Integer(2
5) + Integer(1120)*x**Integer(20) + Integer(448)*x**Integer(15) + Integer(
112)*x**Integer(10) + Integer(16)*x**Integer(5) + Integer(1
))
/usr/lib/python2.7/site-packages/sage/arith/misc.pyc in factor(n, proof,
int_, algorithm, verbose, **kwds)
2476 # Just in case factor method doesn't have a proof option.
2477 try:
-> 2478 return n.factor(**kwds)
2479 except AttributeError:
2480 raise TypeError("unable to factor n")
sage/symbolic/expression.pyx in sage.symbolic.expression.Expression.factor
(/build/sagemath/src/sage-7.2/src/build/cythonized/sage/sy
mbolic/expression.cpp:50369)()
sage/rings/polynomial/polynomial_element.pyx in
sage.rings.polynomial.polynomial_element.Polynomial.factor
(/build/sagemath/src/sage-
7.2/src/build/cythonized/sage/rings/polynomial/polynomial_element.c:33493)()
/usr/lib/python2.7/site-packages/sage/rings/rational_field.pyc in
_factor_univariate_polynomial(self, f)
1275
1276 """
-> 1277 G = list(f._pari_with_name().factor())
1278
1279 # normalize the leading coefficients
sage/libs/pari/gen.pyx in sage.libs.pari.gen.gen.factor
(/build/sagemath/src/sage-7.2/src/build/cythonized/sage/libs/pari/gen.c:12273
8)()
sage/libs/pari/handle_error.pyx in sage.libs.pari.handle_error._pari_err_handle
(/build/sagemath/src/sage-7.2/src/build/cythonized/sa
ge/libs/pari/handle_error.c:3074)()
PariError: bug in gerepile, significant pointers lost, please report
If I run only the second command, sage outputs the correct answer (
(2x^5+1)^8 ).
If I factor these two polynomials directly in the GP/PARI Calculator (using
the same commands as in sage), it works, too.
Best,
Fabian
--
You received this message because you are subscribed to the Google Groups
"sage-devel" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-devel+unsubscr...@googlegroups.com.
To post to this group, send email to sage-devel@googlegroups.com.
Visit this group at https://groups.google.com/group/sage-devel.
For more options, visit https://groups.google.com/d/optout.
factorizations.sage
Description: Binary data
