Community
cancel
Showing results for 
Search instead for 
Did you mean: 
Highlighted
Beginner
72 Views

Speed of DGGHRD and DGGHD3

I noticed something about the routines DGGHRD and DGGHD3. They both reduce a pencil to Hessenberg triangular form for use in the QZ algorithm, with DGGHD3 being the blocked version.

When linked with the reference LAPACK or with OpenBLAS, DGGHD3 is faster than DHHGRD. However, when the test program is linked with MKL, DGGHD3 loses its speed advantage. I also don't see any multicore usage in that routine when linked with MKL. Is it possible that the blocked version is not (fully) implemented in MKL?

Labels (1)
0 Kudos
1 Reply
Highlighted
Beginner
62 Views

Update, it looks like the blocksize returned by ilaenv is always 1 for DGGHD3. Changing that seems to fix the performance issue.

0 Kudos