LAPACK/ScaLAPACK Development

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

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.