|
|
Re: Bug in Jacobian Amplitude?
Posted:
Apr 6, 2013 12:10 PM
|
|
did schrieb:
>
> On Friday, 5 April 2013 03:00:38 UTC+2, Waldek Hebisch wrote:
> >
> > This definition is suboptimal, because it forces discontinuity
> > on real line. The definition I recall were continuous
> > (increasing) on real line.
>
> OK, but maple defines a discontinuous
> function over the real line. Try
> plot(JacobiAM(x,sqrt(3/4)),x=-10..20);
> plot(arctan(JacobiSN(x,sqrt(3/4)),JacobiCN(x,sqrt(3/4)))-JacobiAM(x,sqrt(3/4)),x=-10..20);
>
> MMA, on the other hand, gives a
> continuous increasing function.
>
> Does it mean that maple implement
> the Jacobi Amplitude via the suboptimal
> formula?
>
Having neither Maple nor Mathematica installed, I really have no idea if
the Maple implementation could be improved. Nor do I have the time to
implement the Jacobi amplitude am(u,k) in Derive and play around with
it. Hence an answer in general terms:
According to Gradsteyn & Ryzhik 8.14, am(u,k) possesses an infinite
number of branch points which form a two-dimensional lattice in the
complex plane:
u = 2*m*K + (2*n+1)*K'*#i
where K = K(k) and K' = K(k') = K(sqrt(1-k^2)) involve the complete
elliptic integral of the first kind. These branch points coincide with
poles of the Jacobian elliptic functions sn(u,k) = sin(am(u,k)), cn(u,k)
= cos(am(u,k)) and dn(u,k) = sqrt(1 - k^2*sin(am(u,k))^2). However, the
amplitude am(u,k) is periodic only along the imaginary axis with period
4*K'; along the real axis it has an additional linear component, in
accordance with d/du am(u,k) = dn(u,k), the period here being 2*K.
For numerical evaluation, am(u,k) must be restricted to a single-valued
function by introducing branch cuts. I looks desirable to me that the
pattern of cuts preserves both periods. The straightforward choice is
to connect all the branch points in pairs with the cuts approaching the
two points at u = +- K'*#i from opposite directions, but this still
leaves infinitely many possibilities: cuts aligned with the real axis;
cuts aligned with the imaginary axis; and finally diagonal orientations.
Once the orientation is decided on, periodicity along both axes would
generate all other cuts.
Obviously the Mathematica and Maple programmers have placed their branch
cuts differently, and the actual placement shouldn't be hard to
establish by color-coded 3D plots. I concur with Waldek that continuity
of am(u,k) on the real line is desirable, which would preclude cuts
crossing the real axis. Continuity should probably be abandoned if it
cannot be reconciled with having sn(u,k) = sin(am(u,k)) and cn(u,k) =
cos(am(u,k)) hold for all complex u and all complex k. The details of
all this remain for the Maple and Mathematica programmers to work out -
one would only wish they came up with the same answer!
Martin.
|
|
