• 1 Vote(s) - 5 Average
• 1
• 2
• 3
• 4
• 5
 fast accurate Kneser sexp algorithm sheldonison Long Time Fellow Posts: 664 Threads: 23 Joined: Oct 2008 01/06/2011, 02:53 PM (This post was last modified: 01/06/2011, 04:35 PM by sheldonison.) (01/06/2011, 02:48 AM)JmsNxn Wrote: Thank you for this. This is exactly what I need. I was wondering though, since sexp(5*i) = 0.999 + 4.999i is it safe for me to write sexp(fi) = 1 + fi and assume it's an error in approximation?Perhaps you didn't execute the "loop". Also, the current version gives a warning if you evaluate sexp(z) for imag(z)>I. If you download the kneser.gp program (I would suggest the most recent version, go to this thread), and type in: Code:init(exp(1));loop; .... goes through 13 iterations ..... (07:43) gp > sexp(5*I) !!! WARNING, riemaprx(z) much better than sexp(z) for imag(z)>I !!! %2 = 1.3786691576693131111676650899624 E40 - 3.2923562701722998997666622377240 E40*ITo get the correct result for imag(z)>I, use the riemaprx(z) instead, which is accurate for imag(z)>=0.12*I. Or use splicesexp(z), which is valid everywhere in the complex plane. Code:(07:43) gp > riemaprx(5*I) %3 = 0.32106749434792621140043570196732 + 1.3394349973320786642136709026508*II plan to post again, giving a more clear overview of the sexp and riemaprx routines and how they are generated, and how the different versions of the code improved convergence. The latest version leads to a definition of an infinite sequence of both functions. In particular, the current version allows for a continuous Cauchy integral for the sexp function, as opposed to the earlier version which was required to be strictly a discreet Cauchy integral estimation, which I think had more theoretical problems. I think the current version will help prove convergence, but I'm not there yet. - Shel « 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 2,141 05/05/2021, 04:53 AM Last Post: Gottfried "Kneser"/Riemann mapping method code for *complex* bases mike3 2 9,452 08/15/2011, 03:14 PM Last Post: Gottfried Attempt to make own implementation of "Kneser" algorithm: trouble mike3 9 22,774 06/16/2011, 11:48 AM Last Post: mike3 Numerical algorithm for Fourier continuum sum tetration theory mike3 12 28,493 09/18/2010, 04:12 AM Last Post: mike3 regular sexp: curve near h=-2 (h=-2 + eps*I) Gottfried 2 8,528 03/10/2010, 07:52 AM Last Post: Gottfried Attempting to compute the kslog numerically (i.e., Kneser's construction) jaydfox 11 26,962 10/26/2009, 05:56 PM Last Post: bo198214 regular sexp:different fixpoints Gottfried 6 17,101 08/11/2009, 06:47 PM Last Post: jaydfox sexp(strip) is winding around the fixed points Kouznetsov 8 19,240 06/29/2009, 10:05 AM Last Post: bo198214 sexp and slog at a microcalculator Kouznetsov 0 4,528 01/08/2009, 08:51 AM Last Post: Kouznetsov

Users browsing this thread: 1 Guest(s)