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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
; ]9 C. E7 \/ c0 E; \; t
＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   $ U7 z\" P- }5 v; i4 t8 r/ U
2. tic;! w5 M5 N* i! m
   
3. s=0;
   0 E) ]+ o# `( u. @  r6 A
4. for x=0:10002 O7 v- ~. `* J5 y$ Y
   
5.   for y=0:100
   0 _\" J8 _3 h& W5 ]\" w+ j3 j6 o! ~
6.     for z=0:100
   6 r3 ^: [2 D# }: d
7.       s=s+f(x,y,z);/ r\" T' r$ E8 ~. I, N
   
8.     end
   ! D; O7 u; K0 f) @5 L
9.   end8 d) s  B' v$ o8 O1 A
   
10. end7 S! \) t6 @/ v3 ]5 T
    
11. s
    1 J5 _  v9 A/ ~
12. toc
    ; S$ Q/ n4 w( w& Z: W3 W
13. - }% R- x4 T+ D
    
14. s =2 y* w6 A/ J, g+ ]
    
15.     6.126720600000000e+009
    * E, e) D, I7 w9 ~\" ~) t
16. % r$ s: m\" c0 Y$ D. G
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
, b" ?/ [# V; U( i& `m函数：

1. %file xyz.m# U( ?/ z  I  ~
   
2. function c=xyz(x,y,z)
   ) p4 n. v; J3 n1 |
3. c=x+y+z;
   . _# e: S( C  [  j# t
4. end

复制代码

测试代码：

1. tic;
   / b& U6 s$ D9 @' q1 L/ J& k4 w& g3 V
2. s=0;
   & E7 K' u& h0 C! I% W7 c
3. for x=0:1000
   ; ~. I8 |2 Q' }' L3 O
4.   for y=0:100; \+ r- d8 ~0 g; i: Y
   
5.     for z=0:100
   2 l8 n, W% n3 |; _* U0 K5 e
6.       s=s+xyz(x,y,z);
   ; N7 I+ T0 B/ g% _
7.     end3 t( J- G$ O+ q0 \( t$ d4 o
   
8.   end
   / j2 t) A: N6 g
9. end
   2 ?, a& k5 a8 L& ^& R
10. s
    ; z% }- d: F% X2 T8 o
11. toc
    : L( w( F4 r6 ]6 y* V\" o# ^
12. 
    9 i1 S\" U+ A# L2 G) X. x+ N/ y1 z, }
13. s =
    . A0 E  e) v/ i0 X8 {$ i
14.     6.126720600000000e+009
    - E2 n; g8 f3 o/ r0 r: i5 G2 H
15. 
    $ e5 P0 [4 i; M* o6 |
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
0 T4 G$ ^1 N6 C$ Z- ?/ N7 X; O
5 H/ D" c( C! P  N$ i$ C5 y  s; Y+ ^' wForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   3 Z$ G) G1 w( m* ?/ W
2. mvar:
   6 q9 |8 V8 m\" R' _( K! V1 T
3. t=sys::clock(),3 ~9 i6 k( M, ]! v+ ~9 K
   
4. s=0,
   6 J5 a  |. O' r% y8 a% h, d! N
5. x=0, while{x<=1000,  i- \. i9 S% X2 x0 y& n. |
   
6.   y=0, while{y<=100,
   ( K1 z; I# r* N3 D
7.     z=0, while{z<=100,
   / c2 Y5 P' K+ G4 _$ M
8.       s=s+f(x,y,z),
   - o; u  K* J  q. S8 L
9.       z++: Y4 R/ i; {\" v- R! s
   
10.     },
    ) q. ]6 V1 n  L6 {3 J! [; ~# y
11.     y++
    ) w$ V* A) ?8 U9 d: p
12.   },
    * s) S) A, m5 G
13.   x++/ Y3 K5 K( O, Y7 A8 H# T! G) `
    
14. },
    7 F  R) Y$ G+ w! J9 O
15. s;
    & ?+ ^6 C\" z% ]\" s\" D
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
; i* a9 R7 M& c0 c0 {1.109  秒
8 t9 G9 g# Z! n! H0 Z
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
. i# H& w$ d2 o5 K3 Y8 F1 D; K

hxftntA

5 F4 V: l; T, j! T& s3 O1 Y
4 d$ V% i' ^, g! h4 o- K5 u% D
5 r8 N  B1 d# d1 R, n6 v* O

2 j/ |2 T, `. o- W" f5 o( D" ?0 y
& K# k5 n  x3 W


9 G9 q, ^7 y1 g5 H) [+ V
1 Z, ]; |& L2 v, Q6 |6 b! a
& d1 W* U  @4 l1 N5 l
5 f/ p4 v( W5 O$ W) `6 V


. D$ k7 O: _9 S# _$ Y* q9 ?: Z