Using MAGMA GPU library

Open discussion for MAGMA library (Matrix Algebra on GPU and Multicore Architectures)
Post Reply
mgates3
Posts: 916
Joined: Fri Jan 06, 2012 2:13 pm

Re: Using MAGMA GPU library

Post by mgates3 » Tue Apr 23, 2019 4:51 pm

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.]

Post Reply