效率可以对比matlab
| int m;int i, j, k;double r;double *A, *B, *C;m = 1000;A = new double[m*m];B = new double[m*m];C = new double[m*m];//置初值for (i=0; i<m*m; i++) {A[i] = 1;B[i] = 1;C[i] = 0;}int bf; // blocking factorint jj, kk, im;int minj, mink;bf = 48; // 可以修改for (jj=0; jj<m; jj+=bf)for (kk=0; kk<m; kk+=bf)for ( i=0; i<m; ++i) {minj = (jj+bf)<m ? (jj+bf):m;for (j=jj; j<minj; ++j) {r = 0;im = i*m;mink = (kk+bf)<m ? (kk+bf):m;for (k=kk; k<mink; ++k) {r += A[im+k]*B[k*m+j];}C[im+j] += r;}}delete[] A;delete[] B;delete[] C; |
版权与免责声明
1、本站所发布的文章仅供技术交流参考，本站不主张将其做为决策的依据，浏览者可自愿选择采信与否，本站不对因采信这些信息所产生的任何问题负责。
2、本站部分文章来源于网络，其版权为原权利人所有。由于来源之故，有的文章未能获得作者姓名，署“未知”或“佚名”。对于这些文章，有知悉作者姓名的请告知本站，以便及时署名。如果作者要求删除，我们将予以删除。除此之外本站不再承担其它责任。
3、本站部分文章来源于本站原创，本站拥有所有权利。
4、如对本站发布的信息有异议，请联系我们，经本站确认后，将在三个工作日内做出修改或删除处理。
请参阅权责声明！