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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
+ ?( s1 T. B, {& T: w( ?% C/ R
) `: l* Q! W& F! r6 S＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;8 o. ?* m) F1 R9 j
   
2. tic;
     _* s2 c5 H: a1 h, x
3. s=0;
   * X0 t4 Z; W) ?+ e2 E
4. for x=0:1000
   / G6 {& B/ U/ ?0 {8 h+ V9 U
5.   for y=0:100  U0 l9 k; q% t* w* a( ~
   
6.     for z=0:1006 q' @; Z( b\" w0 V
   
7.       s=s+f(x,y,z);
   $ p: t1 ?; B, O, U4 S0 O7 d8 ?
8.     end/ v9 y; y5 c5 p+ d4 X0 e$ Y
   
9.   end
   ; d! A9 _, l, o
10. end
    # `; k0 M; f$ D6 _7 _
11. s* E1 j3 j/ N& e( H. R4 b
    
12. toc9 }) D( Y* y) O, e' a5 r1 \! C: a( C0 q
    
13. # Z\" P( c* e& O: L1 @\" C\" x% Z. l
    
14. s =
    3 S$ B. M* `& r: x& s
15.     6.126720600000000e+009: ~. x$ s$ W+ S0 ^3 N
    
16. 2 z% M  X, `* \8 i, a) V
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
! W5 E; E3 q$ Z) I' `m函数：

1. %file xyz.m7 L$ w\" Y/ w' g( h. V\" ^9 i
   
2. function c=xyz(x,y,z)
   # u9 }\" I! Y9 p# c0 i: y# V
3. c=x+y+z;& W( ]# \( k+ Z) j1 m8 Y1 _  s8 _
   
4. end

复制代码

测试代码：

1. tic;8 ?  Z8 _* [: k+ T- O( b4 p
   
2. s=0;
   \" ~$ l# U2 M7 R8 K$ n! _4 D
3. for x=0:1000( t; A- u4 G# Q4 O9 c2 j/ e
   
4.   for y=0:100
   ( q, _$ }9 C; F( x
5.     for z=0:100
   + x* F% `; ^# C. }
6.       s=s+xyz(x,y,z);/ V4 P' m  b. \6 q; h1 Y9 q& T
   
7.     end; W3 a/ n( s9 k, f% F. T
   
8.   end$ {# g% ?9 G( F\" b8 c3 v, d
   
9. end7 n. q; z/ S0 t\" q
   
10. s
    1 L& k( \\" _; G  ~! Q
11. toc
    0 ^3 H3 c0 d& a% a! m
12. 
    , d' I5 h8 M1 P, J  n
13. s =
    0 ~# j- n9 i2 [- k
14.     6.126720600000000e+009
    , Q9 ]1 y. C& h# D! z* O6 A! H3 }2 I* U\" @
15. - a5 }* W' g2 I' N( O
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
, z7 W2 l* T, W, t/ L" @: w! Z- z
" |5 c# f: y9 R( f1 d1 G& @2 wForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   1 k\" u9 M% V; N6 b
2. mvar:% c: U/ W1 {  |: o' E
   
3. t=sys::clock(),- }3 Z7 y% P5 C1 H\" k  [) {0 O2 z2 m
   
4. s=0,) |5 x& a6 q# ]. `8 z7 p9 O( @
   
5. x=0, while{x<=1000,
   - [\" M# r$ z9 U  ^) b- n\" j* j. l
6.   y=0, while{y<=100,
   7 k/ L8 x  a+ E8 S4 w/ h# T
7.     z=0, while{z<=100,  `1 Y( s) }* P7 d8 ~! X% c
   
8.       s=s+f(x,y,z),
   ( ^* f+ {0 ~. I2 R& S
9.       z++- v8 I% j2 N( @2 |6 S\" F5 ?
   
10.     },
    # I% W& i8 g8 \$ C4 j
11.     y++3 L3 O! I; T+ q
    
12.   },+ r& m! p7 Y# k% M
    
13.   x++1 k\" Q9 G8 Z2 p4 h- G
    
14. },# h1 U4 n! `+ j. {\" `7 x
    
15. s;) R2 o\" E# Y. m( e0 o8 D
    
16. [sys::clock()-t]/1000;

复制代码

结果：
0 ?$ h  F7 O/ B6 q: q  S& _6 O6126720600.
1.109  秒
/ I& C- z4 V1 p2 @! C  Z, r8 i
6 _5 V: D- C: `( |二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
$ Q7 J3 m* }3 s2 W  ~  `

hxftntA

$ j0 W& S, T# B$ K8 P4 ^% t




8 u6 B$ C  p* h* a5 h

' D- V4 F" c, X- Y
" K+ m% J& _$ v* U- z- ?) ^
9 O8 A: A: q1 }1 C+ `( ~# C

; ]# _& U  F0 A

. E1 Q- p5 N8 o+ y# h
( ]: `) B9 p1 X  a7 b