8Fermer10
PpHdLe 10/05/2008 à 21:45
squalyl (./5) :

It's up to us to create a useable software for calcs, now biggrin.gif


Some news?

By the way, some news. I have optimized a little bit MAY. Previous released branch:
MAY V0.7.0 (GMP V4.2.2 MPFR V2.3.1 CC=gcc CFLAGS=-O3 -fomit-frame-pointer -funroll-loops -ffast-math -march=pentium3 -fexceptions -ffunction-sections -fdata-sections)
eval (sum ai*ai*ai) - quite different - N=100......0.03ms
eval (sum ai*ai*ai) -  quite similar  - N=100......0.03ms
eval (sum ai*ai*ai) - quite different - N=1000......0.50ms
eval (sum ai*ai*ai) -  quite similar  - N=1000......0.50ms
eval (sum ai*ai*ai) - quite different - N=10000......1.71ms
eval (sum ai*ai*ai) -  quite similar  - N=10000......2.00ms
eval (sum ai*ai*ai) - quite different - N=100000......72.00ms
eval (sum ai*ai*ai) -  quite similar  - N=100000......64.00ms
eval (sum ai*ai*ai) - quite different - N=1000000......828.00ms
eval (sum ai*ai*ai) -  quite similar  - N=1000000......884.00ms
eval(x+f(x)+f(f(x))+...+f(5000)(x)), subs f to id...1716ms
expand ((a0+...a500)^2), replace a0, reeval...808ms
expand ((x0+...x2+1)^16*(1+(x0+...x2+1)^16))...1268ms
expand ((1+x)^800*(2+x)^800)...456ms
expand ((17+x)^400*(42+x)^400)...256ms
expand ((1+sqrt(5))^65000)...8ms
expand ((1+x+y)^500)...1645ms
divide ( (1+x)^1000+1 , (1-x)^500)...856ms
divide ( (1+x)^1000+1 , x^3-5*x+17)...256ms
divide ( (1+x+y^2)^50+1 , (1-x)^25+y)...1252ms
divide ( (1+x+y^2)^25+1 , x^3*y-5*x*y^42+17*y+1)...1972ms
gcd ( (1+2*x)^200*(x^3+2*x^2+1) , (1+2*x)^42*(x^3-2*x+42) )...100ms
gcd ( (1+2*x)^200*(x^3+2*x^2+1) , (1+2*x)^42*(x^3-2*x+42)+1 )...728ms
gcd ( (1+2*x+y)^100*(x^3+2*x^2*y+1) , (1+2*x+y)^42*(x^3-2*x+42) )...28030ms
gcd ( (x^2-y^2)*(a+b)^10 , (x-y)*(a-c)^10 )...0ms
gcd ( (x-y)^50+a , (x+y)^50 )...0ms
Total time 41202ms

Current branch:
MAY V0.7.1 (GMP V4.2.2 MPFR V2.3.1 CC=gcc CFLAGS=-O3 -fomit-frame-pointer -funroll-loops -ffast-math -march=pentium3 -fexceptions -ffunction-sections -fdata-sections)
eval (sum ai*ai*ai) - quite different - N=100......0.03ms
eval (sum ai*ai*ai) -  quite similar  - N=100......0.03ms
eval (sum ai*ai*ai) - quite different - N=1000......0.38ms
eval (sum ai*ai*ai) -  quite similar  - N=1000......0.36ms
eval (sum ai*ai*ai) - quite different - N=10000......2.29ms
eval (sum ai*ai*ai) -  quite similar  - N=10000......1.71ms
eval (sum ai*ai*ai) - quite different - N=100000......64.00ms
eval (sum ai*ai*ai) -  quite similar  - N=100000......52.00ms
eval (sum ai*ai*ai) - quite different - N=1000000......628.00ms
eval (sum ai*ai*ai) -  quite similar  - N=1000000......836.00ms
eval(x+f(x)+f(f(x))+...+f(5000)(x)), subs f to id...1620ms
expand ((a0+...a500)^2), replace a0, reeval...748ms
expand ((x0+...x2+1)^16*(1+(x0+...x2+1)^16))...864ms
expand ((1+x)^800*(2+x)^800)...424ms
expand ((17+x)^400*(42+x)^400)...256ms
expand ((1+sqrt(5))^65000)...4ms
expand ((1+x+y)^500)...1657ms
divide ( (1+x)^1000+1 , (1-x)^500)...724ms
divide ( (1+x)^1000+1 , x^3-5*x+17)...120ms
divide ( (1+x+y^2)^50+1 , (1-x)^25+y)...932ms
divide ( (1+x+y^2)^25+1 , x^3*y-5*x*y^42+17*y+1)...1352ms
gcd ( (1+2*x)^200*(x^3+2*x^2+1) , (1+2*x)^42*(x^3-2*x+42) )...108ms
gcd ( (1+2*x)^200*(x^3+2*x^2+1) , (1+2*x)^42*(x^3-2*x+42)+1 )...92ms
gcd ( (1+2*x+y)^100*(x^3+2*x^2*y+1) , (1+2*x+y)^42*(x^3-2*x+42) )...2912ms
gcd ( (x^2-y^2)*(a+b)^10 , (x-y)*(a-c)^10 )...0ms
gcd ( (x-y)^50+a , (x+y)^50 )...0ms
Total time 13396ms

Not too much regressions, so I think I have done a good work smile