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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
0 H; F9 |+ T& m' Q' P
# B5 _- {/ @+ P6 ?/ X＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;( D  O( g: |6 _! s7 `
   
2. tic;0 F* K0 V! L* L0 v. c\" n' v1 n6 t
   
3. s=0;0 Z6 J7 G# z5 Y8 j  G2 k
   
4. for x=0:10004 g\" l' ^1 m, J3 G: Y
   
5.   for y=0:1007 r/ P% V+ ]\" f9 D5 |- [8 C
   
6.     for z=0:100/ E0 S' f$ g+ |3 J4 O* k* s
   
7.       s=s+f(x,y,z);7 S7 W) b3 ~) [/ r
   
8.     end
   2 U6 T* p8 F$ @  @; C/ E( K' j
9.   end
   * D; ~) b) T, I- `. V0 ^
10. end7 ]$ a3 M! ]; n! T% h* j
    
11. s  _5 b4 `4 @! e& Q/ w; U9 W# i$ V
    
12. toc
      n8 _& @0 X( u! Z, a
13. ! H- k9 a6 J5 }5 Z3 Q4 O
    
14. s =' b\" {1 M\" E3 o5 E6 c9 I  _8 Z
    
15.     6.126720600000000e+009' j) ?, f! t9 d$ o
    
16. 
    ( C6 K( `# a+ @
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   * R! \7 H& W5 R  I3 S$ ?) q: q3 Z
2. function c=xyz(x,y,z)
   8 u1 B5 k, b% u( D\" q5 w
3. c=x+y+z;( R0 b. O, I\" {0 z. H
   
4. end

复制代码

测试代码：

1. tic;# A1 X% b8 q; b: F
   
2. s=0;
     R\" j: n( O, v$ N5 ?5 p  d7 W
3. for x=0:1000/ L9 R5 s  d' r. C2 b& ]
   
4.   for y=0:100
   0 I6 p) `6 ?2 b% E, }- c4 K8 l
5.     for z=0:100
   9 i. A7 o4 x5 W4 r
6.       s=s+xyz(x,y,z);
   4 A3 h7 [8 |% v* x; R
7.     end9 P, v6 Y. u; a
   
8.   end
   4 |: c( A\" ~% h: x# `; R
9. end
   4 x- A4 N2 q* A% z. }& S# z
10. s5 b# L% D4 q  Z+ F+ Q! e
    
11. toc' u0 x9 E, G9 _2 a) f3 l+ Y+ a8 P
    
12. 
    & `$ Z  R# B8 p
13. s =
    ! e1 X& f- \1 j* D
14.     6.126720600000000e+0097 r2 O2 F% v( L* v
    
15. 6 M  Z' {: A+ b9 w
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;$ d8 ~4 G! j! I+ J) H4 p( ]
   
2. mvar:
   + }5 n0 p3 K* n8 {
3. t=sys::clock(),
   $ ^8 `  f- G1 T: H: i% I5 t
4. s=0,/ W) I% O7 ?  t. G4 B% K$ M3 b
   
5. x=0, while{x<=1000,
   $ }/ W; o5 R8 a5 V
6.   y=0, while{y<=100,1 M0 J9 C+ i7 {7 ?2 m) n
   
7.     z=0, while{z<=100,& l' T) H, F, S3 [
   
8.       s=s+f(x,y,z),
   7 Y/ Z- O0 X5 k! U/ I% {& L
9.       z++
   : p: y5 v) W0 i6 N& ]\" ?
10.     },) V1 ^6 B: f- Z7 W6 P$ C
    
11.     y++
    ; _; {! P- L1 A+ z; d7 q+ `$ l( R
12.   },
    6 {& u! \: e# Q% P- u
13.   x++* V  Z# @5 s( |/ x* B' U( b5 ], {
    
14. },# I1 r) n4 Z% ^8 Z9 H
    
15. s;
    1 i* [8 u* ~6 z! |, M1 A; E
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
" ?3 O" k6 R) H. i- V, G1.109  秒
* h) ]4 x  w. l4 E8 ?7 T( g
5 [1 r/ e9 Q( a二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
! k/ p2 x2 B' @

hxftntA

5 l" f; q+ l* F' E1 _
9 Q( w, {4 d; w4 k


. M* T$ A9 N0 G* \- Q

" p( Z( d0 g3 J/ T( N

/ \1 }$ ]( G7 C


+ T1 I' H. J+ R2 O: _' ~' q
* N3 A* }8 o6 g2 L4 X8 V5 |2 S3 E
  S+ o8 G' w$ o8 X. f
% i8 o7 Z& e7 O2 b$ ~: v  m/ l