Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
fast accurate Kneser sexp algorithm
(07/15/2011, 11:01 PM)nuninho1980 Wrote:
(07/15/2011, 09:34 PM)sheldonison Wrote: If the precision is manually modified via "\p 134", then "init;loop" for the base e default, it takes about 40 seconds, for 26 iterations, with results accurate to ~225 binary digits.
I tested at same setting and I got 31 seconds on my new AMD X6 1055T@2.8GHz - PC desktop (after old Intel C2D E6600) running Win 7 x64-SP1. Smile

I tried "init(I)" but I got error: "if: incorrect type in comparison."?!?! how? I already read that you can get I^^I. wow! Smile
Thanks. There's a dozen different "if" statements, most of which are in the initialization code, testing the value of B, the sexp base for init(B). If statements <> don't work for complex number. The if statements are required to allow a large range of bases to initialize correctly to guarantee subsequent convergence. Compare statements are required as the base approaches eta, or is less than eta, or becomes larger, to calculate the fixed point, and to setup the initial approximation for sexp(z). Its not practical to have this program work for complex bases, although someday I'd like to calculate a Kneser like mapping for a complex base with two repelling fixed points. I've also written some fixed point scripts, that calculate the superfunction and inverse superfunction for complex bases.
- Shel

Messages In This Thread
The pari-GP code - by sheldonison - 08/07/2010, 09:17 PM
updated code - by sheldonison - 08/19/2010, 02:35 AM
RE: updated code - by nuninho1980 - 08/19/2010, 12:08 PM
RE: updated code - by sheldonison - 08/20/2010, 01:05 AM
update to support B<eta - by sheldonison - 11/15/2010, 02:53 PM
RE: update to support B<eta - by nuninho1980 - 11/15/2010, 03:26 PM
another new version - by sheldonison - 11/17/2010, 06:52 PM
RE: fast accurate Kneser sexp algorithm - by sheldonison - 07/16/2011, 05:55 PM

Possibly Related Threads...
Thread Author Replies Views Last Post
  "Kneser"/Riemann mapping method code for *complex* bases mike3 2 7,465 08/15/2011, 03:14 PM
Last Post: Gottfried
  Attempt to make own implementation of "Kneser" algorithm: trouble mike3 9 17,784 06/16/2011, 11:48 AM
Last Post: mike3
  Numerical algorithm for Fourier continuum sum tetration theory mike3 12 22,222 09/18/2010, 04:12 AM
Last Post: mike3
  regular sexp: curve near h=-2 (h=-2 + eps*I) Gottfried 2 6,639 03/10/2010, 07:52 AM
Last Post: Gottfried
  Attempting to compute the kslog numerically (i.e., Kneser's construction) jaydfox 11 20,826 10/26/2009, 05:56 PM
Last Post: bo198214
  regular sexp:different fixpoints Gottfried 6 13,232 08/11/2009, 06:47 PM
Last Post: jaydfox
  sexp(strip) is winding around the fixed points Kouznetsov 8 14,684 06/29/2009, 10:05 AM
Last Post: bo198214
  sexp and slog at a microcalculator Kouznetsov 0 3,490 01/08/2009, 08:51 AM
Last Post: Kouznetsov

Users browsing this thread: 1 Guest(s)