Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
fast accurate Kneser sexp algorithm
#28
(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
Reply


Messages In This Thread
The pari-GP code - by sheldonison - 08/07/2010, 09:17 PM
updated kneser.gp code - by sheldonison - 08/19/2010, 02:35 AM
RE: updated kneser.gp code - by nuninho1980 - 08/19/2010, 12:08 PM
RE: updated kneser.gp 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 6,948 08/15/2011, 03:14 PM
Last Post: Gottfried
  Attempt to make own implementation of "Kneser" algorithm: trouble mike3 9 16,448 06/16/2011, 11:48 AM
Last Post: mike3
  Numerical algorithm for Fourier continuum sum tetration theory mike3 12 20,479 09/18/2010, 04:12 AM
Last Post: mike3
  regular sexp: curve near h=-2 (h=-2 + eps*I) Gottfried 2 6,189 03/10/2010, 07:52 AM
Last Post: Gottfried
  Attempting to compute the kslog numerically (i.e., Kneser's construction) jaydfox 11 19,183 10/26/2009, 05:56 PM
Last Post: bo198214
  regular sexp:different fixpoints Gottfried 6 12,254 08/11/2009, 06:47 PM
Last Post: jaydfox
  sexp(strip) is winding around the fixed points Kouznetsov 8 13,525 06/29/2009, 10:05 AM
Last Post: bo198214
  sexp and slog at a microcalculator Kouznetsov 0 3,283 01/08/2009, 08:51 AM
Last Post: Kouznetsov



Users browsing this thread: 1 Guest(s)