mm parameter on DHSEIN

Post here if you have a question about LAPACK or ScaLAPACK algorithm or data format

mm parameter on DHSEIN

Postby afani00 » Wed Nov 07, 2012 1:14 pm

Hi, I'm a new user.

I have a question about the subroutine dhsein.

I have a 10x10 real matrix in Hessenberg form. I need to compute all the eigenvectors, so select is a true array of 10 elements. So I choose mm =10.

Sometimes the routine crashes because m exeed 10 (it is 11). I don't understand how this is possible.

Andrea
afani00
 
Posts: 1
Joined: Wed Nov 07, 2012 1:03 pm

Re: mm parameter on DHSEIN

Postby CyLith » Wed Nov 07, 2012 4:22 pm

Is it because you have complex eigenvectors? The comment header says they are stored in two columns, so potentially you would need mm = 2*m.
CyLith
 
Posts: 35
Joined: Sun Feb 08, 2009 7:23 am
Location: Stanford, CA


Return to Algorithm / Data

Who is online

Users browsing this forum: Yahoo [Bot] and 1 guest