Hello,
I am trying to apply Q obtained from LAPACKE_zgelqf to the rectangular matrix from the right using LAPACKE_zunmlq.
#include "mkl.h"
int main() {
int m = 637;
int rows1 = m;
int columns1 = 3 * m;
int lda1 = rows1;
MKL_Complex16* matrix1 = new MKL_Complex16[rows1 * columns1];
for (int i = 0; i < rows1 * columns1; ++i) {
matrix1[i] = { (double)i, (double)i };
}
int rows2 = 3 * m + 1;
int columns2 = 3 * m;
int lda2 = rows2;
MKL_Complex16* matrix2 = new MKL_Complex16[rows2 * columns2];
for (int i = 0; i < rows2 * columns2; ++i) {
matrix2[i] = { (double)i, (double)i };
}
int refl = rows1;
MKL_Complex16* tau = new MKL_Complex16[refl];
int info = 0;
info = LAPACKE_zgelqf(LAPACK_COL_MAJOR, rows1, columns1, matrix1, lda1, tau);
info = LAPACKE_zunmlq(LAPACK_COL_MAJOR, 'R', 'C', rows2, columns2, refl, matrix1, lda1, tau, matrix2, lda2);
delete[] tau;
delete[] matrix1;
delete[] matrix2;
}
I am sure that all parameters have appropriate values but call to the LAPACKE_zunmlq cause access violation.
I run this code using VS 2019 with oneAPI 2022.0 on the machine with Windows 10.
Am I missing something obvious or is there a bug in LAPACKE_zunmlq
