preparing a new release

preparing a new release

[Brian Gough]

I'm preparing release (it will be 1.7) -- I'd appreciate feedback
on this test version:

http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz
http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz.sig

If there are no problems reported then I will release it. Thanks.

-- 
Brian Gough

Here are the changes:

** Switched gsl_randist_binomial to use the faster binomial random
variate TPE algorithm by default.  The previous binomial variate
algorithm is available as gsl_randist_binomial_knuth.  This will
result in a different sequence of binomial variates in programs using
this function.

** Improved the algorithm for gsl_sf_elljac_e to avoid cancellation
errors near quarter periods.

** Fixed the branch selection in gsl_sf_gamma_inc_Q_e to prevent
caused inaccurate results for large a,x where x~=~a.

** The multilinear fitting functions now have forms which accept a
user-specified tolerance for the SVD cutoff and return the
corresponding effective rank of the design matrix.

** The quadratic solvers in poly/ now handle linear equations
gracefully (i.e. quadratrics with a leading coefficient of zero).

** The output of "make check" now only shows test failures by default,
to reduce the amount of output.  Set the environment variable
GSL_TEST_VERBOSE=1 to display all the output.  To assist debugging,
the test number of each failure is shown in square brackets at the
line-end [NNNN].

** Fixed bugs in gsl_linalg_SV_decomp_jacobi which caused
incorrect results for some input matrices.

** Bessl, coulomb, dilogarithm and legendre_H3d functions now use
hypot internally to avoid overflow when computing terms like
sqrt(1+x*x).

** The 'Usage' chapter of the reference manual now explains how to
handle deprecated functions using the GSL_DISABLE_DEPRECATED macro.

** The conflicting enum definitions for 'forward' and 'backward' in
gsl_ftt.h and gsl_wavelet.h are deprecated.  User code should switch
to the new definitions gsl_fft_forward, gsl_fft_backward,
gsl_wavelet_forward and gsl_wavelet_backward. Selectively define
GSL_DISABLE_DEPRECATED before including the headers to use the new
definitions on either or both modules.

** Fixed an error in the the brent minimisation algorithm.  Iterations
should now follow Brent's original description correctly.

** The bound coulomb function gsl_sf_hydrogenicR_e no longer reports
an underflow for exact zeroes of the wavefunction.

** gsl_linalg_SV_decomp_jacobi now reports an error for the
unimplemented case M<N correctly.

** Fixed conformance test for the SYRK and HERK blas functions
gsl_blas_{s,d,c,z}syrk and gsl_blas_{c,z}herk for non-square matrices.

** Configure now checks for presence of ieeefp.h if needed.

** Differential equation solvers now propagate error codes returned
from user-defined functions to the top-level in all cases.

** Sort functions now avoid an infinite loop if Nans are present in
the input vector.  The order of nans in the output is undefined,
although other elements will be sorted correctly.

Re: preparing a new release

[James Theiler]

On Wed, 24 Aug 2005, Brian Gough wrote:

] I'm preparing release (it will be 1.7) -- I'd appreciate feedback
] on this test version:
] 
] http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz
] http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz.sig
] 
] If there are no problems reported then I will release it. Thanks.
] 
] 

configure; make; make check
all worked without a problem on RH Enterprise Linux 3

saw one typo in your list of changes:

** Fixed the branch selection in gsl_sf_gamma_inc_Q_e to prevent
caused inaccurate results for large a,x where x~=~a.

"prevent caused" -> "prevent"

===

By the way, I think there is an error in the documentation (code looks
okay) regarding Gaussian Tail Distribution, section 19.3.

Formula for p(x) should have an extra factor of 1/sqrt(2*pi*sigma^2).

Also, the figure for p(x) is incorrectly normalized -- you can
tell it is incorrect because the area under the curve should be one,
but is obviously much less than that.  It looks like that curve
was generated with a p(x) that *did* have the 1/sqrt(2*pi*sigma^2)
but did *not* have the 1/N(a;sigma) -- that is, it looks like the
tail of the gaussian in the Figure before it in section 19.2, without
any renormalization.

(An example of a correct picture of a tail distribution is in 19.10,
for the Rayleigh tail distribution.)


===

More on documentation, Chapter 7 introduces the Pochhammer symbol 
(a)_x with the function gsl_sf_poch.  Four functions later, in
gsl_sf_pochrel, a new notation is introduced (a,x).  If you want to
include (a,x) as an alternative notation, that could be done with
the gsl_sf_poch function; simpler would be to drop it entirely, and
say for gsl_sf_pochrel, ie:

These routines compute the relative Pochhammer symbol @math{((a,x) -
1)/x} where @math{(a,x) = (a)_x := \Gamma(a + x)/\Gamma(a)}.

  ->

These routines compute the relative Pochhammer symbol @math{((a)_x -
1)/x} where @math{(a)_x := \Gamma(a + x)/\Gamma(a)}.

=== 

