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

forcal

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

. r+ d! j/ A) K8 K( B1 S2 G＝＝＝＝＝＝＝＝＝＝＝＝

- F7 c/ P8 `* _% o! c1 }Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ; ]. f) g: Q& |, \0 o+ T& d! F
2. tic;
   * @0 a3 I, C# r4 u
3. s=0;
   9 ?) w9 I* j* G# x7 R( F0 X
4. for x=0:10001 @0 e\" [  D' C8 o
   
5.   for y=0:100
   / A$ c4 h2 x3 m6 J: B7 d; r\" K
6.     for z=0:1003 j* T' G, q% b
   
7.       s=s+f(x,y,z);. |: A+ |1 u3 a
   
8.     end
   6 }7 [$ k3 f$ a( V+ k! e
9.   end2 o3 @8 H* j; y# D2 ?+ w
   
10. end
    ; o2 U' _5 o3 H4 z+ h
11. s
    \" z0 A- @% b; h. y7 H  Z# |. b
12. toc
      ?* f( Q9 S: M$ }! R0 S) _\" O
13. $ Y! c$ V! Q0 O$ c% v) n5 V
    
14. s =
    & `8 ?- F& S$ D
15.     6.126720600000000e+009
      E) q+ Z# o  |
16. 
    . A( B' |+ ~) Q1 h
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m! X* A, B8 s# m/ \, P
   
2. function c=xyz(x,y,z)
   + ]7 d) u) J4 R+ }) H' A
3. c=x+y+z;! u; J; @* y% m
   
4. end

复制代码

测试代码：

1. tic;
   . I: l% Y\" A7 q\" y$ \- {
2. s=0;
   8 J, P, m$ H7 K\" f2 P) c\" S
3. for x=0:10008 v4 p% q4 X2 x. p: f/ D
   
4.   for y=0:100
   ( w9 P! _& Z5 I& y1 y7 x2 [
5.     for z=0:100  r' l8 g  {* X, i1 m' D
   
6.       s=s+xyz(x,y,z);2 {/ }\" U7 q+ g3 @\" t
   
7.     end3 |+ G, j# }. v$ Q& }
   
8.   end- x2 N6 z\" R9 j
   
9. end
   9 u3 [5 M  p0 s! ~3 G/ t2 a% S
10. s
    3 x' D) Z; d3 l\" Y1 x
11. toc
    $ N7 k+ E  c- ?- O2 j& e8 ?6 q  E. o
12. 
    & u3 v2 Q  L# b. d! `6 v6 m
13. s =8 _' ~8 X6 L- S9 E
    
14.     6.126720600000000e+0092 ?& D. @: h& e# h
    
15. 
    3 j# J) u7 O% M1 X0 C) k
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
2 u/ D1 V4 p! g+ q5 e
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   . h5 O8 J4 V9 {; v% U
2. mvar:7 K- o% K& R4 I\" ]0 g
   
3. t=sys::clock(),
   & q3 S6 o: \+ G+ f
4. s=0,3 R4 N, U4 B$ z1 L' Y
   
5. x=0, while{x<=1000,
   5 ?+ n! L5 ]# Y* r\" l. d. i
6.   y=0, while{y<=100,
   - d3 Q/ L' `\" T. S6 P
7.     z=0, while{z<=100,5 P* Y\" S& U1 [# ~' `
   
8.       s=s+f(x,y,z),
   / x$ s* S6 y3 x% p9 W: s6 V
9.       z++2 O' p% ~4 s) h  e' X
   
10.     },
    7 J' k  g, m! R5 d, Q3 l
11.     y++0 `+ D) ]* N5 @8 w6 [- a
    
12.   },
    2 E9 m3 ]9 n2 X. a9 f
13.   x++' u6 e/ j+ s6 F5 }. A: f5 _
    
14. },( q8 _) P8 G; k# p0 g
    
15. s;
    2 x4 N) I5 u) a; B$ S7 E5 j$ d
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

; f# B/ Q' z2 z+ F( n二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
2 y, \1 x* n6 g5 n

hxftntA

. K& ^. I$ q0 ^. j( {& ]
0 K8 b" \  R0 S* }8 h


  S. N7 B1 ^" }6 K5 k
! D1 ?6 N: D9 }  Q/ ^7 p

6 {, |) q* ~; `

5 S/ X- @( ~$ C8 W0 Y, n
8 b  O2 h5 w1 Y4 D$ D9 o6 ?; i+ ]" f
+ z6 h3 S" F+ g! R' a


+ S& C( n* B% a) a. x

, c  l& p( O6 Q: E- {& H