极限测试之Matlab与Forcal普通函数调用效率

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。

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

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   $ v* N; h! ]# p4 {, i' O
2. tic;( N9 L6 s3 E; l- P5 J
   
3. s=0;
   6 ^\" R& a5 C) c/ C. s' @8 d9 }
4. for x=0:10008 {- c: [2 j* [5 Z! F
   
5.   for y=0:100% x7 Q, m+ j\" i7 R3 b& S2 E; b
   
6.     for z=0:100
   * {9 i7 G% |5 b4 n( k0 \2 o
7.       s=s+f(x,y,z);
   ' [* j. s' c$ C/ W* v- Y/ q
8.     end! c- _0 M% g: `3 M7 S' ?) Q2 v0 k
   
9.   end) C\" P% E- g  V1 E% B: G9 j
   
10. end\" [# p: ~3 l\" T. ]/ g
    
11. s5 q7 \8 _\" ?/ @4 U% I
    
12. toc
    # v/ c2 i1 a& L# P
13. / i2 n# \  N; R. K' h
    
14. s =
    ' q3 z0 u: _) B4 g: j& o, z
15.     6.126720600000000e+009& q/ v7 v% y. D0 n, ?
    
16. 
    0 A. n% C. H1 U* w
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
# Y3 ^% e' C; k1 y! R0 R8 Km函数：

1. %file xyz.m
   4 t9 s- P) i  [0 B2 T# h
2. function c=xyz(x,y,z)' M* n5 {1 X. V& C& U& Y, X
   
3. c=x+y+z;  U3 i! h* t0 ^$ R( `0 N/ y( r1 W
   
4. end

复制代码

测试代码：

1. tic;3 N# u\" n. q  v; ?. i
   
2. s=0;
   ( L* c8 W/ G1 q4 A
3. for x=0:1000\" B: Z* T8 }( S, p1 [3 |# p
   
4.   for y=0:100
   9 d5 Z3 a! Z$ D& {4 O& L
5.     for z=0:100
   : d( A2 v) o) `  V( D& d* c
6.       s=s+xyz(x,y,z);
   ) i) _8 B! n2 ^8 t& @
7.     end
   . n& C$ |2 m2 P' p# h
8.   end
   + c& ], L* u# j; R% m
9. end) _8 y9 O& E$ h4 \$ s( t\" m- D
   
10. s
    $ q, V/ Z: n7 A. b) C
11. toc. j, q( Z+ n, a; y. ]2 d
    
12. 
    ) H  }\" ]# w+ [& e
13. s =2 s. f( P' X! ]
    
14.     6.126720600000000e+0094 L$ n8 P6 U2 _: ~) f: K3 e
    
15. 
    ! W4 S4 `. n\" G* f5 W, o0 U
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   , F0 z' v0 B+ ^
2. mvar:
   5 T5 d7 T0 T, d% c; Q3 N' w3 r
3. t=sys::clock(),
   * ~( p7 s, K% V9 r$ [0 K' @% U' b! X
4. s=0,
   \" f9 I2 X# Z\" U* u9 ]: e) a& P0 R
5. x=0, while{x<=1000,
   $ R! |' F# Y- D
6.   y=0, while{y<=100,+ p: f\" w* W9 n9 ~; }) x
   
7.     z=0, while{z<=100,
   ' ~: Q0 I- C/ N0 [9 }
8.       s=s+f(x,y,z),
   # G0 h6 e, v9 a- R
9.       z++% U7 A) n- v  m0 q
   
10.     },) e$ R. }7 e8 a8 }2 E: P5 f
    
11.     y++
    # M+ |3 E7 c. m( A( w8 \$ \
12.   },  ]0 y/ C9 W3 s5 H
    
13.   x++
    7 O8 O/ H8 z6 A2 L2 [* i) m
14. },* N; u, u1 x9 f+ b
    
15. s;: u( t* h\" j* A* t
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6 E& h& V2 y$ o- e6126720600.
1.109  秒
0 c( i) I. a1 D( |3 C( L
' _. Y: I  @. M  b# e二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

/ t- z6 V# ]" {) Z


+ ~4 X: D( h* l8 s+ Q3 q( f
6 E- ~2 F" ?7 |7 j/ k' w& ~" e: p* O
9 s7 n4 |! F% y7 ~, B3 o


4 {2 J/ g  [0 ^2 u& k2 e
' A" d' g" ]; W7 E* g7 k



* ~- n* N+ O. i
$ ?; k$ t# Q; o5 c- ?" [