C中实现矩阵乘法的一种高效的方法
更新时间：2013年05月04日 17:14:03 作者：
本篇文章介绍了，在C中实现矩阵乘法的一种高效的方法。需要的朋友参考下
如何计算矩阵乘法，这个大家都知道。通常情况下，我们都是用以下代码实现的：
复制代码 代码如下:
for(i=0;i<n;++i)
for(j=0;j<n;++j){
sum=0;
for(k=0;k<n;++k)
sum+=A[i][k]*B[k][j];
C[i][j]+=sum;
}
但是考虑了高速缓存的问题后，其实有一种更好的实现方式：
复制代码 代码如下:
for(i=0;i<n;++i)
for(k=0;k<n;++k){
r=A[i][k];
for(j=0;j<n;++j)
C[i][j]+=r*B[k][j];
}
细看一番就会发现这两种实现语义是等价的，但是后者的实际运行效率却比前者高。
那为什么会如此呢？
那是因为CPU读数据时，并不是直接访问内存，而是先查看缓存中是否有数据，有的话直接从缓存读取。而从缓存读取数据比从内存读数据快很多。
当数据不在缓存中时，CPU会将包含数据在内的一个数据块读到缓存，如果程序具有良好空间局部性，那么第一次cache miss后，之后的几次数据访问就可以直接在缓存中完成。除了空间局部性（程序倾向于引用与当前数据邻近的数据）之外，还有时间局部性（程序倾向于引用最近被引用过的数据）。
回到矩阵乘法。（我们只考虑内循环）
前者对矩阵A，有良好的空间局部性，假设一次能缓存四个元素，则每次迭代对于A只有0.25次miss，但是对于B，则不然，因此B是按列访问的，每次访问都会miss，因此每次迭代总的miss数是1.25。
后者对于矩阵C和矩阵B都有良好的局部性，每次迭代都只有0.25词miss，因此总的miss数是0.5。后者每次迭代多了一次存储（对C[i][j]写入），但是即便如此，后者的运行效率也比前者高。
总而言之，要想程序跑得快，就要在程序中多利用局部性，让缓存hold住你的数据，减少访存次数。要知道CPU可以在3个时钟周期内访问到L1 cache，10个时钟周期左右的时间访问到L2 cache。访问内存却要上百个时钟周期，孰快孰慢，很清楚了吧？
相关文章
这篇文章主要介绍了使用C++绘制GDI位图的基本编写实例,一般来说适用于Windwos下的C++的GUI编程,需要的朋友可以参考下2015-12-12
以下实例是介绍了二分法求多项式在-10 10间值的实现代码。需要的朋友参考下2013-05-05
这篇文章主要介绍了C语言之格式化屏幕输出的相关资料,需要的朋友可以参考下，小编觉得这篇文章写的还不错，希望能够给你带来帮助2021-11-11
本篇文章是对在Linux下使用c语言获取进程退出值的方法进行了详细的分析介绍，需要的朋友参考下2013-05-05
这篇文章主要介绍了C语言使用openSSL库AES模块实现加密功能,详细分析了C语言加密的相关概念、原理及AES模块加密具体实现技巧,需要的朋友可以参考下2017-05-05
对于字符,在计算机内部都是用数字（字符编码）来表示的,而字符串是“字符连续排列”的一种表现,这篇文章主要给大家介绍了关于C语言进阶教程之字符串&内存函数的相关资料,需要的朋友可以参考下2021-09-09
我们使用QT进行界面开发时,可能会遇到需要将窗口置顶的情况,下面这篇文章主要给大家介绍了关于QT窗口/控件置顶方法的相关资料,文中通过代码介绍的非常详细,需要的朋友可以参考下2024-01-01
这篇文章主要介绍了C语言实现在数组A上有序合并数组B的方法,包含了数组操作的完整实现过程以及相应的代码分析与改进,具有不错的借鉴价值,需要的朋友可以参考下2014-09-09
c++三大函数指的是拷贝构造、拷贝赋值、析构函数,下面这篇文章主要给大家介绍了关于C++类中三大函数(构造、析构和拷贝)的相关资料,文中通过实例代码介绍的非常详细,需要的朋友可以参考下2023-03-03
这篇文章主要介绍了C++异步数据交换实现方法，异步数据交换，除了阻塞函数 send() 和 recv() 之外，Boost.MPI 还支持与成员函数 isend() 和 irecv() 的异步数据交换2022-11-11
最新评论