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

forcal

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

! P0 b; Z) @1 V, b+ n1 P1 S1 b＝＝＝＝＝＝＝＝＝＝＝＝
9 I, ~# c$ W. p  w2 [5 g
+ K1 o* B9 y1 p& s- VMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;1 i# r! G/ k9 N9 \0 ]* e& b, x
   
2. tic;; h\" G% S1 N- y
   
3. s=0;
     M& X7 R\" F& ]- |\" l. u& d; R
4. for x=0:1000
   / K+ X, J) M6 d- ^
5.   for y=0:1007 X/ h\" q1 O' N( v. m2 i
   
6.     for z=0:100
   2 b# i, B3 D# i+ u
7.       s=s+f(x,y,z);
   0 f6 Z4 }$ Z: c9 F1 X& e& [
8.     end! S, ]: q7 Z! g9 }5 O% c( ^* _, O2 n\" Y
   
9.   end* H  H0 D  C6 l7 D
   
10. end4 Z* c# c! V- Q4 m# a% X5 v4 C
    
11. s# ^5 m) ?* ~7 M
    
12. toc2 y% {3 j9 D/ ~9 M- ^8 e5 x
    
13. 
    ! ~* v+ g0 b1 L+ o
14. s =1 }+ N) N' X. B0 @0 N
    
15.     6.126720600000000e+009: W! G6 z! v4 }0 J
    
16. & K3 b/ ]) d( m
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
% u: ?* H- I# x2 [m函数：

1. %file xyz.m
   0 y\" l: b  N) I' P# x3 _
2. function c=xyz(x,y,z)
   0 R% T. {& U9 p; o. v9 k
3. c=x+y+z;, Z, @1 X7 N- H  ~8 G% J2 [
   
4. end

复制代码

测试代码：

1. tic;
   7 H; I7 N% O- ]9 Y
2. s=0;
   ( A% [. e8 O, Q! h
3. for x=0:1000
   6 `, O- y- {4 b# Z( [, r
4.   for y=0:100
   0 t' s9 `! F3 n! [
5.     for z=0:100
   # t: o0 O, v( k7 Y+ |* |) L4 j
6.       s=s+xyz(x,y,z);
   ; T/ e, ?* r) ]6 ]- V; h
7.     end; @% W! m* M/ l; U2 J& }* f
   
8.   end
   + x$ i. @+ z, i/ `+ ~8 W, V8 P; {
9. end' k' P: S. S, o
   
10. s
    3 w' m4 M3 m% D
11. toc# \% i6 ]' @1 U5 J2 D/ ]+ o/ M
    
12. 
    2 ]4 w- u+ W* }1 o8 n9 O  B. |& x
13. s =
    0 x7 [( y9 ~\" @1 }5 _* m2 E: r4 @
14.     6.126720600000000e+0099 r% A1 L7 [\" i! z2 J' Q
    
15. 
    5 h( S$ [9 R+ Y! m  b5 c5 G0 g
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

7 h% W; i5 Q0 E% q; e: G5 |Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;9 o% B) y1 i* {0 a' {% L+ m# n) y
   
2. mvar:9 S( o. S& y5 S' v& Z, R3 a9 D* g
   
3. t=sys::clock(),
   6 F5 c' A, U3 b' l0 c( V
4. s=0,
   \" `- {4 Y8 Y) B3 W' ^
5. x=0, while{x<=1000,  [\" o, R9 Z! g+ N
   
6.   y=0, while{y<=100,+ n9 e1 s( u! R' _  e4 Z# f
   
7.     z=0, while{z<=100,- t) ^3 @9 ~3 N
   
8.       s=s+f(x,y,z),
   ; r# `# ^, L1 H\" g3 q1 m
9.       z++: Z  U\" s0 H! D1 u6 ]
   
10.     },: v# b3 j  y/ j! p
    
11.     y++3 C* J; p2 z: H( Q9 e. l3 N
    
12.   },
    ' L2 J! X9 w9 b/ [& K
13.   x++
    4 v& s1 v5 U% T8 n6 @( C
14. },# n' q' u4 Z1 }6 f' V
    
15. s;
      _9 W& W6 H& T$ z# a# X0 M
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
5 |! n. j- t' U% z# h1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
4 |1 m9 e2 @+ D3 E2 V) I

hxftntA

* l. U& p5 ~% [! ^9 `* p1 c$ \8 u' M

/ T: `5 k; s- U" Q
3 f1 E$ [+ s$ Z$ O0 E$ ~- o' q
6 D$ q: ^5 ^; C% S! {' ?) r

0 z  \* H" b1 w# V- S
' N% q+ _  {4 `. N- I, g8 Y! V


+ ]$ D1 u  O) K( v. Y
. F6 `4 Y& I, g* ~% x/ }, D
5 r4 P, t- O8 H4 l- J0 S! X6 d  }4 r. q
0 G# J& V& L6 p
( G7 ^& g6 U& |% W
9 p/ G8 C6 I7 j! M
2 ]* M$ x. G% }! ~9 A% {' R  Z