The documentation in section 7.19 for gsl_sf_beta_inc does not define
B_x(a,b), which is
    B_x(a,b) = \int_0^x t^{a-1} (1-t)^{b-1} dt

===

On a more general note, I'd recommend changing the title of Section 
7.19 from "Gamma Function" to "Gamma and Beta Functions".

And this is more for the ToDo list, but we could use the introduction 
to this section to talk about factorials, combinations, binomial 
coefficients, use of tables for integer values, etc.

As a quickfix, we might just remark on the identity \Gamma(n)=(n-1)!
for integer n>0.  I don't think this is in that section at all.


regards,
jt

-- 
James Theiler            Space and Remote Sensing Sciences
MS-B244, ISR-2, LANL     Los Alamos National Laboratory
Los Alamos, NM 87545     http://nis-www.lanl.gov/~jt

Re: preparing a new release

[Alberto Manuel Brandão Simões]

Make check gives me a lot of errors...

Let me know if I should paste more.
Cheers.
Alberto

FAIL:   SV_decomp_jacobi (5x5) case=33427394 [166953]
call returned status = 14
singular value 2 =    1.41421356237309492 vs previous    1.27608201443340596
FAIL:   SV_decomp_jacobi (5x5) case=33437404 [166963]
call returned status = 14
FAIL:   SV_decomp_jacobi (5x5) case=33438405 [166964]
call returned status = 14
FAIL:   SV_decomp_jacobi (5x5) case=33458425 [166984]
call returned status = 14
FAIL:   SV_decomp_jacobi (5x5) case=33489456 [167015]
FAIL: test
===================
1 of 1 tests failed
===================
make[2]: *** [check-TESTS] Error 1
make[2]: Leaving directory `/home/ambs/Junk/INCOMING/gsl-1.6.91/linalg'
make[1]: *** [check-am] Error 2
make[1]: Leaving directory `/home/ambs/Junk/INCOMING/gsl-1.6.91/linalg'
make: *** [check-recursive] Error 1
[ambs@eremita gsl-1.6.91]$


James Theiler wrote:
> On Wed, 24 Aug 2005, Brian Gough wrote:
> 
> ] I'm preparing release (it will be 1.7) -- I'd appreciate feedback
> ] on this test version:
> ] 
> ] http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz
> ] http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz.sig
> ] 
> ] If there are no problems reported then I will release it. Thanks.
> ] 
> ] 
> 
> configure; make; make check
> all worked without a problem on RH Enterprise Linux 3
> 
> saw one typo in your list of changes:
> 
> ** Fixed the branch selection in gsl_sf_gamma_inc_Q_e to prevent
> caused inaccurate results for large a,x where x~=~a.
> 
> "prevent caused" -> "prevent"
> 
> ===
> 
> By the way, I think there is an error in the documentation (code looks
> okay) regarding Gaussian Tail Distribution, section 19.3.
> 
> Formula for p(x) should have an extra factor of 1/sqrt(2*pi*sigma^2).
> 
> Also, the figure for p(x) is incorrectly normalized -- you can
> tell it is incorrect because the area under the curve should be one,
> but is obviously much less than that.  It looks like that curve
> was generated with a p(x) that *did* have the 1/sqrt(2*pi*sigma^2)
> but did *not* have the 1/N(a;sigma) -- that is, it looks like the
> tail of the gaussian in the Figure before it in section 19.2, without
> any renormalization.
> 
> (An example of a correct picture of a tail distribution is in 19.10,
> for the Rayleigh tail distribution.)
> 
> 
> ===
> 
> More on documentation, Chapter 7 introduces the Pochhammer symbol 
> (a)_x with the function gsl_sf_poch.  Four functions later, in
> gsl_sf_pochrel, a new notation is introduced (a,x).  If you want to
> include (a,x) as an alternative notation, that could be done with
> the gsl_sf_poch function; simpler would be to drop it entirely, and
> say for gsl_sf_pochrel, ie:
> 
> These routines compute the relative Pochhammer symbol @math{((a,x) -
> 1)/x} where @math{(a,x) = (a)_x := \Gamma(a + x)/\Gamma(a)}.
> 
>   ->
> 
> These routines compute the relative Pochhammer symbol @math{((a)_x -
> 1)/x} where @math{(a)_x := \Gamma(a + x)/\Gamma(a)}.
> 
> === 
> 
> The documentation in section 7.19 for gsl_sf_beta_inc does not define
> B_x(a,b), which is
>     B_x(a,b) = \int_0^x t^{a-1} (1-t)^{b-1} dt
> 
> ===
> 
> On a more general note, I'd recommend changing the title of Section 
> 7.19 from "Gamma Function" to "Gamma and Beta Functions".
> 
> And this is more for the ToDo list, but we could use the introduction 
> to this section to talk about factorials, combinations, binomial 
> coefficients, use of tables for integer values, etc.
> 
> As a quickfix, we might just remark on the identity \Gamma(n)=(n-1)!
> for integer n>0.  I don't think this is in that section at all.
> 
> 
> regards,
> jt
> 

