>
> Incidentally I observe that Sympy has the same behavior, so we can't 
> just nick their factoring algorithm -- maybe some other package we can 
> try the same example to see if any of them handle it quickly? 
>
> best 
>
> Robert Dodier 
>
> Hi, I have tried the factorization with giacpy. (cf trac 12375). I had to 
expexpand first before factoring and did this:


sage: from giacpy import libgiac
sage: x=libgiac('x')
sage: s=exp(1024*(x+1))-1
sage: %time s.expexpand().factor()
CPU times: user 0 ns, sys: 0 ns, total: 0 ns
Wall time: 1.32 ms
(exp(x+1)-1)*(exp(x+1)+1)*(exp(x+1)^2+1)*(exp(x+1)^4+1)*(exp(x+1)^8+1)*(exp(x+1)^16+1)*(exp(x+1)^32+1)*(exp(x+1)^64+1)*(exp(x+1)^128+1)*(exp(x+1)^256+1)*(exp(x+1)^512+1)

or

sage: %time s.expexpand().factor().expexpand()
CPU times: user 4 ms, sys: 0 ns, total: 4 ms
Wall time: 1.13 ms
(exp(x)*exp(1)-1)*(exp(x)*exp(1)+1)*((exp(x)*exp(1))^2+1)*((exp(x)*exp(1))^4+1)*((exp(x)*exp(1))^8+1)*((exp(x)*exp(1))^16+1)*((exp(x)*exp(1))^32+1)*((exp(x)*exp(1))^64+1)*((exp(x)*exp(1))^128+1)*((exp(x)*exp(1))^256+1)*((exp(x)*exp(1))^512+1)

best

Frederic

-- 
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 http://groups.google.com/group/sage-devel.
For more options, visit https://groups.google.com/d/optout.

Reply via email to