展开循环来加速计算

最新推荐文章于 2024-01-08 20:00:00 发布

转载最新推荐文章于 2024-01-08 20:00:00 发布 · 557 阅读

1 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/kljfdsa/p/9027991.html

对于循环次数较为有限的for循环，可通过手工展开该循环达到消除此循环的目的，达到加速运行的目的，下面给出一个代码段实例。

1、for循环体

1 // Q为常数
2 for ( k = 0; k < Q; k++ ) {
3     for ( int j = 0; j < Q; j++ ) {
4         m[k] += fp[z][y][x][j] * mat_M[k][j];
5     }
6 }

此段代码由两重循环构成，Q为常数，尝试将内层循环展开来消除内层循环。

2、内存循环展开

1 for ( k = 0; k < Q; k++ ) {
2     m[k] = fp[z][y][x][0] * mat_M[k][0] + fp[z][y][x][1] * mat_M[k][1] + fp[z][y][x][2] * mat_M[k][2] + fp[z][y][x][3] * mat_M[k][3] + fp[z][y][x][4] * mat_M[k][4]
3            + fp[z][y][x][5] * mat_M[k][5] + fp[z][y][x][6] * mat_M[k][6] + fp[z][y][x][7] * mat_M[k][7] + fp[z][y][x][8] * mat_M[k][8] + fp[z][y][x][9] * mat_M[k][9]
4            + fp[z][y][x][10] * mat_M[k][10] + fp[z][y][x][11] * mat_M[k][11] + fp[z][y][x][12] * mat_M[k][12] + fp[z][y][x][13] * mat_M[k][13] + fp[z][y][x][14] * mat_M[k][14]
5            + fp[z][y][x][15] * mat_M[k][15] + fp[z][y][x][16] * mat_M[k][16] + fp[z][y][x][17] * mat_M[k][17] + fp[z][y][x][18] * mat_M[k][18];
6 }