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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

2 M) m- {8 Z( h0 [Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;! F; ]& M% Y$ U
   
2. tic;; b6 ?, ^# s+ ^
   
3. s=0;0 G, p. L: l1 N+ N& _
   
4. for x=0:1000$ F0 T3 C( N) T8 z0 h8 w: o3 Q
   
5.   for y=0:100
   * d5 v4 h  F+ `% U  J' x\" I
6.     for z=0:100
   % t3 H2 B4 n7 u: p
7.       s=s+f(x,y,z);
   ) p# x1 ]- |: I/ j, P1 U8 Y/ R
8.     end1 \- r/ \  `/ Q# x. p
   
9.   end
     y- L\" V; }5 o! y\" }; c\" K5 `
10. end
    3 ~4 ^: x2 m8 `& z9 Z, M
11. s
    + p0 J4 R# h3 Y
12. toc, Q/ U& q# o8 R; L3 \2 t/ D
    
13. 
    8 _8 v7 a; K; y- A# h, z
14. s =4 d' v, }4 @9 m! Z# f: o
    
15.     6.126720600000000e+009
    3 w& I+ D* d3 C8 l6 W& c
16. - g& P$ |0 |  ]
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   # F2 l4 Z8 |' N* x2 E* r
2. function c=xyz(x,y,z)
   9 `* a; @( ~( \3 e' l: n
3. c=x+y+z;
   & m2 {$ |3 S; k4 k# X; t& Y- o' W
4. end

复制代码

测试代码：

1. tic;
   8 W% l& O( z  @4 g
2. s=0;; N- f6 n# Q; p
   
3. for x=0:1000\" R- i1 p% a5 t: T3 j( t
   
4.   for y=0:100
   0 O5 {5 V+ }/ ]) z$ @
5.     for z=0:100( M6 _! x5 O' S( I$ r
   
6.       s=s+xyz(x,y,z);
   8 T4 c$ l1 v  k( c
7.     end# R; ^) |! U' ]% f
   
8.   end+ t0 w- l; F2 h3 \# ~( ^8 r/ m+ }
   
9. end3 u$ K8 I8 ?* q
   
10. s
    4 g8 {* B2 _* j
11. toc  B( v\" m2 N1 M% g! N+ X
    
12. 
    * t! F( ^  r/ D; @
13. s =$ E, X  Y0 H) E& j/ Q+ o
    
14.     6.126720600000000e+009
    ) _: i5 v3 E0 d4 E: Z9 P  c
15. 2 y- w- D$ `0 m2 ]4 f( h. n( \* Q
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
; e1 e+ j4 i, E% E6 o
5 p/ H: T% B$ v$ p* R* HForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ! S2 \- q8 i2 U+ X
2. mvar:3 \% t/ \, Q: |9 p5 Q1 \9 y- d3 o4 i
   
3. t=sys::clock(),
   ! `5 n1 c# O% o, \
4. s=0,: l$ [+ H) H1 x& {) e0 l& l
   
5. x=0, while{x<=1000,1 m. U2 C% V, \; G& B
   
6.   y=0, while{y<=100,
   9 j% g% W- W+ z2 f
7.     z=0, while{z<=100,( A( Q8 n\" l! v/ X) m: s9 }0 f
   
8.       s=s+f(x,y,z),
   ! B4 I& c4 T1 ?
9.       z++2 }# _, |- {( K  K
   
10.     },# C$ J0 ]9 I- L! g  g* b( q7 G
    
11.     y++, k! M; b: A+ e
    
12.   },8 j\" y& q\" W9 c* b2 `
    
13.   x++2 z. g( |5 C- g/ J4 b2 k
    
14. },/ l/ Q0 ^, [% v# v
    
15. s;+ _6 Q( w5 i0 M% H7 b5 Z
    
16. [sys::clock()-t]/1000;

复制代码

结果：
0 ^% [0 ]5 O0 t6 r6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
6 }1 ?; O$ W7 M$ q/ h7 K8 b5 w

hxftntA

& D- @8 X. d) `; t! h1 |
3 Y1 e6 H% X, d" s: j  c/ @# `
! g1 r. U- z# Q4 _$ ~

1 c9 f+ I& I" z& Y8 w





) U2 U6 R4 d3 ~- x. l; h& o# l, r



0 N, p& ?* j8 R, a/ _( D  C
- I$ r5 O3 m  F7 g