This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: help with gsl_matrix_complex_row
- From: Ivo Alxneit <ivo dot alxneit at psi dot ch>
- To: Brian Gough <bjg at network-theory dot co dot uk>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Tue, 20 May 2003 11:09:30 +0200
- Subject: Re: help with gsl_matrix_complex_row
- Organization: PSI / ENE / HTS
- References: <Pine.LNX.4.44.0305152328150.16448-100000@ptt.mi.infn.it> <200305191158.43653.ivo.alxneit@psi.ch> <16073.2415.364552.763652@debian.local>
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
ooops my initial response was not cc'ed to gsl-discuss so i include it here :
> >i get
> >row0(0) = 1, 1.1
> >row1(0) = -3.00652e-182, 7.7651e-62
> >row0(1) = 1.2, 1.3
> >row1(1) = -1.31435e+64, 3.39464e+184
> >on
> >gcc version 2.96 20000731 (Red Hat Linux 7.3 2.96-113)
> >
> >compiler problem?
> >
> >i checked with ddd. gsl_matrix m is correctly initialized but
> >vv1=gsl_matrix_complex_row fails.
> >
> > my initial response ends here
On Monday 19 May 2003 18:42, you wrote:
> Ivo Alxneit writes:
> > gcc version 2.96 20000731 (Red Hat Linux 7.3 2.96-113)
> >
> > compiler problem?
>
> Yes, I don't think the library passes "make check" with that compiler,
> because of this problem with complex numbers.
i just tried make check
make check stops with
make[2]: *** [check-TESTS] Error 1
make[2]: Leaving directory `/usr/local/src/gsl-1.3/eigen'
make[1]: *** [check-am] Error 2
make[1]: Leaving directory `/usr/local/src/gsl-1.3/eigen'
make: *** [check-recursive] Error 1
no FAIL's in gsl/{complex,vector,matrix,block} !
careful: gsl_matrix_complex_row() is not checked but gsl_matrix_complex
get_row PASSes!
i get a lot (1144) of FAIL's in gsl/eigen in the herm(4) and herm(10) part. i
dont include them here.
most of them are that i get completely different numbers than expected but
there are also 'interesting' cases
FAIL: herm(4), eigenvalue(1,0), real, unsorted
(-4.54283729443259248e-305 ob vs -1.77454581813773183e-307 expected)
FAIL: herm(4), eigenvalue(1,2), real, unsorted
(-1.77454581813773144e-307 obs vs -4.54283729443259349e-305 expected)
- -> (1,0) and (1,2) exchanged??
quite some ieee errors like
nan expected versus 0/1 observed
inf expected versus 0 observed
- --
Dr. Ivo Alxneit
Laboratory for Solar Technology phone: +41 56 310 4092
Paul Scherrer Institute fax: +41 56 310 2624
CH-5232 Villigen http://solar.web.psi.ch
Switzerland gnupg key: 0x515E30C7
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.2.0 (GNU/Linux)
iD8DBQE+yfDUAd7CE1FeMMcRAtPkAJ94X7YJQQ59QUpXeKGgD5H7MXP3jQCgvnpz
nFoOutfgDvhCnZ/amrL3HeU=
=zopR
-----END PGP SIGNATURE-----