Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
Dmitrii Kouznetsov's Tetration Extension
#31
Kouznetsov Wrote:Thanks, Henryk. I think, at least for , our tetrations coincide.

I too think so, and above they can not coincide because regular iteration is not defined without fixed point, and development at non-real fixed points does not give real values on the real axis.

Though I have to mention that sometimes differences occur at the level of or even especially if you compare regular iterations at different fixed points.

Quote:I still wonder why do you direct the axis down (at the top it is and at the bottom ),

Because I am just too lazy to properly orientate the drawing Smile. Instead was more focussed to let it look nice, haha.

Quote: but anyway, I believe, you have many correct decimal digits.

I can arbitrarily increase the precision, it just takes longer then.
Reply
#32
bo198214 Wrote:
Kouznetsov Wrote:Henryk, I hope you agree with Figure 2.
As long as you didnt have explained it I can not agree Wink
I have updated the tertation for small b . The minimal of quasi-periods determines the periodicity; the side of longest quasiperiod (negarive real part) has cuts of analyticity. Please, indicate the first statement you consider as "not explained".

Quote:So then give me the next step of your Dmitrii-makes-contourplot-of-tetration-tutorial.
Please, enter into my code, extract few values and compare with yours. I expect you will see disagreement or order of one percent, even at point . This is because I did not adjust well the misplacement along the real axis. Assume, my code evaluates some function . Solve numerically equation and compare function to your tetration.

