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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

  A# V/ W# g& dMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   6 L) S8 f* v1 b7 c2 y  @
2. tic;
   9 s\" N- D# [1 W
3. s=0;
   6 M. D3 S  M3 Z\" Y0 |
4. for x=0:1000  c6 M5 X$ q9 [- T5 B, s
   
5.   for y=0:100
   + W0 [8 Z\" [5 p3 t7 ]7 m2 \- f
6.     for z=0:100
   8 P/ P& f  X3 ?' b9 c
7.       s=s+f(x,y,z);
   : l3 ~3 E\" K' a3 _1 f7 X
8.     end
   9 C; i; V4 H% P2 `) B
9.   end2 z4 m- j0 q5 S( C2 W9 K4 ^3 H
   
10. end
    - X7 }  }! d6 _
11. s
    / M4 G* j7 v+ T' A  \+ @
12. toc
    ) I5 n! N/ y\" [1 [\" A' u% W9 e& K
13. $ z. z1 l. a3 Y- G
    
14. s =
    8 B: e5 n% S! d! U$ O
15.     6.126720600000000e+0098 H; z  T- }5 ^: v4 o
    
16. 
      P! v\" S6 D* g+ v
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m6 H6 |- f; T  Q, F
   
2. function c=xyz(x,y,z)
   ! s: H1 m% g( G1 V, ]; `/ J7 E2 `
3. c=x+y+z;
   0 |' G4 A4 E; d- L4 o5 j
4. end

复制代码

测试代码：

1. tic;4 L9 {/ M0 Q7 c\" `6 J\" E4 g/ T
   
2. s=0;: K& g6 h; G9 f/ u0 }5 N: {# j
   
3. for x=0:1000# k/ L, |( I2 X1 y- I
   
4.   for y=0:100
   : U+ D/ K3 y7 j& y; }; @1 C6 T' a+ C+ Y
5.     for z=0:100  k* M3 H, }0 [4 ?
   
6.       s=s+xyz(x,y,z);
   , C) Q) I+ s. M) N3 v1 @5 b3 I3 F9 k
7.     end0 K; `0 b  p) S/ T5 ?% Q) k
   
8.   end' \% n\" U0 l# m9 L$ L
   
9. end* a0 o/ |6 I+ L2 }\" U
   
10. s
    2 ~7 H1 e+ `/ z! e+ B7 h
11. toc* N! T\" `3 k3 N2 _) D# W9 K
    
12. ! J$ G* K. R$ z0 ~1 T( r\" I+ `4 Q; W
    
13. s =
    / Y2 h2 h5 i2 K+ W- i\" O# N. n) i
14.     6.126720600000000e+009
    ! e% e8 S# W0 ~
15. ) r/ j$ e1 m7 a( r\" a) q& n2 ?2 Z
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

" h# j  L5 U! e  Z3 VForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;$ k% W0 R3 W! o4 W- Q  Q1 i! W
   
2. mvar:
   ( h7 i3 K! f) |/ f. s
3. t=sys::clock(),
   ! R- F5 I% `2 q/ y
4. s=0,4 H1 p( W+ d' R- E% g
   
5. x=0, while{x<=1000,- g' W4 e; j4 W. }
   
6.   y=0, while{y<=100,% I\" T6 S) w0 H2 c% ^0 @
   
7.     z=0, while{z<=100,
   9 I* D0 {2 x0 C\" _; V
8.       s=s+f(x,y,z),2 F$ _( p- Y5 d( q! I0 Q
   
9.       z++$ L9 c: T+ n; |' M. _* L* x( f
   
10.     },: ]$ n' p& l\" ]2 L
    
11.     y++
    : Y, t; h0 g. d7 }3 e\" g
12.   },5 n$ z5 \( P& u5 m# s
    
13.   x++
    5 B7 ~3 J' t7 Y
14. },! {, @3 K3 t4 _2 Q3 Z2 p3 x
    
15. s;; T\" a/ k6 E+ A6 t' d
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

. L/ o& l( G0 h' x3 e9 f4 P; S+ ^二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
2 L4 i1 z5 w. J. q7 k+ z

hxftntA

0 b: g7 `+ Z$ g) q4 `4 z3 c

6 J4 @" {* w- l6 l4 v% X+ \. D

! v) j  ]5 ?- y1 g
* O) n$ r+ l3 n7 z, r
3 \  u6 @# p$ [" r  [/ x' e) Y7 |
* z% y0 r0 |, N8 ~4 z

& O: ~4 z1 ?+ f: c9 q8 C' q


6 i0 \4 F4 ~5 v
0 a6 a' q( x- i: ]
6 T$ Z, s4 }+ O; A. u* @4 r) Z; o, @
8 \8 {! P% b" `: [5 t' V7 @8 g