MKL 进行矩阵向量运算

最新推荐文章于 2024-08-18 21:03:01 发布

AresGod

最新推荐文章于 2024-08-18 21:03:01 发布

阅读量3.3k

点赞数

分类专栏： c/C++编程参考

本文链接：https://blog.csdn.net/AresGod/article/details/51476982

版权

这段代码展示了如何使用Intel MKL库的dgemm函数通过CBLAS接口进行矩阵乘法运算，比较了MKL库的性能与传统三重循环实现的性能差异。初始化数据并进行多次运算，计算每种方法的平均运行时间，从而突出MKL库在矩阵运算中的优势。

摘要由CSDN通过智能技术生成

/* C source code is found in dgemm_example.c */

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

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

int main()
{
float *A, *B, *C;
int m, n, k, i, j;
int l,r;
double s_initial,s_elapsed;
int LOOP_COUNT =100;
float alpha, beta;
double sum;

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

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