Using MAGMA GPU library

[mgates3]

lwork is the length of the work array that is passed to MAGMA's routine. Per the documentation, it should be:

Code: Select all

     -      If JOBZ = MagmaNoVec and N > 1, LWORK >= N + N*NB.
     -      If JOBZ = MagmaVec   and N > 1, LWORK >= max( N + N*NB, 2*N + N**2 ).
            NB can be obtained through magma_get_zhetrd_nb(N).

I guess whoever wrote abi_xhegv.f90 assumed that NB = 33, but that is a bad assumption. It should be queried as described in the documentation above.

-mark

[Sorry for the delay in posting. A misconfiguration left messages in the moderation queue.]