[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Gauss's hypergeometric function
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Gauss's hypergeometric function |
Date: |
Tue, 20 Sep 2011 16:27:18 +0100 |
User-agent: |
Wanderlust/2.15.6 (Almost Unreal) Emacs/23.3 Mule/6.0 (HANACHIRUSATO) |
At Sat, 20 Aug 2011 04:58:42 -0700 (PDT),
Hatef Monajemi wrote:
>
> More specifically, I am trying the following:
> ======================
> for (size_t i=1; i<5 ; ++i){
> double a = gsl_sf_hyperg_2F1 (i-0.5, -0.5, i+1, 0.9);
> std::cout<< "a:"<< a <<std::endl;
> }
>
> ======================
> Note that when 'z < .5', I get the right results. for .5<z<1.0 which is the
> case for me, the results are different than those of Matlab. I am assuming
> Matlab implementation is correct. I do not have access to Mathematica though.
> It seems for 'c-a-b = m' when m is an integer and 0.5<z<1.0 the code is not
> working properly. Any comment?
>
Thanks for the bug report. I agree there is a problem with those
values. I have added a test case for the problem, but I don't have a
fix.
--
Brian Gough
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- Re: [Bug-gsl] Gauss's hypergeometric function,
Brian Gough <=