Revisting my accelerated slog solution using Abel matrix inversion
#14
jaydfox Wrote:the accuracy you mentioned (76 digits)... I assumed this limitation in accuracy was in your solution.   But it finally dawned on me that this limitation was in ... the coefficients I posted.

Correct.  I regenerated my slog, and optimized the sample points to get as much accuracy as possible in the \( \theta_R \) mapping to compare Kneser slog and Jay's slog.

There is the limitation that your 800th term is 10^-126, so we need to keep the sample radius as small as possible.  Here, the \( \theta_\R \) sample points are with slog(z) where |z|<0.706, which is the radius of the red circle.  Sampling on the green dots sexp(0.431i-1.405...0.431i-0.405), we are limited to 35 complex conjugate pairs in the Fourier series, before the Fourier series becomes random noise.
   

With this Fourier series sampling, it is impossible to theoretically do better than 247 decimal digits since 10^-126*0.706^800~=10^-247.  And in fact, for a large chunk of the red circle of radius=0.706, I'm seeing accuracies of about 244 decimal digits with this 35 complex conjugate pair terms in the \( \theta_\R \) mapping.  All points within the red circle are accurate to 1.15E-213 with this theta mapping.  All points within a unit circle are accurate to within 1.5E-126, limited by the 800 terms in Jay's slog series.  It is intuitive to compare the error plot below to the red circle above, and the Fourier series sample points within the red circle.  Here is the error plot equation: precision is limited by the number of terms in Jays slog, and the number of terms in the \( \theta_\R \) mapping.
\( \frac{10^{-240}}{\text{JaySlog}(z)-\text{KneserSlog}(z)-\theta_\R(\text{KneserSlog})} \); 
   

The Kneser slog used was accurate to 276 decimal digits, which is the most accurate computation I've ever done.  A Fourier transform just on the real axis could have 322 decimal digits of consistency, but then all we would have is a line, instead of an enclosed region in the complex plane.
- Sheldon
Reply


Messages In This Thread
Analysis of Jay's slog vs Kneser - by sheldonison - 01/17/2019, 06:44 PM
RE: Analysis of Jay's slog vs Kneser - by jaydfox - 01/18/2019, 06:35 AM
RE: Analysis of Jay's slog vs Kneser - by jaydfox - 01/18/2019, 06:42 AM
RE: Analysis of Jay's slog vs Kneser - by jaydfox - 01/18/2019, 06:17 PM
RE: Revisting my accelerated slog solution using Abel matrix inversion - by sheldonison - 01/19/2019, 02:39 PM

Possibly Related Threads…
Thread Author Replies Views Last Post
  Terse Schroeder & Abel function code Daniel 1 3,792 10/16/2022, 07:03 AM
Last Post: Daniel
  Quickest way to compute the Abel function on the Shell-Thron boundary JmsNxn 0 3,766 04/21/2022, 01:52 AM
Last Post: JmsNxn
  The Promised Matrix Add On; Abel_M.gp JmsNxn 2 5,843 08/21/2021, 03:18 AM
Last Post: JmsNxn
  An incremental method to compute (Abel) matrix inverses bo198214 3 20,217 07/20/2010, 12:13 PM
Last Post: Gottfried
  A note on computation of the slog Gottfried 6 26,101 07/12/2010, 10:24 AM
Last Post: Gottfried
  Improving convergence of Andrew's slog jaydfox 19 67,492 07/02/2010, 06:59 AM
Last Post: bo198214
  intuitive slog base sqrt(2) developed between 2 and 4 bo198214 1 10,067 09/10/2009, 06:47 PM
Last Post: bo198214
  SAGE code for computing flow matrix for exp(z)-1 jaydfox 4 22,086 08/21/2009, 05:32 PM
Last Post: jaydfox
  sexp and slog at a microcalculator Kouznetsov 0 7,239 01/08/2009, 08:51 AM
Last Post: Kouznetsov
  Convergence of matrix solution for base e jaydfox 6 22,902 12/18/2007, 12:14 AM
Last Post: jaydfox



Users browsing this thread: 1 Guest(s)