Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
complex base tetration program
#1
This is the link to tetcomplex.gp pari gp code, which generates tetration for arbitrary bases in the complex plane, using the merged fixed point solution. You'll need to download par-gp.. Then load tetcomplex.gp and init(b), for a complex base, and then sexp(z) gives the resulting tetration function at tet(z). tet(0)=1.
.gp   tetcomplex.gp (Size: 59.81 KB / Downloads: 493) updated March 11th, improved convergence near eta, and works with arbitrary precision settings. Previous updates improved convergence near Shell Thron boundary, and for real bases<eta, and as imag(z) gets bigger The program is now working reasonably well, but is not fully debugged. The program generalizes the ideas Mike and I talked about, and Mike implemented, for generating merged fixed point tetration solutions in the complex plane. Mike has been advocating these merged fixed point solutions as the natural extension of tetration at the real axis, for bases>eta, as well as Dimitrii Kouznetsov. This program goes beyond what Mike did, in that it works for arbitrary bases, inside and outside the Shell Thron region. For bases not too close to one and not too close to the Shell Thron boundary, this program will generate both primary fixed points, using this algorithm, and figures out if they're both repelling, or if one of them is attracting. Then it will generate both Schroder functions and their inverses for both fixed points, and use those to generate superfunctions and superfunction inverses for both fixed points. Then the program iterates, generating two theta(z) mappings, one for the upper fixed point superfunction, and one for the low fixed point superfunction. Between these, is an sexp(z) taylor series. For imaginary(z)>0.175, , and for imaginary(z)<-0.175, . For bases at the real axis > eta, this is mathematically identical to my Kneser tetration algorithm.

The program is fast. The program does have bugs, although it is already working surprisingly well. Also, if the base is too close to the Shell Thron boundary, the program gets really slow, and eventually stops converging. The program comes with a fast complex plane graphing algorithm, MakeGraph. So the default precision is \p 28, which allows very fast accurate results to 12 decimal digits, for making graphs. You can increase the prcecision to \p 67, and results should work for bases between 1.1 and 3. For some bases, higher precision works, up to \p 134, probably depending on the accuracy of the theta, and the superfunction and the initialization for that base. This routine lacks many of the complicated optimizations, that allow kneser.gp to work for a very wide range of real bases. The program will aid in the study of tetration in the complex plane. Below, are a few graphing examples. The program works best with imag(base)>0, and gives novel interesting solutions for some bases with imag(b)<0, where the base is inside the shell thron boundary.
- Sheldon

base=2+0.56i, which is just inside the Shell Thron boundary. This pair of solutions (this one and the one below), seems to indicate that a merged tetration solution on the Shell Thron boundary itself works!
   

base=2+0.51i, which is just outside the Shell Thron boundary, with both fixed points repelling. This is a natural extension of tetration for real bases>eta.
   

base=1.45+0.3i, which is tetration where the two superfunctions are roughly at 45 degree angles, and -135 degree angles, from the sexp(z) at z=0.
   

base=1.45-0.3i. This result totally blew me away, because the upper fixed point, for imag(z)>0, for is 1.2125955-0.6939i, which has imag(z)<0, and the lower fixed point is 1.6459+3.82111i, which has imag(z)>0. Also, the superfunction gets more chaotic as imag(z) decreases for the lower fixed point, which is also interesting. The alternative solution, rotating in the other direction around eta, is the complex conjugate solution of the solution above, which makes much more sense, with the fixed points in the "correct" locations. I wonder if this solution can be extended, as we continue rotating further around eta?
- Sheldon
   
Reply


Messages In This Thread
complex base tetration program - by sheldonison - 02/29/2012, 10:28 PM
RE: complex base tetration program - by mike3 - 03/10/2012, 06:59 AM
RE: complex base tetration program - by mike3 - 03/11/2012, 12:27 AM
RE: complex base tetration program - by Gottfried - 02/06/2016, 01:37 AM
RE: complex base tetration program - by Gottfried - 02/06/2016, 05:36 PM
RE: complex base tetration program - by Gottfried - 02/07/2016, 05:27 AM
RE: complex base tetration program - by Gottfried - 02/07/2016, 01:25 PM
RE: complex base tetration program - by Gottfried - 10/24/2016, 11:50 PM
RE: complex base tetration program - by Gottfried - 10/26/2016, 10:02 AM

Possibly Related Threads...
Thread Author Replies Views Last Post
Wink new fatou.gp program sheldonison 20 18,773 02/11/2020, 04:47 PM
Last Post: sheldonison
  Natural complex tetration program + video MorgothV8 1 1,714 04/27/2018, 07:54 PM
Last Post: MorgothV8
  Mathematica program for tetration based on the series with q-binomial coefficients Vladimir Reshetnikov 0 1,980 01/13/2017, 10:51 PM
Last Post: Vladimir Reshetnikov
  C++ program for generatin complex map in EPS format MorgothV8 0 2,611 09/17/2014, 04:14 PM
Last Post: MorgothV8
  Green Eggs and HAM: Tetration for ALL bases, real and complex, now possible? mike3 27 33,946 07/02/2014, 10:13 PM
Last Post: tommy1729
  "Kneser"/Riemann mapping method code for *complex* bases mike3 2 6,635 08/15/2011, 03:14 PM
Last Post: Gottfried
  fractional iteration with complex bases/a bit of progress Gottfried 1 3,787 07/21/2008, 10:58 PM
Last Post: Gottfried



Users browsing this thread: 1 Guest(s)