From mboxrd@z Thu Jan 1 00:00:00 1970 X-Spam-Checker-Version: SpamAssassin 3.4.4 (2020-01-24) on polar.synack.me X-Spam-Level: X-Spam-Status: No, score=-1.9 required=5.0 tests=BAYES_00,FREEMAIL_FROM autolearn=ham autolearn_force=no version=3.4.4 X-Google-Thread: 103376,35f6cee6f665d64b X-Google-NewGroupId: yes X-Google-Attributes: gida07f3367d7,domainid0,public,usenet X-Google-Language: ENGLISH,ASCII-7-bit Received: by 10.66.85.71 with SMTP id f7mr916727paz.39.1343480860269; Sat, 28 Jul 2012 06:07:40 -0700 (PDT) Path: c10ni45307pbw.0!nntp.google.com!news1.google.com!postnews.google.com!glegroupsg2000goo.googlegroups.com!not-for-mail From: Ada novice Newsgroups: comp.lang.ada Subject: Re: ANN: Ada 2005 Math Extensions 20120712 Date: Sat, 28 Jul 2012 06:01:26 -0700 (PDT) Organization: http://groups.google.com Message-ID: References: <167ecede-3588-45f0-904b-06a8be4cfec7@googlegroups.com> <2ba72d4e-6e88-4900-8232-c075242dec1f@googlegroups.com> <65f76f80-a1e5-4c60-8002-45becc2a1198@googlegroups.com> NNTP-Posting-Host: 193.11.21.204 Mime-Version: 1.0 X-Trace: posting.google.com 1343480860 8199 127.0.0.1 (28 Jul 2012 13:07:40 GMT) X-Complaints-To: groups-abuse@google.com NNTP-Posting-Date: Sat, 28 Jul 2012 13:07:40 +0000 (UTC) In-Reply-To: Complaints-To: groups-abuse@google.com Injection-Info: glegroupsg2000goo.googlegroups.com; posting-host=193.11.21.204; posting-account=Rr9I-QoAAACS-nOzpA-mGxtAlZ46Nb6I User-Agent: G2/1.0 Content-Type: text/plain; charset=ISO-8859-1 Date: 2012-07-28T06:01:26-07:00 List-Id: 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