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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
1 g) ], u* {$ P" j( b5 @
+ B8 E/ t) n! C. B' u! ?9 m＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;0 c9 C$ G/ N1 ?
   
2. tic;
   & s/ p2 Q5 ]2 o3 U
3. s=0;
   % U, {, q2 J$ ^3 _( F
4. for x=0:1000
   \" O1 h, l; U) G% p! \0 a+ K
5.   for y=0:100
   7 O; X& ?- z2 ^9 \. c: C0 V5 B
6.     for z=0:100) w, Y5 G2 }+ m( S
   
7.       s=s+f(x,y,z);
   % ^# h# V3 Y+ i' N
8.     end
   7 R+ }/ L/ S* ^# }6 E) u
9.   end9 e  P. V7 Q; _/ k  a2 ^5 E
   
10. end
    . x8 z& M, |4 Q% o* Y
11. s
    4 E0 {% l+ d8 `; B5 z
12. toc$ T% I0 Q' |9 c( q
    
13. 
    ( K7 F\" ~7 X9 }: s/ N9 C
14. s =
    / E2 x/ R6 L\" V7 B3 d. Z1 ~
15.     6.126720600000000e+009
    ! w* j& t- [( G. j# \. X7 ]. E
16. % e% s+ v; G: [4 W) M2 f
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m: ?* N7 U3 [( W  \* q1 V- ?
   
2. function c=xyz(x,y,z)
   - H0 P8 @7 B& s
3. c=x+y+z;
   9 P9 z1 x9 }0 u4 z2 }
4. end

复制代码

测试代码：

1. tic;
   4 A' }1 J, F/ C
2. s=0;3 ]2 @- ~  d. U
   
3. for x=0:1000
   ( |0 r' C- b/ K- a2 o) w
4.   for y=0:100
   2 y: `4 L0 [( D$ x9 @+ S
5.     for z=0:100- Y2 I: d7 n) ^9 e: n7 H, }0 T
   
6.       s=s+xyz(x,y,z);  a' E5 t: Y) ^: ~& u
   
7.     end
   7 A) D/ [, c+ i& c8 c
8.   end
   ' d0 |+ I2 A# u9 a\" o\" a
9. end/ u! s8 f. {) g
   
10. s+ Y; n8 h+ Z% Z; o' i3 P, N8 s
    
11. toc! z\" F8 J: L) G
    
12. 
    , c0 C; z' J! D7 }
13. s =+ }# G$ @\" w% Z1 f; g) S5 O- J
    
14.     6.126720600000000e+009
    4 X! r3 a( p1 J
15. 
    % D. d7 ]6 E2 V. I- C\" N
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

3 P7 U/ M0 G- O4 ~, e% \Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;4 X7 V, p7 |0 ^
   
2. mvar:
   : W5 f1 H1 a9 T/ S
3. t=sys::clock(),
   ; o& q0 V: s/ h
4. s=0,! s# J+ i/ w8 C/ y1 |  Y! F/ d
   
5. x=0, while{x<=1000,$ L5 v  y\" i8 s# k% j
   
6.   y=0, while{y<=100,
   + B3 c# {# t, Q0 p
7.     z=0, while{z<=100,
   # l1 @8 T9 X% ?' d0 x
8.       s=s+f(x,y,z),  f\" i) o+ Z. Z
   
9.       z++5 U4 s+ ^. J, V, _2 D, V
   
10.     },- x5 h# T  E* H8 K3 c
    
11.     y++
    # ^4 x* A1 v0 n8 e9 ]- P0 ~
12.   },0 ?  @# h0 n7 W6 W0 v. _
    
13.   x++
    5 a- J8 y' [1 p: u* [9 [9 C7 }
14. },- S/ s1 I7 A  }* d8 O
    
15. s;
    * p\" q4 Y  L9 @% ]1 U
16. [sys::clock()-t]/1000;

复制代码

结果：
1 V" e1 G  D) L( r: O4 A: K  K2 t6126720600.
1.109  秒
$ j: d- a! e% `3 Q
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

" A% f, d' U& i' c9 k) P& a6 Y6 e


; w' b; L& G8 d+ L4 V* |


3 P: Q# I2 y" _! F( A

1 u7 e) L, Z( L  f6 _

0 L! f" [) k# s
/ L# R% ^: ]# X


, s' g; H/ v. b5 k7 T- K
, K' O; l' x* @; P5 i4 o
  k/ M  P3 P2 M0 D