+- Tetration Forum (https://tetrationforum.org)
+-- Forum: Tetration and Related Topics (https://tetrationforum.org/forumdisplay.php?fid=1)
+--- Forum: Computation (https://tetrationforum.org/forumdisplay.php?fid=8)
+--- Thread: Designing a Tetration Library (/showthread.php?tid=146)
RE: Designing a Tetration Library - andydude - 05/19/2008
bo198214 Wrote:I described it here. The principle \( f^{\circ t} \circ f=f\circ f^{\circ t} \) with \( {f^{\circ t}}_1={f_1}^t \) is common in the regular iteration literature, however the formula is perhaps too dirtythat I have seen it explicitely mentioned.
Ah, OK. For some reason, when I apply my interpretation of the formula above, I get \( t-1 \) instead of t, perhaps, I'm missing something. I'll take a look at the thread you mentioned.
bo198214 Wrote:Originally I wanted to use "^" via __and__ for the iteration, like I used "**" via __pow__ for the powers. Do you have an idea why this is not recognized? And also the question about how to put the hyperops package into sage that it is automatically recognized is still open for me.
Maybe this is because "^" is __xor__, not __and__.
Andrew Robbins
RE: Designing a Tetration Library - andydude - 05/21/2008
Also, your usage of __call__ goes against Sage usage. For example, PowerSeries implements coefficients through __getitem__ or f[n] while you implement coefficients through __call__ or f(n), which should probably use __getitem__ instead for this. Also, in PowerSeries, __call__ is associated with what you use .compose for, so __call__ should use this in PowerSeriesI as well. If it is OK, I will try and change these in the git repository...
Andrew Robbins
RE: Designing a Tetration Library - bo198214 - 05/21/2008
andydude Wrote:Also, your usage of __call__ goes against Sage usage. For example, PowerSeries implements coefficients through __getitem__ or f[n] while you implement coefficients through __call__ or f(n), which should probably use __getitem__ instead for this. Also, in PowerSeries, __call__ is associated with what you use .compose for, so __call__ should use this in PowerSeriesI as well. If it is OK, I will try and change these in the git repository...
No, I got the same idea and just changed it. See attachment here
But compose is not expressed via __call__. You dont write f(g), either you write f(g(x)) or you write fog.
Accordingly I also changed the compose to o. So that you can now write f.o(g)
RE: Designing a Tetration Library - bo198214 - 05/21/2008
andydude Wrote:If it is OK, I will try and change these in the git repository...
So that means you have your git working?!
Ok, then I will push my actual to the repository.
To push your changes to the repository you need to give me your public ssh key, otherwise you can only read the public repository.
RE: Designing a Tetration Library - Gottfried - 05/28/2008
For the problem of computing the matrix-exponential the following article may be of interest (and may be added to our - to be constructed - virtual library?). I've not read it yet, though.
UWE LUTHER† AND KARLA ROST‡
Abstract. For a given matrix A we compute the matrix exponential e^A
under the assumption that the eigenvalues of A are known, but without
determining the eigenvectors. The presented approach exploits the
connection between matrix exponentials and confluent Vandermonde
matrices V . This approach and the resulting methods are very simple
and can be regarded as an alternative to the Jordan canonical form
methods. The discussed inversion algorithms for V as well as the
matrix representation of V^-1 are of independent interest also in many
other applications.
Key words. matrix exponential, Vandermonde matrix, fast algorithm, inverse.
AMS subject classifications. 34A30, 65F05, 15A09, 15A23.
Electronic Transactions on Numerical Analysis.
Volume 18, pp. 91-100, 2004.
Copyright Ó 2004, Kent State University.
ISSN 1068-9613.
ETNA
Kent State University
etna@mcs.kent.edu
RE: Designing a Tetration Library - bo198214 - 05/29/2008
Gottfried Wrote:For the problem of computing the matrix-exponential
Hm, thanks for the link however do we need the matrix-exponential? We just need \( A^t \), which perhaps can be also expressed as \( \exp(t\log(A)) \) but in this case we also need a fast algorithm for matrix log, or are there other direct uses of the matrix exponential?
RE: Designing a Tetration Library - Gottfried - 05/29/2008
bo198214 Wrote:Gottfried Wrote:For the problem of computing the matrix-exponential
Hm, thanks for the link however do we need the matrix-exponential? We just need \( A^t \), which perhaps can be also expressed as \( \exp(t\log(A)) \) but in this case we also need a fast algorithm for matrix log, or are there other direct uses of the matrix exponential?
Well, Andrew asked for such a routine, if I recall right
RE: Designing a Tetration Library - bo198214 - 06/06/2008
Gottfried Wrote:bo198214 Wrote:Hm, thanks for the link however do we need the matrix-exponential?Well, Andrew asked for such a routine, if I recall right
I now see that the matrix logarithm could be the iterative logarithm, up to multiplicative constant perhaps.
If we define \( (\text{ilog}(f))_n=\log(CM(f))_{1,n} \) then exactly \( \text{ilog}(f^{\circ t})=\log(CM(f)^t)_1=t\text{ilog}(f) \) is satisfied.
RE: Designing a Tetration Library - andydude - 06/06/2008
bo198214 Wrote:If we define \( (\text{ilog}(f))_n=\log(CM(f))_{1,n} \) then exactly \( \text{ilog}(f^{\circ t})=\log(CM(f)^t)=t\text{ilog}(f) \) is satisfied.
Do you mean \( \text{ilog}(f^{\circ t})=(\log(CM(f)^t))_1=t\text{ilog}(f) \)? Then yes, that is very nice!
Andrew Robbins
RE: Designing a Tetration Library - bo198214 - 06/07/2008
andydude Wrote:Do you mean \( \text{ilog}(f^{\circ t})=(\log(CM(f)^t))_1=t\text{ilog}(f) \)? Then yes, that is very nice!
Yes, I add the subscripted 1 in my post.