Quote:your code is in the complete spirit of assembler, as if the world hadnt changed since that time
You have no need to work much inside the "conto" function; the only, may be, suppress the message about "Copyleft"; it will run even faster. (just type // at the beginning of line). Map your tetration with my conto plotter. Ask questions if any doubt.
Reply
#33
Kouznetsov Wrote:
bo198214 Wrote:
Kouznetsov Wrote:Henryk, I hope you agree with Figure 2.
As long as you didnt have explained it I can not agree Wink
I have updated the tertation for small b . The minimal of quasi-periods determines the periodicity; the side of longest quasiperiod (negarive real part) has cuts of analyticity. Please, indicate the first statement you consider as "not explained".

Good with the accompanying text it becomes clear.
One minor correction: In the left part of fig. 2 its not but as the period is purely imaginary. Did you btw play around with other complex fixed points, differently from the one nearest the real axis? For the distribution of fixed points I even created a small "movie" here.


Quote:Please, enter into my code, extract few values and compare with yours.
But you didnt give me your code yet! I have the code for fig. 2 and the example of the contourplot.

Quote:You have no need to work much inside the "conto" function; the only, may be, suppress the message about "Copyleft"; it will run even faster. (just type // at the beginning of line). Map your tetration with my conto plotter.

Attached is the somewhat ameteurish result for base . The red lines are .


Attached Files
.eps   expbs2_20x20_2.eps (Size: 24.44 KB / Downloads: 302)
Reply
#34
bo198214 Wrote:...
In the left part of fig. 2 its not but as the period is purely imaginary.
Yes. Thanks!
Quote:Did you btw play around with other complex fixed points, differently from the one nearest the real axis?
[quote]
Not yet. We shold calculate tetrations with various asymptotic behavior.
[quote]
For the distribution of fixed points I even created a small "movie"
here.
Good! Could you trace them? A movie would be good for evolution of the map of tetration as ln(b) increases. (Yet, I have only 3 stapshots: b=sqrt(2),b=2,b=e).
Quote:But you didnt give me your code yet!
Oops... I was sure I did. Now I see, I forgot. Sorry. Check now, codes for fig.3a and fig.3b should be there.
Quote:Attached is the somewhat ameteurish result for base . The red lines are .
Why do blue lines pass through integer values at the real axis?
(You may reduce a little bit up-to-last argument in the conto function, and increase a little bit the last argument; contour lines will come closer to singulatity).
Reply
#35
Henryk, I have posted the Gauss-legendre quadrature formulas at http://en.citizendium.org/wiki/GauLegExample/code. Could you check it? I hope, it is the last brick you need for evaluation of the contour integral.
By the way, what happened with your application to the http://en.citizendium.org?
I wanted you to post there maps of tetration for different fixed points...
Reply
#36
bo198214 Wrote:.. If you speak French I advice you to read
Jean Écalle, Théorie des invariants holomorphes, Publ. Math. Orsay No. 67–7409
I agree to read French. Could you send please the pdf?
Reply
#37
Just to throw something in, I computed the image of the Cartesian grid for the matrix power super exponential. This can directly be compared with
Figure 7 in D.Kouznetsov. Analytic solution of F(z+1)=exp(F(z)) in complex z-plane. Mathermatics of Computation, in press (200Cool, and it looks quite similar:
   
The picture shows the image of lines with constant real part. The thick red lines are -1,0,1,2.
Reply
#38
bo198214 Wrote:Just to throw something in, I computed the image of the Cartesian grid for the matrix power super exponential. This can directly be compared with
Figure 7 in D.Kouznetsov. Analytic solution of F(z+1)=exp(F(z)) in complex z-plane. Mathermatics of Computation, in press (200Cool...
   
1. I attach the overlap of the figures mentioned. The deviation (if any) is smaller than width of lines; the agreement at least with 2 significant decimal digits.
2. Is the superexponential implemented in the Matrix Power software? Can we cite it?

(for some reason, the image does not load by itself, I have to clisk the url...)

Moderator's Update: Image now visible
Reply
#39
bo198214 Wrote:Thats always the difficulty with those numerical comparison. If the difference is small does that mean that they are equal?
Henryk asked me to calcilate my tetration with displaced basic interval.
In the paper http://www.ils.uec.ac.jp/~dima/PAPERS/2008analuxp99.pdf ,
the basic algorithm evaluates tetration() in the range
; precisely (14 digits) in the range ;
I use it for , and bring to this interval other values using the recurrent relation.
Then, I shifted the interval; the new implementation evaluates the tetration for
.
I run the test; it is precise for .
I use it for , bringing other cases to this intergal.
The plots for old algorithm and for the new one look identical.
Then I plot the difference .
   
The grid covers the range


with unity step.
Levels are shown with thick pink lines
Levels are shown with thin pink lines
Levels are shown with thin green lines
Levels are shown with thck green lines

Levels are shown with thick red lines
Levels are shown with thin red lines
Levels are shown with thin blue lines
Levels are shown with thck blue lines

If variation of function esceed 10 per step of the mesh, then the line is not drawn; so, the right hand side of the figure in vicinity of the real axix left blank; but there, the deviation between these two functions is larger than 1.e-14 and, perhaps, even larger than unity; the only the plotter cannot identify the position of lines, tey are too dence.

I see the jumps of the difference at the integer and at half-integer values of the real part; the jumps are at the level of . Each of the implementations has its own errors, and they can be revealed comparing two finctions. I attribute the humps at the half-integer values to the old algorithm, and jumps at the integer values to the new algorithm. Comparing these two algorithms, I cannot say that one of them is somehow better than another one.

Conclusion: In order to see from numerical evaluation, that two functions are not the same, you should plot them in the complex plane. It is better than to compare the Tailor series.
Reply


Possibly Related Threads...
Thread Author Replies Views Last Post
  Possible continuous extension of tetration to the reals Dasedes 0 1,248 10/10/2016, 04:57 AM
Last Post: Dasedes
  Andrew Robbins' Tetration Extension bo198214 32 41,504 08/22/2016, 04:19 PM
Last Post: Gottfried
  Kouznetsov-Tommy-Cauchy method tommy1729 0 1,946 02/18/2015, 07:05 PM
Last Post: tommy1729
  Non-trivial extension of max(n,1)-1 to the reals and its iteration. MphLee 3 3,706 05/17/2014, 07:10 PM
Last Post: MphLee
  extension of the Ackermann function to operators less than addition JmsNxn 2 3,900 11/06/2011, 08:06 PM
Last Post: JmsNxn
  about kouznetsov again tommy1729 1 3,076 08/13/2010, 05:41 AM
Last Post: BenStandeven
  plz explain kouznetsov slowly tommy1729 7 7,398 07/27/2010, 04:57 AM
Last Post: bo198214
  Tetration Extension to Real Heights chobe 3 5,966 05/15/2010, 01:39 AM
Last Post: bo198214
  Tetration extension for bases between 1 and eta dantheman163 16 18,037 12/19/2009, 10:55 AM
Last Post: bo198214
  Extension of tetration to other branches mike3 15 19,934 10/28/2009, 07:42 AM
Last Post: bo198214



Users browsing this thread: 1 Guest(s)