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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
3 i1 M: w" z4 v$ e2 @" `+ A
7 p$ h6 E; U+ e4 T$ [Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;2 h; E# O4 b/ T1 E4 O6 M$ G& F
   
2. tic;; Q+ T7 H\" `- I& G
   
3. s=0;
   # y5 o. M9 G- p+ _
4. for x=0:1000
   ( h! M$ C( e* S! j2 T9 d
5.   for y=0:1001 m$ G0 X% ?- }3 F0 p/ \; F
   
6.     for z=0:100
   2 C- s% U\" N# k\" g3 D0 i1 l  H
7.       s=s+f(x,y,z);
   ' O1 _! [' ^) r2 t
8.     end1 K* ~2 h7 E/ p1 M
   
9.   end
   ) p6 K+ L- V4 f2 o
10. end; _  z  `6 S$ Q
    
11. s6 k1 ?1 A2 h' s+ O; y
    
12. toc: ~, Q* A1 X; Z: r7 f1 q$ {
    
13. ; v/ o/ K\" P4 z/ v1 J! _
    
14. s =1 `. {; e/ L+ A- o
    
15.     6.126720600000000e+0096 B2 j( ^( s% F% o) i
    
16. 
    $ R( L% E6 F) \7 B% F0 W( G
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
& n& f9 U  K0 R4 Tm函数：

1. %file xyz.m1 I' f# a; ]7 h& D
   
2. function c=xyz(x,y,z)4 W( J0 K\" u$ U
   
3. c=x+y+z;6 E7 U' _\" h/ G$ c3 O
   
4. end

复制代码

测试代码：

1. tic;; t& {! Q, n7 k* U\" P) o) x6 G( v+ x& W2 q
   
2. s=0;! Z% H# ^% W! U
   
3. for x=0:1000! h% M: i' p8 B% u\" E% I
   
4.   for y=0:1009 d9 H\" ?$ X0 v5 p  V  ~2 }& j
   
5.     for z=0:100
   2 f7 L# `9 H( N5 W\" S6 j; Y
6.       s=s+xyz(x,y,z);  E) u8 g& s. W\" J0 J
   
7.     end
   : c\" `# V* z( p; C9 L
8.   end- d, ~\" [, u. `6 g
   
9. end
   \" ^$ l8 {* p' L1 f$ C0 _; i* K) D
10. s+ y1 f& s# k6 L' b
    
11. toc  B) l9 \/ \2 `' U$ R; m- g\" H
    
12. 6 e. H( u7 B6 G3 w3 m( W
    
13. s =
    0 ?( ?' {& S: W% b3 M4 j. i& @
14.     6.126720600000000e+009
    7 f3 x2 s% k$ ]7 }! U
15. 6 D( D\" U7 U. p8 B) a& q
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
& p/ J2 b1 W1 A+ b
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   * k  t/ `/ ~/ [, R$ s7 P0 Z* X
2. mvar:
   $ M* N% W; [: J' P3 r: B2 _
3. t=sys::clock(),) I9 G' t! {3 _1 {
   
4. s=0,* x( Z7 d1 S7 C8 F- X+ O$ D+ z
   
5. x=0, while{x<=1000,2 S/ U- z' R8 q3 i' V
   
6.   y=0, while{y<=100,/ i8 Q1 i# d1 N0 J) Q
   
7.     z=0, while{z<=100,$ v7 g5 z) v7 U0 Y* B$ t
   
8.       s=s+f(x,y,z),! b. }( ~# r) |& T4 a
   
9.       z++1 A9 O- w  R* u
   
10.     },; {* Q( Q2 v- e6 |
    
11.     y++' Q' N  h+ x+ Z
    
12.   },
    4 D9 E\" o  o& H. A
13.   x++
    1 Y% J6 D; a, D% P# x- y\" Q2 p5 ^/ X
14. },\" z3 a; F) D. P, k& I
    
15. s;: |4 A3 Y* ^$ B\" Y8 N5 L) R3 T
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
/ d3 ^7 e% x6 A2 t6 B. g

hxftntA

5 W% g! b9 o% \6 |# `* J


5 F  m9 z3 m( {2 r( ]6 W$ T

4 u# e; ^  H+ s8 s

4 z- a# G' w* O* `

' U( }; C1 s8 l7 C( Z/ s


3 c. ^% Z. F5 j$ u
, G) O% Z3 |0 m7 `' s
* W9 J) H, B6 K5 J: P% u0 v8 o