factorization with GiNaC (update)
Parisse Bernard
Bernard.Parisse at ujf-grenoble.fr
Fri Sep 22 14:36:40 CEST 2000
Hello,
I have updated the polynomial package over GiNaC, especially
factorization:
ftp://fourier.ujf-grenoble.fr/pub/hp48/giac.tgz
Main changes:
* 2 version are available: you can either use giac with the NTL
univariate
factorization library or standalone
* the code is faster.
Installation requires gmp, cln and GiNaC, optionnaly NTL.
Some timimgs of the standalone version on a Celeron 533 Mhz (# Time 0
means <=0.01 sec)
(c+x+b)*(b^2-c*x+x^2-c*b+c^2-x*b) # Time0
(x+y+z)*(x-y+z)*(x-y-z) # Time0.01
(-5+3*x+y-y^2)*(-1+3*x+x*y-y^2) # Time0.01
(x-3*y+y*z)*(2*x-z+x*y)*(1+x*z+x*y+6*x^3) # Time0.02
(x*y+y^2+x^2)*(-x*y+y^2+x^2)*(x-y)*(-x-y) # Time0.01
(2+x^3+13*x*y*z^2-21*x*y^3*z)*(1+x*z+x*y+6*x^3) # Time0.01
2+x^3+13*x*y*z^2-21*x*y^3*z # Time0.01
x^2*y^2*(-x+y)*z^5*(-x-y-z) # Time0.01
-x^3*y^2*z*(x-y-z)*(-x+y-z)*(-x-y-z) # Time0.01
(8+15*x)*(1+8*x)*(11+5*x)*(-31+30*x)*(-13+12*x)*(-209+120*x) # Time0.01
(1+x^5+x^20+x^10+x^15)*(1-x^5+x^20+x^10-x^15) # Time0.02
(1+x^35+x^5-x^20+x^40-x^25-x^15)*(1-x^35-x^5-x^20+x^40+x^25+x^15)*(1+x-x^4+x^8-x^5-x^3+x^7)*(1-x-x^4+x^8+x^5+x^3-x^7)*(1+x+x^2)*(1-x+x^2)
# Time0.13
If you are factoring big polynomials, the NTL version should be used
instead of the
standalone version (since it does not include e.g. FFT polynomial
multiplication ...)
Bernard Parisse
More information about the GiNaC-devel
mailing list