1000*1000的矩阵A、B、C，C=A*B，用C/C++实现 参考：blocking factor，http://home.iae.nl/users/mhx/mm.c 效率可以对比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 factor int 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;

这个程序比matlab要慢1倍，我曾到网上搜索到一个fortran程序(分块dgemm)， 其速度比matlab慢0.5倍左右，不知道有谁能写一个能与matlab媲美的程序？

MATLAB 语言是当今国际上科学界 (尤其是自动控制领域) 最具影响力、也是最有活力的软件。它起源于矩阵运算，并已经发展成一种高度集成的计算机语言。它提供了强大的科学运算、灵活的程序设计流程、高质量的图形可视化与界面设计、便捷的与其他程序和语言接口的功能。MATLAB 语言在各国高校与研究单位起着重大的作用。 MATLAB 语言由美国 The MathWorks 开发，2003 年推出了其全新的 MATLAB 6.5.1 正式版。目前最新版本 Release 14 (MATLAB 7.0) 的测试版，2004 年 7 月 (?) 将正式推出。

呵呵 不好意思 我只知道一些关于matlab的简介``

[QUOTE]这个程序比matlab要慢1倍，我曾到网上搜索到一个fortran程序(分块dgemm)， 其速度比matlab慢0.5倍左右，不知道有谁能写一个能与matlab媲美的程序？[/QUOTE]

我这点水平的确很难搞定 我也等待高手的佳音`````

嗯

这个很难，连fortran都难以比拟。要知道fortran在矩阵运算中的能力是其他语言没法比的。