08/29/2007, 06:15 PM

Gottfried Wrote:jaydfox Wrote:Using gp, I computed Andrew's solution for base e, using a 50x50 matrix. (Side question: Andrew, how much faster are other libraries at solving these large matrices?)Jay,

could you provide me your gp-code for comparision? You may put it here or send me an email.

Also I would like to check, whether your last computations can be smoothed by accelerating convergence using Euler-summation of terms. But I've only Pari/GP and maxima(the latter with extremely little experience) and I don't know how complex your computations were and how much work it would be to get it implemented.

What do you think?

Gottfried

I posted my code here:

http://math.eretrandre.org/tetrationforu...213#pid213

I gave up on PARI/gp, the solve_right() function in SAGE (using maxima) was about 50-100 times faster. The matsolve function in gp must not be a low-level function, and I didn't bother looking for one once I got it working with maxima.

~ Jay Daniel Fox