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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
0 I. z6 {% X# \. x! ?! D
＝＝＝＝＝＝＝＝＝＝＝＝
, X' A- {: x! e" T
+ C9 Y% i8 @1 }, _0 x, aMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;9 C( ]8 N7 ^1 ]- }3 h: C0 r) k
   
2. tic;
   9 i' M0 L' F. ?3 E) f( e' _. I
3. s=0;
   % G+ E. B6 ^9 u: K% y4 M! u0 ^9 [
4. for x=0:10006 E6 Q* s+ {6 g1 _- i) W2 s
   
5.   for y=0:100  e, T/ a- ]9 u4 @$ m5 W
   
6.     for z=0:100
   7 S\" P! G7 q) V* `. l$ P3 J
7.       s=s+f(x,y,z);
   / D; R1 o/ ~2 c, U* r/ y$ M2 U! a# a) n
8.     end! M, _% C; g3 {' N1 V% }- A
   
9.   end
   * [: t7 S  O& h, @& \- T6 k9 g
10. end! |. v+ O! l+ j. r3 `) G3 m7 \+ i
    
11. s
    ' i9 w1 R2 a/ g; A- I* I' F, T
12. toc
    9 |/ g4 s$ I$ Y6 d\" \/ v# a) J
13. 
    ) F# z  c5 K4 K3 g% L+ T
14. s =; {5 d0 t) r$ ~# G2 }9 T$ l
    
15.     6.126720600000000e+009
    5 g9 z7 y( T8 q8 G2 I
16. 
    : Z3 @7 i4 q; r- Z/ M2 k. }\" D
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
* I& ~  }& R- v# ^$ Y+ gm函数：

1. %file xyz.m, ?: _! p9 ?% X
   
2. function c=xyz(x,y,z)  Q- `% `, a( O) B
   
3. c=x+y+z;
   ! e7 l\" V. b/ Y
4. end

复制代码

测试代码：

1. tic;8 S' w) _% I! A
   
2. s=0;7 v) ~! Q- E0 y% s\" W7 z
   
3. for x=0:1000
   - J* y6 Z1 f7 G  h- n4 t2 Z\" l$ s& i
4.   for y=0:1001 \7 ~5 w0 `+ J# ]6 H1 U
   
5.     for z=0:100$ X2 C' d, A, t; T% I
   
6.       s=s+xyz(x,y,z);' @( F  |8 t3 f6 _7 G! B# }
   
7.     end
   7 p) _- X4 l: M
8.   end# z  F: W8 O; \2 F4 U
   
9. end
   * K* p0 f4 G5 n4 f: S4 \0 N4 V
10. s
      B3 _5 G$ M' G5 {
11. toc
    2 B! L( w2 x- n7 J, f% ]
12. ( p' U. f& d# \
    
13. s =: d; q( U' R6 h* Z( ~+ I! J
    
14.     6.126720600000000e+009
    / ]* X  N1 v2 n% {# M8 h+ R2 ^
15.   S' h5 r4 |3 {. q. ~0 k
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
/ C1 F; S4 Q& `
6 A2 r1 }: B3 ?7 ZForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
     |) L7 p8 F1 W$ G' D
2. mvar:
   9 ]) g\" `: w; A- o2 C
3. t=sys::clock(),* i! g* v6 W; L; o
   
4. s=0,) w+ R3 R6 B# O3 ^
   
5. x=0, while{x<=1000,
   6 x' K  r: n6 q9 O5 m
6.   y=0, while{y<=100,2 y4 O- c/ |) P2 s+ o# {: E
   
7.     z=0, while{z<=100,
   ! n* z0 G3 p\" ]# G' u- w4 Z4 }6 t+ h. m
8.       s=s+f(x,y,z),+ e! d, h, S8 g; h6 `2 b
   
9.       z++6 O+ w3 y: ^5 M6 s7 l1 e/ }
   
10.     },
    & t  N4 |, ]3 i- b6 ?8 L5 \2 v% \
11.     y++
    ' K- |& l0 O1 D
12.   },
    8 U) W- r; M; b- o2 Z: N* t
13.   x++
    \" ^& n4 E; c' I, v9 [* J/ z6 U
14. },- q. @9 G: T- F3 Q0 N
    
15. s;+ L, V2 Z1 c5 s4 f4 O/ Q) V0 |
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
! s3 E$ v2 P, |+ i+ J8 K0 W; H
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
1 v  `( h1 }8 ^) d' w/ d+ M

hxftntA

  T: ~; l0 M7 w5 K* B7 _9 S5 P
4 S& W* _9 T7 v
, q9 n. z0 \: E! P" e* U/ M( ?

: j" q$ v4 q8 Z2 A  w! P


; H! B5 \8 H3 R* s, q- h
; k) i" z0 u$ T: X8 j



* t( ~8 T$ @4 M" X  R5 l

* R2 p; I8 U; U7 Q. a