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

forcal

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

# a5 c4 }  l" W. y4 N＝＝＝＝＝＝＝＝＝＝＝＝

3 ^/ P1 q  X1 X6 u1 _) cMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;. V# T) n$ ?5 E
   
2. tic;
   & [2 m# T2 p5 n) b9 [
3. s=0;( r. J) b\" R2 ~8 \6 O
   
4. for x=0:1000
   2 h\" J  b3 H# w& T\" q( Q- K% Q
5.   for y=0:100$ ]; l. `$ O. _
   
6.     for z=0:100* x8 b7 I9 h! T
   
7.       s=s+f(x,y,z);7 ~6 h0 p% j. y& E$ @
   
8.     end/ Q( }- F\" u9 \: V! H
   
9.   end
   , v& b/ ~+ t) F: o3 |' z' O  D9 X
10. end+ Y# N1 W7 a3 w
    
11. s
    ) Q3 e1 _9 Q' M: d6 [
12. toc
    # b% {% J8 i2 M
13. 2 V; o+ X  R; ~2 M
    
14. s =
    % G5 a) j; |6 |9 p
15.     6.126720600000000e+009
    9 }\" @* a) `) Y5 s+ u. O* J
16. 
    ; B8 W$ r' a. M: E' a3 I+ Y
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m- j9 k# V8 [% q7 o% M4 Z
   
2. function c=xyz(x,y,z)
   6 d# r4 n8 t/ _; t/ T* @; c5 s
3. c=x+y+z;5 Y% U. q0 x2 b% P3 J8 q! C' y  \) A
   
4. end

复制代码

测试代码：

1. tic;
   2 b0 H/ ]. _) I
2. s=0;6 @' [, I) p  L, I1 u6 u/ Y* N
   
3. for x=0:1000
   0 I) w: y: P1 I: i1 n1 N3 M
4.   for y=0:1001 H7 R( `2 _& Z. }% Q
   
5.     for z=0:100; b# D% o0 w5 {, h) v* U! g, D
   
6.       s=s+xyz(x,y,z);$ X3 m$ r, @( n  F, s% x/ o8 q
   
7.     end
   & V% S# V! h6 @% ~: K
8.   end
   ) ]\" B( F8 N  F
9. end
   9 h8 e/ e; ^  V) ^7 R
10. s; r1 A* ]5 b9 Z$ l5 ?
    
11. toc\" L1 j* Y\" ]- M# ]9 I! Y
    
12. 
    ) }1 ]' e! x# m( ]
13. s =* u' q7 P. y) _7 q' J
    
14.     6.126720600000000e+009
    & t( S8 L* |2 C3 M! A
15. 
    3 m# b0 j3 h5 b) R$ _2 P
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;, @) Z6 x! ]7 L  f5 y& j( D
   
2. mvar:
   . e* U: Y3 }* ]
3. t=sys::clock(),' z% k3 y: i/ ^9 v
   
4. s=0,
   & p\" {6 r5 Y! U* k* ~
5. x=0, while{x<=1000,
   1 ~9 k7 E( U2 b3 X+ S3 h& b9 f
6.   y=0, while{y<=100,1 v) e, b2 V5 s+ k1 |
   
7.     z=0, while{z<=100,
   8 j' y' R; T5 {
8.       s=s+f(x,y,z),
   % l/ L9 i$ D\" t- q# ~' C
9.       z+++ M$ ^1 O8 D$ E$ g
   
10.     },
    * p4 J$ I- z! C, C2 Y  a1 n
11.     y++, \2 e: L/ @0 C
    
12.   },6 U1 k5 T5 r- w  i/ E2 @# c# {
    
13.   x++
    ( E. r1 Z8 o/ y: \9 U1 n1 @
14. },+ A+ T9 ]+ ~* ~4 P4 U% I
    
15. s;4 T, s0 O% n8 G9 k
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
' ?+ ?/ }0 T3 A- v* e' w! c2 s1.109  秒

! _2 b! {3 @* }5 D二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
8 s" c' `3 j; A7 Z3 y8 v" b

hxftntA

* K# |. n! j' ]7 S: P
, s0 ^8 T: Q& |8 S


+ a) @/ n. p# H0 K3 r4 L- k

! W; B) x0 k/ I0 @) V$ ?




& l' y  {5 u& R/ W+ `' Q
0 H: C6 o4 w% o# v9 b
- P4 Q& i  g5 Z