repro.c

// cc -o repro repro.c libopenblas.a
// or
// cc -o repro repro.c -L. -lopenblas
// then
// ./repro
#include <stdio.h>
#include "cblas.h"

int main(void)
{
    int m = 2;
    int n = 2;
    int k = 2;
    double alpha = 1.0;
    double lhs[4] = {1.0, 0.0, 0.0, 1.0};
    int lhs_stride = 2;
    double rhs[4] = {1.0, 1.0, 1.0, 1.0};
    int rhs_stride = 2;
    double beta = 0.0;
    double c[4] = {0.0, 0.0, 0.0, 0.0};
    int c_stride = 2;
    cblas_dgemm(
        CblasRowMajor,
        CblasNoTrans,
        CblasNoTrans,
        m,
        n,
        k,
        alpha,
        lhs,
        lhs_stride,
        rhs,
        rhs_stride,
        beta,
        c,
        c_stride);
    for (int i = 0; i < 4; i++)
    {
        printf("%f\n", c[i]);
    }
    return 0;
}