[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-gsl] Different value for mathieu_ce in Mathematica and GSL

From: Patrick Alken
Subject: Re: [Help-gsl] Different value for mathieu_ce in Mathematica and GSL
Date: Fri, 17 Feb 2017 15:17:57 -0700
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:45.0) Gecko/20100101 Thunderbird/45.7.0


I'm not familiar with the mathieu functions, but I entered your report into the bug tracker so hopefully someone with more knowledge than I can take a look and fix it.

On 02/17/2017 08:53 AM, Phyks wrote:

I have some code that I prototyped in Mathematica and am now writing in C using GSL, that makes use of Mathieu functions. I have different results between the two of them, and I cannot figure out whether this is a bug in GSL, Mathematica or simply some misunderstanding from my part.

I am using `MathieuC` function in latest Mathematica ( which should be the same function as `gsl_sf_mathieu_ce` ( except that the former one takes a single `a` argument being the characteristic value whereas the GSL 2.3 implementation takes the order `n` and the `q` parameter directly.

So, I guess,
N[MathieuC[MathieuCharacteristicA[0, -1], -1, 2*Pi/180]]

should be equivalent to
gsl_sf_mathieu_ce(0, -1.0, 2.0 * M_PI / 180.0)
which gives a totally different value: 0.99751942347886335.

I tried to debug with different values, and the discrepancies between Mathematica and GSL seems to appear only when the `q` parameter (-1.0 here) is negative. If I take 1.0 instead, I get values in agreement. I tried to find yet another implementation to debug it, and found Scipy ( which relies on Fortran SPECFUN library apparently, and is in agreement with GSL.

I am missing something? Thanks!

reply via email to

[Prev in Thread] Current Thread [Next in Thread]