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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
# P7 e8 h' g1 @6 @0 g! O, P1 j. c
＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   # c! t5 p6 t5 j+ Z! {, |  K0 }
2. tic;( t3 ^6 s* d, Q4 l9 d% |\" k5 F
   
3. s=0;2 g# J+ _! Y* R! S% c% n: d$ O
   
4. for x=0:10003 q, B0 ]8 }) z  e
   
5.   for y=0:100
   ; e7 A* R0 `9 T
6.     for z=0:1003 p' X5 D/ m3 U& a
   
7.       s=s+f(x,y,z);
   1 E3 U5 l5 e) w' o% s0 b. y4 O+ x
8.     end9 f  ?5 U8 \6 Z9 B) ]6 i+ p
   
9.   end$ q8 y2 w. g  `. u7 X0 W
   
10. end
    , l, E- |% z- I; A# y' h6 j
11. s6 d) v3 s& i5 N) K* L* Q& Z
    
12. toc) Y5 o  n% d4 u7 C
    
13. 6 E9 _# j! X, u7 o& u% a! K- b
    
14. s =# Y2 f8 I* m5 \, _+ @
    
15.     6.126720600000000e+009. s9 }( I% x* Z2 ?, Y
    
16. 0 t, i3 i$ T\" |* V  B3 Q
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
0 x4 E9 _7 w9 {. `: b/ pm函数：

1. %file xyz.m
   2 R  C2 A0 m5 I; \0 \
2. function c=xyz(x,y,z)- C' R, R- Y# g0 [
   
3. c=x+y+z;. n; T% ^2 d0 a6 L% j
   
4. end

复制代码

测试代码：

1. tic;- W& J6 j0 F' M
   
2. s=0;/ q+ Q2 [- s5 ?\" {3 ?  L, B
   
3. for x=0:10009 I& R5 n! z0 d, h5 T% ?/ O
   
4.   for y=0:100
   \" G. c  R5 z2 `2 p2 }  g
5.     for z=0:100% T) k  }, W; N3 g* ^3 ?# a
   
6.       s=s+xyz(x,y,z);( a\" j: |* d* W4 U: E
   
7.     end, K; P* E\" H  `: e  o1 z0 e+ Y) n
   
8.   end
   ( V/ _3 c$ y1 |
9. end
   3 V7 |8 Q- `1 B; r- k' v3 m
10. s
    8 w5 H! O& s+ Z& ^1 ?
11. toc4 J( x* V7 x% l- j! x/ Y7 }
    
12. ! R9 d! {1 x) ?1 R( S
    
13. s =- p) ?& a+ r5 U- v8 O* A; s  x8 @. d* W
    
14.     6.126720600000000e+009
    9 P9 g! b; H; e* L: }8 U0 s& v
15. - R( e1 a; f- ^& O* a8 J  n
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
; e" ^! C% E8 k. X- R7 u
( X9 V+ Q( C9 D0 c  Y% YForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;1 w: f) G! I8 G3 K
   
2. mvar:
     P6 I% a7 _5 e( x
3. t=sys::clock(),8 f/ H% e8 W# Z2 P
   
4. s=0,& y# o\" @% n8 t; n
   
5. x=0, while{x<=1000,
   # E  R; x  W) ]5 c6 ]2 |. ?
6.   y=0, while{y<=100,) X9 o% I; m0 ^# q7 P5 A
   
7.     z=0, while{z<=100,' d& x$ g  b9 k: U
   
8.       s=s+f(x,y,z),3 @  T2 d# Y6 e2 q. w\" o- w5 L6 n
   
9.       z++
   4 d* W6 I1 o5 c9 J- E8 t2 ?
10.     },8 r1 }\" L5 W- f) s, z$ O
    
11.     y++
    - C% ]1 E8 ?# |* ~4 S: j# s% e
12.   },. R3 N& ]4 ^7 G0 T: u. n
    
13.   x++
    : D* X, S- g3 q, [7 K$ G5 {( L
14. },  G& j1 O$ g\" g# U# D
    
15. s;
    8 Q' N! S  j0 M! f3 ~6 ?( d5 O; ^
16. [sys::clock()-t]/1000;

复制代码

结果：
7 `5 ~4 b5 \7 K8 M7 o6126720600.
2 R  ]8 \- K: D5 F1 J- Z1.109  秒
" G% h, }, Z2 S' Z+ K
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
, o0 \( o/ y* A; ?3 |" Z% N7 u

hxftntA

* {4 q6 \/ h, ~  g. G

0 z6 \3 I( Q. Q3 j& d* n

& J; J4 t' [5 N0 L8 O
: P, J& I$ p; A$ n
2 H$ T& m) B/ A7 A4 y
! u3 S7 K4 S' s1 k! ]



4 _$ `  e; l" o, w4 _
! r6 w& S* d( O: k0 A

! w- g. k8 y# U+ _' V. ?1 J