When I call (mkl from fortran)
NB = ILAENV( 1, 'DGETRF', ' ', 4000, 4000, -1, -1 )
it always gives NB=1. For other methods it seems to work ok (DGEQRF gives NB=128).
I don't want to use mkl's getrf. I'm running Intel(R) Math Kernel Library 11.0 Update 5 for Linux* OS.
check with the latest version on MKL 11.2
here what I see on my side. win64.
dsytrd blocksize Output DATA: 16
dormtr blocksize Output DATA: 128
DGETRF blocksize Output DATA: 164
!Intel(R) Math Kernel Library Version 11.2.0 Product Build 20140723 for Intel(R) 64 architecture applications