在非平方Q中使用LAPACK的DORMQR

#include <stdio.h>
#include <stdlib.h>
#include <lapacke.h>

void multiply_q_c(int m, int n, int k, double *q, int ldq, double *c, int ldc) {
    lapack_int info;
    char side = 'L'; // 左乘
    char trans = 'N'; // 不转置
    double *work = (double *)malloc(k * sizeof(double));
    lapack_int lwork = k;

    info = LAPACKE_dormqr(LAPACK_COL_MAJOR, side, trans, m, n, k, q, ldq, work, c, ldc);

    if (info != 0) {
        printf("Error in DORMQR: %d\n", info);
    } else {
        printf("Matrix multiplication successful.\n");
    }

    free(work);
}

int main() {
    int m = 3, n = 2, k = 3;
    double q[k * k] = {
        1, 0, 0,
        0, 1, 0,
        0, 0, 1
    };
    int ldq = k;
    double c[n * m] = {
        1, 2, 3,
        4, 5, 6
    };
    int ldc = m;

    multiply_q_c(m, n, k, q, ldq, c, ldc);

    return 0;
}