[Help-gsl] Problem with gsl_cdf_gaussian_X

[Frédéric Schoenahl]

Hello,
I posted a similar message on gsl-discuss by mistake.. please apologize! I had a
problem with the gsl_cdf_gaussian_P function on two stations with a simple
program  :

----8<------------------------------------------------
#include <gsl/gsl_randist.h>
int
main(int argc, char *argv[]) {
  double i = -10.;
  while(i<10.) {
    printf("%f %f\n", i, gsl_cdf_gaussian_P(i, 2.));
    i+=0.1;
  }
  return 1;
}
----8<------------------------------------------------

With gsl-1.6 on linux with gcc 3.3.4 I get a buggy stdout (stripped for display)
here is a sample :


...
8.400000 -9255875844002258...351113699328.000000
8.500000 nan
8.600000 0.000000
8.700000 19035747065432754...926242906112.000000
8.800000 -0.000000
8.900000 -9255870940004835...966634991616.000000
9.000000 nan
9.100000 0.000000
9.200000 19035734098215151...708181815296.000000
9.300000 -0.000000
9.400000 -9255866036026119...541521784832.000000
9.500000 nan
9.600000 0.000000
...

on sparc-solaris with gsl-1.5 and gcc-3.4.1 I have

...
8.400000 4.2000000
8.500000 4.2500000
8.600000 4.3000000
8.700000 4.3500000
8.800000 4.4000000
8.900000 4.4500000
9.000000 4.5000000
9.100000 4.5500000
9.200000 4.6000000
9.300000 4.6500000
9.400000 4.7000000
9.500000 4.7500000
9.600000 4.8000000
...

i.e. gsl_cdf_gaussian(x, s) = x/s !

Do you have the same problem?
Thanks for your time,

Fred

-- 
Use the list email to answer !