极限测试之Matlab与Forcal编译运行大型程序的效率

forcal

这里没有实际的大型程序，仅对大型程序进行了模拟。

首先写5个几乎相同结构的m函数：

1. %file aa.m; m& p8 |0 ~+ b, X6 [* c( U3 k( T8 C
   
2. function c=aa(x,y)
   * K1 g8 v; u$ g# H5 o
3. c= a000(x,y)+a001(x,y)+a002(x,y)+a003(x,y)+a004(x,y)+a005(x,y)+a006(x,y)+a007(x,y)+a008(x,y)+a009(x,y)...
   : M  d, q/ `$ ?; T9 }
4.   +a010(x,y)+a011(x,y)+a012(x,y)+a013(x,y)+a014(x,y)+a015(x,y)+a016(x,y)+a017(x,y)+a018(x,y)+a019(x,y)...
   $ \1 T\" x3 L$ p8 r  ]# n! [
5.   ... ...
     Y# I  L) a* [
6.   +a990(x,y)+a991(x,y)+a992(x,y)+a993(x,y)+a994(x,y)+a995(x,y)+a996(x,y)+a997(x,y)+a998(x,y)+a999(x,y);' V; l, K( Q' y+ Y0 p# M
   
7. end% H5 p# A6 |: Z( U, I
   
8. / c! t\" P, {! J- p; C9 c1 ]' m
   
9. function c=a000(x,y)
   6 U* T% H6 ]5 q( i
10. c=x+y;6 ?! D, V3 z- x
    
11. end
    & V# O+ h; [1 @+ x* k0 ^
12. function c=a001(x,y)
    * d9 k: P  C7 Y, \4 D
13. c=x+y;& u, p+ R  y+ ?
    
14. end! n& }1 h! J, |\" F5 P1 l/ ~; D9 Q. {
    
15. ... ...
    7 B8 k6 k\" o  r* V+ D( S# y
16. function c=a999(x,y)
    2 Y7 \4 q* }5 t) _$ P% k
17. c=x+y;
    . p9 p$ x: h( T. j
18. end

复制代码

1. %file bb.m! x* {8 w# h9 d
   
2. function c=bb(x,y)
   & q1 T7 R# ]% z& T7 I3 {. Z5 @
3. c= a000(x,y)+a001(x,y)+a002(x,y)+a003(x,y)+a004(x,y)+a005(x,y)+a006(x,y)+a007(x,y)+a008(x,y)+a009(x,y)...! F# G\" c* v; l) m' N; ]) i: {
   
4.   +a010(x,y)+a011(x,y)+a012(x,y)+a013(x,y)+a014(x,y)+a015(x,y)+a016(x,y)+a017(x,y)+a018(x,y)+a019(x,y)...
   8 n\" p$ z& n' t& h
5.   ... ...
   % I\" g% K\" O: |+ j+ t
6.   +a990(x,y)+a991(x,y)+a992(x,y)+a993(x,y)+a994(x,y)+a995(x,y)+a996(x,y)+a997(x,y)+a998(x,y)+a999(x,y);
   ; w0 x% n* D1 X\" v6 Y1 k
7. end# T! D8 b4 B9 c7 j, `2 t% `\" B
   
8.   k  k) P1 A4 r7 ]8 K
   
9. function c=a000(x,y)- V' F* q/ ~: R3 S! n# _# B( b
   
10. c=x+y;
    6 i1 d# @' \. I$ M- N9 y$ ~& i4 ~# ?& P
11. end
    - f7 X6 ]3 \8 m+ M+ I( U. s
12. function c=a001(x,y). Y0 X$ h$ B3 @5 \( U6 Z- D( T% R
    
13. c=x+y;
    $ U! L+ z- M( d1 J1 I
14. end
    1 t* t: {% I) h) {# V9 D
15. ... ...
    7 N5 ^+ Q5 K. Z8 e
16. function c=a999(x,y)
    . N5 J3 r2 h9 w. J% S( s
17. c=x+y;
    / y) M/ K! ^3 O: @
18. end

复制代码

还有类似的m文件：cc.m、dd.m、ee.m

在matlab 2009a中运行以下代码3次：

1. >> tic;
   % N1 W9 H3 A: i
2. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1)
   $ e6 W. v8 ~1 s6 h% {/ P, T) }0 L
3. toc
   3 {- k* U' Y\" W
4. ! e* G5 V% r' v! X- G
   
5. ans =
   9 b) K: L, \! X3 c( n# u' L9 a' o
6. 
   / @! F; I0 G) W0 A$ E% U! B# b& W
7.        10000, J+ ]8 L$ J* t% u6 l% I
   
8. 7 C- a: z) d/ F( z! C* ^; j
   
9. Elapsed time is 1.853170 seconds.
   2 x8 l, @% }  H+ |+ E0 y) g7 g& ~  }
10. >>
    4 i5 P% B* Q\" X, m$ U
11. >> tic;) Z4 X5 |( |6 D: U/ ~' K+ y5 I
    
12. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1)
    9 l, q$ d# \  @; T( t7 _1 {3 n) a, ~: Q
13. toc) U5 [- E# ?1 _) v3 X& a: R4 a: {. s
    
14. 5 U; J0 v' n! c6 q( \! @7 Q% x1 w4 I6 \
    
15. ans =
    3 R3 J\" ]7 [/ @8 P) K% M/ f1 {
16. 
    ! R% P6 Q7 F5 r# K# R; \) E
17.        100004 d4 o  D' N; a, g2 B
    
18. $ H- [  F1 r1 F& }' p
    
19. Elapsed time is 0.051328 seconds.
    # z5 l# A0 Y& {$ L
20. >> tic;5 n0 J: N7 u- s1 l# \+ ]1 f) m  S
    
21. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1)0 |: _8 |; B: E$ q# E8 g
    
22. toc- v# P& s\" g- J$ x* I) e% m( u: l
    
23. \" v$ M, S% L4 L  u7 o. N( k
    
24. ans =+ C4 F3 z6 w# `( g
    
25. 
    & m4 V, F/ s  }
26.        10000' h; I2 }+ `7 K* n
    
27. 1 `. Q$ q/ H. {* L, y; I0 ~
    
28. Elapsed time is 0.058188 seconds.

复制代码

可以看出，首次运行时进行了编译，耗时1.853170秒，以后直接运行，耗时约0.051328 秒。

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
: ~, J3 B$ ~$ u( h9 p* M
写5个与matlab相同的Forcal模块：

1. //aa.txt
   . Q, r7 g& @+ Y& A/ N4 l
2. a000(x,y)=x+y;
   ) W2 {9 S2 U8 u' `  t- l) D
3. a001(x,y)=x+y;, ~3 ~/ R8 \, G& L
   
4. a002(x,y)=x+y;& }' n/ ^3 A8 G% V\" H$ h+ e$ A
   
5. ... ...& v  k4 N3 H! u9 n) T6 R
   
