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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
' r$ S' p8 n4 l; T
6 x! n+ P& T  G" G$ n6 ~7 w＝＝＝＝＝＝＝＝＝＝＝＝
' d" w4 h. I3 J6 B+ A: c
. V4 a! O/ Z- ~  \Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;5 O; ~1 Q- E/ K: \$ P0 E
   
2. tic;2 |1 k' L& m( i( y
   
3. s=0;\" i: }  n! s$ @8 V$ j7 @9 H
   
4. for x=0:10000 U! @, \# x( t; I$ v+ N
   
5.   for y=0:100: S5 z, d: C# A
   
6.     for z=0:100
   - {( X* v. c/ H' T7 t; C; B. A6 |
7.       s=s+f(x,y,z);
   + h- c( q; O3 s! P\" I! U
8.     end! {7 z: g\" N) h, U( Q
   
9.   end
   6 m# z) H  M4 ^+ F- E3 |
10. end) v1 |6 |5 d8 {# o6 ~/ k6 B  R
    
11. s5 ^9 Z3 l- {! O, l% Z
    
12. toc, U6 [' H# M5 o4 O( Y8 e8 c* A
    
13. 
    - _5 ~5 R, U# x9 K
14. s =
    : b& |8 T, n* b/ j\" l8 U5 J
15.     6.126720600000000e+009* Q6 j9 B( y% l- @9 f% O
    
16. ; I2 g0 l2 k5 U
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   ( Z: b+ r: P, w* C# H
2. function c=xyz(x,y,z)
   * U% h  s' m4 t! m0 j
3. c=x+y+z;
   ) h2 R* u, f4 R, C/ ^$ M
4. end

复制代码

测试代码：

1. tic;. h9 _6 j2 w4 p( a) F$ W. A# [5 U
   
2. s=0;2 u$ f* @0 m; m8 m4 m- P$ B/ `
   
3. for x=0:1000* y% q; b) d. K' q/ m5 s' }0 {1 N, n
   
4.   for y=0:100; c/ s+ Y$ @; p
   
5.     for z=0:100
   9 y0 A\" H) W( {9 C. O
6.       s=s+xyz(x,y,z);
   ' {1 j6 N, ~# p7 c- k9 C/ a+ D
7.     end
   9 G8 h% Q0 V7 I/ B9 q6 ^
8.   end: k; N3 l% D9 {* {/ P, x9 g
   
9. end' Y' c4 ]* x* @% r& j
   
10. s
    % k: Z  G. k5 ~, v3 x
11. toc: m2 G$ G' ?4 e4 e! L
    
12. 4 P\" q! ^3 O% G4 A; i! x: q2 ^
    
13. s =% G) F. R/ [+ k& M0 i% v- m
    
14.     6.126720600000000e+009
    ( F+ p( S4 j\" y0 T1 f1 l
15. 
    , a6 U. W: r6 K2 {* K, |/ X9 N
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   9 H' _* Z, i0 W* i6 v$ h8 E2 Y: V
2. mvar:* r- l% T: ]$ N# }0 t
   
3. t=sys::clock(),# k' ^; x7 V6 _+ k
   
4. s=0,
   0 {- T* _. l5 ]3 x# z9 x
5. x=0, while{x<=1000,8 J3 E/ ?% f2 A
   
6.   y=0, while{y<=100,
   \" x2 b% y* ]* A\" u* |9 W5 L2 X
7.     z=0, while{z<=100,, f+ I5 ^3 K; }3 X
   
8.       s=s+f(x,y,z),
   . ~. E+ x: Q, u
9.       z++% u6 G5 M, _9 N* Q8 v! G' |
   
10.     },: i3 x( r$ a/ ^( g7 f
    
11.     y++# M, L  T# }7 E8 y5 N8 g0 T
    
12.   },
    - K7 o; r; Y1 m* x& k
13.   x++
    ! D. a% P4 f9 i5 l6 R- A2 z
14. },
    3 s1 T1 n! L6 h/ r; c# V. u8 V
15. s;
    , n5 @. Z1 X6 b5 H9 w
16. [sys::clock()-t]/1000;

复制代码

结果：
( Z& V% e" q! Z5 v# W; Q; @) I6126720600.
, p/ O0 S2 m0 O) a" @+ B: {1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

4 K/ K5 D* u; W7 m( Z

" `9 ^2 \3 E& d2 F$ y
* g0 o) m5 U: r" e. d6 s+ |+ B6 \+ x7 E


4 E: ~/ @1 u0 p3 l5 h0 k


) y# i) {- n* ]1 O



$ I. `5 {1 k* s, n, A5 u