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

forcal

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

4 r( j5 z% O  g9 U0 {9 y＝＝＝＝＝＝＝＝＝＝＝＝
; x; S5 W' f4 x
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;# T) a( L! X7 X3 n& h, c9 X+ A
   
2. tic;6 l- _1 T. u3 D) Y# l% X5 I4 S
   
3. s=0;
   6 V- c. k  d( I2 V
4. for x=0:10001 S\" e) `7 ^7 p0 L' Q
   
5.   for y=0:100
   & c0 H  R- C% O, ?
6.     for z=0:100
   ! i2 ]# n$ h, J& M
7.       s=s+f(x,y,z);7 B  D# S2 {, c& v2 W* l7 _7 |
   
8.     end
     n6 }3 L7 o2 [, F
9.   end' n, H, Y  H- o7 Q7 }
   
10. end
    4 m+ T0 |  n6 Y8 X- R! \
11. s1 c$ k% o- j! m5 L% ?7 Y
    
12. toc$ I( A4 S3 t  B5 Q2 o
    
13. 
    * V8 G. A* |) I& v. {* N5 _  U- `$ p
14. s =) c3 m4 l/ d% V3 d4 D) [3 O
    
15.     6.126720600000000e+009
    : a& B. {2 ?8 A! ^& ~  q$ i6 A/ U
16. 
    # Q% S! _& L: \+ z; b\" H: H% t6 o
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
3 ^: {6 x% L8 F, T1 M, Em函数：

1. %file xyz.m& ?- |9 {. w4 M
   
2. function c=xyz(x,y,z)
   * H% p/ K6 w% u* t$ n
3. c=x+y+z;. V! G) C2 g. {+ P3 V# ~! `& x5 b
   
4. end

复制代码

测试代码：

1. tic;: f/ V) u4 r( v, I2 W6 Y0 w/ I1 m
   
2. s=0;; |& j1 W% W! _' K7 s  ?
   
3. for x=0:1000; t+ R& l0 D. [! U, S4 g2 O
   
4.   for y=0:100
   2 M7 C* w, B/ @2 ^, e4 g
5.     for z=0:100
   ( X0 K# d1 i9 }, @( P  F1 W) t\" K
6.       s=s+xyz(x,y,z);& d$ P( q, Y\" U' W9 x+ q
   
7.     end
   : ~. c\" r, E- }) O; a5 @( j$ d9 d$ j
8.   end
     Q9 F0 Y) A3 `- \( y) z4 X/ Y# w
9. end
   4 {9 }* d9 I  U0 G
10. s
    # I* z# f8 N0 ?9 t& d
11. toc7 P1 b$ I$ V6 P, {/ K
    
12. $ M; V4 ?9 b7 i
    
13. s =
    ' L1 z: C0 \- l. x) {8 P
14.     6.126720600000000e+009
    3 B3 T+ B) G0 D- `% L4 J: s
15. ; s+ m, L  L1 A2 m9 V
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
  P% A) J& T* s( q
7 X) {* t# Y. j6 U& r- GForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;# i* }+ [; J+ t, s% `9 T. i
   
2. mvar:% |* \1 m+ D/ ^2 n& Z
   
3. t=sys::clock(),
   & |: b+ m+ {' D  C- R3 ~, j
4. s=0,' v3 i8 ?, w! }/ Z
   
5. x=0, while{x<=1000,2 y% g: H8 _! X4 I\" J\" T
   
6.   y=0, while{y<=100,/ Q( s! Y( f7 S+ ?3 a$ S3 E6 D, l
   
7.     z=0, while{z<=100,
   ! V4 p% l5 U; {8 C0 @
8.       s=s+f(x,y,z),. G- [0 X' d$ e, B
   
9.       z++
   , M3 @$ w% u3 c( }3 W5 }  B
10.     },
    9 M+ K\" F4 [  ^1 h\" e
11.     y++
    3 `: B7 r% [  F
12.   },+ R) }2 C\" u8 N  l5 o3 B
    
13.   x++# i1 K; @) T; O+ ]
    
14. },
      s7 H- {7 f9 Z
15. s;9 I& l. f6 K( ]( ]5 @8 t: d
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
% N/ W4 Q9 g' }& _  H) i1.109  秒

. D& E) Z3 h+ I$ J, S( b二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

5 b& c9 K$ M) w5 G* b- }
: |9 G$ i# ^4 W1 ]



, }: O; t* t9 t* O+ M9 h* p+ k


! M) B2 K* B- Z" R

4 Q: g. Y! n5 h* D6 I

5 @9 B( {7 r. U4 Y8 l
  h' d0 f8 @* `2 w; R/ Q: V
: h' y) q% n! [. o

% p: ?0 c2 I: h3 l3 ]