Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
Attempt to make own implementation of "Kneser" algorithm: trouble
(06/15/2011, 03:18 PM)sheldonison Wrote: The first thing I notice is that your "step 2" is trying to compute the coefficients for theta, via the Fourier integral at the real axis, from -1 to zero. Theta(z) has a really nasty singularity at integer coefficient. Analytic Fourier analysis works really great with very minimal number of terms if you're generating the Fourier analysis of an analytic function, but the closer to the singularity you are, the more of a problem you have. As noted in the post above, I got around this by computing the integral at imag(z)=0.12i. I actually computed theta(z) from -0.5+0.12i to +0.5+0.12i, because that was centered very close to the middle of my sexp(z) approximation, where it was most accurate. 0.12i is far enough away from the singularity that the theta(z) function is very well behaved, and requires minimal terms in the Fourier analysis. fyi, my sexp(z) function is sampled in a unit circle around z=0. But the price you pay is that theta(z) now only converges for imag(z)>=0.12i. I went through many many iterations about how to deal with that, originally using a discreet Cauchy sampling, before I got to the current implementation, that splices together the Kneser/Riemann approximation with log(z)/exp(z), from the previous iterations sexp(z) approximation, which works really nicely. So that would be a change to your step (3) for your Cauchy integral.

So you would use a countour like that of Kouznetsov's Cauchy method, right (i.e. half-circle connected to another with straight lines to form a pill- or oval-like shape), and when , you switch to integrating over the Taylor appproximation from the last run? I'm not sure what you mean by "log/exp" of it, though -- if you're integrating with zero at the center, the Taylor already converges best there, so I don't see why you need to take log/exp.

And I guess that why it doesn't converge below is that the function is not perfectly periodic initially, so when the area integrated over is repeated, it forms something not analytic along the given parallel. But I suppose this goes away as the algorithm progresses, and so eventually the theta mapping converges right up to near the real axis.

Messages In This Thread
RE: Attempt to make own implementation of "Kneser" algorithm: trouble - by mike3 - 06/15/2011, 10:22 PM

Possibly Related Threads...
Thread Author Replies Views Last Post
  Kneser-iteration on n-periodic-points (base say \sqrt(2)) Gottfried 11 5,107 05/05/2021, 04:53 AM
Last Post: Gottfried
  fast accurate Kneser sexp algorithm sheldonison 38 117,949 01/14/2016, 05:05 AM
Last Post: sheldonison
  Attempt to find a limit point but each step needs doubling the precision... Gottfried 15 34,430 11/09/2014, 10:25 PM
Last Post: tommy1729
  "Kneser"/Riemann mapping method code for *complex* bases mike3 2 10,289 08/15/2011, 03:14 PM
Last Post: Gottfried
  Numerical algorithm for Fourier continuum sum tetration theory mike3 12 31,288 09/18/2010, 04:12 AM
Last Post: mike3
  Attempting to compute the kslog numerically (i.e., Kneser's construction) jaydfox 11 29,523 10/26/2009, 05:56 PM
Last Post: bo198214

Users browsing this thread: 1 Guest(s)