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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
& i' N6 u8 P) v' @
3 ], g( \* O. {+ X0 q＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;+ o9 C  x6 K& I+ g- s. M+ h& T
   
2. tic;0 Y! b+ N0 W$ s
   
3. s=0;9 K; c- b5 \- u\" w
   
4. for x=0:10007 v' i: ~8 N8 Q- F5 ]. d
   
5.   for y=0:100
   4 C\" R3 g- o( S& `) o6 C' ?/ r
6.     for z=0:1004 C! m9 m0 T, C
   
7.       s=s+f(x,y,z);: F6 Q/ N. b, q( [  W
   
8.     end
   ( o9 u, u/ @  u  {$ a\" x
9.   end& K7 S\" b! i# ^& R- y3 ~, v
   
10. end
    1 |6 n$ R' s% E7 Q% Y1 H( o( }1 H
11. s
    : W2 R0 f, Q; ~0 V- X6 s3 z  M8 ^
12. toc2 O! x  L8 S8 \: m. i5 }
    
13. * j2 n4 t+ G, v
    
14. s =
    & }7 S' J, m( N6 a6 U& O
15.     6.126720600000000e+009* t( K4 p, `$ F0 V3 X) j
    
16. ) D; o* b' z. w/ ?0 G6 `4 G
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
+ y$ f/ K4 I: K! G  p3 z/ L6 ym函数：

1. %file xyz.m3 o# P/ w6 ]6 ?
   
2. function c=xyz(x,y,z)
   ; l\" m) G) }; p4 S& N$ h
3. c=x+y+z;4 @  B. `8 s/ {  S% C
   
4. end

复制代码

测试代码：

1. tic;
   5 }' h$ V/ h7 X6 J
2. s=0;
   7 ]! I0 @1 A& C& Q; u6 L2 q+ E2 e1 @
3. for x=0:1000
   ! C5 V1 K. O- S; \; x' `4 M  p
4.   for y=0:100
   / b' V4 u. t0 C% \, v1 _$ Z
5.     for z=0:1005 h4 w+ T: Z( _' A6 _7 d7 L# N& V+ A
   
6.       s=s+xyz(x,y,z);
   8 d# r7 U. _& ]+ U
7.     end3 u( R6 P. X- X\" K2 _- `1 L
   
8.   end9 ?4 g5 U! I. U+ Z: u0 O
   
9. end
   + p+ t& e/ ^0 r4 ?8 D
10. s
    - L7 b\" ~& p) k8 u
11. toc
    . ]1 n: \! Y: K3 }* T
12. 
    % J/ W- w3 b1 u3 d0 _4 b
13. s =
    5 L/ Y1 f) o2 }) T5 ^\" L
14.     6.126720600000000e+009
    * G0 {7 Q9 j+ v# C
15. 
    ; d. V7 z  ?% Y6 `, u0 H3 s
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
     M& k, Q3 E, L5 f
2. mvar:; [5 V2 P+ y) b& D* d- W
   
3. t=sys::clock(),
   . p. m, D( N! L
4. s=0,
   $ y5 O& u! L. w
5. x=0, while{x<=1000,+ K: Z2 n! l# M\" v
   
6.   y=0, while{y<=100,2 r2 v+ \9 W$ ?$ n
   
7.     z=0, while{z<=100,& G( T6 Q* i& R4 V  E5 z2 |/ \
   
8.       s=s+f(x,y,z),) ]: r2 U; g! _\" q, b: m% o
   
9.       z++. `; M; q  K+ q4 B; f/ y
   
10.     },; _  G9 \% S! }. |( H, S
    
11.     y++
    $ i# D; w/ d9 x2 Y/ a
12.   },* t+ I6 ~1 ]1 {8 @
    
13.   x++3 }. t- A( V6 U. _2 N1 {
    
14. },
    / n' P+ Z' M$ C% H
15. s;$ d0 Z1 X$ O: f* Q
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
1 C* o, U. [( A, z8 ~% w- l, D6 Q
2 a4 Y4 O: d" {: ~二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

7 w4 h8 p( T0 b7 _5 b% h/ d

, D  f8 U' b& ~4 I

2 u" ^. o' Z* D1 Y) l$ R


) c) l& ^* h4 J% P

  n$ l- |' U% a; o7 c( o

7 `; b! b' f; B8 a3 A  f
* H* |4 q, j* \% M, H/ r8 J9 P