矩阵乘法的分块优化

普通矩阵乘法一般采用3重循环完成。  void dgemm (int n, double* A, double* B, double* C)  {   for (int i = 0; i < n; ++i)     for (int j = 0; j < n; ++j)     {      double cij = C[i+j*n]; /* cij = C[i][j] */      for(i

>>阅读原文<<