[SciPy-dev] hyper1f1, ticket #659

Moritz Helias helias@bccn.uni-freiburg...
Wed May 28 07:22:34 CDT 2008


Hello,

Johann Cohen-Tanugi wrote:

> hello,
> I confirm this on my linux box with a very recent version of scipy  
> from
> the bzr trunk. It seems somewhat related to
> http://www.scipy.org/scipy/scipy/ticket/419 ....

true. The problem seems to be the error estimation of Kahan's sum  
algorithm which is too optimistic in this case.

>
> Anyway, this should be fixed....
> Johann
>

special/specfun/specfun.f contains an implementation of M(a,b,x) which  
is called chgm.
This works fine for me. It just has to be wrapped (which I did on my  
machine).

Greetings,

Moritz

> Moritz Helias wrote:
>> Hi,
>>
>> I found some strange behavior of the hypergeometric function with  
>> real
>> arguments (hyp1f1). I do not know, who is the specialist on this
>> subject, but
>> could somebody please try to reproduce it (see ticket #659,
>> attached .py file): For large negative x, the function becomes quite
>> jagged.
>>
>> I also described a possible reason and a workaround in the ticket. My
>> work relies on this function, so currently I have to patch each of my
>> scipy installations to continue working.
>> If the bug can be reproduced by others, is there the possibility that
>> I provide a patch?
>>
>> Thanks in advance,
>>
>> greetings,
>>
>> Moritz
>> _______________________________________________
>> Scipy-dev mailing list
>> Scipy-dev@scipy.org
>> http://projects.scipy.org/mailman/listinfo/scipy-dev
>>
> _______________________________________________
> Scipy-dev mailing list
> Scipy-dev@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-dev



More information about the Scipy-dev mailing list