• 0 Vote(s) - 0 Average
• 1
• 2
• 3
• 4
• 5
 SAGE code for computing flow matrix for exp(z)-1 bo198214 Administrator Posts: 1,389 Threads: 90 Joined: Aug 2007 08/21/2009, 12:38 PM You may also use the powerseries package that was an outcome/branch of some effort of Andrew and me to put together sage routines for hyperoperations. I recently adapted the code to work with sage 4.1.1: http://github.com/bo198214/hyperops/raw/...rseries.py You would use it for the purpose of iterating exp(x)-1 like: Code:```---------------------------------------------------------------------- | Sage Version 4.1.1, Release Date: 2009-08-14                       | | Type notebook() for the GUI, and license() for information.        | ---------------------------------------------------------------------- In [2]: load formal_powerseries.py In [3]: P = FormalPowerSeriesRing(QQ) In [4]: w = PolynomialRing(QQ,'w')([0,1]) In [5]: time sum(P.Dec_exp.it(w)[0:100]); CPU times: user 326.53 s, sys: 1.49 s, total: 328.02 s Wall time: 345.83 s In [7]: P.Dec_exp.it(w) Out[7]: [0, 1, 1/2*w, 1/4*w^2 - 1/12*w, 1/8*w^3 - 5/48*w^2 + 1/48*w, ...]```(This is on a Athlon XP 2500, it consumes roughly 200-300 MB memory.) "it(w)" does just what one would expect: it iterates the given powerseries by w. « Next Oldest | Next Newest »

 Messages In This Thread SAGE code for computing flow matrix for exp(z)-1 - by jaydfox - 08/19/2009, 08:39 PM RE: SAGE code for computing flow matrix for exp(z)-1 - by jaydfox - 08/19/2009, 09:17 PM RE: SAGE code for computing flow matrix for exp(z)-1 - by jaydfox - 08/19/2009, 11:47 PM RE: SAGE code for computing flow matrix for exp(z)-1 - by bo198214 - 08/21/2009, 12:38 PM RE: SAGE code for computing flow matrix for exp(z)-1 - by jaydfox - 08/21/2009, 05:32 PM

 Possibly Related Threads... Thread Author Replies Views Last Post Revisting my accelerated slog solution using Abel matrix inversion jaydfox 21 12,277 02/09/2019, 02:25 PM Last Post: sheldonison C++ code for tet, ate and hexp MorgothV8 0 3,255 07/10/2014, 04:24 PM Last Post: MorgothV8 "Kneser"/Riemann mapping method code for *complex* bases mike3 2 7,198 08/15/2011, 03:14 PM Last Post: Gottfried An incremental method to compute (Abel) matrix inverses bo198214 3 9,740 07/20/2010, 12:13 PM Last Post: Gottfried Sage Question? rsgerard 1 4,312 05/09/2010, 11:40 AM Last Post: bo198214 Single-exp series computation code mike3 0 3,145 04/20/2010, 08:59 PM Last Post: mike3 computing teh last digits without computing the number deepinlife 3 6,248 02/24/2009, 09:09 AM Last Post: deepinlife How to free memory in SAGE? jaydfox 2 5,893 12/21/2007, 06:38 PM Last Post: andydude Convergence of matrix solution for base e jaydfox 6 9,675 12/18/2007, 12:14 AM Last Post: jaydfox Computing Abel function at a given center jaydfox 10 13,387 11/30/2007, 06:44 PM Last Post: andydude

Users browsing this thread: 1 Guest(s)