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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
2 e" f- z8 s. [/ q
. w/ j) y0 ?8 u  H% ]) f＝＝＝＝＝＝＝＝＝＝＝＝
% g1 [( E! A7 E' C0 z4 @) V
- _# \/ y1 U& z% M( T+ QMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;7 \9 H& y4 y6 p1 u( h
   
2. tic;
   0 y; q- K2 t7 F\" x8 r) ~0 W
3. s=0;
   / e, v9 A- ?& W3 r. v
4. for x=0:1000& T5 t+ H+ ?9 X7 N8 }
   
5.   for y=0:100
   / e1 L: T/ `; g: @0 i2 C
6.     for z=0:1004 O5 E* _# X2 h6 m% ]
   
7.       s=s+f(x,y,z);
   ) `: T; H0 m) b) {$ F
8.     end
   * P& \$ H* A6 V/ |4 `2 N
9.   end
   . }( G2 d1 h* @
10. end
    6 I2 K& R\" F! A. c8 }2 `* T8 q\" x0 L
11. s
    ) P\" N6 e; `9 X- i2 t! a! k
12. toc8 s* z$ H4 o5 l6 U: ~
    
13. 7 ?* l- |# k& P( J. ]2 T- {+ _
    
14. s =
    1 G% }3 q, k+ A+ q6 P9 E# J7 N& E
15.     6.126720600000000e+009! S4 Q1 G! `+ j3 `9 \* I3 ?' X
    
16. 
    ( q2 \# m) ^4 ^% @# D
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
m函数：

1. %file xyz.m
   % I' t; S5 P# G6 B: P6 [+ R8 G
2. function c=xyz(x,y,z)% c6 ]2 Q\" w$ K5 m
   
3. c=x+y+z;* C2 a; q; Z4 C
   
4. end

复制代码

测试代码：

1. tic;( i2 ]; l# K& Y4 L2 k2 T* h5 N
   
2. s=0;
   9 y5 {  v# y; G, N
3. for x=0:1000: r# s) M+ S* p* `
   
4.   for y=0:100
   ; a' D. d4 w/ D# s
5.     for z=0:100
   : v, n) J4 Y( P. r+ C6 K( b; N7 W5 W* X
6.       s=s+xyz(x,y,z);
   / y' k7 U& g; p
7.     end1 x, `* z7 v, ^1 q3 J
   
8.   end0 q& M9 p' Y& y0 |+ s8 c( u4 d
   
9. end
   : g) R* ]/ G( d  @( j
10. s8 b4 I+ J# X: r3 {, h; q
    
11. toc
    3 U% z) S\" i* N  V9 r
12. 
    . @* o& C% K) P0 L# @# ]: w& T3 K- Z& ?
13. s =
    : O& F1 ^, L' Y% T! E  p
14.     6.126720600000000e+009& i  N. c7 U! D' x  p8 }\" F
    
15. 
    - d/ T! l$ _1 M# ]! _5 P
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

) W2 S$ G5 F) S' M! E+ Q/ DForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;6 L2 @: x. v) Q
   
2. mvar:
   * T0 c1 `+ t/ h6 p\" I0 t; w
3. t=sys::clock(),: G$ I' U5 @! R# _8 A( G- d\" a
   
4. s=0,
   , v/ M# ^* p% N+ V6 }
5. x=0, while{x<=1000,6 x* m' R+ w& W; C! S2 m
   
6.   y=0, while{y<=100,
   5 d  w# x) t4 H6 K
7.     z=0, while{z<=100,' g- G$ G) F' l3 a0 o4 x
   
8.       s=s+f(x,y,z),; H( y( N% `1 b
   
9.       z++
   / G, y8 F# L, g
10.     },' y7 ?$ v  o/ t/ W/ J( Z
    
11.     y++) H; d! u! U9 u' ~( R2 Y
    
12.   },
    # h( w2 x0 p& K& M* i
13.   x++, J3 P0 x3 k. b, V2 k
    
14. },
    9 s* v. N9 F& ?, n
15. s;
    % @% B9 ^7 R, \' o, u5 H
16. [sys::clock()-t]/1000;

复制代码

结果：
) M7 H+ o8 m' y3 C, O6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

1 Y! v. ^7 _7 v


) c/ q% A0 z, ]4 {( k3 m6 I

7 G4 B2 F  F4 M$ }* H  I
# L! b8 U$ B4 _  T. y( v



. q1 Q! k$ v( k. T7 ?





6 u# O" B7 s3 ~# }