6. a999(x,y)=x+y;
   4 R- t& y% _1 B: v) h
7. / d\" Q: [) I, S2 S
   
8. ~aa(x,y)=& `; N$ E4 b7 ^* I4 W' k; ^
   
9.    a000(x,y)+a001(x,y)+a002(x,y)+a003(x,y)+a004(x,y)+a005(x,y)+a006(x,y)+a007(x,y)+a008(x,y)+a009(x,y)
   1 h' D1 [  d1 O/ k
10.   +a010(x,y)+a011(x,y)+a012(x,y)+a013(x,y)+a014(x,y)+a015(x,y)+a016(x,y)+a017(x,y)+a018(x,y)+a019(x,y): k1 J' j( t2 L+ b
    
11.   ... ...# D) M( [5 S; k+ d: n  _% k( B# K
    
12.   +a990(x,y)+a991(x,y)+a992(x,y)+a993(x,y)+a994(x,y)+a995(x,y)+a996(x,y)+a997(x,y)+a998(x,y)+a999(x,y);

复制代码

1. //bb.txt
   4 U# c- ~0 a9 q
2. a000(x,y)=x+y;* Z6 y+ f- _& y& h
   
3. a001(x,y)=x+y;2 J$ ~) g$ b2 k9 q+ A- e\" P
   
4. a002(x,y)=x+y;
   , W) ?# E# T: m2 c
5. ... ...9 U1 p\" w$ X9 l! H\" n& Q# b$ O
   
6. a999(x,y)=x+y;8 i2 n  T# D$ Y! ?# X\" P
   
7. 
   ( f' H2 v5 R+ l$ W9 U' `
