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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
) v* n7 a' ]) \+ m8 r
; C( }0 {, d$ |＝＝＝＝＝＝＝＝＝＝＝＝
: S/ M- [( F4 [+ z. L
  K5 H1 b9 u: }1 z+ x0 aMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   / U# u) q; D4 f5 \
2. tic;7 l$ f: a  V1 u/ P. x
   
3. s=0;
     o) ?# j! [% w8 k2 w& O
4. for x=0:1000* n1 ^& w: _' s: O2 Q+ y4 p! H9 G
   
5.   for y=0:100
   3 B9 Z4 {$ v$ I- E& N# x  f: Y
6.     for z=0:1004 g3 C& P, i3 ^) n
   
7.       s=s+f(x,y,z);& P) z; i  i( K% n
   
8.     end
   8 X6 |1 H! k: b
9.   end' K8 {% P  ^4 L/ x
   
10. end$ v+ j0 W5 L/ k3 H
    
11. s& A0 _  r( a* ~
    
12. toc
    $ ~' \* f3 i+ y% n& _7 N
13. 
    ; D; t! j4 E: `* A- e$ _
14. s =& ^1 T\" a2 g/ b
    
15.     6.126720600000000e+009
    7 H$ d+ j9 y# W7 C5 ]* A
16. 
    2 P! Q, v& i5 o# b) o1 S' w  \
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
% x' X* B& _' W: e" jm函数：

1. %file xyz.m, k6 K1 V+ r\" k3 [
   
2. function c=xyz(x,y,z)
   9 {: l$ u) }* _) o4 Z
3. c=x+y+z;
   8 E5 M- l* j3 i4 _
4. end

复制代码

测试代码：

1. tic;( U# T- \' Y$ B4 E9 X. O2 L7 f
   
2. s=0;
   2 `: d& ]6 a+ h# ?  s
3. for x=0:1000
   3 {- I% I' L9 |4 U
4.   for y=0:1004 M, Q% n! ~\" d3 o; R, H) t
   
5.     for z=0:1002 r' A% h+ E6 _
   
6.       s=s+xyz(x,y,z);
   ) \$ }9 I2 _! W6 f
7.     end; k+ J. A9 d- q! ^0 Y+ e$ y
   
8.   end
   4 x1 k. ?/ ]# @5 X; f, q! j
9. end
     o7 j# s7 z; p9 o( |
10. s  q+ A. C% ?. {\" U& f3 Z% U+ ?
    
11. toc, f& U( ~2 v' `8 Z6 l+ z6 y
    
12. 7 F  d/ E% R+ |- ~, p7 F
    
13. s =
    2 {$ _! C8 h( i) P+ B* ]
14.     6.126720600000000e+009
    * y5 f; `' B& m1 w
15. ' v$ e% d) o2 ?1 ?& Y# Z/ \
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

: a' D' X3 L$ ]# ~* rForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   + i\" H2 R  O) w6 B: g
2. mvar:9 e1 o# m9 m$ f( T) `% B+ z! G( d+ j
   
3. t=sys::clock(),% P/ g$ C$ T/ w' l
   
4. s=0,/ O  ?. F6 [% N  E: [
   
5. x=0, while{x<=1000,
   1 ^) e* {' [; q! G/ {3 `
6.   y=0, while{y<=100,9 l, G4 L2 U, E6 `! Y
   
7.     z=0, while{z<=100,
   , o7 W9 o4 _: C, Q0 t, S# M
8.       s=s+f(x,y,z),! P* P; e& p9 ~' n
   
9.       z++
   \" ]8 w8 \! G$ z+ y2 M
10.     },
    ' A; g\" p+ T2 X7 H
11.     y++' `# R5 e7 L! p' ]\" ~8 a' r  T: E
    
12.   },
    % c9 a3 I  d7 B/ F1 A2 q4 I
13.   x++
    5 h' A2 L7 Y\" }( c7 J; {% L' L( O
14. },; K9 C3 q1 p6 k+ z! q
    
15. s;
    7 N\" y2 b3 p) {& F
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
( W6 X) y! @. e0 Y# T1.109  秒

3 D; d" L. y+ v) Z7 \! M. Y9 \二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

  Y! v( o( ?1 L' k, u8 u, A( T3 M
, o6 ^7 M, I; u9 Y

7 v8 K$ Q/ u- Z. `( [, [! S: [
8 l9 ~  a4 T1 Y1 S" z4 d
1 a4 N- ~& [" v1 a



2 V2 Q9 c. [- v/ n( M2 Y3 h( M( f, t! d
: z/ M$ t; o6 n; V- g$ w- y
6 `* @8 U/ U9 m8 n  g# A

1 f4 a4 y7 J5 Q

, q! P4 c/ L$ t  _3 r9 t( k4 Z