极限测试之Matlab与Forcal编译计算含自变量的简单长表达式

forcal

Matlab 2009a编译了如下简单长函数，多次编译，耗时均大约为3分钟：
. L9 [4 T" k6 k5 t
4 s0 Y5 Y- Z3 J7 rf=@(x)x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
% }- ?8 n6 E; ~* B7 @; Ox+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
8 K: m0 z  d! `' T! g7 h... ...（省略）
% M( z6 g, ~) J9 i. P7 cx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
1;

完整的程序见附件： 多个x加1的Matlab代码.txt (136.83 KB, 下载次数: 0)

2011-7-30 12:06 上传

点击文件名下载附件
下载积分: 体力 -2 点

测试如下代码速度：

1. tic;
   1 n; x\" d; l  v
2. s=0;+ R) [$ K, K2 V& j  g* i* G
   
3. for i=0:1000& E' A: l- J* C5 P% W- y( p
   
4. s=s+f(1);7 A. W4 g: e# f6 R6 k7 o
   
5. end& Y# Y6 |8 s8 F5 _$ ]5 b2 `
   
6. s
   2 t0 w+ H) o4 f% f$ R# p; ]
7. toc
   % C5 w( ]' |7 V: ~, R4 v; R3 H
8. 
   . v, U4 W& v  v\" f
9. s =
   ) P) K2 F% T, |3 r7 Y5 [
10. 
    . g' v\" ]9 R3 @* r5 U- s) d
11.     67100033
    ! r! b) ~; g  n. Q2 T5 \\" D. B
12. : e7 F  u8 m) k, A. T7 t
    
13. Elapsed time is 11.902260 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）编译了如下简单长函数，第一次编译，耗时约0.688秒，以后编译耗时约6.3e-002秒：

4 }  v6 {  C- ]. M0 r( F/ E6 P; zf(x)=
# w' N8 l2 a4 rx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
" h* }1 o5 {/ f( a5 E  w  }: zx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
... ...（省略）
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
1;

* j& T" X5 R# }: Y: q! C$ v完整的程序见附件： 多个x加1的Forcal代码.txt (133.36 KB, 下载次数: 0)

2011-7-30 12:06 上传

点击文件名下载附件
下载积分: 体力 -2 点

9 e/ D. d# _+ y; m4 V6 O
测试如下代码速度：

1. mvar:
   ; v) x; [7 F9 t- O
2. t0=sys::clock(),
   - R, U' D. k8 T0 M
3. s=0, i=0, while{i<=1000, s=s+f(1), i++}, s;
   - H* I3 o, U% z/ \+ C  V
4. [sys::clock()-t0]/1000;

复制代码

结果：
67100033.
0.64  秒
$ v5 f% t) ?& G" }$ b( q6 ^1 M
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

可以看出，Matlab运行大的程序，仅编译就会耗费不少时间。

# s: i' B1 @2 C  Q2 s" c9 ?. d6 o＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

- e1 p: y7 M$ l% N- W! e  q8 K" P为什么Forcal第一次编译耗时较长（约0.688秒），而以后编译耗时较少（约6.3e-002秒）？

Forcal编译器在初次使用时，给编译模块分配的内存较少，此时若所编译的表达式较短，则编译所用内存足够用，在编译时就不用再分配内存；若遇到一个较长的表达式，编译模块就要重新分配内存，直到能编译该表达式为止，故耗费了一些时间；以后，若遇到的表达式都短于以前曾编译的长表达式，就不再重新分配内存，故不再降低编译速度。

* `8 z! p# X: K+ {2 n3 y! m! E目前，需重新初始化Forcal，才释放给编译模块分配的较大内存，但一般不需要这样做。
7 c- C3 Z  @6 U7 y
- U5 m/ W) X5 Z9 @6 A+ k对于脚本，编译效率和运行效率同样重要。

wangchen881202

“拆东墙补西墙”之升级版“拆别人墙补自己墙”

支持一下
+ z8 g6 o+ L! x2 d" W

yunbuhuiku

如果跟导师讲不清楚，那么就把他搞胡涂吧!

6 }4 x( x, @. `* T) E% Q$ m3 R2 U. C来电来函..欢迎洽购...
4 U3 R$ P2 q$ R1 L/ c; S! G

安平池伟

没来得急看，应该不错，先帮你顶

jmedwardliu1205

geilivable