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

forcal

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

$ e# z8 `6 m0 ^+ m8 T＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   7 @( Y\" g6 V/ y0 V: w# ^8 C7 c8 x
2. tic;
   \" J( i+ f  R8 |+ h5 A1 E
3. s=0;. K, l2 _$ q6 l
   
4. for x=0:1000; O5 U% j7 V5 G5 y$ Q% b' Q
   
5.   for y=0:100, D1 ?4 G# _' P* O
   
6.     for z=0:100
   ; }, n  l/ K3 b6 p( E6 V9 e
7.       s=s+f(x,y,z);
   0 `1 D9 v% g5 Z\" ]
8.     end# c7 i- \2 q9 G0 e0 s* g: w
   
9.   end& }& w2 }  h1 `3 C3 l
   
10. end3 t3 {4 d6 g6 |5 F. |! L
    
11. s
    3 Z6 z9 V8 d3 y  |0 I: D
12. toc
    5 }\" I. I5 b& ]1 B
13. 
    # h; Y  f8 u5 I4 P% o
14. s =) O6 X8 |, }9 t3 D
    
15.     6.126720600000000e+009
    # Z! a( X+ x$ y- P
16. 7 p) \( Z, s* A* {) b) l, y
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m( U( y2 N5 I' r0 x: ~. ~% g& _7 U
   
2. function c=xyz(x,y,z)
   3 d7 J; s% q. \' Q5 t
3. c=x+y+z;
   0 C! }. S. o, _$ c( U: w; e  Q
4. end

复制代码

测试代码：

1. tic;1 c# \! m\" K5 N+ U/ A& W0 z7 X9 }
   
2. s=0;
   5 J! u) E. s) H
3. for x=0:1000) k. i- v3 R: a1 w: H; R, s
   
4.   for y=0:100
   ( o9 k: p' ]. N3 q0 Q$ O& ^: j
5.     for z=0:1008 N7 y5 c\" A. x! A
   
6.       s=s+xyz(x,y,z);
   \" z4 g; }  V9 Q/ g& Z
7.     end
   $ h! I4 v. ~! W/ s. h8 W3 r
8.   end- l9 ~\" _6 Z5 k% ]
   
9. end
   * j! ]7 c( s; k  R
10. s
    - ]/ J; u: e! U/ Q
11. toc
    2 X8 @# |, p$ r
12. 5 U\" w$ R2 g: a& M$ X6 B
    
13. s =
    9 O' I- m' q+ G0 u5 U0 y7 d8 \
14.     6.126720600000000e+009/ B# [$ J, J. g5 z& d$ b% F% i
    
15. 
    \" \: a7 D. j3 s
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
3 D: g6 Z+ Y5 _3 K
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;4 G+ J3 [5 m0 G
   
2. mvar:3 \. r6 \! w$ }. ~) N6 X9 P
   
3. t=sys::clock(),( v! s3 @. {- U* ]# N  ]/ j6 }
   
4. s=0,4 A* j! D. q\" s\" x% G; }
   
5. x=0, while{x<=1000,6 i- `: m5 K# E, A8 x- c
   
6.   y=0, while{y<=100,$ x4 l  E; o# Y
   
7.     z=0, while{z<=100,$ J\" G0 W3 y! r0 z% l
   
8.       s=s+f(x,y,z),5 {# ]3 o! n  d, V  r2 N
   
9.       z++0 Y/ u8 }: C9 J5 R0 c3 R' W
   
10.     },
    $ ~0 P3 [3 ^# S, n( B8 T4 o
11.     y++, }\" r$ \7 T4 T( I4 s
    
12.   },$ m. W' \/ P& i! Z  \' F2 @# E$ r1 u
    
13.   x++
    ) |/ d. L6 j- u; \- _
14. },  l; M3 K5 ]9 g9 N& J
    
15. s;
    6 b: R$ Z! J% }% t- `
16. [sys::clock()-t]/1000;

复制代码

结果：
3 i) \' N! Y. d6126720600.
1.109  秒
% |' I9 J' m6 Z
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
0 [' t# A0 e. P  n0 G% \

hxftntA

9 e5 G6 J: R7 a6 e3 E


" l, j, J. v  _  r; J
2 }/ ]5 U9 v$ p0 Z, }' H
/ x9 T1 I7 z7 m( P, K2 y