寻找最高效的矢量化算法

forcal

要不大家练练手？ 两个例子： 1、在1＃的例子：

1. //用C++代码描述为：
2. s=0.0;
3. for(x=0.0;x<=1.0;x=x+0.0011)
4. {
5. for(y=1.0;y<=2.0;y=y+0.0009)
6. {
7. s=s+cos(1-sin(1.2*x^(y/2)+cos(1-sin(1.2*y^(x/2)))));
8. }
9. }

复制代码

2、矩阵乘 matlab代码：

1. tic
2. a = rand(1000,1000);
3. b = rand(1000,1000);
4. k = a*b;
5. k(1:3,5:9)
6. toc

复制代码

Forcal代码：

1. !using["math","sys"];
2. (:a,b,k,t0)=
3. t0=clock(),
4. oo{
5. a=rand[1000,1000], b=rand[1000,1000],
6. k=a*b, //矩阵乘
7. k[1,3:5,9].outm()
8. },
9. [clock()-t0]/1000;

复制代码

无论用什么方法，只要速度快就行。例子1要超过C/C++的速度，例子2超过matlab的速度似乎较难，先超过Forcal的速度即可（这也是C/C++的速度）。 如果各位不愿给出源代码，生成exe文件比较一下即可（只要知道大家可以得到更高的速度就心满意足了，呵呵）。 Forcal程序可下载OpenFC（绿色免安装，2M）演示：http://www.forcal.net/xiazai/forcal9/openfc32w.rar