bad work array sizes in QRP about matrix-toolkits-java HOT 1 CLOSED

/*
  gcc -O3 mkl-bug.c common.c -o mkl-bug -I. -L. -lblas -lgfortran -lm

  Then choosing the backend is a runtime choice:

  LD_LIBRARY_PATH=. ./mkl-bug
  LD_LIBRARY_PATH=/usr/lib/atlas-base:. ./mkl-bug
  LD_LIBRARY_PATH=/usr/lib/openblas-base:. ./mkl-bug
  LD_LIBRARY_PATH=/opt/intel/mkl/lib/intel64:. ./mkl-bug
*/

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <cblas.h>
#include "common.h"

int main() {
    int m = 5;
    int n = 5;
    double a[25] = {
        1.0, 0.0, 0.0, 0.0, 0.0,
        0.0, 1.0, 0.0, 0.0, 0.0,
        0.0, 0.0, 1.0, 0.0, 0.0,
        0.0, 0.0, 0.0, 1.0, 0.0,
        0.0, 0.0, 0.0, 0.0, 1.0
    };
    int lda = 5;
    int jpvt[5] = { 0, 0, 0, 0, 0 };
    double tau[5] = { 0.0, 0.0, 0.0, 0.0, 0.0 };
    double *work;
    int lwork = -1;
    int info = 0;
    double optimal = 0.0;
    int i, j;

    dgeqp3( &m, &n, a, &lda, jpvt, tau, &optimal, &lwork, &info );

    lwork = (int)optimal;
    printf("optimal = %i\n", lwork);

    work = malloc(lwork * sizeof(*work));

    dgeqp3( &m, &n, a, &lda, jpvt, tau, work, &lwork, &info );

    printf("result = %i\n", info);
    for(i = 0; i < 5; i++) {
        for(j = 0; j < 5; j++) {
            printf("%f ", a[j * 5 + i]);
        }
        printf("\n");
    }

    free(work);

    return(0);
}