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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
0 Y- U  X7 _. X2 X
＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;5 y3 x$ O* U4 o. q5 M
   
2. tic;
   2 M8 k4 O: c# v! v( g1 k& I, c
3. s=0;
   ' c1 U1 _  C& O8 Q% n! G1 Q: T
4. for x=0:1000
   4 G, J( |$ G7 T# y1 x5 R
5.   for y=0:100/ h( G1 X  `2 N6 }: ~
   
6.     for z=0:100
   & `0 ]+ `, p: ?  A
7.       s=s+f(x,y,z);
   / D' O  p* ~0 I  N3 B3 O( ?
8.     end& N0 _, ~+ n# f2 N4 k& D\" w7 V
   
9.   end& b# }2 H& T0 C8 X, U* X) s
   
10. end
    $ k  H% ]: n+ V3 O* `2 W$ N6 g5 g3 N
11. s
    ) ^* H0 H\" C3 J
12. toc
    8 c+ L; Q1 d$ O
13. 
    % q. _: q- Q% Y
14. s =
    # O# ~6 c% \& a: W9 {
15.     6.126720600000000e+009
    ( e/ [3 K% k7 l; ]. U9 g# P+ a
16. 9 o) ]* b5 O( g: y7 R  d% B& r
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m9 X6 {8 n1 N) z+ m8 b1 n  s; f! I* n
   
2. function c=xyz(x,y,z)
   6 E4 N- I- i  O2 y
3. c=x+y+z;/ _# C  T9 ?* c' Y8 X
   
4. end

复制代码

测试代码：

1. tic;
   , S  \! y& ~+ w: ^
2. s=0;+ x& J9 J' P* ^6 r
   
3. for x=0:1000
   - r( O; N2 D  D& y
4.   for y=0:1009 ~9 n( `5 H+ m# n. i5 ^# K
   
5.     for z=0:100
   . ]' J) Y6 Y( u0 F8 i
6.       s=s+xyz(x,y,z);
   . K% @2 [6 C\" r7 m9 f1 N/ A
7.     end
   ' b8 }' C8 ]% L  }3 k) w% I( o
8.   end
   0 ?/ [: ?5 C, C\" S& ]
9. end
   % H# V\" B6 m) \, P- p. B$ ^( J
10. s
    / F' w& [9 l0 g. K% ~$ q; [* y
11. toc) z  G) k' ]4 z: C) M& h, V
    
12.   w\" k' s# i' `) K+ D7 x
    
13. s =
    ; {: j3 @, M6 [6 m\" s0 [! _
14.     6.126720600000000e+0096 y+ p, h3 v/ U8 Z\" i# D+ T0 f
    
15.   m! O5 z/ l) A7 C- ~
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
! q& h& c4 E; J5 c% e# W
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;/ K9 \9 a+ t/ p7 v7 t
   
2. mvar:6 U! I- l) J. \# m5 l# @
   
3. t=sys::clock(),# I1 V+ T) I0 K( R
   
4. s=0,
   ) m' i* D\" Y( z
5. x=0, while{x<=1000,
   % U\" ~( G3 T6 t' b
6.   y=0, while{y<=100,: e+ T, R9 r$ t
   
7.     z=0, while{z<=100,
   \" a  `5 T2 X& d- D; r. T
8.       s=s+f(x,y,z),
   , x3 k, f4 m9 H2 k
9.       z++# [\" F$ L8 r& }# q6 e& i
   
10.     },
    9 ?. Z- S7 w  n- _! H7 K7 C
11.     y++
    ' Q9 |( s\" [9 L4 O1 E+ N6 E
12.   },2 r/ f5 B7 x\" }$ |; ^+ x- v; X
    
13.   x++, L) U* y& ~2 e) p
    
14. },: P6 |- [8 y& _8 F\" y
    
15. s;$ w% x/ C7 f3 s) k, S) k0 a$ J
    
16. [sys::clock()-t]/1000;

复制代码

结果：
% j, Q6 H$ w0 a! X% e6126720600.
1.109  秒

! ^0 @8 E: j: J' D二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

4 D4 {& @3 q  V$ `* h( K6 q4 L. l+ H1 S
* z% C! j3 C6 e' q0 ?

& i2 p5 A# `! }' B# K1 T; X

. R* h  C) _% N8 x2 `

# @' j1 S7 W0 C2 P4 i


  e) l+ ]7 X6 q( z! u