comp.lang.ada
 help / color / mirror / Atom feed
From: Simon Wright <simon@pushface.org>
Subject: Re: ANN: Ada 2005 Math Extensions 20120712
Date: Sat, 28 Jul 2012 20:33:32 +0100
Date: 2012-07-28T20:33:32+01:00	[thread overview]
Message-ID: <m2ehnvpto3.fsf@nidhoggr.home> (raw)
In-Reply-To: e1d5d31c-f873-4910-8bf3-6700a17007d1@googlegroups.com

Ada novice <shai.lesh@gmx.com> writes:

> When running "demo_extensions", I saw some differences in the outputs
> (in the last example for instance).

The last example comes from
http://www.nag.co.uk/lapack-ex/node122.html. The results
(http://www.nag.co.uk/lapack-ex/examples/results/zggev-ex.r) look a lot
more like the Debian output you list than the ACML output.

FWIW the Mac OS X 10.7.4 output is

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.44663E-02)( 7.44662E-02, 1.53495E-01)(-1.53495E-01, 7.44663E-02)
Eigenvalue( 2) = ( 2.00000E+00,-5.00000E+00)
Eigenvector( 2) = (-6.29583E-01,-3.70417E-01)(-4.14824E-03, 4.65175E-04)(-3.95739E-02,-2.32833E-02)( 2.32834E-02,-3.95739E-02)
Eigenvalue( 3) = ( 3.00000E+00,-1.00001E+00)
Eigenvector( 3) = ( 9.77535E-01, 2.24653E-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.37645E-02)( 7.43045E-03,-6.87512E-03)(-3.02078E-02, 3.12552E-03)( 1.45859E-02, 1.40970E-01)

which has similar absolute values but the signs of all except (3) are
reversed; and are the same as the outputs of the NAG test program/data.

I'm not a numerical analyst, so I have no idea whether these differences
represent bugs or features.



  parent reply	other threads:[~2012-08-01  3:36 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
2012-07-28 19:10             ` Simon Wright
2012-07-28 19:33             ` Simon Wright [this message]
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