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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
2 |1 P$ ?( e- T( A2 g% n: U/ o4 \
＝＝＝＝＝＝＝＝＝＝＝＝
# @1 s. [6 B+ v" {
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;- A. I2 E& T  j7 r6 _
   
2. tic;7 G3 m) X! b% d4 f& ^) r
   
3. s=0;
   / T/ }5 p0 }/ D. G% T% r
4. for x=0:1000
   + P$ w/ p% z. [1 N+ K
5.   for y=0:100& t\" H% x( K1 F6 W0 x  t
   
6.     for z=0:100
   ' n1 V! ]+ B! Z
7.       s=s+f(x,y,z);
   1 a: N/ u7 C; d  r
8.     end
   7 L% _* O, l$ F5 o2 n! n! O
9.   end7 s3 q, y- A; n- x& x& h9 w
   
10. end
    9 `- H$ _% A0 `& G- ~
11. s
    # y8 N+ x/ o; K7 m
12. toc
    + l\" }  h0 |/ E( r+ R! q: F$ d0 i. \
13. 
    \" n5 K( _$ k6 S, [$ A! O1 O* f+ o
14. s =
    . A  u( \  q5 s) k. F
15.     6.126720600000000e+009
    0 X+ l* p  d+ g  k$ B\" ^
16. ( Z+ |( O+ i3 o9 z0 p# S4 t
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
- l  M1 J; W9 t! F( D- rm函数：

1. %file xyz.m
   7 M\" d& L! z- n. `% V2 c) C$ D( Y
2. function c=xyz(x,y,z)9 _& q/ d; A0 @7 x2 G6 b- P
   
3. c=x+y+z;
   # s( L/ x3 R5 ?: I' v
4. end

复制代码

测试代码：

1. tic;
   8 j9 b  e$ v' W, S1 u  k
2. s=0;- R9 u; D6 W; |' Z1 F
   
3. for x=0:1000- d% W+ `# R2 S' c
   
4.   for y=0:1008 X; y2 W# u9 w6 }0 b
   
5.     for z=0:100
   , A+ Q0 V9 f' d0 T. L( K
6.       s=s+xyz(x,y,z);& V/ G' [5 c. s4 B/ R8 N( i
   
7.     end! l\" s% y3 W6 Q% M# H' Z
   
8.   end
   / x7 \- b5 P\" O) |
9. end5 ?7 W- H8 ]$ o\" V, J0 ^% `. l
   
10. s& ]: E0 O# b7 t
    
11. toc; ^' M5 Z  K! O/ b* J
    
12. 
    $ N) U7 I' m9 k: y$ T
13. s =
    \" D# G) E6 l: F* Q
14.     6.126720600000000e+009) J0 i& `+ \8 m, z& [
    
15. 
    4 ], r1 v6 P2 N& p; z3 }$ v$ S
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ' O( i+ L+ ^! l\" L9 r
2. mvar:  b' Q* h. g+ [2 d8 m
   
3. t=sys::clock(),% {, G! Q7 f$ B8 q$ s8 w* X
   
4. s=0,8 Q\" @& G0 z7 a) M/ d! K( {
   
5. x=0, while{x<=1000,/ n2 s: O3 ?/ a' M
   
6.   y=0, while{y<=100,$ G; E6 z8 f! P
   
7.     z=0, while{z<=100,
     d5 q( I; w8 w5 g8 V+ J
8.       s=s+f(x,y,z),
   , X5 W& w7 ~1 j% D+ N2 V
9.       z++3 }) ]9 M& R! T% m* K5 }! E\" F
   
10.     },1 X: O. _9 d7 v9 g4 K  r8 k: P/ u  p+ f
    
11.     y++
    ! y1 R/ s, n- S  z7 J
12.   },
    \" c+ ^0 a& N4 d! w9 R
13.   x++4 k% Z$ ^! W& H! X# V5 T
    
14. },
      m/ O! T, n# U' o3 o% p; u9 o8 i
15. s;
    ! U: Y2 b4 A' N; W+ e' @, g
16. [sys::clock()-t]/1000;

复制代码

结果：
) j8 v1 Q6 [% D$ U1 x+ f3 o; y6126720600.
1.109  秒
) v% h8 ?) n! [1 P
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

( _5 L3 G: [/ h' s" k
+ r, j" @7 a  j9 m3 B) l7 I, e! C
& e7 V. }- j/ {% h* T5 U
4 O7 C( W& s$ S# r: n" \8 L

/ T/ O: s! T1 y. Z


5 h5 I+ G( u! w8 I- ^




" X  d( [  F4 N- \# K8 ?! n' ~) u6 z
& T" |  v: S! H9 w
# |7 C) t6 m: j