-- 
Alberto Simões - Departamento de Informática - Universidade do Minho
                  Campus de Gualtar - 4710-057 Braga - Portugal

Re: preparing a new release

[Brian Gough]

Alberto Manuel Brandão Simões writes:
 > FAIL:   SV_decomp_jacobi (5x5) case=33427394 [166953]
 > call returned status = 14
 > singular value 2 =    1.41421356237309492 vs previous    1.27608201443340596

Can you post your compiler version, platform details etc.  Thanks.

-- 
Brian Gough

Re: preparing a new release

[Aaron Schweiger]

Running make check I noticed that the script frequently was interrupted by
permission problems:

make[2]: Entering directory `/c/tmp/gsl-1.6.91/histogram'
/bin/sh: ./test.exe: Permission denied
FAIL: test.exe

Re-running make check allowed me to complete the sequence without incident.
This looks like a problem with MSYS - but I wanted to report it here in case
others run into it.

==================
All 1 tests passed
==================
make[2]: Leaving directory `/c/tmp/gsl-1.6.91'
make[1]: Leaving directory `/c/tmp/gsl-1.6.91'


I am running MSYS/MINGW:

Thread model: win32
gcc version 3.2.3 (mingw special 20030504-1).

On a side note, I am curious why we don't include a single "gsl.h"
that includes all the gsl headers - with the expection of gsl_cblas.h
(because it's exceptionally large)?  Most modern computers can handle 
the large
number of files in a trivial amount of time?

Aaron

Brian Gough wrote:

>I'm preparing release (it will be 1.7) -- I'd appreciate feedback
>on this test version:
>
>http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz
>http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz.sig
>
>If there are no problems reported then I will release it. Thanks.
>
>  
>

Re: preparing a new release

[albie]

Erm, 

 Slackware 10.1, gcc 3.2.2., glibc 3.2.3

Cheers
Alberrto

On Wed, Aug 24, 2005 at 06:06:28PM +0100, Brian Gough wrote:
| Alberto Manuel Brandão Simões writes:
|  > FAIL:   SV_decomp_jacobi (5x5) case=33427394 [166953]
|  > call returned status = 14
|  > singular value 2 =    1.41421356237309492 vs previous    1.27608201443340596
| 
| Can you post your compiler version, platform details etc.  Thanks.
| 
| -- 
| Brian Gough

-- 
Departamento de Informatica - Universidade do Minho

"The C Programming Language -- A language which combines the
flexibility of assembly language with the power of assembly language."

Re: preparing a new release

[Linas Vepstas]

> On Wed, 24 Aug 2005, Brian Gough wrote:
> 
> ] I'm preparing release (it will be 1.7) -- I'd appreciate feedback
> ] on this test version:
> ] 
> ] http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz
> ] http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz.sig
>
> configure; make; make check

Hi,

Seemed to work just fine on an IBM PowerPC (power5) running 
SLES9 (SuSE/Novell Linux Enterprise Server version 9) 
and 
RHEL4 (Red Hat Enterprise Linux version 4), 
as below:

================================================
The SLES9 machine: 

gcc --version:  gcc (GCC) 3.3.3 (SuSE Linux)

RPM files: rpm -qa |grep glibc

gcc-3.3.3-43.26.1
glibc-2.3.3-98.35.3
glibc-devel-2.3.3-98.35.3

cat /proc/cpuinfo
processor       : 0
cpu             : POWER5 (gr)
clock           : 1656.408000MHz
revision        : 2.1

processor       : 1
cpu             : POWER5 (gr)
clock           : 1656.408000MHz
revision        : 2.1

processor       : 2
cpu             : POWER5 (gr)
clock           : 1656.408000MHz
revision        : 2.1

processor       : 3
cpu             : POWER5 (gr)
clock           : 1656.408000MHz
revision        : 2.1

timebase        : 207051000
machine         : CHRP IBM,9111-520


=================================================
The RHEL4 machine:

gcc --version
gcc (GCC) 3.4.3 20050227 (Red Hat 3.4.3-22.1)

glibc-devel-2.3.4-2.9
glibc-2.3.4-2.9

cat /proc/cpuinfo
processor       : 0
cpu             : POWER5 (gs)
clock           : 1898.100000MHz
revision        : 2.0

processor       : 1
cpu             : POWER5 (gs)
clock           : 1898.100000MHz
revision        : 2.0

timebase        : 237601000
machine         : CHRP IBM,9111-520

=================================================

--linas

Re: preparing a new release

[Andrej Prsa]

Hi,

> I'm preparing release (it will be 1.7) -- I'd appreciate feedback
> on this test version:
> 
> http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz
> http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz.sig
> 
> If there are no problems reported then I will release it. Thanks.

No problems on i686, Debian testing (gcc 4.0.1). :)

Best wishes,

Andrej

Re: preparing a new release: 1 test failed MacOSX 10.4

[yoshiki9]

  One test failed on MacOS X 10.4 (Tiger).

yoshiki@hastings[6] uname -a
Darwin hastings 8.2.0 Darwin Kernel Version 8.2.0: Fri Jun 24  
17:46:54 PDT 2005; root:xnu-792.2.4.obj~3/RELEASE_PPC Power Macintosh  
powerpc
yoshiki@hastings[7] gcc -v
Reading specs from /usr/libexec/gcc/darwin/ppc/3.3/specs
Thread model: posix
gcc version 3.3 20030304 (Apple Computer, Inc. build 1809)

--
./configure ---> OK
make          ---> OK
make test   ---> 1 test failed at the wavelet section.

Making check in wavelet
make  test
gcc -DHAVE_CONFIG_H -I. -I. -I.. -I..    -g -O2 -c test.c
In file included from ../gsl/gsl_vector_complex_long_double.h:26,
                  from ../gsl/gsl_vector.h:4,
                  from test.c:6:
../gsl/gsl_complex.h:39: warning: use of `long double' type; its size  
may change
in a future release
../gsl/gsl_complex.h:39: warning: (Long double usage is reported only  
once for e
ach file.
../gsl/gsl_complex.h:39: warning: To disable this warning, use -Wno- 
long-double.
)
/bin/sh ../libtool --mode=link gcc  -g -O2   -o test  test.o  
libgslwavelet.la ..
/blas/libgslblas.la ../cblas/libgslcblas.la ../matrix/ 
libgslmatrix.la ../vector/
libgslvector.la ../block/libgslblock.la ../ieee-utils/ 
libgslieeeutils.la  ../err
/libgslerr.la ../test/libgsltest.la ../sys/libgslsys.la ../utils/ 
libutils.la -lm
gcc -g -O2 -o .libs/test test.o  ./.libs/libgslwavelet.al ../ 
blas/.libs/libgslbl
as.al -L../cblas/.libs -lgslcblas ../matrix/.libs/libgslmatrix.al ../ 
vector/.lib
...
creating test
make  check-TESTS
FAIL: bspline(309)-2d standard other data untouched, n = 2, tda = 3  
[6680]
FAIL: bspline-centered(309)-2d standard other data untouched, n = 2,  
tda = 3 [66
82]
FAIL: bspline(309)-2d nonstd other data untouched, n = 2, tda = 3 [6684]
FAIL: bspline-centered(309)-2d nonstd other data untouched, n = 2,  
tda = 3 [6686
]
FAIL: bspline(307)-2d standard other data untouched, n = 2, tda = 3  
[6688]
FAIL: bspline-centered(307)-2d standard other data untouched, n = 2,  
tda = 3 [66
90]
FAIL: bspline(307)-2d nonstd other data untouched, n = 2, tda = 3 [6692]
FAIL: bspline-centered(307)-2d nonstd other data untouched, n = 2,  
tda = 3 [6694
]
...
FAIL: haar(2)-2d standard other data untouched, n = 64, tda = 69 [12132]
FAIL: haar-centered(2)-2d standard other data untouched, n = 64, tda  
= 69 [12134
]
FAIL: haar(2)-2d nonstd other data untouched, n = 64, tda = 69 [12136]
FAIL: haar-centered(2)-2d nonstd other data untouched, n = 64, tda =  
69 [12138]
FAIL: test
===================
1 of 1 tests failed
===================
make[2]: *** [check-TESTS] Error 1
make[1]: *** [check-am] Error 2
make: *** [check-recursive] Error 1
(END)

