comp.lang.ada
 help / color / mirror / Atom feed
From: marius63 <amado.alves@gmail.com>
Subject: Re: Eigenvalues to find roots of polynomials
Date: Tue, 26 Jul 2011 07:07:53 -0700 (PDT)
Date: 2011-07-26T07:07:53-07:00	[thread overview]
Message-ID: <ba4f6508-9969-4649-9181-31c37cea11a0@x7g2000vbk.googlegroups.com> (raw)
In-Reply-To: 12rqoo5jd0lqy$.1kgmuqfqym6z2$.dlg@40tude.net

> > --    |   -A(1)/A(0)   -A(2)/A(0)   -A(3)/A(0)   ...   -A(N)/A(0)   |
> > --    |        1            0            0       ...        0       |
> > --    |        0            1            0       ...        0       |
> > --    |       ...          ...           1       ...        0       |
> > --    |_       0            0            0       ...        0      _|
>
> The matrix passed to Eigenvalues (ARM G.3.2) must be Hermitian. The above
> matrix is not Hermitian.
>
> P.S. A matrix M is Hermitian if M (I,J) if the complex conjugate of M
> (J,I)), i.e. Re (M(I,J)) = Re (M(J,I)) and Im (M(I,J)) = -Im (M(J,I))
> (Dmitry A. Kazakov)

I suspected that. I couldn't tell an hermitian matrix from a
hamiltonian gargle blaster so when the program did not raise
Argument_Error I assumed the matrix was kosher enough. It seems not.
(So, Ada Eigenvalues is not a very good way to find roots, is it?)
Thanks.



  reply	other threads:[~2011-07-26 14:07 UTC|newest]

Thread overview: 17+ messages / expand[flat|nested]  mbox.gz  Atom feed  top
2011-07-25 17:36 Eigenvalues to find roots of polynomials marius63
2011-07-25 19:04 ` Dmitry A. Kazakov
2011-07-26 14:07   ` marius63 [this message]
2011-07-25 20:22 ` Simon Wright
2011-07-26  7:31   ` Ada novice
2011-07-26 10:44     ` Simon Wright
2011-07-27  4:54       ` Ada novice
2011-07-26 14:01   ` marius63
2011-07-26 14:04     ` Simon Wright
2011-07-26 14:47       ` marius63
2011-07-26 17:33         ` Simon Wright
2011-07-26 18:04         ` Georg Bauhaus
2011-07-26 19:45           ` marius63
2011-07-26  6:21 ` John B. Matthews
2011-07-26 15:40   ` marius63
2011-07-27  3:21     ` John B. Matthews
2011-07-27  9:33       ` marius63
replies disabled

This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox