针对以下矩阵相乘的C循环语句.docVIP

针对以下矩阵相乘的C循环语句，直接在其源程序基础上做出循环优化（代码外提、强度消弱等） 设矩阵如下： int a[20][20],b[20][20],c[20][20];//width of an integer is 4 // a = b * c; for(i=0;i20;i++) for(j=0;j20;j++) for(k=0;k20;k++) a[i][j] = a[i][j] + b[i][k]*c[k][j]; 首先考虑代码外提。 对于k－循环而言，a[i][j]和b[i]为循环不变式，可以外提至k－循环外（j－循环内），程序变换如下： for(i=0;i20;i++) { for(j=0;j20;j++) { t1 = addr(a[i][j]); t2 = addr(b[i]); for(k=0;k20;k++) { *t1 = *t1 + t2[k]*c[k][j]; } // k-loop } // j-loop } // i-loop 考虑a[i]和b[i]在中间层循环j-loop中仍然是循环不变式，可以进一步外提出j-loop(提到i-loop中)，程序变换为： for(i=0;i20;i++) { t3 = addr(a[i]); t4 = addr(b[i]); for(j=0;j20;j++) { t1 = addr(t3[j]); t2 = t4; for(k=0;k20;k++) { *t1 = *t1 + t2[k]*c[k][j]; } // k-loop } // j-loop } // i-loop 考虑复写传播，可以删除变量t2。程序变换如下： for(i=0;i20;i++) { t3 = addr(a[i]); t4 = addr(b[i]); for(j=0;j20;j++) { t1 = addr(t3[j]); for(k=0;k20;k++) { *t1 = *t1 + t4[k]*c[k][j]; } // k-loop } // j-loop } // i-loop 然后考虑数组下标的情况（设数组元素按行优先排列），即将其变换为指针形式，即地址计算的展开形式。 t3 = addr(a[i]); 其地址表述为：i*20*4+A-0 = i*80+A //A 为数组a的起始地址； //C语言数组下标从0开始。这是一个基本归纳变量i族的归纳表达式 t4 = addr(b[i]); 其地址表述为：i*20*4+B-0 = i*80+B //B 为数组b的起始地址； //C语言数组下标从0开始。这是一个基本归纳变量i族的归纳表达式 t1 = addr(t3[j]); 其地址表述为：j*4+t3-0 = j*4+t3 //C语言数组下标从0开始。这是一个基本归纳变量j族的归纳表达式 t4[k]的地址表述为： k*4+t4-0 = k*4+t4 //C语言数组下标从0开始。这是一个基本归纳变量k族的归纳表达式 c[k][j]的地址表述为： k*20*4+j*4+C-0=k*80+j*4+C //C 为数组c的起始地址； //C语言数组下标从0开始。这是一个基本归纳变量k族的归纳表达式 //因为上面式子中的j*4+C在k-loop中不变 现在进行强度消弱。注意每次的强度消弱分三个步骤完成，即， 引入新的变量s，并将原归纳表达式替换为变量s； 在循环中基本归纳变量增值语句后，添加变量s的增值语句； 在循环外基本归纳变量的初值后（一般在前置块末尾）添加变量s的初值计算语句； t4[k]: 针对上面的归纳表达式k*4+t4引入变量t5 c[k][j]: 针对上面的归纳表达式k*80+j*4+C引入变量t6; addr(t3[j]): 针对上面的归纳表达式j*4+t3引入变量t7; addr(b[i]): 针对上面的归纳表达式i*80+B引入变量t8; addr(a[i]): 针对上面的归纳表达式i*80+A引入变量t9; 程序继续变换为： t8 = B;// addr(b[i])的初值 t9 = A; // addr(a[i])的初值 for(i=0;i20;i++) { t3 = t9; t4 = t8; t7 = t3; // addr(t3[j])的初值； for(j=0;j20;j++) { t1 = t7; t5 = t4; //t4[k]的初值 t6 = j * 4