Bill Allombert on Wed, 13 Jan 2010 01:30:16 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: error finding eigenspaces 

On Tue, Jan 12, 2010 at 04:18:41PM -0600, Eliot Brenner wrote:
> When trying to find eigenspaces of a matrix with mateigen() command, I get
> the error message
> 
>   ***   missing eigenspace. Compute the matrix to higher accuracy, then
>   ***   restart eigen at the current precision.
> 
> can you tell me what this means and how to address it?

It means exactly what it is written. The algorithm used to compute eigenspace
is numerically unstable and sometimes will fail with this error which suggests
to increase the precision (use \p).

Cheers,
Bill.