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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
5 L6 X1 b+ J2 n& U
% N% z. P- w" ~% F6 f＝＝＝＝＝＝＝＝＝＝＝＝

6 h3 ~# ]& I# d" ^  kMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;- j; a6 ?1 c8 V\" f8 \( m
   
2. tic;
   / ?; O& A, D1 |' P
3. s=0;
   ( d- k) L9 k- o# N7 W1 ]/ v! K\" i$ i
4. for x=0:1000: H& Y8 q; R( v% ]
   
5.   for y=0:100+ K$ V+ @7 x! c) m/ p, p# b
   
6.     for z=0:1002 Z8 V6 r4 C) @! \) F
   
7.       s=s+f(x,y,z);
   ( `# t: ?9 Q. A. v( Y3 C+ Z# q
8.     end& W  d+ z! g$ e& ^' q$ t# U
   
9.   end- k  |3 y8 S* o& _5 c, L% B* w
   
10. end
    - R/ w# G6 ?; k/ Y2 m  m
11. s
    $ Q7 ?( p: ?\" L7 \, t
12. toc
    , W( P. o8 b1 W& q& l* ]
13. 
    8 p7 ?+ C. Y, D8 n
14. s =
    8 S9 ~! ^  ^3 x: ~6 \# }
15.     6.126720600000000e+009
    $ H\" D3 ~/ k, |- \
16. 9 V% j6 @- W1 I& H! \: J
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m: K# ?# |0 O' D4 z
   
2. function c=xyz(x,y,z)
   % o* V8 z3 {2 Y: \
3. c=x+y+z;
   . F8 [$ L: I2 Z8 }! {* `2 K; M! V+ I+ \
4. end

复制代码

测试代码：

1. tic;
     I5 ?, I1 _9 A) ~& r, \- R
2. s=0;% N! \# W\" B5 l: d- R
   
3. for x=0:10002 T\" Y, l5 H+ B) H) L; o$ ^/ \
   
4.   for y=0:100# F; _5 w0 I- F3 V
   
5.     for z=0:1001 C. B  R! o1 m+ ?1 `+ Q2 i
   
6.       s=s+xyz(x,y,z);7 S) n6 E, D9 n$ a7 W9 x
   
7.     end) o, K. \6 H# b, i
   
8.   end
   . r/ a- X7 w# S- o\" Q, D# Y8 g
9. end
   ! l& P/ T: H- Z/ y, F' ?% \
10. s9 p1 Y* @( X. _\" q# F' @3 p* D
    
11. toc
    $ @3 ^$ @. T  R7 k- M2 Q% x
12. 
    - A% n0 Z1 X6 J; k$ e! t5 |8 {
13. s =0 y  U\" v) p* t7 ]8 o6 F
    
14.     6.126720600000000e+009
    $ c! B/ t7 _6 p; `8 e
15. 
    , }\" K) q, {- V
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
3 R1 c( u% D! y" {# p
/ ^1 Q" A/ W* h! FForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;5 h) ^- y2 T( v1 Y/ S
   
2. mvar:, h) M: K0 j  k$ F# @6 f) G
   
3. t=sys::clock(),
   & f3 h) G8 r; n' f! m' X, j
4. s=0,6 ?7 x$ k\" [  o6 n
   
5. x=0, while{x<=1000,
   1 b5 |% O( S; k3 o
6.   y=0, while{y<=100,% H0 J) ~\" d/ H- n4 u
   
7.     z=0, while{z<=100,
   % A\" [% E4 |  O0 v/ @- J
8.       s=s+f(x,y,z),
     N& E$ Q1 l! N) i2 y' ~4 P
9.       z++* \( A4 t4 B, q+ h
   
10.     },
    - P% \( S- q( |1 U
11.     y++
    ; i5 t) |! T  W
12.   },
    ' \1 R\" d+ H7 M+ X+ ~
13.   x++
    # [/ n& ]# {2 T7 v# L  v/ s( ^
14. },
    . o! ^( u3 d8 Z( E5 b! S. h
15. s;) n+ g& K3 Z3 W
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

/ y, `% r8 h1 R" g$ v- [二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

* k. @6 C1 L& E* ?


/ ~" ]) Y3 C$ t! ?

& @6 L8 O6 K  H
' T: |9 Y" G5 C7 G

+ N" W+ r5 A/ A4 x' X" H$ k0 B

& k% E: ^4 V4 {9 `" a
; g6 [# U* I" C8 V, b5 R# F

8 a9 g0 t" u7 h2 g! L
9 w. n1 J# W( P

9 f$ R9 J8 g7 ?) a  X, d, `