I work on a program in which I'd like to substitute ZHETRD function with MAGMAF_ZHETRD. After substitution I got error message with information that parameter 9 is invalid (length of work array). LAPACK documentation says

- Code: Select all
`LWORK is INTEGER`

The dimension of the array WORK. LWORK >= 1.

For optimum performance LWORK >= N*NB, where NB is the

optimal blocksize.

in magma_zhetrd you have the sane comment on lwork size, but strictly check if lwork>=n*nb

- Code: Select all
`} else if (lwork < nb*n && ! lquery) {`

*info = -9;

}

such behaviour seems not to be good.