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

forcal

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

' |0 c, z6 d% U* v首先写5个几乎相同结构的m函数：

1. %file aa.m
   - c4 ]8 t# Y8 |# w# G
2. function c=aa(x,y)8 z3 D9 ?0 @  w2 F9 E
   
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)...
   $ \: g2 e4 @2 g  T6 v1 p7 z
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)...
   7 v4 \6 N/ x% H0 Y5 O
5.   ... ...  o8 U9 w7 Y  A0 z: X0 b8 [# q; _
   
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);$ f; N# z\" q! @' M
   
7. end
   7 V* Z2 o- z7 V: |% R* N5 h
8. & v' X7 n8 F+ w4 b
   
9. function c=a000(x,y)
     H- S) [0 ~1 J! G! }
10. c=x+y;
    3 S8 y4 {- j0 c5 `% E4 N. [! ^! t
11. end
    \" v7 B% w: _( v* F- ]\" f
12. function c=a001(x,y)
    : j/ C! u$ a3 }/ C7 E* g\" r
13. c=x+y;
    3 I# u+ P\" C2 U, G
14. end
    & p7 [3 M; F, X$ M: M% ]
15. ... ...- z  l7 Z) {! E& ^1 d8 U
    
16. function c=a999(x,y); ?3 T# o) d$ ^: A( @
    
17. c=x+y;
    / w0 [% s4 E$ G
18. end

复制代码

1. %file bb.m
     O9 B7 H; `$ {  `
2. function c=bb(x,y)
   3 R/ E: }- k1 B\" {- U, Y
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)...% A4 Q: Z+ i/ a; i: k
   
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)...4 e9 Y  g4 r1 P3 x
   
5.   ... .../ }9 ]# H- l' W# z  S1 l: @
   
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);
   8 }: g/ {7 U5 d2 N/ B! a) G
7. end
   4 N, t% _+ p, R0 v) S
8. / `7 ?5 ^\" n- R8 `& d7 G. L: r; a
   
9. function c=a000(x,y); i* a0 I. @7 }2 g. [# w+ X+ A
   
10. c=x+y;
    2 X' O4 g/ o0 r& t' M0 p# o7 t
11. end\" a8 I7 Y* e- t/ A9 h, G% G4 ^
    
12. function c=a001(x,y)! P6 U- C# b8 I4 c# A0 F
    
13. c=x+y;  `/ \, g! t, A
    
14. end  q- h/ s  Y+ ]\" e/ {
    
15. ... .... F% ^  v/ l; w9 S
    
16. function c=a999(x,y)! \\" B/ S5 c& s2 p; W
    
17. c=x+y;
    5 h9 A( H  a& D8 l( b
18. end

复制代码

还有类似的m文件：cc.m、dd.m、ee.m
/ ~: O+ ?# O/ Q
在matlab 2009a中运行以下代码3次：

1. >> tic;
   ! @' E% D: ]5 ]; X% C5 m9 n. m
2. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1)
   - M4 y# C. Y, y8 e
3. toc
   \" D\" ~& N. u3 D3 T# u  E5 b0 j
4. 
   2 H* U+ v4 P# h
5. ans =  X7 R& b% t+ N5 a) c  q0 I( H
   
6. 4 b2 r5 i& C3 P. P/ O+ ?
   
7.        10000
   ! \' A& P+ K( j5 X' X/ y2 g6 Z8 q
8. 
   , |, e+ Q7 S( r) w1 ?2 i
9. Elapsed time is 1.853170 seconds.! Q  p: K: i, J0 w\" F6 y0 H
   
10. >>
    * C5 \  v, S9 C( i6 B
11. >> tic;1 Y4 l- N7 y$ _. L4 \
    
12. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1). G; x: \4 X! [; k
    
13. toc9 ~6 h  `: w3 f. n$ Y8 Y4 Y
    
14.   }' Y) W2 M) I) w3 O1 g
    
15. ans =3 v5 n  E4 X1 G5 H8 J  K
    
16. 2 y, c, F3 K' u! E
    
17.        10000
    3 j7 V! z# J/ [) X- ]9 E7 d# a9 [) G1 x
18. 
    . D4 A* P$ w) K) X- s
