Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
fixed point formula
#5
(10/01/2014, 04:43 AM)jaydfox Wrote: ...
Another benefit of a DFT (compared to, say, Lagrange or Newton interpolation), is that it uses considerably less memory, and it's quite a bit faster if you use an FFT (Fast Fourier Transform). I'll see what I can put together over the next few days.
....

Okay, here's a somewhat stable version of my DFT/FFT library, a slightly updated version of the fixed point finder, and a piece of code that calculates the fixed point Taylor series using a DFT (i.e., a "Cauchy integral").

I was able to calculate 512 terms to 96 digits of precision, in just a couple seconds. It took my other code (closed-form calculation of terms, using reversion of a power series) over 2 minutes to get the same result. I then calculated 1600 terms, again in only a few seconds. It took nearly two days using my other code.

Needless to say, in this specific scenario, the DFT method is extremely fast. And just as accurate.
.....
Note the performance of the FFT. I can use thousands of samples and still get a result in a few seconds. Heck, I can do tens of thousands of samples and still get a result in a reasonable timeframe. I did 8! (40320) samples in 11 seconds at 67 digits of precision, or 30 seconds at 144 digits.

I missed the posts in the computation section! This looks like great fun. I've never written an FFT, but I would like to. Its just that a generic DFT works pretty good within reasonable radius of convergence scenarios.

Back to the original topic, generating the fixed points from the . The form is also very interesting because it means the tetration itreration function switches from iterating to iterating . Here,

Another comment. Sometime after this post, I did succeed in getting a complex base tetration algorithm to work; its posted online here http://math.eretrandre.org/tetrationforu...hp?tid=729. But it doesn't converge for bases near the Shell Thron boundary. Its also dog slow for bases close to . So I'm thinking about a tetration algorithm that iterates generating the slog, centered exactly between the two fixed points, using the "k" form ... and using the very same fixed point Taylor series as from above; well actually the . The algorithm I have in mind should have convergences of , where I think r should be the real(pseudo period+1), so a 500 Taylor series terms for base e would be accurate to 15 decimal digits or so. Smaller values of k (base closer to eta) have larger real periods.... An FFT would help because it allow for thousands of terms instead of hundreds of terms.... So I may eventually need your FFT.
- Sheldon
Reply


Messages In This Thread
fixed point formula - by sheldonison - 02/25/2012, 06:35 AM
RE: fixed point formula - by jaydfox - 09/23/2014, 11:01 PM
RE: fixed point formula - by jaydfox - 09/28/2014, 09:22 PM
RE: fixed point formula - by jaydfox - 10/01/2014, 04:43 AM
RE: fixed point formula - by sheldonison - 03/11/2015, 01:24 PM
RE: fixed point formula - by sheldonison - 03/30/2015, 05:25 AM
RE: fixed point formula - by mike3 - 05/23/2015, 04:32 AM

Possibly Related Threads...
Thread Author Replies Views Last Post
  Attempt to find a limit point but each step needs doubling the precision... Gottfried 15 29,451 11/09/2014, 10:25 PM
Last Post: tommy1729
  Find all fixed points of exp[b] MorgothV8 10 23,738 10/07/2014, 11:00 AM
Last Post: Gottfried
  Road testing Ansus' continuum product formula mike3 33 48,306 09/23/2009, 08:26 PM
Last Post: mike3
  sexp(strip) is winding around the fixed points Kouznetsov 8 18,517 06/29/2009, 10:05 AM
Last Post: bo198214
  An error estimate for fixed point computation of b^x bo198214 0 3,643 05/31/2008, 04:11 PM
Last Post: bo198214



Users browsing this thread: 1 Guest(s)