Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
computing the iterated exp(x)-1
#12
Hi, this is actuallly in reference to jaydfox's posts about Sage and the like. I rewrote the Mathematica code for the custom IDM method, and it works better than what you would usually do in Sage to get the same result, but it still not as fast as Mathematica (which doesn't surprise me, given that most of the work is being done behind the scenes by Maxima, which Sage gives a front-end to). So jaydfox, I would also like to correct you on your PARI vs. Sage question. Sage provides a front-end for PARI as well, so it's not a question of which is better. Its like saying "ROM vs. Computer" which is better? obviously they're not quite comparable.

Anyways, without further delay I suppose I should give you the Sage code:

Code:
# Define the custom IDM:
def parabolic_idm(expr, x, size):
  ret = []
  ser = x
  for i in xrange(size):
    ser = taylor(ser(expr), x, 0, size)
    cof = ser.coeffs(x)
    ret.append([0] + [cof[int(k)][0] for k in xrange(len(cof))])
  top = [0, 1] + [0 for i in xrange(size-1)]
  return [top] + ret

# Get the flow from custom IDM:
def parabolic_flow_coeff(idm_matrix, t):
  size = len(idm_matrix) - 1
  ret = [sum([(-1)^(n-k-1)*idm_matrix[k][n]*
    binomial(t, k)*binomial(t-k-1, n-k-1)
    for k in xrange(n)]) for n in xrange(size)]
  return ret

# Get just the coefficients of (x,t):
def parabolic_flow_matrix(flow_coeff, t):
  size = len(flow_coeff)
  ret = []
  for i in xrange(size):
    if (i < 2):
      ret.append([flow_coeff[i]])
    else:
      cof = flow_coeff[i].expand().coeffs(t)
      ret.append([0] + [cof[int(k)][0] for k in xrange(len(cof))])
  return ret

# Use the custom IDM:
from sage.calculus.calculus import SymbolicVariable
x = SymbolicVariable('x')
t = SymbolicVariable('t')
size = 20
mtx = parabolic_idm(exp(x)-1, x, size)
fco = parabolic_flow_coeff(mtx, t)
fma = parabolic_flow_matrix(fco, t)

The only way I can see making this even faster is to implement it in Maxima directly, since thats the CAS that Sage uses behind the scenes to calculate this kind of thing. Maybe my next post will include a Maxima version. For more information on the CASs that Sage uses, see http://sagemath.org.

Andrew Robbins
Reply


Messages In This Thread
computing the iterated exp(x)-1 - by Daniel - 08/13/2007, 10:47 PM
RE: computing the iterated exp(x)-1 - by andydude - 08/16/2007, 01:28 AM
RE: computing the iterated exp(x)-1 - by jaydfox - 08/16/2007, 06:51 AM
RE: computing the iterated exp(x)-1 - by bo198214 - 08/16/2007, 07:48 AM
RE: computing the iterated exp(x)-1 - by andydude - 08/17/2007, 08:44 PM
RE: Iterability of exp(x)-1 - by bo198214 - 08/13/2007, 10:50 PM
RE: Iterability of exp(x)-1 - by Daniel - 08/14/2007, 06:51 PM
RE: Iterability of exp(x)-1 - by jaydfox - 08/14/2007, 01:01 AM
RE: Iterability of exp(x)-1 - by Gottfried - 08/14/2007, 12:45 PM
RE: Iterability of exp(x)-1 - by bo198214 - 08/14/2007, 04:11 PM
RE: Iterability of exp(x)-1 - by Gottfried - 08/14/2007, 04:35 PM
RE: Iterability of exp(x)-1 - by jaydfox - 08/14/2007, 02:42 AM
RE: Iterability of exp(x)-1 - by Gottfried - 08/14/2007, 03:08 AM
RE: Iterability of exp(x)-1 - by jaydfox - 08/14/2007, 05:09 AM
RE: Iterability of exp(x)-1 - by Gottfried - 08/14/2007, 05:09 AM

Possibly Related Threads...
Thread Author Replies Views Last Post
  iterated derivation Xorter 0 865 06/09/2019, 09:43 PM
Last Post: Xorter
  1st iterated derivatives and the tetration of 0 Xorter 0 1,666 05/12/2018, 12:34 PM
Last Post: Xorter
  Iterated nand Xorter 2 4,310 03/27/2017, 06:51 PM
Last Post: Xorter
  Iterated compositions Xorter 0 1,899 08/20/2016, 01:19 PM
Last Post: Xorter
  Grzegorczyk hierarchy vs Iterated differential equations? MphLee 0 2,407 01/03/2015, 11:02 PM
Last Post: MphLee
  Iterated polynomials JmsNxn 4 8,198 12/16/2010, 09:00 PM
Last Post: JmsNxn
  computing f(f(x)) = exp(x) with g(g(x)) = - exp(x) tommy1729 2 4,815 06/06/2009, 10:28 PM
Last Post: tommy1729
  The fractal nature of iterated ln(x) [Bandwidth warning: lots of images!] jaydfox 16 20,525 09/09/2007, 01:21 AM
Last Post: jaydfox
  Bell formula for iterated exponentiation bo198214 1 4,461 08/26/2007, 12:27 PM
Last Post: Gottfried



Users browsing this thread: 1 Guest(s)