Difference between revisions of "AuTra"
m (Text replacement  "\$([^\$]+)\$" to "\\(\1\\)") 

(One intermediate revision by one other user not shown)  
Line 1:  Line 1:  
−  [[File:AuTraMapT.jpg400pxthumb 
+  [[File:AuTraMapT.jpg400pxthumb\(u\!+\!\mathrm i v= \mathrm{AuTra}(x\!+\!\mathrm i y)\)]] 
−  [[File:AuTraPlotT.jpg200pxthumb 
+  [[File:AuTraPlotT.jpg200pxthumb\(y= \mathrm{AuTra}(x)\)]] 
−  [[AuTra]] 
+  [[AuTra]] is the [[inverse function]] of [[SuTra]]. 
−  [[AuTra]] is [[Abel function]] of the [[Trappmann function]], 
+  [[AuTra]] is [[Abel function]] of the [[Trappmann function]], \(\mathrm{tra}(z)=z+\exp(z)\). 
[[AuTra]] satisfies the Abel equation 
[[AuTra]] satisfies the Abel equation 

−  +  \( \mathrm{AuTra} \Big( \mathrm{tra}(z) \Big)= \mathrm{AuTra}(z)+1\) 

and the additional condition 
and the additional condition 

−  +  \( \mathrm{AuTra}(1)=0~\) . 

−  This condition corresponds to relation 
+  This condition corresponds to relation \(~\mathrm{SuTra}(0)=1~\). The same choice of the initial value is used also for many other [[superfunction]]s, including [[tetration]], [[SuZex]] function, [[Tania function]] and [[Shoka function]]. 
−  [[AuTra]] is realholomorphic, 
+  [[AuTra]] is realholomorphic, \(\mathrm{AuTra}(z^*)=\mathrm{AuTra}(z)^*\) 
−  ==Branch points 
+  ==Branch points, cut lines and growth== 
−  [[AuTra]] has two branch points 
+  [[AuTra]] has two branch points \(~B_{\pm}=1\! \pm\!\mathrm i \pi ~\); 
−  <! The Abel function of for the Trappmann function with another choice of the cut lines. let it be called AuTraM, can be expressed with 

−  $\mathrm{AuTra}(z)=\mathrm{AuTraM}(z)~$ for $~\Im(z) < \pi~$ 

+  \(\mathrm{Au}_{\pm} := \mathrm{AuTra}(B_{\pm}) \approx 3.4101257504807645 \pm 1.4101841452081931 \, \mathrm i\) 

−  Along the real axis and its vicinity, AuTra shows fast growth 

−  $\mathrm{AuTraM}(z)=\mathrm{AuZex}( \mathrm e^z )$ 

+  These branch points corresponds to the zero derivative of function [[SuTra]]; 

+  \(\mathrm{SuTra}(\mathrm{Au}_{\pm}) = 1\! \pm\!\mathrm i \pi ~\) 

−  Between the cut llines, AuTra shows almost linear decrease, that is represented with almost uniform rectangular grid at the complex map. 

+  
−  Outside the cut lines, AuTra shows slow growth, similar to that of the [[ArcTetration]]. 

+  \(\mathrm{SuTra}^\prime (\mathrm{Au}_{\pm}) =0 ~\) 

−  !> 

+  
+  The cut lines 

+  
+  \(\{ z\in \mathbb C : ~ \Im(z)\!=\!\pm i, ~ \Re(z)\le 1 \}\) 

+  
+  are directed to the left hand side of the complex plane, parallel to the real axis. In [[TORI]], this is default choice of the cut lines. 

+  
+  Between the cut llines, [[AuTra]] grows exponentially in the direction of decrease of the real part of the argument. 

+  
+  Outside the cut lines, AuTra]] grows slowly, similarly to the [[ArcTetration]]. 

+  In particular, [[AuTra]] shows this slow growth along the positive direction of the real axis. 

+  This growth is a little bit faster than that of function [[AuZex]]: 

+  
+  \(\mathrm{AuTraM}(z)=\mathrm{AuZex}( \mathrm e^z )\) 

+  
+  Once [[AuZex]] is implemented, this expression can be used for the numerical implementation of [[AuTra]] at moderate values of the real part of the argument, while its exponential still fits the range of values allowed for the variable used in the computation. 

+  
+  ==Asymptotic expansion== 

+  The asymptotic expansion for AuTra for large negative values of the argument can be obtained, inverting series for function [[SuTra]]. 

+  The inverse of \(F\), id est, \(G=F^{1}\) has the following asymptotic: 

+  
+  \(\displaystyle G(z)=\frac{z}{2}\mathrm e^{z}\) \(\displaystyle 

+  \frac{\mathrm e^z}{6}\) \(\displaystyle 

+  +\frac{\mathrm e^{2 z}}{16}\) \(\displaystyle 

+  \frac{19 \mathrm e^{3 z}}{540}\) \(\displaystyle 

+  +\frac{\mathrm e^{4 z}}{48}\) \(\displaystyle 

+  \frac{41 \mathrm e^{5 z}}{4200}\) \(\displaystyle 

+  + \mathcal O( \mathrm e^{6 z}) 

+  \) 

+  
+  The coefficients of the expansion above are calculated and evaluated with the [[Mathematica]] code below: 

+  <poem> 

+  g0[z_] = z/2  Exp[z] + Sum[c[n] Exp[n z], {n, 1, 20}] 

+  n = 1; s[n] = Series[g0[Log[t]] + 1  g0[tra[Log[t]]], {t, 0, n + 1}] 

+  u[n] = Extract[Solve[Coefficient[s[n], t^(n + 1)] == 0, c[n]], 1] 

+  g[n, z_] = ReplaceAll[g0[z], t[n]] 

+  For[n = 1, n < 20, n++; 

+  s[n] = Series[ g[n  1, Log[t]] + 1  g[n  1, tra[Log[t]]], {t, 0, n + 1}]; 

+  u[n] = Extract[Solve[Coefficient[s[n], t^(n + 1)] == 0, c[n]], 1]; 

+  g[n, z_] = ReplaceAll[g[n  1, z], u[n]]; 

+  ] 

+  g[n, z] 

+  Table[Coefficient[g[n, z], Exp[n z]], {n, 1, 20}] 

+  N[Table[Coefficient[g[n, z], Exp[n z]], {n, 1, 20}], 18] 

+  </poem> 

+  For some fixed integer \(M\), define the primary approximation as truncation of the series above: 

+  
+  \(\displaystyle 

+  G_M(z)=\frac{z}{2}\mathrm e^{z}\) \(\displaystyle +\sum_{m=1}^{M} c_m \mathrm e^{mz} \) 

+  
+  Function \(G\) can be defined as limit 

+  
+  \(\displaystyle 

+  G(z)=\lim_{n\rightarrow \infty} (G_M(\mathrm{ArcTra}^n(z))+n)\) 

+  
+  AuTra can be expressed thirough \(G\) with 

+  
+  \( \mathrm{AuTra}(z)=G(z)G(0) \approx G(z)+1.1259817765745026\) 

+  
+  The constant \(G(0)\) can be interpreted also as coefficient \(c_0\), id est, \(c_0=G(0)\); then, 

+  finction \(g\) can be defined with 

+  
+  \(\displaystyle 

+  g_M(z)=\frac{z}{2}\mathrm e^{z}\) \(\displaystyle +\sum_{m=0}^{M} c_m \mathrm e^{mz} \) 

+  
+  and 

+  
+  \( \mathrm{AuTra}(z)=\lim_{n\rightarrow \infty} g_M(\mathrm{ArcTra}^n(z))\) 

+  
+  This representation is used for the numericcal implementation of AuTra described below. 

+  
+  ==Numerical implementation of AuTra== 

+  [[File:Autran0m9tes64t.jpg160pxthumbAgreement \(\mathcal A\!=\!A_9(x\!+\!\mathrm i y)~\) ]] 

+  
+  This section describe the [[complex double]] numerical implementation of AuTra. For \(M=9\), the figure at right shows the agreement 

+  
+  \(\displaystyle 

+  A_M(z)=\lg\left( \frac{ \mathrm{SuTra}(g_M(z))z}{ \mathrm{SuTra}(g_M(z))+z} \right)\) 

+  
+  of the primary approximation \(g_M\). 

+  
+  The thick outher loop indicates \(\mathcal A=1\), function \(g_M\) does not approximate AuTra. This cannot be improved increasing \(M\). 

+  The inner loop corresponds to \(\mathcal A=16\), and it looks a little bit irregular, that can be attributed to the rounding errors at the use of the [[complex double]] arithmetics. 

+  For the same reason, the level \(\mathcal A=15\) also looks scratched. 

+  
+  The shaded region refers to \(y<3\), \(y<2.51.5y\); in this region, the primary approximation provides of order of 16 decimal digits. This precision seems 

+  to be best that can be done for complex double implementation. Values of \(g\) from the shaded region are qualified as "precise" and used for evaluation of AuTra also at other values of the argument, applying the transfer equation: 

+  Function [[ArcTra]] is applied to the argument so many times as necessary in order to bring it to the shaded region. 

+  Note that the shaded region happens to be almost twice wider, that the initial strip \(y<\pi/2\approx 1.57\), where the asymptotic expansion is valid. 

+  
+  [[File:Autran9tes20t1.jpg300pxthumbAgreement for the final implementation of [[AuTra]]]] 

+  The agreement for the resulting secondary (and, perhaps, final) approximation is shown in figure at right. The figure indicates, that the implementation returns of order of 14 decimal digits. At least partially, the deviation can be attributed to the implementation of function [[SuTra]], that is also implemented with 14 decimal digits. 

+  These 14 digits of precision are close to the maximum that can be realised with complex double arithmetics. This precision greatly exceeds the precision, required to plot the cameraready pictures of the complex maps of the related functions. 

==References== 
==References== 

Line 36:  Line 127:  
==Keywords== 
==Keywords== 

[[Abel function]], 
[[Abel function]], 

+  [[AuZex]] 

−  [[Trappmann function]] 

+  [[Numeric implementation]], 

+  [[Superfunction]]. 

+  [[SuZex]] 

+  [[Trappmann function]], 

+  
[[Category:AuTra]] 
[[Category:AuTra]] 

[[Category:Abel function]] 
[[Category:Abel function]] 

+  [[Category:SuTra]] 

[[Category:Trappmann function]] 
[[Category:Trappmann function]] 
Latest revision as of 18:48, 30 July 2019
AuTra is the inverse function of SuTra.
AuTra is Abel function of the Trappmann function, \(\mathrm{tra}(z)=z+\exp(z)\).
AuTra satisfies the Abel equation
\( \mathrm{AuTra} \Big( \mathrm{tra}(z) \Big)= \mathrm{AuTra}(z)+1\)
and the additional condition
\( \mathrm{AuTra}(1)=0~\) .
This condition corresponds to relation \(~\mathrm{SuTra}(0)=1~\). The same choice of the initial value is used also for many other superfunctions, including tetration, SuZex function, Tania function and Shoka function.
AuTra is realholomorphic, \(\mathrm{AuTra}(z^*)=\mathrm{AuTra}(z)^*\)
Contents
Branch points, cut lines and growth
AuTra has two branch points \(~B_{\pm}=1\! \pm\!\mathrm i \pi ~\);
\(\mathrm{Au}_{\pm} := \mathrm{AuTra}(B_{\pm}) \approx 3.4101257504807645 \pm 1.4101841452081931 \, \mathrm i\)
These branch points corresponds to the zero derivative of function SuTra;
\(\mathrm{SuTra}(\mathrm{Au}_{\pm}) = 1\! \pm\!\mathrm i \pi ~\)
\(\mathrm{SuTra}^\prime (\mathrm{Au}_{\pm}) =0 ~\)
The cut lines
\(\{ z\in \mathbb C : ~ \Im(z)\!=\!\pm i, ~ \Re(z)\le 1 \}\)
are directed to the left hand side of the complex plane, parallel to the real axis. In TORI, this is default choice of the cut lines.
Between the cut llines, AuTra grows exponentially in the direction of decrease of the real part of the argument.
Outside the cut lines, AuTra]] grows slowly, similarly to the ArcTetration. In particular, AuTra shows this slow growth along the positive direction of the real axis. This growth is a little bit faster than that of function AuZex:
\(\mathrm{AuTraM}(z)=\mathrm{AuZex}( \mathrm e^z )\)
Once AuZex is implemented, this expression can be used for the numerical implementation of AuTra at moderate values of the real part of the argument, while its exponential still fits the range of values allowed for the variable used in the computation.
Asymptotic expansion
The asymptotic expansion for AuTra for large negative values of the argument can be obtained, inverting series for function SuTra. The inverse of \(F\), id est, \(G=F^{1}\) has the following asymptotic:
\(\displaystyle G(z)=\frac{z}{2}\mathrm e^{z}\) \(\displaystyle \frac{\mathrm e^z}{6}\) \(\displaystyle +\frac{\mathrm e^{2 z}}{16}\) \(\displaystyle \frac{19 \mathrm e^{3 z}}{540}\) \(\displaystyle +\frac{\mathrm e^{4 z}}{48}\) \(\displaystyle \frac{41 \mathrm e^{5 z}}{4200}\) \(\displaystyle + \mathcal O( \mathrm e^{6 z}) \)
The coefficients of the expansion above are calculated and evaluated with the Mathematica code below:
g0[z_] = z/2  Exp[z] + Sum[c[n] Exp[n z], {n, 1, 20}]
n = 1; s[n] = Series[g0[Log[t]] + 1  g0[tra[Log[t]]], {t, 0, n + 1}]
u[n] = Extract[Solve[Coefficient[s[n], t^(n + 1)] == 0, c[n]], 1]
g[n, z_] = ReplaceAll[g0[z], t[n]]
For[n = 1, n < 20, n++;
s[n] = Series[ g[n  1, Log[t]] + 1  g[n  1, tra[Log[t]]], {t, 0, n + 1}];
u[n] = Extract[Solve[Coefficient[s[n], t^(n + 1)] == 0, c[n]], 1];
g[n, z_] = ReplaceAll[g[n  1, z], u[n]];
]
g[n, z]
Table[Coefficient[g[n, z], Exp[n z]], {n, 1, 20}]
N[Table[Coefficient[g[n, z], Exp[n z]], {n, 1, 20}], 18]
For some fixed integer \(M\), define the primary approximation as truncation of the series above:
\(\displaystyle G_M(z)=\frac{z}{2}\mathrm e^{z}\) \(\displaystyle +\sum_{m=1}^{M} c_m \mathrm e^{mz} \)
Function \(G\) can be defined as limit
\(\displaystyle G(z)=\lim_{n\rightarrow \infty} (G_M(\mathrm{ArcTra}^n(z))+n)\)
AuTra can be expressed thirough \(G\) with
\( \mathrm{AuTra}(z)=G(z)G(0) \approx G(z)+1.1259817765745026\)
The constant \(G(0)\) can be interpreted also as coefficient \(c_0\), id est, \(c_0=G(0)\); then, finction \(g\) can be defined with
\(\displaystyle g_M(z)=\frac{z}{2}\mathrm e^{z}\) \(\displaystyle +\sum_{m=0}^{M} c_m \mathrm e^{mz} \)
and
\( \mathrm{AuTra}(z)=\lim_{n\rightarrow \infty} g_M(\mathrm{ArcTra}^n(z))\)
This representation is used for the numericcal implementation of AuTra described below.
Numerical implementation of AuTra
This section describe the complex double numerical implementation of AuTra. For \(M=9\), the figure at right shows the agreement
\(\displaystyle A_M(z)=\lg\left( \frac{ \mathrm{SuTra}(g_M(z))z}{ \mathrm{SuTra}(g_M(z))+z} \right)\)
of the primary approximation \(g_M\).
The thick outher loop indicates \(\mathcal A=1\), function \(g_M\) does not approximate AuTra. This cannot be improved increasing \(M\). The inner loop corresponds to \(\mathcal A=16\), and it looks a little bit irregular, that can be attributed to the rounding errors at the use of the complex double arithmetics. For the same reason, the level \(\mathcal A=15\) also looks scratched.
The shaded region refers to \(y<3\), \(y<2.51.5y\); in this region, the primary approximation provides of order of 16 decimal digits. This precision seems to be best that can be done for complex double implementation. Values of \(g\) from the shaded region are qualified as "precise" and used for evaluation of AuTra also at other values of the argument, applying the transfer equation: Function ArcTra is applied to the argument so many times as necessary in order to bring it to the shaded region. Note that the shaded region happens to be almost twice wider, that the initial strip \(y<\pi/2\approx 1.57\), where the asymptotic expansion is valid.
The agreement for the resulting secondary (and, perhaps, final) approximation is shown in figure at right. The figure indicates, that the implementation returns of order of 14 decimal digits. At least partially, the deviation can be attributed to the implementation of function SuTra, that is also implemented with 14 decimal digits. These 14 digits of precision are close to the maximum that can be realised with complex double arithmetics. This precision greatly exceeds the precision, required to plot the cameraready pictures of the complex maps of the related functions.
References
Keywords
Abel function, AuZex Numeric implementation, Superfunction. SuZex Trappmann function,