blas_misunderstanding_amp.cpp

#include<stdio.h>

int main(void){

    double a[4][4];
    double b[4][4];
    double c[4][4];

    int i, j;
    int  n = 4;
    double alpha = 1.0, beta = 0.0;

    //aとbの初期化

    DGEMM('T', 'T', &n, &n, &n, &alpha, a, &n, b, &n, &beta, c, &n);

    for(i = 0; i < n; i++){
        for(j = 0; j < n; j++){
            printf("c[%d][%d] =  %5.2f; ", i, j, c[j][i]);
        }
        printf("\n");
    }

    return 0;
}