06/15/2011, 10:22 PM
(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 \( |\Im(z)| < 0.12 \), 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 \( |\Im(z)| = 0.12 \) 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.
