[Bug-gsl] [bug #45925] Incomplete Gamma Functions flipped?

[Patrick Alken]

URL:
  <http://savannah.gnu.org/bugs/?45925>

                 Summary: Incomplete Gamma Functions flipped?
                 Project: GNU Scientific Library
            Submitted by: psa
            Submitted on: Fri 11 Sep 2015 10:07:55 AM GMT
                Category: Runtime error
                Severity: 3 - Normal
        Operating System: 
                  Status: None
             Assigned to: None
             Open/Closed: Open
                 Release: 
         Discussion Lock: Any

    _______________________________________________________

Details:

from vsakharuk =at= hbk =dot= com

Hello all,
Running 32 bit gsl on windows.
Got flipped results for incomplete gamma Functions.

#include <iomanip>
#include <iostream>
#include <gsl\gsl_specfunc.h>

using namespace std;

int main( int argc, char* argv[])
{
    cout<<"GSL incomplete gamma=" << gsl_sf_gamma_inc(1.000000, 3.000000)<<
endl;
    cout<<"GSL incomplete gamma P=" << gsl_sf_gamma_inc_P(1.000000,
3.000000)<< endl;
    return 0;
}

Output:
GSL incomplete gamma=0.0497871
GSL incomplete gamma P=0.950213

IMHO, GSL incomplete gamma should be equal to 0.950213
This is output of incomplete gamma P

Regards V.





    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?45925>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/