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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
9 T5 R0 R9 o8 X, Y; I
. r" R* h; S4 n  N  kMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;) V$ [  b* L\" j& G
   
2. tic;1 D6 R* O9 y) w( ~' @
   
3. s=0;/ u0 [$ G1 Z2 i. W: X: h/ W$ [
   
4. for x=0:1000/ q3 m7 ^( M  K) ~. {, ~' `/ Y/ F
   
5.   for y=0:100! w5 E9 n7 ]5 x+ X
   
6.     for z=0:100
   9 O  v. z! C8 \
7.       s=s+f(x,y,z);0 j. W( r\" X3 S
   
8.     end. m0 ]# V5 `- o. v8 f5 k9 v% {
   
9.   end\" v: h! c( E& [7 G7 C; x
   
10. end) ^\" ^2 F4 y\" @) E8 t, H2 _
    
11. s1 y% R' k0 [3 y
    
12. toc, Q/ @* ~& K+ [\" i\" m
    
13. 
    ! [' I# o3 S& `+ r5 ~' j\" k
14. s =
    , m/ a; ?- d  F; g( P' Y( p
15.     6.126720600000000e+0097 X9 U; c: A( x- c# w$ U; e; X
    
16. \" d% |. f2 G: {! H) u
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
# @1 q4 r- E$ }2 t7 Am函数：

1. %file xyz.m9 @# L\" z\" K9 D% S
   
2. function c=xyz(x,y,z)
   # m8 |' v; F0 e6 n# l\" ^
3. c=x+y+z;  T# c5 Q8 h+ G: V; c, f  [
   
4. end

复制代码

测试代码：

1. tic;
   9 A) v9 F4 N% e, A  Z! c- s
2. s=0;
   : d\" O+ ~/ R: x+ q' y+ C) J
3. for x=0:1000
   , g1 m# \* ^) Y4 U( Q
4.   for y=0:100
   ( Z: m6 A' X0 s
5.     for z=0:100
   % s: P$ z! @1 c
6.       s=s+xyz(x,y,z);$ ]! }' c  V( `6 S! ^
   
7.     end7 ~$ L5 u3 T9 o  Z
   
8.   end
   ) Q8 Q' a& {: i\" f
9. end* l) }+ o5 F) y4 I+ S3 p3 ?
   
10. s1 y9 M9 O\" L8 a
    
11. toc
    9 L+ s. E5 `) M, c1 V8 [
12. 
    / R) Q+ Y% D9 _3 N& n6 |
13. s =1 i5 Q6 G/ X; Z3 d* d% V! m% J' a# g
    
14.     6.126720600000000e+009
    & |! P% j( z6 y/ K) @) w8 S
15. 
    ( D. T& N! B\" m* L) k# ~
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
" E% q; p% `0 f( ^) n' ~  ^, K8 t$ q
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;- Y* R4 G$ @7 O
   
2. mvar:/ k, _! a  g6 n\" p) o$ W* |% B$ b
   
3. t=sys::clock(),
   ; f5 N: ]! J\" C0 L3 t! Q
4. s=0,
   1 `: G/ j/ n% I
5. x=0, while{x<=1000,
   2 m9 x6 b- r& L# h: t, m
6.   y=0, while{y<=100,& q, }) G' ~2 G\" c+ a\" Q
   
7.     z=0, while{z<=100,+ X/ T$ Q+ i, j: a
   
8.       s=s+f(x,y,z),$ T0 |8 L\" b1 K, G* F8 P
   
9.       z++
   6 H4 N8 `6 J  F& k7 |* L& C
10.     },& w$ O, W2 h, C. S; H0 N
    
11.     y++
    2 E2 n% K( ]; I; q/ ~3 \
12.   },8 n0 ~: v. E  h
    
13.   x++
    , J! E1 O& ?\" }6 @
14. },
    ; Y% x' Z. \: Z/ l% d9 V: f) K
15. s;9 L  ?, ?! l# b: a, }2 I: ~
    
16. [sys::clock()-t]/1000;

复制代码

结果：
8 `& M9 Y/ K: m1 Y$ P6126720600.
1.109  秒

- r' M2 e: J' O二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
  o7 v: t. c- M

hxftntA

! S# q4 U, G, r+ J8 Q/ d
8 D! [8 q1 S& i$ G

  [2 C2 M. c6 @8 \3 U



) @! U% u1 `5 |
1 n6 f2 I% f8 x: ~$ O! c" f% p


( y% Z( t2 o6 X/ H8 e
$ S! u$ ~- B. @, {8 q

6 E3 R& Q6 V  T- O! Q. |/ ^
" H/ p0 G  {* r1 f- U4 h/ ^" d