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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
5 O9 \* a' H0 e2 R* A8 B2 X* x  P0 `9 K+ n
) h7 B7 f6 {5 O# {0 }＝＝＝＝＝＝＝＝＝＝＝＝

2 Z5 F4 E# t) G# uMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   * Z) O4 C; ]6 n; d% i/ q
2. tic;
   ! q1 u. u3 u* v. N5 H
3. s=0;
   + b& ~3 _/ p' B/ T
4. for x=0:1000/ P: e( S# q3 R8 ^  v' v& E
   
5.   for y=0:100: p1 _' l! |2 B  N5 c+ o
   
6.     for z=0:100& U\" f7 @* ^, m
   
7.       s=s+f(x,y,z);
   ' n2 X8 E1 P/ U  U
8.     end
   - Y, K+ g1 o8 I& @( x6 P
9.   end# J; L+ L2 Z4 X, y- r5 _
   
10. end
    - u3 }2 j, @: P2 U3 z7 I* q6 k; p
11. s4 Z  x0 o0 b: p2 y* j' t3 ~
    
12. toc
    - n5 S2 M8 {  o3 M$ Z. \) D( D
13. 
    & Z! u+ a6 a* h  Y
14. s =( m1 A% p' g/ j. c, I; V
    
15.     6.126720600000000e+009( |5 X2 \/ C- V+ P6 s
    
16. / G! E5 X* r, G: X' j% N' O- T
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
- L/ C7 i/ ^" k+ v* G% }/ Tm函数：

1. %file xyz.m& E% S! Z1 ]\" [+ y
   
2. function c=xyz(x,y,z), m7 F$ u8 y' d\" v- L0 E8 @. u; O\" |
   
3. c=x+y+z;4 ?0 J+ r5 b: ]- _; ]5 q* o
   
4. end

复制代码

测试代码：

1. tic;
   8 C8 R: ~, r+ `# J
2. s=0;
   \" c. |& ^% n3 M; Q* Z
3. for x=0:1000
   ( E2 @2 q6 l/ p# J2 \' n
4.   for y=0:1007 W4 t. y' B2 x7 G) O( F  P% m
   
5.     for z=0:100. d0 F% G# F+ a8 o: \4 d
   
6.       s=s+xyz(x,y,z);
   : r6 s. Q# d7 Z% c, |) {
7.     end
   - d1 ^& `  E+ {
8.   end
   % g$ S/ I. x3 v. y% w0 @/ S6 ]
9. end\" k' j! |' C1 ~$ k6 Q# x\" b
   
10. s) Q0 N9 o- B- A
    
11. toc- A# @& x* o7 i! h& \' W
    
12. 
    , \5 B! Y* g# H# L9 ^7 I! @: c
13. s =
    7 K0 [' x& U$ H3 r\" c& s% t
14.     6.126720600000000e+009+ O& t8 T9 W8 E- f& p
    
15. 3 R9 \% o) C! X2 ~
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

! T+ x; Q, i' ]Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;9 a1 b  \  S! s+ x- @
   
2. mvar:
     F  n/ j% W, S7 }9 k8 u( D
3. t=sys::clock(),\" J, t3 J, B/ _! X\" f6 R% k( X
   
4. s=0,8 g; I, H# m4 {6 T3 ?- H% D6 Q
   
5. x=0, while{x<=1000,
   : f: M+ U' b  Q! _1 B# s, T
6.   y=0, while{y<=100,
   \" L: A; L+ K) l+ `
7.     z=0, while{z<=100,
   * m; k' L- a( q
8.       s=s+f(x,y,z),
   ' L7 v: M  F/ t4 t# }
9.       z++, [/ j\" V% g6 J6 k! [& O9 M
   
10.     },
    , S% f; j  U2 E2 _0 G6 Y( U) C/ D' _
11.     y+++ s( a! U0 T( w* C8 _; e; p
    
12.   },
    ! X* D# U3 r7 G6 a; b/ p) ?8 A7 u
13.   x++5 N/ ~/ B& r$ D3 A
    
14. },- _# H5 L\" o, H3 J$ q+ e3 Y, b8 M4 y
    
15. s;; j5 I+ x3 {1 q
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
( d. Q+ \" n+ D" u, a! m
* w) x3 z" ?! w! i$ L1 _; m二者效率有几倍的差距。

hxftntA

. c6 Q4 W* n% s1 e7 J# r6 I% J. s
! _/ b: Y. r. x  d




7 }8 }/ Z, A$ D

2 A4 h, X- Q1 w4 U
( A) w, ?% i, w) c/ |* T8 Y+ T
; {" k; {6 q4 n0 T" \1 X
9 K% Z7 q# r* I* ^7 R

% Z; Z- y6 i( }
, Z, W: ~( |9 F/ ~4 l
+ X; v+ J7 W; ^% c& y
' P/ N: |8 _8 T3 x( n# n. b

ehi28

嗯,不错,支持一下.