Open discussion for MAGMA library (Matrix Algebra on GPU and Multicore Architectures)
-
mgates3
- Posts: 918
- 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.]