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

forcal

Matlab 2009a编译了如下简单长函数，多次编译，耗时均大约为3分钟：

f=@(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+...
+ B8 N/ r, {9 m" P/ q* q2 Kx+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+...
... ...（省略）
* S$ D! r. \0 _" g: M3 a5 j2 \1 Ix+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;
) R1 O/ s1 D1 S, M& s
* W: P0 z. U! `; I" Q' v7 z" M完整的程序见附件： 多个x加1的Matlab代码.txt (136.83 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

测试如下代码速度：

1. tic;
   ( D: j/ x* ^4 ]- f2 q/ t
2. s=0;
   1 q/ M1 y6 G  l) I; W
3. for i=0:1000* D& [9 ]9 [9 Z5 g& p: u' J
   
4. s=s+f(1);( A5 S; t' B5 V
   
5. end0 I7 T' x- w, J* H7 N6 G' S& J1 E
   
6. s
   0 Q0 M8 A7 @) F- ?4 a
7. toc
   . N5 N1 r1 z& f, _1 ^
8. - q6 m\" Q; E' D/ k$ Y( z. O8 ?
   
9. s =7 U4 |1 Y. Q% o: ~\" u1 }/ l
   
10. 
    ) ?, }2 R* y* W1 @+ V$ G4 v\" w
11.     67100033
    ' h, ]) x) A4 a. E
12. 
    0 }, B9 d6 Q\" j  P
13. Elapsed time is 11.902260 seconds.

复制代码

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

Forcal 9（OpenFC演示）编译了如下简单长函数，第一次编译，耗时约0.688秒，以后编译耗时约6.3e-002秒：
9 C% Y) U. k6 n
5 ?; l7 o7 C5 @2 B0 K% q- v9 q& Tf(x)=
" n/ Z  C2 A% h$ ^+ gx+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+
... ...（省略）
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+
9 z0 U" D; |. j1;

' O, x7 C' @8 B( S完整的程序见附件： 多个x加1的Forcal代码.txt (133.36 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

, G/ h" R- i# |' C" l- `测试如下代码速度：

1. mvar:
   $ Y: G. j! {6 `) I# b8 `5 z\" l
2. t0=sys::clock(),! R% z9 V! N5 Y8 T: I$ V\" l
   
3. s=0, i=0, while{i<=1000, s=s+f(1), i++}, s;+ @1 R8 M  C: j\" {, d6 b) l  t
   
4. [sys::clock()-t0]/1000;

复制代码

结果：
67100033.
* c( r: q9 T: U& M  k) Q0.64  秒

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
. l8 X/ `# t( Z
! r! K( _0 v) L' b可以看出，Matlab运行大的程序，仅编译就会耗费不少时间。
+ g/ [/ @5 A" ]% I2 p& [& ^$ R
5 t/ c8 f: f2 W8 G" K＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
4 `* t% e0 Q+ K  a7 T
为什么Forcal第一次编译耗时较长（约0.688秒），而以后编译耗时较少（约6.3e-002秒）？

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

目前，需重新初始化Forcal，才释放给编译模块分配的较大内存，但一般不需要这样做。
; z9 S" R9 v% O% M/ i4 m
7 C! y% R! @" b, h对于脚本，编译效率和运行效率同样重要。

jmedwardliu1205

geilivable

安平池伟

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

yunbuhuiku

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

来电来函..欢迎洽购...

wangchen881202

“拆东墙补西墙”之升级版“拆别人墙补自己墙”
0 H; ~9 n- E( c' U) u& T
支持一下
' w' ~# M6 H: x- k$ m