From: Ada novice <shai.lesh@gmx.com>
Subject: Re: ANN: Ada 2005 Math Extensions 20120712
Date: Sat, 28 Jul 2012 06:01:26 -0700 (PDT)
Date: 2012-07-28T06:01:26-07:00 [thread overview]
Message-ID: <e1d5d31c-f873-4910-8bf3-6700a17007d1@googlegroups.com> (raw)
In-Reply-To: <m2mx2kq5me.fsf@nidhoggr.home>
I apologise for the late reply. Thanks for the information. ACML works fine also by just setting the Linker_Option as "-lacml". GREAT!
To compare ACML and the "Debian" lapack:
With Debian lapack:
$ ldd demo_extensions
linux-gate.so.1 => (0xb7756000)
liblapack.so.3gf => /usr/lib/liblapack.so.3gf (0xb6ef2000)
libblas.so.3gf => /usr/lib/libblas.so.3gf (0xb6be9000)
libc.so.6 => /lib/i386-linux-gnu/i686/cmov/libc.so.6 (0xb6a8b000)
libpthread.so.0 => /lib/i386-linux-gnu/i686/cmov/libpthread.so.0 (0xb6a72000)
libgfortran.so.3 => /usr/lib/i386-linux-gnu/libgfortran.so.3 (0xb696f000)
libgcc_s.so.1 => /lib/i386-linux-gnu/libgcc_s.so.1 (0xb6952000)
libm.so.6 => /lib/i386-linux-gnu/i686/cmov/libm.so.6 (0xb692c000)
/lib/ld-linux.so.2 (0xb7757000)
libquadmath.so.0 => /usr/lib/i386-linux-gnu/libquadmath.so.0 (0xb68b8000)
and with ACML:
$ ldd demo_extensions
linux-gate.so.1 => (0xb777c000)
libacml.so => /opt/acml4.4.0/gfortran32/lib/libacml.so (0xb6d96000)
libc.so.6 => /lib/i386-linux-gnu/i686/cmov/libc.so.6 (0xb6c0a000)
libpthread.so.0 => /lib/i386-linux-gnu/i686/cmov/libpthread.so.0 (0xb6bf1000)
librt.so.1 => /lib/i386-linux-gnu/i686/cmov/librt.so.1 (0xb6be8000)
libgfortran.so.3 => /usr/lib/i386-linux-gnu/libgfortran.so.3 (0xb6ae5000)
libm.so.6 => /lib/i386-linux-gnu/i686/cmov/libm.so.6 (0xb6abf000)
/lib/ld-linux.so.2 (0xb777d000)
libquadmath.so.0 => /usr/lib/i386-linux-gnu/libquadmath.so.0 (0xb6a4b000)
I'm curious to know this "libquadmath" that appears in both of the above outputs at the end. What is it?
When running "demo_extensions", I saw some differences in the outputs (in the last example for instance). It is hard to compare on the screen; I will have to print and see the differences well. The outputs are:
(a) With Debian lapack:
$ ./demo_extensions
--------------------------------
Values from <143ef70b-7e74-426b-a621-a5fd157849be@x21g2000yqa.googlegroups.com>
42 => ( 2.00000, 4.00000)
43 => ( 2.00000,-4.00000)
44 => ( 1.00000,-0.00000)
--------------------------------
Values in Test16 of http://people.sc.fsu.edu/~jburkardt/f_src/lapack/lapack_OSX_prb_output.txt
using Complex_Arrays.Eigenvalues
6.00000
4.00000
2.00000
-0.00000
-2.00000
-4.00000
-6.00000
using Extensions.Eigenvalues
( 6.00000, 0.00000)
( 4.00000, 0.00000)
(-6.00000, 0.00000)
( 2.00000, 0.00000)
(-0.00000, 0.00000)
(-4.00000, 0.00000)
(-2.00000, 0.00000)
--------------------------------
Values from http://en.wikipedia.org/wiki/Skew-symmetric_matrix
( 0.00000, 4.58258)
(-0.00000, 0.00000)
( 0.00000,-4.58258)
--------------------------------
Results from http://en.wikipedia.org/wiki/Orthogonal_matrix
Eigenvalues:
( 0.00000, 1.00000)
(-0.00000,-1.00000)
(-1.00000, 0.00000)
Eignesystem Values:
( 0.00000, 1.00000)
(-0.00000,-1.00000)
(-1.00000, 0.00000)
Eigensystem Vectors:
( 0.70711, 0.00000) ( 0.70711, 0.00000) (-0.00000, 0.00000)
( 0.00000,-0.56569) (-0.00000, 0.56569) (-0.60000,-0.00000)
( 0.00000,-0.42426) (-0.00000, 0.42426) ( 0.80000, 0.00000)
--------------------------------
Generalized eigensystem of real non-symmetric matrix.
The solutions are such that beta*a - alpha*b is singular, ie
its determinant is zero. We'll show that it's small for
a selection of randomly-benerated matrices.
j: 1 determinant: 0.00000E+00
j: 2 determinant: 6.41364E-07
j: 3 determinant: 1.09388E-09
j: 4 determinant: 3.61017E-08
j: 5 determinant: 1.99069E-08
j: 6 determinant: 1.59409E-08
j: 1 determinant: 2.73842E-08
j: 2 determinant: 1.97363E-09
j: 3 determinant: 2.46558E-09
j: 4 determinant: 1.55393E-05
j: 5 determinant: 4.60858E-08
j: 6 determinant: 3.19744E-09
j: 1 determinant: 2.77505E-08
j: 2 determinant: 1.87194E-08
j: 3 determinant: 6.03069E-09
j: 4 determinant: 1.63023E-09
j: 5 determinant: 8.51147E-08
j: 6 determinant: 4.31280E-09
j: 1 determinant: 7.87980E-10
j: 2 determinant: 2.10041E-10
j: 3 determinant: 3.72882E-11
j: 4 determinant: 1.11441E-10
j: 5 determinant: 3.37806E-08
j: 6 determinant: 6.60784E-12
j: 1 determinant: 2.75188E-10
j: 2 determinant: 0.00000E+00
j: 3 determinant: 1.48421E-13
j: 4 determinant: 1.43467E-09
j: 5 determinant: 7.14908E-10
j: 6 determinant: 3.36777E-10
j: 1 determinant: 4.88385E-09
j: 2 determinant: 4.37879E-07
j: 3 determinant: 2.77873E-08
j: 4 determinant: 9.66340E-09
j: 5 determinant: 2.69842E-08
j: 6 determinant: 4.55744E-08
--------------------------------
ZGGEV example at http://www.nag.co.uk/lapack-ex/node122.html
Eigenvalue( 1) = ( 3.00000E+00,-9.00000E+00)
Eigenvector( 1) = (-8.37379E-01,-1.62621E-01)(-1.53495E-01, 7.44662E-02)(-7.44662E-02,-1.53495E-01)( 1.53495E-01,-7.44662E-02)
Eigenvalue( 2) = ( 2.00000E+00,-5.00000E+00)
Eigenvector( 2) = ( 6.29583E-01, 3.70417E-01)( 4.14831E-03,-4.65165E-04)( 3.95740E-02, 2.32833E-02)(-2.32833E-02, 3.95740E-02)
Eigenvalue( 3) = ( 3.00000E+00,-1.00000E+00)
Eigenvector( 3) = ( 9.77535E-01, 2.24654E-02)( 1.59101E-01,-1.13710E-01)( 1.20898E-01,-1.53710E-01)( 1.53710E-01, 1.20898E-01)
Eigenvalue( 4) = ( 4.00000E+00,-5.00000E+00)
Eigenvector( 4) = (-9.06236E-01, 9.37639E-02)(-7.43055E-03, 6.87512E-03)( 3.02078E-02,-3.12554E-03)(-1.45858E-02,-1.40970E-01)
(b) with ACML
$ ./demo_extensions
--------------------------------
Values from <143ef70b-7e74-426b-a621-a5fd157849be@x21g2000yqa.googlegroups.com>
42 => ( 2.00000, 4.00000)
43 => ( 2.00000,-4.00000)
44 => ( 1.00000,-0.00000)
--------------------------------
Values in Test16 of http://people.sc.fsu.edu/~jburkardt/f_src/lapack/lapack_OSX_prb_output.txt
using Complex_Arrays.Eigenvalues
6.00000
4.00000
2.00000
-0.00000
-2.00000
-4.00000
-6.00000
using Extensions.Eigenvalues
( 6.00000, 0.00000)
( 4.00000, 0.00000)
(-6.00000, 0.00000)
( 2.00000, 0.00000)
( 0.00000, 0.00000)
(-4.00000, 0.00000)
(-2.00000, 0.00000)
--------------------------------
Values from http://en.wikipedia.org/wiki/Skew-symmetric_matrix
( 0.00000, 4.58258)
( 0.00000, 0.00000)
(-0.00000,-4.58258)
--------------------------------
Results from http://en.wikipedia.org/wiki/Orthogonal_matrix
Eigenvalues:
( 0.00000, 1.00000)
( 0.00000,-1.00000)
(-1.00000, 0.00000)
Eignesystem Values:
( 0.00000, 1.00000)
( 0.00000,-1.00000)
(-1.00000, 0.00000)
Eigensystem Vectors:
( 0.70711, 0.00000) ( 0.70711, 0.00000) (-0.00000,-0.00000)
( 0.00000,-0.56569) ( 0.00000, 0.56569) (-0.60000,-0.00000)
( 0.00000,-0.42426) ( 0.00000, 0.42426) ( 0.80000, 0.00000)
--------------------------------
Generalized eigensystem of real non-symmetric matrix.
The solutions are such that beta*a - alpha*b is singular, ie
its determinant is zero. We'll show that it's small for
a selection of randomly-benerated matrices.
j: 1 determinant: 3.38316E-08
j: 2 determinant: 3.48750E-06
j: 3 determinant: 7.27814E-10
j: 4 determinant: 5.19989E-08
j: 5 determinant: 1.88069E-08
j: 6 determinant: 1.16859E-08
j: 1 determinant: 3.75126E-09
j: 2 determinant: 9.37482E-09
j: 3 determinant: 7.31456E-08
j: 4 determinant: 7.91349E-06
j: 5 determinant: 3.26590E-09
j: 6 determinant: 1.97660E-08
j: 1 determinant: 2.24647E-08
j: 2 determinant: 9.01645E-08
j: 3 determinant: 1.23797E-08
j: 4 determinant: 4.89069E-09
j: 5 determinant: 7.94987E-08
j: 6 determinant: 6.34529E-09
j: 1 determinant: 3.00981E-09
j: 2 determinant: 7.41406E-10
j: 3 determinant: 2.78002E-10
j: 4 determinant: 2.11384E-10
j: 5 determinant: 6.75613E-08
j: 6 determinant: 3.30392E-11
j: 1 determinant: 1.10075E-09
j: 2 determinant: 1.52874E-08
j: 3 determinant: 3.25786E-11
j: 4 determinant: 5.33539E-09
j: 5 determinant: 2.72676E-09
j: 6 determinant: 1.34711E-09
j: 1 determinant: 8.40869E-09
j: 2 determinant: 6.92644E-07
j: 3 determinant: 3.07913E-08
j: 4 determinant: 8.94615E-10
j: 5 determinant: 1.75984E-08
j: 6 determinant: 5.06383E-08
--------------------------------
ZGGEV example at http://www.nag.co.uk/lapack-ex/node122.html
Eigenvalue( 1) = ( 3.00000E+00,-9.00000E+00)
Eigenvector( 1) = ( 5.22150E-01, 4.77850E-01)( 1.40886E-01, 1.37983E-02)(-1.37980E-02, 1.40886E-01)(-1.40886E-01,-1.37980E-02)
Eigenvalue( 2) = ( 2.00000E+00,-5.00000E+00)
Eigenvector( 2) = (-9.41937E-01, 5.80634E-02)(-4.10722E-03, 3.49492E-03)(-5.92076E-02, 3.64964E-03)(-3.64968E-03,-5.92078E-02)
Eigenvalue( 3) = ( 3.00001E+00,-1.00000E+00)
Eigenvector( 3) = ( 9.37074E-01, 6.29255E-02)( 1.57483E-01,-1.02381E-01)( 1.22517E-01,-1.42381E-01)( 1.42381E-01, 1.22517E-01)
Eigenvalue( 4) = ( 4.00000E+00,-5.00000E+00)
Eigenvector( 4) = (-9.06228E-01, 9.37721E-02)(-7.43061E-03, 6.87517E-03)( 3.02073E-02,-3.12593E-03)(-1.45869E-02,-1.40969E-01)
YC
next prev parent reply other threads:[~2012-07-28 13:07 UTC|newest]
Thread overview: 17+ messages / expand[flat|nested] mbox.gz Atom feed top
2012-07-12 19:44 ANN: Ada 2005 Math Extensions 20120712 Simon Wright
2012-07-13 11:12 ` Ken Thomas
2012-07-14 8:13 ` Ada novice
2012-07-14 10:43 ` Simon Wright
2012-07-14 14:18 ` Ada novice
2012-07-27 19:23 ` Ada novice
2012-07-27 21:03 ` Simon Wright
2012-07-28 13:01 ` Ada novice [this message]
2012-07-28 19:10 ` Simon Wright
2012-07-28 19:33 ` Simon Wright
2012-07-29 14:05 ` Ada novice
2012-07-29 14:22 ` Nasser M. Abbasi
2012-07-29 14:34 ` Ada novice
2012-07-29 14:52 ` Nasser M. Abbasi
2012-07-29 15:02 ` Ada novice
2012-07-29 15:31 ` Ada novice
2012-07-30 18:59 ` Ada novice
replies disabled
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox