VS2019配置MKL教程(Windows)

#include <stdio.h>#include <stdlib.h>
#include "mkl.h"

#define min(x,y) (((x) < (y)) ? (x) : (y))

intmain()
{
    double* A, * B, *C;
    intm, n, k, i, j;
    doublealpha, beta;

    printf("This example computes real matrix C=alpha*A*B+beta*C using 
"
        "Intel(R) MKL function dgemm, where A, B, and  C are matrices and 
"
        "alpha and beta are double precision scalars

");

    m = 2000, k = 200, n = 1000;
    printf("Initializing data for matrix multiplication C=A*B for matrix 
"
        "A(%ix%i) and matrix B(%ix%i)

", m, k, k, n);
    alpha = 1.0; beta = 0.0;

    printf("Allocating memory for matrices aligned on 64-byte boundary for better 
"
        "performance 

");
    A = (double*)mkl_malloc(m * k * sizeof(double), 64);
    B = (double*)mkl_malloc(k * n * sizeof(double), 64);
    C = (double*)mkl_malloc(m * n * sizeof(double), 64);
    if (A == NULL || B == NULL || C ==NULL) {
        printf("ERROR: Can't allocate memory for matrices. Aborting... 

");
        mkl_free(A);
        mkl_free(B);
        mkl_free(C);
        return 1;
    }

    printf("Intializing matrix data 

");
    for (i = 0; i < (m * k); i++) {
        A[i] = (double)(i + 1);
    }

    for (i = 0; i < (k * n); i++) {
        B[i] = (double)(-i - 1);
    }

    for (i = 0; i < (m * n); i++) {
        C[i] = 0.0;
    }

    printf("Computing matrix product using Intel(R) MKL dgemm function via CBLAS interface 

");
    cblas_dgemm(CblasRowMajor, CblasNoTrans, CblasNoTrans,
        m, n, k, alpha, A, k, B, n, beta, C, n);
    printf("Computations completed.

");

    printf("Top left corner of matrix A: 
");
    for (i = 0; i < min(m, 6); i++) {
        for (j = 0; j < min(k, 6); j++) {
            printf("%12.0f", A[j + i *k]);
        }
        printf("");
    }

    printf("Top left corner of matrix B: 
");
    for (i = 0; i < min(k, 6); i++) {
        for (j = 0; j < min(n, 6); j++) {
            printf("%12.0f", B[j + i *n]);
        }
        printf("");
    }

    printf("Top left corner of matrix C: 
");
    for (i = 0; i < min(m, 6); i++) {
        for (j = 0; j < min(n, 6); j++) {
            printf("%12.5G", C[j + i *n]);
        }
        printf("");
    }

    printf("Deallocating memory 

");
    mkl_free(A);
    mkl_free(B);
    mkl_free(C);

    printf("Example completed. 

");

    system("PAUSE");
    return 0;
}