-----
On 2005/08/24, at 23:12, Brian Gough wrote:

> I'm preparing release (it will be 1.7) -- I'd appreciate feedback
> on this test version:
>
> http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz
> http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz.sig
>
> If there are no problems reported then I will release it. Thanks.
>

--------------------------------------------------
  Yoshiki TSUNESADA

Re: preparing a new release

[Dirk Eddelbuettel]

Brian,

It built just fine in a Debian unstable chroot -- the resulting Debian
packages have been uploaded to Debian and reports of builds on the 11
non-i?86 platforms should trickle in over the next few days and will become
visible at  http://buildd.debian.org/build.php?&pkg=gsl

Regards, Dirk

-- 
Statistics: The (futile) attempt to offer certainty about uncertainty.
         -- Roger Koenker, 'Dictionary of Received Ideas of Statistics'

Re: preparing a new release

[Kevin Bube]

[-- Attachment #1: Type: text/plain, Size: 10508 bytes --]

Hello,

first the good news: './configure && make && make check' successful on
AMD64:

,----
| kevin@wien:~/tmp/2/gsl-1.6.91 $ uname -a; /lib64/libc.so.6; gcc -v
| Linux wien 2.6.5-7.147-default #1 Thu Jan 27 09:19:29 UTC 2005 x86_64 x86_64 x86_64 GNU/Linux
| GNU C Library stable release version 2.3.3 (20040325), by Roland McGrath et al.
| Copyright (C) 2004 Free Software Foundation, Inc.
| This is free software; see the source for copying conditions.
| There is NO warranty; not even for MERCHANTABILITY or FITNESS FOR A
| PARTICULAR PURPOSE.
| Configured for x86_64-suse-linux.
| Compiled by GNU CC version 3.3.3 (SuSE Linux).
| Compiled on a Linux 2.6.4 system on 2004-04-05.
| Available extensions:
|         GNU libio by Per Bothner
|         crypt add-on version 2.1 by Michael Glad and others
|         linuxthreads-0.10 by Xavier Leroy
|         GNU Libidn by Simon Josefsson
|         BIND-8.2.3-T5B
|         libthread_db work sponsored by Alpha Processor Inc
|         NIS(YP)/NIS+ NSS modules 0.19 by Thorsten Kukuk
| Thread-local storage support included.
| pthread library is compiled with floating stack support enabled.
| Report bugs using the `glibcbug' script to <bugs@gnu.org>.
| Reading specs from /usr/lib64/gcc-lib/x86_64-suse-linux/3.3.3/specs
| Configured with: ../configure --enable-threads=posix --prefix=/usr --with-local-prefix=/usr/local --infodir=/usr/share/info --mandir=/usr/share/man --enable-languages=c,c++,f77,objc,java,ada --disable-checking --libdir=/usr/lib64 --enable-libgcj --with-gxx-include-dir=/usr/include/g++ --with-slibdir=/lib64 --with-system-zlib --enable-shared --enable-__cxa_atexit x86_64-suse-linux
| Thread model: posix
| gcc version 3.3.3 (SuSE Linux)
`----

The bad news: make check fails on dual Alpha system:

,----
| kevin@alphaserver:/tmp/gsl-1.6.91 $ uname -a ; /lib/libc.so.6.1 ; gcc -v
| Linux alphaserver 2.4.4-SMP #2 SMP Sun May 6 08:08:09 GMT 2001 alpha unknown
| GNU C Library stable release version 2.2, by Roland McGrath et al.
| Copyright (C) 1992,93,94,95,96,97,98,99,2000 Free Software Foundation, Inc.
| This is free software; see the source for copying conditions.
| There is NO warranty; not even for MERCHANTABILITY or FITNESS FOR A
| PARTICULAR PURPOSE.
| Compiled by GNU CC version 2.96 20000731 (SuSE Linux 7.1/Alpha).
| Compiled on a Linux 2.4.0 system on 2001-04-23.
| Available extensions:
|         GNU libio by Per Bothner
|         crypt add-on version 2.1 by Michael Glad and others
|         Berkeley DB glibc 2.1 compat library by Thorsten Kukuk
|         linuxthreads-0.9 by Xavier Leroy
|         NoVersion patch for broken glibc 2.0 binaries
|         BIND-8.2.3-T5B
|         NIS(YP)/NIS+ NSS modules 0.19 by Thorsten Kukuk
|         NSS V1 modules 2.1 by Thorsten Kukuk
|         libthread_db work sponsored by Alpha Processor Inc
| Report bugs using the `glibcbug' script to <bugs@gnu.org>.
| Reading specs from /usr/lib/gcc-lib/alpha-suse-linux/2.96/specs
| gcc version 2.96 20000731 (SuSE Linux 7.1/Alpha)
`----

'./configure && make -j 2 && make check' fails with:

,----
| make  check-TESTS
| make[2]: Entering directory `/tmp/gsl-1.6.91/specfunc'
| FAIL: gsl_sf_coupling_9j_e(4, 2, 4, 3, 3, 2, 1, 1, 2, &r) [463]
|   expected:  -0.0408248290463863
|   obtained: -0.05103103630798286   1.133116629592098e-16  2.22045e-15
|   fracdiff:   0.1111111111111109
|   value/expected not consistent within reported error
|   value not within tolerance of expected value
|   -0.0510310363079828633  1.13311662959209844e-16
| FAIL: gsl_sf_coupling_9j_e(8, 4, 10, 7, 3, 8, 1, 1, 2, &r) [464]
|   expected:  0.02545875386086578
|   obtained:                    0                      0  nan
|   fracdiff:                    1
|   value/expected not consistent within reported error
|   value not within tolerance of expected value
|                        0                       0
| FAIL: Coupling Coefficients [483]
| FAIL: gsl_sf_legendre_sphPlm_e(10, 5, -0.5, &r) [1880]
|   expected:  -0.3725739049803294
|   obtained:                  nan                    nan  nan
|   fracdiff:                    0
|   value/expected not consistent within reported error
|                      nan                     nan
| FAIL: gsl_sf_legendre_sphPlm_e(10, 5, 1.0e-08, &r) [1881]
|   expected: -3.137723358937679e-08
|   obtained:                  nan                    nan  nan
|   fracdiff:                    0
|   value/expected not consistent within reported error
|                      nan                     nan
| FAIL: gsl_sf_legendre_sphPlm_e(10, 5, 0.5, &r) [1882]
|   expected:   0.3725739049803294
|   obtained:                  nan                    nan  nan
|   fracdiff:                    0
|   value/expected not consistent within reported error
|                      nan                     nan
| FAIL: gsl_sf_legendre_sphPlm_e(10, 5, 0.999, &r) [1883]
|   expected: -6.23687067472737e-06
|   obtained:                  nan                    nan  nan
|   fracdiff:                    0
|   value/expected not consistent within reported error
|                      nan                     nan
| FAIL: gsl_sf_legendre_sphPlm_e(200, 1, -0.5, &r) [1887]
|   expected:   0.3302975570099493
|   obtained:                  nan                    nan  nan
|   fracdiff:                    0
|   value/expected not consistent within reported error
|                      nan                     nan
| FAIL: gsl_sf_legendre_sphPlm_e(200, 1, 0.5, &r) [1888]
|   expected:  -0.3302975570099493
|   obtained:                  nan                    nan  nan
|   fracdiff:                    0
|   value/expected not consistent within reported error
|                      nan                     nan
| FAIL: gsl_sf_legendre_sphPlm_e(200, 1, 0.999, &r) [1889]
|   expected:   -1.406979205554626
|   obtained:                  nan                    nan  nan
|   fracdiff:                    0
|   value/expected not consistent within reported error
|                      nan                     nan
| FAIL: L[10]+0.0 [1891]
|   expected:  0.07452710323813559
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: L[95]+0.0 [1892]
|   expected:   0.2586535599088016
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: gsl_sf_legendre_sphPlm_array(100, 5, 0.5) [1893]
| FAIL: L[10]+0.0 [1895]
|   expected: 4.854315031308678e-14
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: L[95]+0.0 [1896]
|   expected: 8.313898496365084e-12
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: gsl_sf_legendre_sphPlm_array(100, 2, 1.0 - 2^(-50)) [1897]
| FAIL: L[95]+0.0 [1899]
|   expected: -8.313898496365084e-12
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: gsl_sf_legendre_sphPlm_array(100, 2, -1.0 + 2^(-50)) [1900]
| FAIL: DL[2]+0.0 [1917]
|   expected:    1.505354723003901
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: DL[10]+0.0 [1918]
|   expected:   0.7357655966864824
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: DL[98]+0.0 [1919]
|   expected:    28.44458995026438
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: gsl_sf_legendre_sphPlm_deriv_array(100, 2, 0.5) [1920]
| FAIL: DL[2]+0.0 [1923]
|   expected:  -0.9830996902900139
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: DL[3]+0.0 [1924]
|   expected:   -2.772827098895454
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: DL[10]+0.0 [1925]
|   expected:   -5.740713331544348
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: DL[95]+0.0 [1926]
|   expected:   -25.89393462474739
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: gsl_sf_legendre_sphPlm_deriv_array(100, 5, 0.5) [1927]
| FAIL: DL[2]+0.0 [1930]
|   expected: 1.648269720073467e-21
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: DL[95]+0.0 [1931]
|   expected: 3.989054946607135e-15
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: gsl_sf_legendre_sphPlm_deriv_array(100, 5, 1.0 - 2^(-50)) [1932]
| FAIL: DL[2]+0.0 [1935]
|   expected: -1.648269720073467e-21
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: DL[95]+0.0 [1936]
|   expected: 3.989054946607135e-15
|   obtained:                  nan
|   fracdiff:                    1
|   value not within tolerance of expected value
|                      nan
| FAIL: gsl_sf_legendre_sphPlm_deriv_array(100, 5, -1.0 + 2^(-50)) [1937]
| FAIL: Legendre Functions [2099]
| FAIL: test
| ===================
| 1 of 1 tests failed
| ===================
| make[2]: *** [check-TESTS] Error 1
| make[2]: Leaving directory `/tmp/gsl-1.6.91/specfunc'
| make[1]: *** [check-am] Error 2
| make[1]: Leaving directory `/tmp/gsl-1.6.91/specfunc'
| make: *** [check-recursive] Error 1
`----


Regards,

Kevin

-- 
publickey 2048R/0AFDFB19: http://www.icbm.de/~bube/publickey.asc
fingerprint: 542B 1378 04AA AF1F 572E  78BF 1BF5 5C71 0AFD FB19

[-- Attachment #2: Type: application/pgp-signature, Size: 480 bytes --]

Re: preparing a new release

[Ivo Alxneit]

[-- Attachment #1: Type: text/plain, Size: 1022 bytes --]

>On Wednesday 24 August 2005 16:12, Brian Gough wrote:
> I'm preparing release (it will be 1.7) -- I'd appreciate feedback
> on this test version:
>
> http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz
> http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz.sig
>
> If there are no problems reported then I will release it. Thanks.
everything ok on scientific linux
gcc-v
Reading specs from /usr/lib/gcc-lib/i386-redhat-linux/3.2.3/specs
Configured with: ../configure --prefix=/usr --mandir=/usr/share/man 
--infodir=/usr/share/info --enable-shared --enable-threads=posix 
--disable-checking --with-system-zlib --enable-__cxa_atexit 
--host=i386-redhat-linux
Thread model: posix
gcc version 3.2.3 20030502 (Red Hat Linux 3.2.3-52)

-- 
Dr. Ivo Alxneit
Laboratory for Solar Technology   phone: +41 56 310 4092
Paul Scherrer Institute             fax: +41 56 310 2688
CH-5232 Villigen                   http://solar.web.psi.ch
Switzerland                   gnupg key: 0x515E30C7

[-- Attachment #2: Type: application/pgp-signature, Size: 189 bytes --]

Re: preparing a new release

[Brian Gough]

albie@alfarrabio.di.uminho.pt writes:
 > Erm, 
 > 
 >  Slackware 10.1, gcc 3.2.2., glibc 3.2.3

Can you confirm the version of glibc, the current glibc series is 2.x
Thanks.

-- 
Brian Gough

Re: preparing a new release

[Alberto Manuel Brandão Simões]

Brian, sorry..

glibc 2.3.2 :) not 3.2.3 :)

Brian Gough wrote:
> albie@alfarrabio.di.uminho.pt writes:
>  > Erm, 
>  > 
>  >  Slackware 10.1, gcc 3.2.2., glibc 3.2.3
> 
> Can you confirm the version of glibc, the current glibc series is 2.x
> Thanks.
> 

-- 
Alberto Simões - Departamento de Informática - Universidade do Minho
                  Campus de Gualtar - 4710-057 Braga - Portugal

Re: preparing a new release: 1 test failed MacOSX 10.4

[Brian Gough]

yoshiki9@mac.com writes:
 >   One test failed on MacOS X 10.4 (Tiger).
 > 
 > yoshiki@hastings[6] uname -a
 > Darwin hastings 8.2.0 Darwin Kernel Version 8.2.0: Fri Jun 24  
 > 17:46:54 PDT 2005; root:xnu-792.2.4.obj~3/RELEASE_PPC Power Macintosh  
 > powerpc
 > yoshiki@hastings[7] gcc -v
 > Reading specs from /usr/libexec/gcc/darwin/ppc/3.3/specs
 > Thread model: posix
 > gcc version 3.3 20030304 (Apple Computer, Inc. build 1809)

Thanks. I think that is due to an optimisation bug in
the compiler (see the INSTALL file for more details).

-- 
Brian Gough

Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/

Re: preparing a new release

[Brian Gough]

Alberto Manuel Brandão Simões writes:
 > Brian, sorry..
 > 
 > glibc 2.3.2 :) not 3.2.3 :)
 > 
 > Brian Gough wrote:
 > > albie@alfarrabio.di.uminho.pt writes:
 > >  > Erm, 
 > >  > 
 > >  >  Slackware 10.1, gcc 3.2.2., glibc 3.2.3
 > > 
 > > Can you confirm the version of glibc, the current glibc series is 2.x

