• 1 Vote(s) - 5 Average
• 1
• 2
• 3
• 4
• 5
 fast accurate Kneser sexp algorithm sheldonison Long Time Fellow Posts: 683 Threads: 24 Joined: Oct 2008 10/15/2010, 09:20 PM (This post was last modified: 11/21/2011, 09:37 PM by sheldonison.) (10/15/2010, 04:03 PM)nuninho1980 Wrote: lastest your code for faster. I tried bases 1.5), 1.68, 1.7 and 100 and worked (not work on your older code). it's better. but I get errors after try init(1.6) and init(200). errors here: base 1.6 base 1.5 ....Hi Nuninho, I finally got around to debugging this (you reported it last time too). What's happening, is that the isuperf/superf equations starts misbehaving as the base approaches $\eta$. Below, there is a a graph for base 1.7, which just barely works. You can see the isuperf "jumps" when the sexp(z) starts getting slightly less than zero; sexp(z)<0, for z<-1. For the riemup code to work, the superf/isuperf needs to be consistent, over at least the unit circle. The problem is occurring in the last logarithm of the isuperf equation, after subtracting the fixed point "L". There is a term trying to counteract this multi-valued logarithm problem, by multiplying by c^n, before the last log_c, and this term works for bases>=1.7, but not for smaller bases, approaching eta. I'm checking in a "patch", that adds a new "xterminc", to push the function back into the well behaved region for bases B>=1.49. It still fails for B=1.48, but for other reasons, since the algorithm to generate "L" failed. The number of log/exp iterations required for the superf/isuperf for bases near $\eta$ starts to grow very large. In my program, this is the variable "scnt", and the algorithm for generating "scnt" does not work for bases, B<1.49. - Sheldon, small typo fixed in new code, uploaded a new version. Also, fixed so that it works down to B=1.47, but B<1.47 still fails. These bases close to eta run very slow. Code:xterminc = 0; /* xterminc used to stabilize isuperf for bases<1.7 */ if ((B<1.7),  xterminc=1); if ((B<1.55), xterminc=2); .... Code:init(1.7) ploth(t=-2,3,y=sexp(t+0.1*I);z=isuperf(y)-t;[real(y),real(z)]);     [attachment=780] For the most recent code version: go to the Nov 21st, 2011 thread. Attached Files   kneser.gp (Size: 15.78 KB / Downloads: 719) « Next Oldest | Next Newest »

 Messages In This Thread fast accurate Kneser sexp algorithm - by sheldonison - 08/07/2010, 06:53 AM The pari-GP code - by sheldonison - 08/07/2010, 09:17 PM RE: fast accurate Knesser sexp algorithm - by bo198214 - 08/08/2010, 03:54 PM RE: fast accurate Knesser sexp algorithm - by sheldonison - 08/08/2010, 06:46 PM RE: fast accurate Kneser sexp algorithm - by nuninho1980 - 08/15/2010, 12:09 AM 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 RE: fast accurate Kneser sexp algorithm - by sheldonison - 10/14/2010, 10:00 PM RE: fast accurate Kneser sexp algorithm - by nuninho1980 - 10/15/2010, 04:03 PM RE: fast accurate Kneser sexp algorithm - by sheldonison - 10/15/2010, 09:20 PM kneser.gp modified for large bases - by sheldonison - 10/19/2010, 03:33 AM small update to allow more flexibility, faster - by sheldonison - 10/30/2010, 09:47 PM update to support B

 Possibly Related Threads... Thread Author Replies Views Last Post Kneser-iteration on n-periodic-points (base say \sqrt(2)) Gottfried 11 4,300 05/05/2021, 04:53 AM Last Post: Gottfried "Kneser"/Riemann mapping method code for *complex* bases mike3 2 10,081 08/15/2011, 03:14 PM Last Post: Gottfried Attempt to make own implementation of "Kneser" algorithm: trouble mike3 9 24,359 06/16/2011, 11:48 AM Last Post: mike3 Numerical algorithm for Fourier continuum sum tetration theory mike3 12 30,594 09/18/2010, 04:12 AM Last Post: mike3 regular sexp: curve near h=-2 (h=-2 + eps*I) Gottfried 2 9,120 03/10/2010, 07:52 AM Last Post: Gottfried Attempting to compute the kslog numerically (i.e., Kneser's construction) jaydfox 11 28,874 10/26/2009, 05:56 PM Last Post: bo198214 regular sexp:different fixpoints Gottfried 6 18,379 08/11/2009, 06:47 PM Last Post: jaydfox sexp(strip) is winding around the fixed points Kouznetsov 8 20,723 06/29/2009, 10:05 AM Last Post: bo198214 sexp and slog at a microcalculator Kouznetsov 0 4,826 01/08/2009, 08:51 AM Last Post: Kouznetsov

Users browsing this thread: 1 Guest(s)