19. Elapsed time is 0.051328 seconds.
    ; L$ V$ i0 |# t9 _
20. >> tic;
      D! k7 a( d- \7 V! o5 l8 T
21. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1)' O3 U0 b$ |) \! F5 e% {
    
22. toc
      o7 M\" y( R0 y% r( \3 @
23. 5 g5 L* d' q. q# _/ y
    
24. ans =' Y- @; j6 M1 {: }' s( Z% j+ ]0 w$ l0 J
    
25. 
    1 Y* ~* K2 O7 U- m; b3 {4 e+ f+ x0 B& K
26.        10000
    ! J+ j& h( A  t6 a7 i: {
27. ; `8 ^1 L( K9 e$ G- N- m) x! }, Q
    
28. Elapsed time is 0.058188 seconds.

复制代码

可以看出，首次运行时进行了编译，耗时1.853170秒，以后直接运行，耗时约0.051328 秒。
0 c4 a6 {# G: W7 R3 `; `
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

写5个与matlab相同的Forcal模块：

1. //aa.txt. ^9 F* \\" P  T- L
   
2. a000(x,y)=x+y;
   ( N- n0 _: h' K8 l/ ~6 c
3. a001(x,y)=x+y;7 V) `! ~% r  V
   
4. a002(x,y)=x+y;
   4 f0 W$ d! X9 Z1 l7 M
5. ... ...
   0 J! T6 a2 c. d4 b
6. a999(x,y)=x+y;. _6 F& H8 P+ B$ D
   
7. 
   1 V' v5 T! ~9 G
8. ~aa(x,y)=7 p/ N% [$ `0 _3 @* X
   
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)
   / N) Z  r3 y  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)4 J9 }' q( j7 \7 _
    
11.   ... ..., `* J, D% c( g- {; S, ]0 {5 v
    
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
   1 w4 P# V' W5 X\" ]) s' r+ s! q
2. a000(x,y)=x+y;+ `8 R/ F: w5 s$ L6 A/ ?9 v9 k! ^
   
3. a001(x,y)=x+y;
   2 }5 x( _7 r/ X* d
4. a002(x,y)=x+y;# ^( [\" e8 @/ ~5 _& p
   
5. ... ...+ @; J# d# V; o
   
6. a999(x,y)=x+y;
   1 ^) {! Y/ |9 d; Z2 q% |
7. 7 R2 k1 |1 U$ W( \, Q% N
   
8. ~bb(x,y)=
   0 \4 M0 J0 d5 O+ z
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)0 T+ {: m$ v' \* r: l8 }$ s
   
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): [\" d# O6 s8 {\" m/ |% n: T4 b$ @
    
11.   ... ...
    ( d' U: t7 N2 d9 h1 y/ d2 w
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:
   - O5 Z8 A  @& n( E
2. !t=sys::clock();
   % M5 Y6 E$ H! y$ {7 E5 J0 G% J6 H
3. #USE# D:\Forcal\Forcal32W\matlab\aa.txt;1 p' n& e+ O; p: ~! T2 B6 N. |
   
4. #USE# D:\Forcal\Forcal32W\matlab\bb.txt;+ t5 L; S2 H5 U\" x6 g, E4 d% d
   
5. #USE# D:\Forcal\Forcal32W\matlab\cc.txt;5 \\" r$ P/ {( u- X$ A3 u
   
6. #USE# D:\Forcal\Forcal32W\matlab\dd.txt;
   % D3 g# D, N  x! v
7. #USE# D:\Forcal\Forcal32W\matlab\ee.txt;
   . w) e% |8 x; T9 y8 t) \\" L$ I
8. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1);( w- T' Q; m6 c7 x
   
9. [sys::clock()-t]/1000;

复制代码

结果：
10000.
+ Z; j; L" w) t& u0.141  秒
10000.
0.  秒
10000.
0.  秒
4 l/ W9 u& {# h& n; F
, M9 ?) T( x/ A0 V' k, M可以看出，Forcal首次运行时进行了编译，耗时0.141秒，以后直接运行，耗时约0秒（耗时少，测不出）。
/ r& t( V8 L6 ^: }) g3 e+ S9 |
; D# R$ N0 w+ s% s* A7 W- I) G: H＝＝＝＝＝＝＝＝＝＝＝
8 p  _0 m) H2 s, t
2 Z* \5 @& |9 i2 ~本例模拟测试Forcal编译运行效率大致是matlab 的10倍多。
2 f# a- E# T/ S本例的源代码下载： matlab Forcal模块.rar (18.76 KB, 下载次数: 0)

2011-7-31 15:27 上传

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

sxjm567

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

forcal

justftdsky 发表于 2011-8-1 09:51
我去。。。。。这个还是和电脑性能有关的吧

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

justftdsky

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

jmedwardliu1205

好啊谢谢啊

forcal

为了使测试更符合实际，改写函数：

1. function c=a000(x,y): B& W4 p# A9 x\" s* i
   
2. c=x+y;' a9 c- M* C% j6 G2 |$ I6 o
   
3. end

复制代码

为下面更一般的形式：

1. function c=a000(x,y): z( w, [+ @( c, R& y3 {
   
2. s=0;
   5 W( L; g$ {6 p
3. if x>y
   ) \, x% G% q( s/ c( v6 O
4.   s=x+y;+ K% y! V! _3 B
   
5. else
     Z5 R; b$ A) H: M
6.   s=x-y;& \: j\" W& T. P
   
7. end
   * T, {2 ]* ]3 e4 P
8. for i=0:2: B. V# x7 `. `- {\" q\" y
   
9.   s=s+i;# h$ \# U% g6 i. Q' H6 _, p
   
10. end
    5 z& m2 P; G% ^3 h' Y
11. c=s+x+y-x*(exp(y-x)+sin(x+cos(y))-0.5);+ C& B( p& ?. U- V\" v) w1 t
    
12. end

复制代码

对aa.m、bb.m、cc.m、dd.m、ee.m中所有的类似函数都进行改写。
9 B  B4 |/ u  f* T- Y
在Matlab 2009a中的运行结果为：

1. >> clear all
   \" u( P  W2 r- h: n, u# r$ p
2. >> tic;/ a8 t/ \% ]1 B+ B$ ^0 z0 m
   
3. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1)  H7 s- r- w2 n9 P! O
   
4. toc+ P8 o, L& i7 |
   
5. + M( E  o; o9 X; ]
   
6. ans =
   $ M6 z\" q; C' c8 f; C* {) s
7. : _* u5 Z7 t$ P; s
   
8.   1.7502e+004
   : F# R, I/ r1 j4 m/ Y& j
9. 6 W- i3 x  ~& U* b) A
   
10. Elapsed time is 8.762889 seconds.0 j. N# Q( X2 Y' [6 M
    
11. >> tic;9 _/ ~) B  _9 D
    
12. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1)
    1 k) L' V0 Z1 W9 `* v) \# r- x
13. toc
    ! `6 T2 b! j5 Z3 H: n& D& M+ z0 N
14. 
    , S4 ?$ @, j3 s1 C7 U8 e* H
15. ans =
    ( V) Z9 Q- X4 X/ F\" l5 t8 g; O
16. 
    $ o, T9 U! e) |6 w
17.   1.7502e+004
    8 f; C! A% O* P; J9 r- E2 a+ p
18. 
      u  U& \8 D) p! m
19. Elapsed time is 0.069484 seconds.5 F! @\" N) X3 v: m5 p! d! b1 r
    
20. >> tic;1 S\" M5 b  E  y/ @
    
21. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1)
    ; T* `# N- p) I) _  c) b
22. toc. F! m' Z6 U8 D, H0 N
    
23. 
    : }; }$ a, G( p\" v4 A6 K
24. ans =
    $ r# S6 I# [& P' A0 R! i
25. 
    : g3 @0 d9 e* ^\" m, Q/ ?
26.   1.7502e+0042 f8 M$ Q0 M7 K# D! s) D\" v
    
27. 
    0 M7 i% w0 R9 @, k+ |8 H/ Q6 h
28. Elapsed time is 0.068631 seconds.

复制代码

和以前一样，第一次运行进行了编译，耗时较长，其余2次不再编译，故耗时较短。
& D8 Y' x5 _. J4 i) y/ `
' [& y0 A; Q& c  c7 t* O' o2 L＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
( K4 r3 G4 X- Y7 @6 Z0 n) A! T
5 F" |4 {. T' I9 _3 A$ t对Forcal程序做同样的更改，即将函数：

1. a000(x,y)=x+y;

复制代码

改为：

1. a000(x,y:s,i)=
   * ?6 n' e5 i! {! K6 i) c
2.   s=0,
   9 a) j, v' J, q4 A) H: y& r
3.   s=which{ x>y,1 |  h% d3 _0 H) R7 j
   
4.     x+y,+ E) p( i3 a: C( V* B, c
   
5.     x-y% h+ B: r: j, L$ H2 |. N1 R
   
6.   },
     C$ [. r, u$ k' Q' @& p
7.   i=0, while{i<=2,
   ) w\" `: p+ @% N  t8 k$ ~1 T6 W
8.     s=s+i, ++i
   + {0 b7 {6 ?' c; l\" B3 f1 g. g+ C3 f
9.   },
   # G\" @9 m9 m0 ^$ D\" R. I\" F3 U
10.   s+x+y-x*(exp(y-x)+sin(x+cos(y))-0.5);

复制代码

其实改成以下更好，但为了使代码看起来比较长，还是采用上面的写法。

1. a000(x,y:s,i)=
   $ b2 z( v9 [4 }4 `, _
2.   s=0,
     @$ D$ P3 `2 x4 I0 y\" K
3.   s=which{ x>y,   x+y,   x-y},; y! M8 m1 r9 N\" O4 y
   
4.   i=0, while{i<=2,   s=s+i++ },
   4 Q- h/ K; ^+ C! N; O4 ^
5.   s+x+y-x*(exp(y-x)+sin(x+cos(y))-0.5);

复制代码

对aa.txt、bb.txt、cc.txt、dd.txt、ee.txt中所有的类似函数都进行改写。
( H- l4 D+ w+ {. S, L- `1 e$ _
3 v, n- p$ m! M& N2 |在OpenFC中的运行以下代码：

1. mvar:
   1 Z3 d# |0 u( J$ U. j& ~0 o! [
2. !t=sys::clock();
   0 x8 j& h8 q  G
3. #USE# D:\Forcal\Forcal32W\matlab\aa.txt;
   \" f/ _: F# X9 b* H\" i3 y
4. #USE# D:\Forcal\Forcal32W\matlab\bb.txt;# |5 g9 _/ L' G7 n( k) d: W
   
5. #USE# D:\Forcal\Forcal32W\matlab\cc.txt;
   \" M9 K9 s: H6 q: P
6. #USE# D:\Forcal\Forcal32W\matlab\dd.txt;
   , T; E: G4 h& q( e: ]' F3 a' X
7. #USE# D:\Forcal\Forcal32W\matlab\ee.txt;  j0 L1 w6 W. q( e4 K- W
   
8. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1);& x5 l4 U( P- I/ x: Z2 A
   
9. [sys::clock()-t]/1000;

复制代码

结果为：
17502.32453314312
0.281  秒
17502.32453314312
0.  秒
17502.32453314312
# o- m- j% u4 Y- S1.5e-002  秒
! c2 t; y: v* v
) a& i$ l* W) o  R1 N( N  z5 l＝＝＝＝＝＝＝＝＝＝＝
0 B+ {1 D0 A2 ^. I: U
可以看出，随着函数变得复杂，Matlab编译耗时增加较快，本例中，Matlab与Forcal耗时之比大约为31:1。
6 V( R7 T0 Q6 ?3 T- n6 i
( n: [3 f! a# {% J＝＝＝＝＝＝＝＝＝＝＝
! u7 `: l! b8 Y# c
; ^4 `+ B8 D$ o% x, ]' H- P) ]4 W为了比较Matlab与Forcal运行耗时，多次运行以下Matlab代码，时间基本相同：

1. >> tic;8 w5 o8 S1 }( m7 K0 v2 o) V
   
2. s=0;
   \" \3 z. c, C' m# M  x
3. for i=1:108 X! S/ n! g2 `5 G& [: U
   
4.   s=s+aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1);
   ! [: M$ A1 K' @' n) K+ S
5. end4 w' v\" F3 C. P3 U5 W5 T  ?. v\" Q
   
6. s: V/ A5 E* k) S: `+ p* {
   
7. toc- M( L! V' w6 i- {; s
   
8. 9 x# S/ o6 s% Q) L7 a! A4 ~
   
9. s =
   0 D) I: t: I4 ~' G
10. , t\" m. }0 D2 p  \! n; ~
    
11.   1.7502e+0055 @, O4 ^. G: x4 D6 v
    
12. 
    6 O5 w- u% M\" d9 }( c0 m) n
13. Elapsed time is 0.426280 seconds.

复制代码

多次运行以下Forcal代码：

1. mvar:
   - `. x8 }\" T# N' l! ~
2. !t=sys::clock();9 C- g) U% `8 a( S8 z& Z  p
   
3. #USE# D:\Forcal\Forcal32W\matlab\aa.txt;
   ' U0 A; r6 {& ^2 _% b0 m) H/ o$ @' s
4. #USE# D:\Forcal\Forcal32W\matlab\bb.txt;2 Y$ N9 a6 z% |! P9 Q
   
5. #USE# D:\Forcal\Forcal32W\matlab\cc.txt;
   # b( b: O( i; g) ]# r0 N2 u, N6 h
6. #USE# D:\Forcal\Forcal32W\matlab\dd.txt;) C3 \7 {+ f8 V1 j( J& Z5 ^
   
7. #USE# D:\Forcal\Forcal32W\matlab\ee.txt;. N) `5 t, f3 y. N( I* ], q
   
8. i=0, s=0, while{i<10, s=s+aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1), i++}, s;
   ) c4 q& H( d% X# Q
9. [sys::clock()-t]/1000;

复制代码

结果：
175023.2453314311
6.3e-002  秒

故本例中，Matlab与Forcal运行耗时之比大约为6:1。
本例源代码下载： matlab Forcal复杂模块.rar (22.03 KB, 下载次数: 0)

2011-7-31 15:35 上传

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

forcal

实际上，Forcal的以上5个模块可以放在一个模块文件中，或者直接写在OpenFC的代码窗口：

1. mvar:
   ' @0 a, K) `3 h$ O5 t$ S
2. !t0=sys::clock();    //用t0记住此时的时间。该表达式编译后立即执行。
   7 p2 O( h! h( @3 c8 u
3. ; C4 o* L\" }& v. y9 E4 s7 J+ \
   
4. #MODULE#
   : M; O! S3 Q$ }
5. a000(x,y)=x+y;
   7 r) v( U4 H+ ]; U2 F
6. a001(x,y)=x+y;
   & n  Z9 @; q9 e3 N! R
7. a002(x,y)=x+y;
     {, ?& w6 W% B8 G6 Y
8. ... ...7 V+ o+ z: ?6 L1 U' D6 t
   
9. a999(x,y)=x+y;# ]$ M; t- y; h/ V: O) m
   
10. 
    ) z' X4 W* p4 h5 R
11. ~aa(x,y)=( i5 ~4 J; l* f  N* K1 R
    
12.    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)
      b2 O0 z- A6 _. ?' I
13.   +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)
    3 |9 B+ ]! O( p# g. O) ?
14.   +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);
    6 n5 i1 e( p. T) [  c' Z
15. #END#1 q% S0 h! h6 N* z  M6 [  e! T
    
16. * o$ [9 j4 b2 ~$ C
    
17. #MODULE#; L2 }6 }/ g( o
    
18. a000(x,y)=x+y;+ X\" d0 Q7 S! i! m6 ^
    
19. a001(x,y)=x+y;0 V  l2 q) j  q5 e\" {% C
    
20. a002(x,y)=x+y;
      Y* m' T# G7 g# H/ |5 K; @
21. ... ..., R\" }6 u2 H5 U0 @) w6 S& s
    
22. a999(x,y)=x+y;
    8 |\" @* K5 `, x
23. ) k& B9 a* A' z7 U9 w
    
24. ~bb(x,y)=7 d' I# S: Y! A8 {
    
25.    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)- u( ?& i5 C1 y- K0 r
    
26.   +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)
    7 T0 S8 b& B0 H
27.   +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);
    8 d) }8 ~5 v2 U5 y/ B% a1 E
