Open discussion for MAGMA library (Matrix Algebra on GPU and Multicore Architectures)

mgates3
 Posts: 916
 Joined: Fri Jan 06, 2012 2:13 pm
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.]