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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
# w( k( g7 A( l$ r7 p# ?
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   3 c, [9 T0 r\" T7 c4 Z; ?5 o
2. tic;
   % e' x4 z% N) h) w- A0 a/ B
3. s=0;
   9 g! E0 U\" b0 V# Q1 |+ ~
4. for x=0:1000
   ' H0 w* N, j0 O9 K) x' M, r% P
5.   for y=0:100! i\" p  x9 c, b1 ]# ?- z8 J
   
6.     for z=0:100( @! s$ [3 s2 x  U
   
7.       s=s+f(x,y,z);
   ( h+ S5 h0 d6 Q5 c; B
8.     end: y% k1 k8 r: V; y
   
9.   end. b! `  p$ C8 k) b6 Q* s5 o
   
10. end
    2 V; F% H. G+ Y& p
11. s
    - F\" ]& @6 I* I\" l1 R7 [5 \
12. toc
    8 s8 h. T6 \  |9 C  a  P
13. 
    * E( _+ i5 B$ W& ]) U8 G
14. s =
    5 h/ |( K* T9 \' X\" e
15.     6.126720600000000e+009
    0 {% T/ ~% E: X. ^
16. 
    - G% P2 S6 W9 m- A, n4 o
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
0 B3 ?) |& r7 u! r" G5 k) F" `m函数：

1. %file xyz.m
   3 ], S6 S* x  L9 ?, i, J, a: Q( Y
2. function c=xyz(x,y,z)
   % E5 Z. U2 c  Y* k  n* s2 Y
3. c=x+y+z;
   % k8 ~( C2 m+ K* S
4. end

复制代码

测试代码：

1. tic;, C& z6 D+ `$ g) E( K1 x
   
2. s=0;
   ( \1 O# S% Y5 I) E  t2 J6 z
3. for x=0:1000# b) C$ ^2 g' o9 @9 m
   
4.   for y=0:100
   # _, }. p9 Q* x7 C* u+ ~% d
5.     for z=0:100* z6 z' J8 |3 R* m( r+ m( K
   
6.       s=s+xyz(x,y,z);
   % w8 a4 Z3 d. C) z# J- Y' x. h9 B3 }
7.     end: Y0 n' M6 a  b5 E6 @% R
   
8.   end
   2 q7 z0 u& E. I8 W8 h3 g& z\" N
9. end% r% m+ g. r& p2 M
   
10. s
    ) p! @, F* m5 d4 Z6 ^7 R: |
11. toc
    : M- e1 j1 D9 V! m; V& J# t
12. 
    9 {4 h/ K/ E6 D9 b$ K
13. s =2 B0 l, ?1 t9 F. k6 a' t5 d/ |
    
14.     6.126720600000000e+009) G& }- F  ~) t7 J
    
15. 
    % ~+ B\" m5 [) A
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
. I# @* y8 `3 n8 G( U" G% N
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   & a) g2 n- j+ R+ m! F
2. mvar:
   6 F' Y1 p/ {) s, o' e+ j
3. t=sys::clock(),- T( T2 b2 y: W1 n1 h; r* s6 O
   
4. s=0,
   2 U* K! T/ ^) Q0 Q0 p
5. x=0, while{x<=1000,6 v2 Q5 ~: o\" X7 V: W7 L
   
6.   y=0, while{y<=100,6 p/ H9 @5 Z1 A0 v
   
7.     z=0, while{z<=100,
   6 [4 ]5 a\" V& S, \# N
8.       s=s+f(x,y,z),3 D' \6 A! Q7 n( @3 `% `: f
   
9.       z++2 c5 F# T3 X6 S- [
   
10.     },7 g) P- p* d# r, i+ t
    
11.     y++$ k7 W1 ~9 S2 \/ n6 _7 N\" n
    
12.   },1 V& h- E4 L% g
    
13.   x++! r6 Z( _- U5 [) S2 `( q  L' R; e
    
14. },4 l. G$ e* e; `4 @( B8 P1 A8 Y
    
15. s;
    & D4 C% V. y5 c: K
16. [sys::clock()-t]/1000;

复制代码

结果：
1 w& ^) S* m( {6 ]$ A! d  i- c% }  ]6126720600.
3 z4 y1 X* e; j& m8 z  r8 _1.109  秒
: `+ c: x: E$ H* l7 S7 d: _" q
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
) @5 j( P. o& O2 X0 f; d; ]" P! v

hxftntA

: v% _1 F% l, ]. M7 s7 J
. i( _6 [% D( P1 s6 h8 W
# @& R( {& q, c
) l# E. b7 O* S& A5 x* V







9 r2 d9 U8 }) Z


/ k2 Q; B7 Y0 J& y/ p  n7 B