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