complex base tetration program
#21
(02/07/2016, 12:28 PM)sheldonison Wrote: In the case at hand, the other fixed point, for the lower half of the complex plane for sexp(z), is -1.862-0.411i, which Gottfried has listed strangely as 3-periodic, where as it is a simple repelling fixed point for \( i^z \). (...)
Upps, that was an error in the manual transfer of the values into the table, sorry. I've updated the picture.

Gottfried
Gottfried Helms, Kassel
#22
Hi Sheldon -

inspired by the bugs, rooted in interferences between the procedures in fatou.gp and my std-user-library I was also reviewing the older tetcomplex.gp procedure which you've posted earlier (and which worked in my environment with very little renaming-conventions), and did also some minor edits for numerical improvements (replaced "eval"-procedures by "subst") and deleted all "kill(x)"-calls.


Using the base \( \small B=2 \) I could nicely reproduce the problem in MSE: does exist and then what is \( \small x \gt 2 \) in \( \;^2x = \;^x2 \)

By some other idea I investigated the behave near above the real axis, and specifically I looked at the point \( \small z_0= \text{Re} (L) + \varepsilon I \) and found, that the required height \( h \) to iterate from \( z_0 \) to\( L \) is nearly \( \small 2 \pi I \) , so

\( \small \text{sexp}(2 \pi I+ \text{slog}(z_0 ) )\approx L \qquad \qquad \text{ using }z_0 \text{ with } \varepsilon=1e-6,1e-12,1e-24,... \)

That this is very near suggests, that here equality is intended (and it is also different from the Schröder-method). However, the error is in the near of 0.003 which looks large in comparision to the numerical accuracy for other computations (20 to 30 digits correct by my initialization) .

Could that result be improved by some finetuning? (Or did I speculate wrongly?)

Gottfried
Gottfried Helms, Kassel
#23
(10/24/2016, 11:50 PM)Gottfried Wrote: By some other idea I investigated the behave near above the real axis, and specifically I looked at the point \( \small z_0= \text{Re} (L) + \varepsilon I \) and found, that the required height \( h \) to iterate from \( z_0 \) to\( L \) is nearly \( \small 2 \pi I \) , so

\( \small \text{sexp}(2 \pi I+ \text{slog}(z_0 ) )\approx L \qquad \qquad \text{ using }z_0 \text{ with } \varepsilon=1e-6,1e-12,1e-24,... \)

That this is very near suggests, that here equality is intended....

Gottfried,

You're using a real base>exp(1/e). Either base(2) or base(e), so you don't really need the complex base support in tetcomplex. Kneser.gp would work just fine, and as you noted, tetcomplex won't work with the latest versions of pari-gp.

For any real valued sexp base as imag(z) increases, sexp(z) will go to the fixed point, since the Kneser tetration function is a 1-cyclic mapping of S(z) where S(z) is the Schröder solution superfunction from the complex fixed point. By definition, \( S(z+1)=\exp_b(S(z))\;\; \). In the upper half of the complex plane, we use \( S_U(z) \) and in the lower half of the complex plane we use \( S_L(z) \) generated from the complex conjugate fixed point.

\( \text{If}\; \Im(z)>0 \;\;\text{sexp}(z)=\text{S}_U(z+\theta_U(z)) \)

\( \text{If} \; \Im(z)<0 \;\;\text{sexp}(z)=\text{S}_L(z+\theta_L(z)) \)

The solution in the upper half of the complex plane approaches \( S_U(z+k) \) as \( \Im(z) \) increases. This is because the periodic terms in the 1-cyclic \( \theta_U(z) \) mapping quickly decay to zero as imag(z) increases leaving only the constant term of k. Try plugging in z=2*I into this equation for theta, to see how quickly the periodic terms in \( \theta_U(z) \) decay.

\( \theta_U(z) = k + \sum_{n=1}^{\infty} \; a_n \exp(2n\pi i z) \)

