Legendre 'make check' failure

Jari Häkkinen jari@chiralcomp.com
Fri Feb 20 10:16:00 GMT 2004

Hi everyone.

Just reporting a 'make check' error. 'make check' fails for 
'gsl_sf_conicalP_cyl_reg_e' when using -O2 optimization, but passes with 
no optimization. This issue has been reported some time ago (June/July 
2002) in this mailing list, but still exists (CVS February 19, 2004). In 
a previous mail dated '12 Nov 2002 13:02:29 -0700' Gerard Jungman wrote 
something about loosening tolerance for the test, and I followed these 
ideas for my local CVS copy for gsl.

The error reported from 'make check' is very small:

FAIL: gsl_sf_conicalP_cyl_reg_e(20, 1.0, 0.5, &r)
  expected: 7.066408031229072e-24
  obtained: 7.066408031229079e-24   1.880369398635769e-36  2.661e-13
  fracdiff: 5.198425890615777e-16
  value not within tolerance of expected value
  7.06640803122907946e-24  1.88036939863576925e-36

I've change the test tolerance on my CVS copy of gsl, and would like 
someone to include something similar to the original CVS tree:

#> cvs diff test_legendre.c
Index: test_legendre.c
RCS file: /cvs/gsl/gsl/specfunc/test_legendre.c,v
retrieving revision 1.29
diff -r1.29 test_legendre.c
<   TEST_SF(s, gsl_sf_conicalP_cyl_reg_e, (20, 1.0, 0.5, &r), 
7.066408031229072207e-24,     TEST_TOL0, GSL_SUCCESS);
 >   TEST_SF(s, gsl_sf_conicalP_cyl_reg_e, (20, 1.0, 0.5, &r),  
7.066408031229072207e-24,     TEST_TOL1, GSL_SUCCESS);

Thank you for listening,


I am compiling on a Linux SUSE 9 box. Compiler and processor information:

#> gcc --version
gcc (GCC) 3.3.2
Copyright (C) 2003 Free Software Foundation, Inc.
This is free software; see the source for copying conditions.  There is NO

#> cat /proc/cpuinfo
processor       : 0
vendor_id       : AuthenticAMD
cpu family      : 6
model           : 6
model name      : AMD Athlon(TM) XP2000+
stepping        : 2
cpu MHz         : 1679.056
cache size      : 256 KB
fdiv_bug        : no
hlt_bug         : no
f00f_bug        : no
coma_bug        : no
fpu             : yes
fpu_exception   : yes
cpuid level     : 1
wp              : yes
flags           : fpu vme de pse tsc msr pae mce cx8 sep mtrr pge mca 
cmov pat pse36 mmx fxsr sse syscall mmxext 3dnowext 3dnow
bogomips        : 3309.56

More information about the Gsl-discuss mailing list