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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;/ a\" Z  Y3 W8 V( q! E, @
   
2. tic;
     b! T\" ]& D4 R/ T9 Y) z4 Q: g
3. s=0;! W* i* @- n3 T% T. q+ m+ _* L
   
4. for x=0:10006 h3 N$ B3 g: h
   
5.   for y=0:100
   9 A' ^+ x; p' X; M9 m
6.     for z=0:1002 T: a8 D3 b: J. Z  q  J/ n; b
   
7.       s=s+f(x,y,z);
   & G# b& e$ j, b
8.     end! A( a0 W1 b  B& q\" o6 R
   
9.   end( u3 ]2 b% P+ E6 r8 `1 A
   
10. end
      L\" {\" z( }2 t1 P4 P( g1 E
11. s$ f8 j- }/ N$ o: M2 k7 R
    
12. toc
    9 o6 d4 @4 K  C9 Y2 A, F
13. 
    5 |% Y0 `) b3 L5 c, v) s( h
14. s =/ J7 [. {8 @/ ?\" h
    
15.     6.126720600000000e+009
    6 e0 X- I! g  D9 E
16. 
    - W. A( A; F( I2 `6 v
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
) I/ G; U/ C$ H6 w$ e( m% jm函数：

1. %file xyz.m- X5 F+ }7 ?* c- h7 c6 L
   
2. function c=xyz(x,y,z)  a# ]# R0 z' w) r) Z
   
3. c=x+y+z;) J3 ^1 ?\" b* w; @9 q/ X$ _
   
4. end

复制代码

测试代码：

1. tic;
   ! a\" o, Z6 _9 A( L8 y. J8 P5 A
2. s=0;
   ; W* V8 K9 I4 H! {
3. for x=0:10001 k: t/ e9 C0 ?
   
4.   for y=0:1004 ~* J' ]2 K# ?8 c: q
   
5.     for z=0:100
   / @( z$ A, W1 q2 W2 _4 e
6.       s=s+xyz(x,y,z);  G0 w, b\" b; j% O. i
   
7.     end
   3 g) B1 M9 S! ^+ `: V: N
8.   end
   * o- l$ M4 P- T0 e0 ~' B
9. end/ m, u6 Q2 ?7 T7 C% m- r! z! ^
   
10. s
    8 I! ^  F' @7 v6 U* _. R
11. toc
    8 \6 [5 R) H6 a; n3 d
12. 4 E- B4 E3 \+ M$ `# J
    
13. s =2 V3 \! z4 h% j
    
14.     6.126720600000000e+009# w8 j4 f4 {/ j3 u8 h4 I, L. }0 e
    
15. & r& c% P\" x6 u2 }4 d7 H5 E
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
) p; ~6 E4 v6 z8 Y1 J; ]
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   5 l9 z+ v/ T\" |# d: V. p
2. mvar:
   4 _: X\" x% d8 k. W$ O: M0 v- b$ j
3. t=sys::clock(),
   \" v5 T/ P& [6 n8 U, y
4. s=0,2 ~$ c8 n( X7 i$ M9 N/ M
   
5. x=0, while{x<=1000,
     ?+ r- Q) p7 ?  M9 k0 K\" a
6.   y=0, while{y<=100,
   7 c' r\" R! ~3 S( H9 J# Z
7.     z=0, while{z<=100,7 |& w) p3 O: _* g+ R
   
8.       s=s+f(x,y,z),\" n% {9 H0 h! C# D
   
9.       z++
   ! E- o2 H; L8 o: U  v. p
10.     },
    ' D5 H7 Q# S  }( I
11.     y++3 K; O/ \0 e4 V, O* U
    
12.   },; h: B) x, m6 R1 p, w: d
    
13.   x++* L! ?8 s\" J2 E: P3 P9 j* v5 _1 @1 o
    
14. },/ \( M# e; W1 U9 p* Q' M
    
15. s;2 _6 Y4 w' X7 B9 \\" n  O5 l6 A
    
16. [sys::clock()-t]/1000;

复制代码

结果：
8 R% }! c8 C" R2 \6126720600.
  m5 b0 O, E4 c: ^1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

5 z5 g1 s1 U; v1 u% N
( t9 Y0 E6 b/ u' L
! T: y. m; z7 U$ H4 p" U
' P% U! p8 l: V
+ ^2 m2 b* U' Y- `% o


8 k8 t; M) C  |
. C. a8 u" M( \7 D


  s& Y9 a7 Y$ l! `; J



2 `) _# k% Q' X