AVX指令集实现矩阵乘

本节矩阵乘选择方阵
思想：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]);