The \( S_U(z)\;S_L(z) \) functions are also periodic with a different period that decays more slowly. For base2 the Su(z) Schröder solution superfunction has a period of 5.584 + 1.0542*I, and Su(z) decays to the fixed point as \( \Im(z) \) increases, but the periodic terms in Su(z) decay much more slowly than the \( \theta_U(z) \) function.

At the real axis \( \theta_U(z) \) has a really complicated singularity at integer values of z, since the Schröder superfunction solution never takes on the values of 0 anywhere in the complex plane, and yet sexp(-1)=0. So \( \theta_U(z) \) is only used in the upper half of the complex plane; and its complex conjugate is used in the lower half of the complex plane.

So sexp_2(2*Pi*I) = 0.827783486030388 + 1.56642681578434*I; and the fixed point is 0.824678546142074 + 1.56743212384965*I.
abs(sexp_2(2*Pi*I)-L)=0.00326
abs(sexp_2(10*I)-L)=0.0000572
abs(sexp_2(20*I)-L)=1.101E-9
abs(sexp_2(30*I)-L)=2.105E-14
I hope this helps.
- Sheldon
#24
(10/26/2016, 08:01 AM)sheldonison Wrote: So sexp_2(2*Pi*I) = 0.827783486030388 + 1.56642681578434*I; and the fixed point is 0.824678546142074 + 1.56743212384965*I.
abs(sexp_2(2*Pi*I)-L)=0.00326
abs(sexp_2(10*I)-L)=0.0000572
abs(sexp_2(20*I)-L)=1.101E-9
abs(sexp_2(30*I)-L)=2.105E-14
I hope this helps.
Hi Sheldon -

yes, this helps, thank you. My speculation was unfounded.

Gottfried

<hr>
P.s. I have still to use Pari/GP in the 2.2.11 version because in the next version they removed the feature, switched on be option -emacs, to give a textual "prompt" after each output (*) when Pari/GP was used in as subprocess in a GUI. I could use that "prompt" to detect end-of-answer and construct the conversation with the GUI around this. I don't know why they did this and after a query at Karim Belabas, I'm left with little hope that this behaviour would have attention on it. So as long as I use my nice GUI I've accordingly little hope to be able to upgrade some time...

(*) upps, it was modified in the next version and switched off in the >2.3 versions. First I thought it is a modification by the "readline.dll" people, but it seems to be in the sources of Pari/GP itself. Only if you use Pari/GP in a DOS-box/as a Windows-"console"-program, that "prompt" is still produced. I do not really understand what is going on here...
Gottfried Helms, Kassel


Possibly Related Threads…
Thread Author Replies Views Last Post
Wink new fatou.gp program sheldonison 33 68,552 07/10/2022, 02:23 AM
Last Post: Catullus
  The beta method program JmsNxn 0 1,398 02/25/2022, 03:05 AM
Last Post: JmsNxn
  My new ABEL_L.gp program JmsNxn 13 11,317 10/06/2021, 07:18 PM
Last Post: Ember Edison
  Natural complex tetration program + video MorgothV8 1 6,796 04/27/2018, 07:54 PM
Last Post: MorgothV8
  Mathematica program for tetration based on the series with q-binomial coefficients Vladimir Reshetnikov 0 5,466 01/13/2017, 10:51 PM
Last Post: Vladimir Reshetnikov
  C++ program for generatin complex map in EPS format MorgothV8 0 5,737 09/17/2014, 04:14 PM
Last Post: MorgothV8
  Green Eggs and HAM: Tetration for ALL bases, real and complex, now possible? mike3 27 76,688 07/02/2014, 10:13 PM
Last Post: tommy1729
  "Kneser"/Riemann mapping method code for *complex* bases mike3 2 12,672 08/15/2011, 03:14 PM
Last Post: Gottfried
  Expansion of base-e pentation andydude 13 48,187 07/02/2011, 01:40 AM
Last Post: Cherrina_Pixie
  fractional iteration with complex bases/a bit of progress Gottfried 1 7,558 07/21/2008, 10:58 PM
Last Post: Gottfried



Users browsing this thread: 1 Guest(s)