       Behavior of Eigenvalues and Eigensystem

• To: mathgroup at smc.vnet.net
• Subject: [mg115072] Behavior of Eigenvalues and Eigensystem
• From: "Eric Michielssen" <emichiel at eecs.umich.edu>
• Date: Thu, 30 Dec 2010 04:11:10 -0500 (EST)

```Hi,

I am running into some trouble using Eigenvalues and Eigensystem.

I run

Eigensystem[mat1,mat2,3, Method -> Arnoldi]

where mat1 and mat2 are both sparse, Hermitian, positive definite matrices.
Mathematica confirms this: when I run HermitianMatrixQ[mat1] or
PositiveDefiniteMatrixQ[mat1] both return True, same for mat2. I also
explicitly compute these matrices' eigenvalues, and they all are positive.

Yet, when I run the Eigensystem command then one of two things happens:

1. Mathematica complains that the matrices are not Hermitian. OK, they
sometimes are ever so slightly non Hermitian, say by having an imaginary
component 10^(-21) on the main diagonal. I can eliminate that by
Chopping/symmetrizing, but that comes at the cost of additional memory, and
mat1 and mat2 can be huge.

2. Even when I do so, Mathematica complains that the second matrix is not
positive definite:

Eigensystem::chnpdef: The second matrix
SparseArray[Automatic,{4,4},0.\[VeryThinSpace]+0.
I,{1,{{0,3,6,9,12},{{1},{2},{3},{1},{2},{4},{1},{3},{4},{2},{3},{4}}},{4.333
33\[VeryThinSpace]+0. I,-0.666667+0. I,-2.81637+0.719138 I,-0.666667+0.
I,4.33333\[VeryThinSpace]+0. I,-2.81637+0.719138 I,-2.81637-0.719138
I,4.33333\[VeryThinSpace]+0. I,-0.666667+0. I,-2.81637-0.719138
I,-0.666667+0. I,4.33333\[VeryThinSpace]+0. I}}]
in the first argument seems not to be positive definite, which is required
for the Arnoldi method. >>

Stuck...

Note: I use Method -> Arnoldi as Mathematica suggests I do so as my matrices
are highly sparse, and huge. I cannot find any documentation on Method ->
Arnoldi, though. Are there other dials that can be set when using this
option? Are there other options available for Eigenvalues/system?

Thanks,
Eric Michielssen

```

• Prev by Date: Re: Font problem under Linux
• Next by Date: Re: Mathematica daily WTF
• Previous by thread: Re: Wolfram free input shortcut is "cntrl ="
• Next by thread: Re: Behavior of Eigenvalues and Eigensystem