8. ~bb(x,y)=
   . W, f/ B5 A+ O
9.    a000(x,y)+a001(x,y)+a002(x,y)+a003(x,y)+a004(x,y)+a005(x,y)+a006(x,y)+a007(x,y)+a008(x,y)+a009(x,y)6 J. @/ }0 |& {8 x. t; g
   
10.   +a010(x,y)+a011(x,y)+a012(x,y)+a013(x,y)+a014(x,y)+a015(x,y)+a016(x,y)+a017(x,y)+a018(x,y)+a019(x,y)
    8 b. ^: p0 j* f7 {8 F\" `! B# U! `( ~
11.   ... ...! ?5 X9 Z4 a% {, M
    
12.   +a990(x,y)+a991(x,y)+a992(x,y)+a993(x,y)+a994(x,y)+a995(x,y)+a996(x,y)+a997(x,y)+a998(x,y)+a999(x,y);

复制代码

还有类似的Forcal模块文件：cc.txt、dd.txt、ee.txt

在OpenFC中运行以下代码3次：

1. mvar:
   # R, Q1 D2 R0 M! X$ t1 ]' N6 G
2. !t=sys::clock();
   5 B\" e. o, I- O: `( X, I- u, i
3. #USE# D:\Forcal\Forcal32W\matlab\aa.txt;8 h' k# X7 _- a: Q
   
4. #USE# D:\Forcal\Forcal32W\matlab\bb.txt;  r! t- [# `\" g8 ?% t, k* K
   
5. #USE# D:\Forcal\Forcal32W\matlab\cc.txt;. u! f+ P( w5 a8 m& }. X& G, C3 W
   
6. #USE# D:\Forcal\Forcal32W\matlab\dd.txt;% `8 a  a0 u7 M+ x
   
7. #USE# D:\Forcal\Forcal32W\matlab\ee.txt;
   4 {! r' ]. R% z! l3 ^, R1 T* \+ Y
8. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1);
   3 k* c- w* P- u; `: E1 F7 \1 g# v/ I$ k( ~
9. [sys::clock()-t]/1000;

复制代码

结果：
" I1 x1 y0 y. ?10000.
5 h& ~& G+ v9 ^) U9 |0.141  秒
# S. h7 Y& s9 P( z/ F' s( w10000.
0.  秒
10000.
0.  秒
" Q# C2 i. r; j
可以看出，Forcal首次运行时进行了编译，耗时0.141秒，以后直接运行，耗时约0秒（耗时少，测不出）。
6 \5 s6 u6 a" _0 u: L
- V3 @1 a+ s8 s8 g! h+ f, q: {  J＝＝＝＝＝＝＝＝＝＝＝

% R* J+ w9 w) o$ I! m( G+ B- l本例模拟测试Forcal编译运行效率大致是matlab 的10倍多。
- I1 E1 j& i% n+ I本例的源代码下载： matlab Forcal模块.rar (18.76 KB, 下载次数: 0)

2011-7-31 15:27 上传

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

sxjm567

好东西哦，大家不要光看不顶

forcal

justftdsky 发表于 2011-8-1 09:51
1 r% I1 P* o; j% Z/ B: f我去。。。。。这个还是和电脑性能有关的吧

应该和电脑性能有关。
我是在同一台电脑上进行测试的，但结果也仅限于此电脑的此种配置。
若有不同电脑或matlab不同版本的测试结果，可在这里贴出，以便于比较。

justftdsky

我去。。。。。这个还是和电脑性能有关的吧

jmedwardliu1205

好啊谢谢啊