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

forcal

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

9 O$ s* W  w& N2 m1 l+ G' [＝＝＝＝＝＝＝＝＝＝＝＝
1 }( E7 j, A6 x
' D1 |) B! p% |' b" C! _Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;9 o  c1 N$ M$ C\" p: X
   
2. tic;- z/ G\" `2 j9 b1 Q9 Z  g
   
3. s=0;/ D: e6 q; U, @5 Y% i
   
4. for x=0:1000
   & u2 T- i/ J\" c
5.   for y=0:100\" C+ s3 G. v  l% p8 f
   
6.     for z=0:1005 s! v3 ?2 l% ?
   
7.       s=s+f(x,y,z);
   : |' N% K. r% [+ Q$ {! X% S
8.     end
   \" O4 X7 y* w# D% x6 B. t3 y
9.   end
   . ^$ H- O$ F\" s3 Y, U
10. end
    : F+ j- e! T\" F8 X* H- a  n
11. s4 x- [7 L. z7 c  R6 M- E
    
12. toc
    # C# h: ^& w- t* Z* x# n. G\" a
13. \" Z6 t) {+ Y, f' O* ^
    
14. s =
    8 f2 R! `; O) `; s- O! n; m\" \; A
15.     6.126720600000000e+009
    # F4 ?& i7 A! h1 _
16. ' L& S' ^! F, C6 @8 _4 x1 W0 D9 k
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   & I5 k9 K* r$ Y: h: A& E
2. function c=xyz(x,y,z)
   ; ^3 h3 K8 t# _3 M; P, {2 v  z
3. c=x+y+z;; v) z+ y- |; \# G
   
4. end

复制代码

测试代码：

1. tic;
   : S( h7 s1 g9 ^\" }0 p( {
2. s=0;+ L+ x3 |& ^$ R* _( r& `0 n
   
3. for x=0:1000
   / p9 O5 O. V* r. T
4.   for y=0:100
   ! E1 d+ G) G# l9 ]( Q
5.     for z=0:100
   : f: H# ^\" h; P) G! g2 K$ z
6.       s=s+xyz(x,y,z);
   + ~\" {0 F4 x0 [) h8 e2 z8 R& n
7.     end
   6 d3 @\" g, K/ a( M\" L
8.   end
   # s\" E+ F7 i: p  v9 U3 p# B2 B
9. end3 P) W3 q. X2 O, `  a1 b! `9 `  r1 S
   
10. s, A' c; W\" l/ c, `& @( j9 e
    
11. toc9 I+ G4 k. I# N/ ^3 b$ e8 |
    
12. 
    & E7 W# h* U& V9 t
13. s =! l5 U' f2 j\" D: c8 s+ m8 |
    
14.     6.126720600000000e+009
    2 J5 \/ `. S( D6 J
15. 
    : C1 l2 e5 T3 u- j
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
     f. ?* v& f/ L* w. h
2. mvar:
   4 A$ x1 C. |$ T! t0 \# ?% B  K( r/ ?
3. t=sys::clock(),5 p\" a5 l9 G0 c
   
4. s=0,* C- D; `2 F' Z' s: l
   
5. x=0, while{x<=1000,
   % G1 Z1 U( Z6 |3 o. c0 [6 s9 R8 x2 V# W
6.   y=0, while{y<=100,  o4 _3 k' D6 E
   
7.     z=0, while{z<=100,
   : {) S8 e  M. P$ V6 z# M
8.       s=s+f(x,y,z),1 P* H# h2 X* z) c! h
   
9.       z++4 W' s0 ]2 n* w
   
10.     },
    8 h; v- Q- M! E1 i/ _\" L7 G5 I
11.     y++
    % k4 J% ]7 W: b9 S: u3 z( [  W
12.   },6 o: \( I. `8 {8 _6 K+ Q  T
    
13.   x++  n8 E# p( F3 ]  `
    
14. },: N\" J/ ^\" D  f% [
    
15. s;& w, q2 W9 [7 o/ m
    
16. [sys::clock()-t]/1000;

复制代码

结果：
, F* U- ]0 s# d* T& b# d6126720600.
1.109  秒
& |$ ^# q% @% ~8 S# J4 B& ^) Q
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

0 g7 O" L5 I4 o0 b+ ~& h1 @

3 _5 N, w9 k: d/ l4 l2 Q0 t
7 `  m9 a, y5 e4 q2 D" a) s



5 v/ g+ l' `# w/ V% P7 `" b



5 g3 r4 v7 C2 p



9 E: N2 }5 q8 V0 s5 @/ B, |