This is the mail archive of the gsl-discuss@sources.redhat.com mailing list for the GSL project.


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

Re: help with gsl_matrix_complex_row


-----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-----


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]