28. #END#4 J: ]4 e- d: O* r0 E+ |& a# b3 f
    
29. ) [% ?- [& F9 Z0 h
    
30. #MODULE#1 ~# F# ~/ A9 V1 z. T
    
31. a000(x,y)=x+y;
    0 A, t% W3 T+ m
32. a001(x,y)=x+y;
    ( R8 u: {\" ~1 |2 b' U8 G
33. a002(x,y)=x+y;
    9 d  |1 Y( _) c- z1 d0 D
34. ... ...  W; @/ l5 i. C\" M% ~& Y6 k\" n
    
35. a999(x,y)=x+y;
    & W# J1 L# h, Y
36. 
    + r8 `+ @) j7 X1 V* W& M7 g
37. ~cc(x,y)=
    5 b3 j8 T6 o: z; h  m! ~
38.    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)* Z- _* Q( G2 Y* c$ n, U
    
39.   +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)! S+ |6 c6 q* u2 |8 u& s
    
40.   +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);7 p2 x- a$ L- Q- H9 L  m* ], {
    
41. #END#2 D; g% W, v/ L2 _# }' M  a
    
42. 3 Y; j2 X- J+ U
    
43. #MODULE#; d\" Q. |+ [2 {' R/ l8 v, l; `
    
44. a000(x,y)=x+y;
    % W1 u\" w9 J! v9 n6 ]3 q* E6 g
45. a001(x,y)=x+y;
    4 ~& j6 o' ~9 P# ?; z' B
46. a002(x,y)=x+y;
    # w$ m0 _5 a8 g4 a+ D
47. ... ...
    5 E3 P3 N  W\" l. |2 u! C* ?! n
48. a999(x,y)=x+y;0 R\" M5 T. r+ z# P9 P& c
    
49. 
    6 ~: G1 t0 q* ]
50. ~dd(x,y)=
    8 X0 \# l; q. @! H/ O5 F
51.    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)
    ; O/ ^* k0 O1 f& w% D0 e/ l
52.   +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)
    ' K! f; Q1 ~' Y7 Q! p5 C$ w
53.   +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);
    0 n- Z( b* R/ _3 T
54. #END#
    # G( F4 B$ O/ a7 \3 ?6 ^- F! \
55. 
    * X, Z: f9 @\" i. o
56. #MODULE#4 g2 n0 A3 F: U4 Y0 G8 f
    
57. a000(x,y)=x+y;0 O, o2 b3 h( F6 j5 f. q$ J6 S
    
58. a001(x,y)=x+y;3 O, @) P+ J8 o$ k, M. ^
    
59. a002(x,y)=x+y;3 R\" G5 J( S. {. q* ]- I( ]0 |; Y
    
60. ... ...
    3 Q, s5 Y) D+ D5 ~2 D