Ok, I'm still trying to reproduce this.  I downloaded Slackware 10.1
and installed it on a spare partition.  The versions of gcc and glibc
were 3.3.4 and 2.3.4 and they passed make check.  How exactly did you
put your system together?

-- 
Brian Gough

Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/

Re: preparing a new release

[Brian Gough]

Aaron Schweiger writes:
 > make[2]: Entering directory `/c/tmp/gsl-1.6.91/histogram'
 > /bin/sh: ./test.exe: Permission denied
 > FAIL: test.exe
 > 
 > Re-running make check allowed me to complete the sequence without incident.
 > This looks like a problem with MSYS - but I wanted to report it here in case
 > others run into it.

Strange that re-running the tests works -- if the file permissions are
wrong I wouldn't expect them to change.

 > On a side note, I am curious why we don't include a single "gsl.h"
 > that includes all the gsl headers - with the expection of gsl_cblas.h
 > (because it's exceptionally large)?  Most modern computers can handle 
 > the large
 > number of files in a trivial amount of time?

Some of us don't have modern computers :-(
Also it just seems better style (to me) to keep things modular.

-- 
Brian Gough

Re: preparing a new release

[Brian Gough]

Kevin Bube writes:
 > The bad news: make check fails on dual Alpha system:

Since there are some nans generated I'd suspect it might be a problem
with optimisation---you could try compiling with -O1 or -O0 to see if
it makes any difference.

-- 
Brian Gough

Re: preparing a new release

[Claudio Attaccalite]

configure; make; make check
all worked perfectly without any problem on Mandrake 10.2

with gcc 3.4.3 and CFLAGS="-O3 -march=athlon-xp -msse -mmmx -m3dnow -mfpmath=sse,387"

I tried  also with
icc version 9.0 and CFLAGS="-O3 -xK -ip -mcpu=pentiumpro"
I obtained a series of warnings like this:

vector_source.c(55): warning #290: function "gsl_vector_char_ptr" has 
already been defined
 FUNCTION (gsl_vector, ptr) (TYPE (gsl_vector) * v, const size_t i)

test_complex_source.c(479): warning #267: the format string requires 
additionalarguments
      FUNCTION (gsl_vector, fprintf) (f, v, OUT_FORMAT);

plus these two:

coulomb.c(994): warning #592: variable "G_lam_G" is used before its 
value is set
        G->val  = G_lam_G;
                  ^
coulomb.c(996): warning #592: variable "Gp_lam_G" is used before its 
value is set
        Gp->val = Gp_lam_G;
                  ^

and then make check remains stack here:

==================
All 1 tests passed
==================
make[2]: Leaving directory `/home/attacc/gsl-1.6.91/sort'
make[1]: Leaving directory `/home/attacc/gsl-1.6.91/sort'
Making check in ieee-utils
make[1]: Entering directory `/home/attacc/gsl-1.6.91/ieee-utils'
make  test
make[2]: Entering directory `/home/attacc/gsl-1.6.91/ieee-utils'
icc -DHAVE_CONFIG_H -I. -I. -I.. -I..    -O3 -xK -ip -mcpu=pentiumpro -c test.c
/bin/sh ../libtool --mode=link icc  -O3 -xK -ip -mcpu=pentiumpro   -o test  test.o libgslieeeutils.la ../err/libgslerr.la ../test/libgsltest.la ../sys/libgslsys.la ../utils/libutils.la -lm
icc -O3 -xK -ip -mcpu=pentiumpro -o test test.o  ./.libs/libgslieeeutils.al ../err/.libs/libgslerr.al ../test/.libs/libgsltest.al ../sys/.libs/libgslsys.al ../utils/.libs/libutils.al -lm
make[2]: Leaving directory `/home/attacc/gsl-1.6.91/ieee-utils'
make  check-TESTS
make[2]: Entering directory `/home/attacc/gsl-1.6.91/ieee-utils'

Ciao


Brian Gough ha scritto:

>I'm preparing release (it will be 1.7) -- I'd appreciate feedback
>on this test version:
>
>http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz
>http://www.network-theory.co.uk/download/gsl/gsl-1.6.91.tar.gz.sig
>
>If there are no problems reported then I will release it. Thanks.
>
>  
>


-- 
==========================================================
Claudio Attaccalite

S.I.S.S.A.
Via Beirut 2-4 34014 Grignano (TS)
email: attacc@sissa.it

ICQ number: 23970278
Skype: claudioattaccalite
web site: www.attaccalite.altervista.org
===========================================================
A quart of ale is a dish for a king.  
                                      -William Shakespeare                                           			 Winter's Tale
	

Re: preparing a new release

[Brian Gough]

keith.briggs@bt.com writes:
 > gsl-1.6.92.tar.gz passed all tests on AMD64 with gcc-4.0.1.
 > 
 > Keith
 > 
 > PS: perhaps the grammar in "All 1 tests passed" could be improved.

One for the Automake maintainers I think -- I've cc'ed this to
bug-automake@gnu.org.

The "make check" messages are from automake (maybe "All tests passed
(1/1)" would do the job).

-- 
Brian Gough

Re: preparing a new release

[keith.briggs]

gsl-1.6.92.tar.gz passed all tests on AMD64 with gcc-4.0.1.

Keith

PS: perhaps the grammar in "All 1 tests passed" could be improved.

Re: preparing a new release

[Brian Gough]

Reinh.Bader@t-online.de writes:
 >  the Intel compiler is another matter. If it works, the code
 >  performs much better on Itanium, but I need to specifically sort
 >  out which routines to compile with low optimization. Is there any
 >  further interest in icc from the side of the GSL developers?

You can send me the output of make check if you like.

-- 
Brian Gough

Re: preparing a new release

[Reinh.Bader]

  gcc -v
Reading specs from /usr/lib/gcc-lib/ia64-suse-linux/3.3.3/specs
Configured with: ../configure --enable-threads=posix --prefix=/usr --with-local-prefix=/usr/local --infodir=/usr/share/info
--mandir=/usr/share/man --enable-languages=c,c++,f77,objc,java,ada --disable-checking --libdir=/usr/lib --enable-libgcj
--with-gxx-include-dir=/usr/include/g++ --with-slibdir=/lib --with-system-zlib --enable-shared --enable-__cxa_atexit ia64-suse-linux
Thread model: posix
gcc version 3.3.3 (SuSE Linux)

 on SUSE Enterprise Linux 9 (glibc 2.3.3) for Itanium is OK (no failures).

 the Intel compiler is another matter. If it works, the code performs much better on Itanium, but I need to
 specifically sort out which routines to compile with low optimization. Is there any further interest in icc
 from the side of the GSL developers?


Best regards


-- 
 Dr. Reinhold Bader

 Leibniz-Rechenzentrum, Abt. Hochleistungssysteme | Tel. +49 89 289 28825
 Barerstr. 21, 80333 Muenchen                     | email Bader@lrz.de