Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
Revisting my accelerated slog solution using Abel matrix inversion
#1
Happy New Year!

A few weeks ago, I started getting interested in revisiting my accelerated slog solution.  The basics of my method are discussed in this thread:
Improving convergence of Andrew's slog

A note on that thread.  The inverse Abel matrix approach was not discovered by me.  I don't know for sure who first used and published the method, but I'm reasonably certain that my first introduction to the method was through Andrew Robbins, so I give him the credit.  As far as I know, however, I was the first to try to accelerate convergence with the method described in the linked thread.  Indeed, I don't think I've seen it described anywhere else.  If someone knows of another person using this method, please let me know!

Anyway, I bring this all up, because I am once again playing with this solution.  I've modified my matrix solver to use Gaussian elimination in order to perform an LU decomposition.  Due to memory constraints, I don't actually save the L matrix, and I've modified the code to swap rows of the U matrix to disk.  This has allowed me to solve larger and larger systems, without RAM becoming a bottleneck.  The downside is a 20%-40% reduction in speed.  (I'm fortunate to have a SSD disk.  The performance might not be as tolerable on a traditional hard drive.)

I've managed to solve a 4096x4096 system with 7168 bits of floating point precision.  It took about six days.  Based on previous estimates of the convergence and system precision loss, the final result should have over 1500 bits of precision, which is far more than needed (accuracy is probably on the order of 80-100 bits for practical applications).

I will post the SAGE code soon.  The forum isn't letting me attach it, so I might need to change the file name/format.  At any rate, it needs to be cleaned up a bit.

However, I will attach the first 1,000 coefficients of my accelerated 4096x4096 matrix solution.  Why only a thousand?  Well, two reasons.  First of all, based on some meta-analysis I've done (which I'll post in a few weeks when I have the time to write it up properly), I believe only the first 600-700 terms are accurate anyway.  For example, if you were wanting to perform a series reversion to get the sexp function, you would pretty much get "noise" after about the 600th term.

The second reason is that the forum is limiting me to 200kB file size for attachments.  The full 4096 terms was well over the limit.  So I had to truncate the file anyway.  As long as I was truncating it, I figured I'd limit it to something reasonable.

Please note, these coefficients are for the power series for the slog residue (i.e., the slog function minus the logarithms at the primary fixed points).  I know, "residue" is already a well-defined mathematical concept.  I still don't have a better name for it yet.  It's what I called it in my original posts on this subject over 11 years ago, so there's no sense changing the name flippantly.  I will work on finding a better name.  Anyway, calculate the taylor series for the logarithms at the primary fixed points, and add it to this series.  For "practical" calculations, just use the logarithms, and add this function.  I'll maybe devote a forum post later to demonstrate.


Attached Files
.txt   slog_coeffs_4096_acc_1000terms.txt (Size: 81.91 KB / Downloads: 88)
~ Jay Daniel Fox
Reply


Messages In This Thread
Revisting my accelerated slog solution using Abel matrix inversion - by jaydfox - 01/06/2019, 12:54 AM
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

Possibly Related Threads...
Thread Author Replies Views Last Post
  An incremental method to compute (Abel) matrix inverses bo198214 3 8,686 07/20/2010, 12:13 PM
Last Post: Gottfried
  A note on computation of the slog Gottfried 6 8,995 07/12/2010, 10:24 AM
Last Post: Gottfried
  Improving convergence of Andrew's slog jaydfox 19 23,848 07/02/2010, 06:59 AM
Last Post: bo198214
  intuitive slog base sqrt(2) developed between 2 and 4 bo198214 1 3,357 09/10/2009, 06:47 PM
Last Post: bo198214
  SAGE code for computing flow matrix for exp(z)-1 jaydfox 4 7,699 08/21/2009, 05:32 PM
Last Post: jaydfox
  sexp and slog at a microcalculator Kouznetsov 0 2,973 01/08/2009, 08:51 AM
Last Post: Kouznetsov
  Convergence of matrix solution for base e jaydfox 6 8,249 12/18/2007, 12:14 AM
Last Post: jaydfox
  Computing Abel function at a given center jaydfox 10 11,299 11/30/2007, 06:44 PM
Last Post: andydude
  Matrix-method: compare use of different fixpoints Gottfried 23 22,151 11/30/2007, 05:24 PM
Last Post: andydude
  SAGE code implementing slog with acceleration jaydfox 4 6,695 10/22/2007, 12:59 AM
Last Post: jaydfox



Users browsing this thread: 1 Guest(s)