本节矩阵乘选择方阵
思想:c语言默认按行优先存储,矩阵a * b,a的行连续,可以连续访存,大大提高效率;但是b要按列取数,所以去b的列向量浪费时间,解决办法是:将b转置存储,这样b就可以按行进行连续访问。但是要牺牲空间去转存,还浪费了转存的计算时间。但实际仍然有加速效果。
另外还有很多算法可以提升效率。
#include <stdio.h>
#include <time.h>
#include <stdlib.h>
#include <x86intrin.h>
const int M = 1000;
//串行矩阵乘
void matrix(float **a, float **b, float **c){
int i, j, k;
for(i = 0; i < M; i++){
for(k = 0; k < M; k++){
for(j = 0; j < M; j = j+1){
c[i][j] += a[i][k] * b[k][j];
//c[i][j+1] += a[i][k] * b[k][j+1];
//c[i][j+2] += a[i][k] * b[k][j+2];
//c[i][j+3] += a[i][k] * b[k][j+3];
}
}
}
}
//两层循环展开
void matrix_loop_two(float **a, float **b, float **c){
int i, j, k;
for(i = 0; i < M; i++){
for(k = 0; k < M; k++){
for(j = 0; j < M; j = j+2){
c[i][j] += a[i][k] * b[k][j];
c[i][j+1] += a[i][k] * b[k][j+1];
//c[i][j+2] += a[i][k] * b[k][j+2];
//c[i][j+3] += a[i][k] * b[k][j+3];
}
}
}
}
//数组赋值
void value(float **a){
int i, j, t = 1.0;
for(i = 0; i < M; i++){
for(j = 0; j < M; j++){
a[i][j] = t;
}
t++;
}
}
//打印数组
void print(float **a){
int i, j;
for(i = 0; i < M; i++){
for(j = 0; j < M; j++){
printf("%.2f ", a[i][j]);