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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
/ P9 \0 ^8 ?3 G  m  b7 b* C6 w
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;. c/ ^+ v/ l$ Z& P& W# e8 I
   
2. tic;
   . A6 Y' g* w6 v* e0 n; r
3. s=0;7 U. j9 P8 [  n9 ]9 T4 @
   
4. for x=0:1000
   ! t3 V+ I1 o! a- W5 S' R* h/ O8 A
5.   for y=0:100
   / n/ @1 o7 T4 Y; |: a4 ?0 ~
6.     for z=0:1005 O, q: L  y' h
   
7.       s=s+f(x,y,z);
   $ e$ I\" Y% Z8 m. n, Q
8.     end8 a- W! o/ y4 d3 ]* d* B\" X
   
9.   end
   1 h0 Q6 _4 t3 I# o1 m& g
10. end
    : f7 [4 Q- X' r, j
11. s
    7 R% \* y7 D% c5 o% Z; I: c; u
12. toc: b\" h' B- K9 b5 v, g/ Z4 n, I
    
13. 2 C  G- K: c- E' |: P+ I8 r
    
14. s =; x: E$ O: q( O
    
15.     6.126720600000000e+009
    1 h' i. h) |' u2 O, W4 C
16. 
    : _2 ]7 V+ w: u# L3 ?
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m( k, I9 |8 i( @, h/ N
   
2. function c=xyz(x,y,z)
   * H6 x7 ^7 e) p+ C% k
3. c=x+y+z;
   4 G* C/ V! h7 m8 d
4. end

复制代码

测试代码：

1. tic;
   % u4 g5 |. r; L: _5 _  Z$ \5 i
2. s=0;4 f2 t3 K! z% D, _7 t' P6 e, Q' C
   
3. for x=0:10007 ^% {. {6 t+ e0 e' z* o
   
4.   for y=0:100. D0 U& m. v$ A; w) N! g2 H\" u; J# n
   
5.     for z=0:100
   - ]( x* \$ l& c# Z, T
6.       s=s+xyz(x,y,z);5 i, f$ N3 d& A! C2 P6 G8 e
   
7.     end# {9 w: |( ^* ~3 k- h
   
8.   end
   6 i  Z+ _+ b  b+ R, P
9. end! }& q, U4 s2 b3 B) P. v, ~9 C
   
10. s+ L1 W8 n* Z6 ~) ^
    
11. toc
    ; g- k1 r7 @' E8 D' z
12. 
    & k) `4 ?9 Y3 p; U5 T
13. s =
    3 ^- c' g% T0 T' z
14.     6.126720600000000e+009/ H$ M( J4 @- H( y1 P& o% }
    
15. / n$ Z# @4 N4 w: H' I0 D+ |+ w, j
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
7 D: X& n2 ^: o( c! {+ c
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;( I) ^5 z5 p( J9 Y
   
2. mvar:/ c* t' g; {% [: z9 |* G' s
   
3. t=sys::clock(),
   6 y4 H& O8 V+ ]' Y7 d+ T. V
4. s=0,' Z2 h( [6 D. P
   
5. x=0, while{x<=1000,9 d6 v& {$ i. J' B
   
6.   y=0, while{y<=100,
   : h, T# m6 x. i0 U1 Z
7.     z=0, while{z<=100,
   . }6 F/ R  h, s8 h5 H- g; [
8.       s=s+f(x,y,z),
   . L3 `8 }. P: |, E\" w
9.       z++
   ; L4 J6 E8 d# W& b$ m  U\" _& J
10.     },
    5 h  R9 a, z6 R9 r
11.     y++  C- o; F+ T6 Q. v8 w
    
12.   },
    1 ?- C1 q* U1 y2 J. n6 ]5 j& Z
13.   x++7 V3 p# x2 ?+ j9 T4 J
    
14. },
    , A7 J) o/ d' T
15. s;$ O4 i0 K6 W, r\" q
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
# e. P) K, S/ M+ f4 F6 J# }' d1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
( C% F/ ~! \  D) J6 B: {" \# ]- ]: w

hxftntA

6 G( K: n* K! A" V* M' B) I

2 v  _' V( |& d0 I8 `

- e, M6 ]% d1 Y7 A8 Z3 ^2 g/ X0 T


% {$ d$ N( g8 _
# T  V/ i8 {* ~2 k8 a7 s  C$ c5 t
; D3 @7 h. T0 ]( f& ~# S2 |
+ e7 U% ]' Y) E( e. ^) ~

9 h+ w1 G( N( N6 Q" m' H2 x
) V& h0 B8 A' W+ f/ l