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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
) a5 ^2 r8 Y2 F7 R
＝＝＝＝＝＝＝＝＝＝＝＝
/ G* u# G% ^$ {( m# y$ I6 ]) V
5 r* }; y) o! Y6 G- SMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;) f( @# c* V$ w5 u/ m% \: R; n
   
2. tic;
   # O3 ]1 s3 B7 k* m\" K
3. s=0;
   6 K6 K( L$ m3 c) z. r; S* N& k  u\" q
4. for x=0:1000
   4 e; n8 Z! ?' s9 q% ]. X9 g5 W\" y
5.   for y=0:1009 T4 D5 T3 |* `, a: g
   
6.     for z=0:100
   * H! x7 b3 o3 _$ |4 n
7.       s=s+f(x,y,z);7 m) H; S! I7 ]5 i
   
8.     end, g- E* F( U1 N5 |+ H; Q  s
   
9.   end
     J# T8 f0 J% z/ q0 `
10. end& B4 O# ]9 ^6 K/ T. W% l: q9 R
    
11. s
    ( @5 M\" u+ M4 E) Y
12. toc
    : d# [! j# _# B3 U: s& y5 {
13. ; X; f+ W6 q\" E3 S. W% @/ S( ?& {
    
14. s =
    8 D; K/ E* [. H) |( C1 W
15.     6.126720600000000e+009# z1 e. }+ H* g, M
    
16. 3 j# ]  j% |- M5 [3 n5 O( L' l
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   / Z6 C2 O. u; Y2 e8 ]4 s
2. function c=xyz(x,y,z)) b1 T7 i  f/ N4 O1 j; y  ]
   
3. c=x+y+z;( E  |& c% H# ]8 ^: C# h6 G; m
   
4. end

复制代码

测试代码：

1. tic;# M; h# Q9 k) O  W
   
2. s=0;
   $ {4 ]2 H( ^$ z9 g
3. for x=0:1000( ^' _' W$ w$ B, C3 J. N
   
4.   for y=0:100
   7 N; `- a% k& V: ?\" e5 N. |
5.     for z=0:100
   % a  `) ^9 b1 q( J2 \/ c; O+ D
6.       s=s+xyz(x,y,z);- b; @: b7 F0 g. _: s, g
   
7.     end8 D! ]$ S& m9 L, e! _
   
8.   end
   ) q- F# S& o  v2 e% m3 t* u
9. end5 N& ~9 A) x3 v4 z
   
10. s
    : v\" R% ]# q( _% v2 S% N; g/ N\" l
11. toc
    5 C; K/ N- _, g  c2 F1 b- H4 s
12. 4 R: e! a0 ]; l/ J8 w
    
13. s =
    \" G9 e- e& D2 J) J/ R) ?\" S$ b, i6 R
14.     6.126720600000000e+009
    ; o! ~; a0 ?1 ?' t/ W
15. 
    & V' {. [- ?6 T, d
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
2 u' e. T0 s" n6 Y
+ P9 \6 J* ?7 I9 V- O: ]# T6 vForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;% X% k* E( b0 a/ ?' j' P; Q; J% K
   
2. mvar:
   & V. v: ~7 a! b  n
3. t=sys::clock(),
   $ K% C1 n# v4 D
4. s=0,8 o( t9 j# G+ B( _) l1 y
   
5. x=0, while{x<=1000,
   9 c2 B( X- Z\" B6 X  W5 F0 s
6.   y=0, while{y<=100,$ _7 x' r$ ~1 r' j, g8 p& p
   
7.     z=0, while{z<=100,
   6 x9 ?\" A7 _0 X9 ~( r
8.       s=s+f(x,y,z),+ r3 v8 P, o- ?& M
   
9.       z++6 W1 T9 ~& p) W2 @! i
   
10.     },
    & `0 W, r- `# g) r$ S
11.     y++) U8 ^: e; ^4 t* Q; T- c
    
12.   },
    % ^% i8 v\" s- M3 \% u
13.   x++) q, h, _* c$ e. D) j+ S4 t2 i
    
14. },* Y3 L+ R& C; R
    
15. s;8 ^6 B; m  t) X3 k% F9 @! k
    
16. [sys::clock()-t]/1000;

复制代码

结果：
$ ~4 k  O, k+ ~. R3 `, H6126720600.
1.109  秒
/ j( [& t6 p0 `, @; {" U+ P
* B- A3 l2 S- @, f二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
$ F" e* w$ f' g/ p

hxftntA

; z/ X9 a- C9 ~& G) S






, i, _  G4 u  V; i- r; a" [
2 [5 J, |8 V/ x7 ^5 z

2 ~/ M2 Q; H& C, _% i; I, U2 `



) u3 z% R$ D* H2 B. P1 a/ {