02/06/2016, 04:04 PM
(This post was last modified: 02/06/2016, 04:35 PM by sheldonison.)
(02/06/2016, 01:37 AM)Gottfried Wrote: Hmm, today I tried tetcomplex.gp with init(I). Unfortunately the program hangs/loops infinitely after the message
generating Schroder2 taylor series for isuperf2 function, scnt2 27
Using a real base, say init(1.44) works immediately.
What to do?
update: I found that it hangs in the routine "loop" when called from the "init(I)"-procedure. (Surprisingly "init()" doesn't provide an argument to the routine loop while that has a formal parameter "t")
update2: it enters "thetaup" and does not come back...
Gottfried
hmmm, I haven't used tetcomplex in awhile. It is much much more limited in terms of what bases it will converge for than the newer fatou.gp program. For example, tetcomplex has no hope of converging anywhere near base(i), whereas fatou.gp works just fine. Also, the old program has no "fallback" algorithm to use for rationally indifferent fixed points, that are on the ShellThron boundary, whereas as the newer algorithm can work (although with less precision), without any Schroeder function whatsoever. The newer algorithm is much better, except for bases<eta, where someday I may allow for rotation angles >180 degrees for fatou.gp, but not yet. By the way, I finally had some time to post my answer to sexp base(i) on mathstack; http://math.stackexchange.com/questions/...35#1643235. I wojuld like to post more about fatou.gp here on the tetration forum as well.
Code:
\r fatou.gp
setmaxconvergence(); /* base i is hard to compute */
sexpinit(i);
sexp(0.5)
1.07571355731392 + 0.873217399108003*I
- Sheldon

