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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
* P% E1 t% F/ U0 G4 l
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;\" ]% Y- L4 q4 Z' y' a' `! F
   
2. tic;
   7 G) J+ ]: i4 ^! N/ r
3. s=0;
   3 ^2 U% ?' t; ^1 s- x
4. for x=0:1000
   , }' ?$ E' T% x: i: S8 x, }/ W2 ^; D
5.   for y=0:100: r! @\" J9 |' V5 S, r- p; ~
   
6.     for z=0:1006 C& x1 v& }' w' N1 d( _
   
7.       s=s+f(x,y,z);
   ( {2 W/ X; ~0 B
8.     end
   3 W! V5 S) w\" y
9.   end
     o+ A( q' I7 ?. L
10. end
    2 r0 K0 M4 O+ ?  N
11. s' _* G+ c+ i\" ?
    
12. toc6 x5 T( n2 \9 r( x
    
13. 
    5 o2 Z5 c& ?7 V2 t( o
14. s =, e! P; R! _) v; k( p' ^
    
15.     6.126720600000000e+009
    \" b& n8 O% ?5 H$ z, j. _# v- P
16. 8 r% s9 H\" h  z# S- ?
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
. P/ {: T& I; T  A% b+ S2 l) @m函数：

1. %file xyz.m3 F8 @( \! i) B+ v* l
   
2. function c=xyz(x,y,z)
   ) _! r\" K; M% U5 H, `
3. c=x+y+z;
   5 d, B0 _4 |) P( C# O7 ?8 ^
4. end

复制代码

测试代码：

1. tic;4 }# e% v  Y( \: L8 {( r
   
2. s=0;
   3 ?  O! J5 z. g8 R! ~, B; E
3. for x=0:1000
   : V( u5 w& U# v; {$ U\" {
4.   for y=0:100. B3 X0 W! M/ l
   
5.     for z=0:1002 o8 P8 }- j3 q\" y
   
6.       s=s+xyz(x,y,z);) A: t7 E+ e6 ]% k; b
   
7.     end- m1 C8 ]5 w1 e) |- E, v: J
   
8.   end
   $ R0 d! ]3 Q; N' v  r; Z
9. end
   / S) D) M6 n3 n7 g7 p3 ?
10. s
    8 N) C0 B1 k4 f- [2 [- ~3 I
11. toc
    ( a- u+ a$ [4 y9 Y- X7 ?9 [
12. 
    / F) I. Z$ e- Z0 m' v1 t6 s
13. s =
      b, H: a0 x2 |& C. Q# q5 T2 U/ E, \
14.     6.126720600000000e+009% U0 Y( J# ~* [5 W$ }3 Z# ]
    
15. 
    9 M: _( r' g6 u4 e$ ~) U
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
6 ^9 e# V7 @. @; M% C- f( i
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   \" U. Q: j, @) p$ ]& n0 s
2. mvar:
   : e3 \0 ]4 ~8 q' S: V
3. t=sys::clock(),0 b9 I- Y; \' g' _4 G
   
4. s=0,8 O# ]% F2 p2 T- z$ t5 e) B3 Q
   
5. x=0, while{x<=1000,* x6 S* U; s0 M' E- @; R/ A
   
6.   y=0, while{y<=100,, ^' P3 s\" W- R, T( D: T
   
7.     z=0, while{z<=100,
   ! Q\" C4 q1 G, s+ B8 R$ ^  C
8.       s=s+f(x,y,z),
   % y: D- t& z* N; C
9.       z++
   , B6 B\" s1 K\" c( Z6 w\" M
10.     },
    \" C7 _0 }# }6 M+ Q\" B1 t
11.     y++
    4 ^- T& z' a% `6 M+ v: W' z! x2 X
12.   },% W: k\" S. f6 N0 |8 M( o8 K
    
13.   x++
    : g1 X' c- w/ k* H' k- x
14. },
    - C5 R7 e1 N% t) s8 n7 r6 R
15. s;2 o: l9 y( ?5 @/ w, u* u' [
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
2 k2 ~/ k. m! q7 {& _2 G6 J5 o
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
8 A9 R# z, s) f" ~

hxftntA

0 }0 U8 `7 I) g
% ^. o2 P# U' P  G! s2 E* ?$ C; F
, ~+ t9 N9 r& i1 @" G! I* p



6 h7 f9 A2 K* `) y

% x/ Q, J, {4 |# f& \
5 |# v) E2 s* G9 f! D/ h
; u. V+ D" e: _3 M6 [6 H

% A7 G* l2 {+ S

( W# n/ }/ h: @! p+ C

% ^# B+ C6 D) @