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

forcal

Matlab 2009a编译了如下简单长函数，多次编译，耗时均大约为3分钟：
& s" G$ v( d0 f$ F# G7 N/ k6 w
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+...
x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
- F% S" D5 z/ r( i5 @3 g, tx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+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 W2 F; o; Z" Yx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+...
  O1 X3 r! v( I& f! y1;

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

2011-7-30 12:06 上传

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

测试如下代码速度：

1. tic;
   * \/ g& ?% U* ~' n; ]/ Q, O3 t* W
2. s=0;+ F# y; M6 s( O/ }
   
3. for i=0:1000
   ( `% [$ c$ |+ b5 j# f\" Z
4. s=s+f(1);
   ( ]/ X8 r1 G, a/ @$ P% V7 ?
5. end
   ( J3 [$ ]7 I/ \* Y2 F
6. s
   9 L9 X* F$ {6 o; U
7. toc
   6 l9 v( m' ^8 ~5 n0 J\" |2 z
8. 
   \" ^2 Z- |\" w* C% K  t5 N
9. s =
   ; Z$ U  L% Z\" o1 r: V
10. 1 x/ Y- F7 j8 t\" e, f8 x
    
11.     67100033
    ! w5 N  o1 V# e, P- Y9 c! }* r
12. 2 G5 l6 Z& l/ H% T3 f1 a& q# [
    
13. Elapsed time is 11.902260 seconds.

复制代码

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

1 C+ I6 O: i: lForcal 9（OpenFC演示）编译了如下简单长函数，第一次编译，耗时约0.688秒，以后编译耗时约6.3e-002秒：
# E! g& T, Y8 ]: l% w- m) m
5 {. t2 E6 f/ ?; {) af(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+
6 r4 G3 g6 m; N4 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+
, _* _1 e+ m/ p) @4 a0 ?" u... ...（省略）
$ k# _) z) S: {9 p( fx+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+x+1+
, c; h, {3 f" t- H9 ^0 U$ @1;

$ Y+ s/ a7 g: I$ [/ ?- P5 G完整的程序见附件： 多个x加1的Forcal代码.txt (133.36 KB, 下载次数: 0)

2011-7-30 12:06 上传

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

/ S: t3 Q$ B) G) j1 Q
测试如下代码速度：

1. mvar:
   7 Q& R0 s0 d8 _: L* ]5 L
2. t0=sys::clock(),% @; Q1 u/ h; r! Q, j( M- U3 L3 K
   
3. s=0, i=0, while{i<=1000, s=s+f(1), i++}, s;
   ( v1 R! L7 p+ b% {
4. [sys::clock()-t0]/1000;

复制代码

结果：
: ]+ \2 n8 |2 y% a1 |9 N67100033.
0.64  秒
5 f, i0 e  Y9 K8 F( m$ H4 v2 C% K
( t, z0 c3 p1 a8 [＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
7 `1 N2 u+ u. S5 l- E
可以看出，Matlab运行大的程序，仅编译就会耗费不少时间。

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

为什么Forcal第一次编译耗时较长（约0.688秒），而以后编译耗时较少（约6.3e-002秒）？

Forcal编译器在初次使用时，给编译模块分配的内存较少，此时若所编译的表达式较短，则编译所用内存足够用，在编译时就不用再分配内存；若遇到一个较长的表达式，编译模块就要重新分配内存，直到能编译该表达式为止，故耗费了一些时间；以后，若遇到的表达式都短于以前曾编译的长表达式，就不再重新分配内存，故不再降低编译速度。
# R3 ^5 b7 r- @1 s4 f
目前，需重新初始化Forcal，才释放给编译模块分配的较大内存，但一般不需要这样做。
4 F  o  g7 k$ f! f
/ o9 K- Y, x$ s3 l1 _3 A% @9 R对于脚本，编译效率和运行效率同样重要。

jmedwardliu1205

geilivable

安平池伟

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

yunbuhuiku

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

来电来函..欢迎洽购...
# o: O+ U! e, W& B$ X

wangchen881202

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

支持一下