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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
% {; l' ^5 ~( [
% q2 b/ u2 w' L% Y＝＝＝＝＝＝＝＝＝＝＝＝
$ P& @9 }& M' l/ A
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;) l/ {/ `7 [6 b\" B2 _
   
2. tic;
   8 |; p4 V$ t, V\" w5 `
3. s=0;
   1 _+ t3 q# }8 d  d5 k! Y7 }
4. for x=0:1000* s/ p\" T1 r4 c4 a
   
5.   for y=0:100% f5 `/ W/ Y' D2 l
   
6.     for z=0:100
   3 k2 u; n\" a+ j: r8 {# `
7.       s=s+f(x,y,z);) i  S2 F2 F$ b/ i, F
   
8.     end
   9 }! E$ D: \8 m4 w
9.   end
   ( {+ D4 e\" w6 n\" X* r
10. end
    + F+ m' ^\" S3 `& x
11. s1 p4 e  P7 G' I- R! V0 {
    
12. toc
    ! |' J9 E3 x/ C/ C
13.   L1 [) a! ?. f; u; l) j! g+ d, E
    
14. s =& Z) V0 `4 K! n$ m! L4 u
    
15.     6.126720600000000e+009
    7 I: g1 j1 c- z+ }5 @) U: R( K
16. 
    1 C1 L' |9 }! @0 J2 M/ k\" `3 C
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
, ]+ Z9 m0 Y, I8 T3 Q/ w* Y" _: r( Pm函数：

1. %file xyz.m
   1 B: Q7 o; V- T* T
2. function c=xyz(x,y,z)+ d! U4 m7 v$ ?& U$ p
   
3. c=x+y+z;
   3 R' U+ P( T: M. C
4. end

复制代码

测试代码：

1. tic;
   $ U3 D- H. m, l- j! V4 t% l: m
2. s=0;6 Q' y+ J$ ]) m* x( x! k
   
3. for x=0:1000
   0 {$ }& ~2 i- ?  }7 R% S. h6 H
4.   for y=0:100
   6 |& q) q. ]% ~8 N* g' l
5.     for z=0:1008 ^  E% O; Z$ K4 X
   
6.       s=s+xyz(x,y,z);
   6 y' J- V! }7 ?: _, [# p) M
7.     end0 @& a) k6 ~0 f- ?  y
   
8.   end' Y; [5 g4 O# L3 s: t0 Z& |
   
9. end% Q* l0 q8 @. H2 x
   
10. s
    2 `( n( x\" ^  j
11. toc$ m8 i$ x6 p9 E. j+ W, G
    
12. . ]5 S3 J7 i8 N# @  G3 R% F, F
    
13. s =
    4 I$ H, h- ]$ r3 t6 s
14.     6.126720600000000e+009- V( X; k+ m) j1 r
    
15. 
    7 @! Y$ \( [7 _
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   5 V- i8 h6 d\" W2 B\" `4 |
2. mvar:8 b  N2 h( W! F3 t2 a3 H3 J
   
3. t=sys::clock(),
   & W( L+ x4 e' `8 Q$ Z
4. s=0,; t/ }8 c9 h1 ^
   
5. x=0, while{x<=1000,6 ?6 U2 _1 T* i\" v+ m6 O) a, ^- E
   
6.   y=0, while{y<=100,
   , `  q5 R8 @& n/ `4 N\" ^! q
7.     z=0, while{z<=100,% f6 O  V- D) V! p% {, {+ E
   
8.       s=s+f(x,y,z),
   + O) |6 C1 \0 s6 M- v; [$ q
9.       z++
   2 D0 W' V* ^) d+ H: y
10.     },
    9 O5 s# Q: j- p1 [3 U3 p
11.     y++. Q& `: u5 j( e  ^3 `/ |6 b
    
12.   },
    # l' p+ w( J6 o; N# d8 ~: N9 {8 W$ x
13.   x++
    2 c9 Z8 h! y( l4 C! H6 t: a/ d$ n! k
14. },! w$ z! U5 b9 \, M
    
15. s;7 ?4 a6 e: _& ~
    
16. [sys::clock()-t]/1000;

复制代码

结果：
9 ]6 T, @1 O, i% |5 p6126720600.
. Q& X9 o4 }: z0 T1.109  秒
& A2 H% E3 Z' t& R6 G' M
5 x' m: ^( _: Y- }( Q) k9 ?5 P) d二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

" r# G% x9 ]) A, ^2 y! [9 l





" y1 }5 W3 N6 |' E8 ?
" n! _; N' [, X0 K, {7 v
4 e8 w1 L7 J6 ~
/ x! N9 p  J' Z8 m

8 `  l, h( I" f  T" x: E/ o
; t& Z5 ^& p0 W, e3 j) H$ P

; G( W: F* `& L# r4 X
( N" Y) ^" o  j& ]