Re: Eigenvalues Okay, Eigenvectors Fail

*To*: mathgroup at smc.vnet.net*Subject*: [mg88742] Re: Eigenvalues Okay, Eigenvectors Fail*From*: jfuite at phys.ualberta.ca*Date*: Thu, 15 May 2008 06:52:34 -0400 (EDT)*References*: <fupm4m$s4r$1@smc.vnet.net>

Okay, a week later, I can answer my own question for the benefit of someone else someday, perhaps. To recap, I had a large (n = 733), square (n x n), real matrix, mSub, and asked why, if I could find eigenvalues in short order, my computer could not find the eigenvectors? Such that, A. Timing@Eigenvalues[mSub] = {2 Second, l}, where l is a list of eigenvalues of length n, but, B. Eigenvectors[mSub] does not complete. Anyway, I still do not completely understand what is happening, but I get the sense that Mathematica was simply not converging to an answer quickly with the default settings for the function Eigenvectors. Reading Trott, "Mathematica Guidebook for Numerics", chapter 1.4, and later Mathematica Help, "Advanced Documentation - Linear Algebra", I learned of some obscure options for the function Eigenvectors, and set up the following function which runs quickly and provides the kind of answers I need: Eigensystem[mSub, 4, Method->{"Arnoldi", Criteria->"BothEnds", MaxIterations->10^3, Tolerance->10^-3}]