61. a999(x,y)=x+y;
    1 Y4 M! U- F4 M. H
62. 
    + D* H2 G; r+ p
63. ~ee(x,y)=7 K: I3 U2 z) e6 s, }) B6 b
    
64.    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)
    0 Z5 h0 ^2 ]* S4 {9 A6 J5 q2 s
65.   +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); z2 H. n' {! p; R1 }( b4 e& q. o
    
66.   +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);
    0 E/ ?; Y. n7 J7 j4 E
67. #END#
    4 z2 D  g: y+ U4 e
68. ) Z. G1 b! |2 H3 y
    
69. mvar:
    4 S/ d& n1 D2 `3 @9 J7 L
70. aa(1,1)+bb(1,1)+cc(1,1)+dd(1,1)+ee(1,1);/ W% F$ r8 M8 q: p
    
71. [sys::clock()-t0]/1000;

复制代码

结果：
10000.
9.4e-002  秒
3 b' F* V6 C; a8 y, H" ?
( W! |3 q/ ?/ C+ c: d( U: q多次运行的结果相同，每次都是先编译后计算。编译和计算共耗时9.4e-002  秒，由于不读文件，故耗时减少了。
2 i4 M8 ]  t; G0 s+ L2 d
1 O% s* ?; \  S6 z9 C" a7 n每个模块文件可包含许多子模块，每个子模块可输出任意多个全局函数，这是Forcal的优点。模块及函数不必存为磁盘文件也可编译运行，也是Forcal的优点。