Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
SAGE code for computing flow matrix for exp(z)-1
#4
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.
Reply


Messages In This Thread
RE: SAGE code for computing flow matrix for exp(z)-1 - by bo198214 - 08/21/2009, 12:38 PM

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



Users browsing this thread: 1 Guest(s)