矩阵乘法优化之分块矩阵

当矩阵尺寸过大时,数据的大小将超过缓存的大小,这是容易出现满不命中现象。

将矩阵进行分块可以解决这个问题,以下是完整的矩阵乘法代码:

vord brck(array A, array B, array C, int n, int bsize) 
{
    int r, c, k, kk, cc;
    double sum;
    int en = bsize * (n/bsize); /* Amount that frts evenly into blocks */
  
    for (r = 0; r < n; r++)
        for (c = 0; c < n; c++)
            C[r][c] = 0.0;

    for (kk = 0; kk < en; kk += bsize) { 
        for (cc = 0; cc < en; cc += bsize) {
            for (r = 0; r < n; r++) {
                for (c = cc; c < cc + bsize; c++) {
                    sum = C[r][c];
                    for (k = kk; k < kk + bsize; k++) {
                         sum += A[r][k]*B[k][c];
                    }
                    C[r][c] = sum;
                }
            }
        }

  分析思路:

  1. 矩阵分块前后的乘法计算总数恒定不变,分块前是n^3 。

  2. 现将矩阵按mxm进行分块,整个矩阵被分成n^2/m^2 个子矩阵,乘法计算总是 (n^2/m^2) x n x m^2  。

  3. 由2可知,分块后,矩阵以mxm为单位进行乘法运算,它被嵌套在三层循环内。

  ps: 本算法分块后的并不是正方形矩阵,而是在n行或者n列上的矩阵乘积和。优点: 每次都在相邻位置上进行读写,提高了访问性能。

 

转载于:https://www.cnblogs.com/wlzy/p/8997483.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值