Changeset 1287 for trunk

Timestamp:

Apr 24, 2008, 11:24:12 PM (15 years ago)

Author:

Peter

Message:

Fixes #87 - replace cumulative hypergeometric with GSL implementation
NOTE: This means that yat 0.5 will require GSL 1.8 or
later. configure.ac is updated to check for that.

Location:

trunk

Files:

• NEWS (modified) (1 diff)
• README (modified) (1 diff)
• configure.ac (modified) (1 diff)
• yat/statistics/utility.cc (modified) (1 diff)

trunk/NEWS

@@ -4 +4 @@
 
 Version 0.5 (released DATE)
+
+  - yat now requires GSL 1.8 or later
 
 Version 0.4 (released 11 April 2008)

trunk/README

@@ -35 +35 @@
 
 GNU Scientific Library, [http://www.gnu.org/software/gsl/ GSL] version
-1.6 or later.
+1.8 or later.
 
 === BLAS ===

trunk/configure.ac

@@ -84 +84 @@
 # Including AM_PATH_GSL macro from gsl.m4 distributed by GSL
 gsl_found="no"
-gsl_version="1.6"
+gsl_version="1.8"
 sinclude(./build_support/gsl.m4)
 AX_PATH_GSL(${gsl_version},gsl_found="yes",gsl_version_check="no")

trunk/yat/statistics/utility.cc

@@ -40 +40 @@
                               unsigned int n2, unsigned int t)
   {
-    double p=0;
-    for (size_t i=0; i<=k; i++)
-      p+= gsl_ran_hypergeometric_pdf(i, n1, n2, t);
-    return p;
+    return gsl_cdf_hypergeometric_P(k, n1, n2, t);
   }