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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
* q( G/ ]0 m0 \2 P7 Z: {
＝＝＝＝＝＝＝＝＝＝＝＝

. x2 i8 a( G: C: TMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;* j7 [; `$ P8 [\" ?: v  a% m5 q
   
2. tic;
   $ _1 x  D7 b! d
3. s=0;
   \" D; ?5 |* q0 w' N+ |- R5 B8 ~
4. for x=0:1000
     R) M2 Q% n: a2 c) k3 c
5.   for y=0:100
   : W6 H* J6 j& `% w
6.     for z=0:100' h: ?/ G1 S% \/ [4 [! h/ U$ ]5 s
   
7.       s=s+f(x,y,z);/ ~7 T& c$ U# `) D  F6 v
   
8.     end# E* m8 i: |0 n' [, h
   
9.   end
   $ F* h0 u$ R3 o\" {! S
10. end
    \" V# R0 L# P* q- }' U; I
11. s5 ~% w/ T- B! i$ W
    
12. toc( d9 K8 D. S' v\" t' m6 U
    
13. 
    ( b& j0 g- z3 d# t% @( ~% F1 l
14. s =
    - R' B0 M- N% {1 ^- z/ Q0 E( T
15.     6.126720600000000e+009
      }8 V2 N1 k' }7 X5 v4 b, h
16. * m. r! E9 I6 n( a5 t2 ^- _# g
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m* X) H' e5 m+ S* @1 Y
   
2. function c=xyz(x,y,z)) |\" I! [. h3 A: z. {- ^
   
3. c=x+y+z;( c3 Z( ~) y' o8 k- b
   
4. end

复制代码

测试代码：

1. tic;
   ( b& O! T% z0 i# e
2. s=0;7 Q5 V- y9 A/ |, J
   
3. for x=0:1000
   * ^) w- d6 M2 S% w# G' f
4.   for y=0:100/ G. e4 V: ?5 z2 O9 i& k/ d
   
5.     for z=0:100; j& G; @' f& k\" \: D( A8 J6 D
   
6.       s=s+xyz(x,y,z);6 h0 R+ X\" R! R9 ^# s$ R
   
7.     end4 J% F) d  |  [4 m8 J+ F+ d0 Y  Y
   
8.   end
   ' q& E* ^& n* j+ W
9. end' r3 w, P! V+ ]6 S$ B) R
   
10. s: E. c6 m( R8 @6 o4 A$ L6 Y' n$ g& T
    
11. toc
    * `/ N. O* N( O/ m9 J; Z
12. 
    ' w! _( J; ?% b) ]
13. s =
    3 M0 U' {9 i% m
14.     6.126720600000000e+009- s4 v7 Y& i, l$ H2 y2 M3 x
    
15. 
    ; P\" @+ ?/ r; B) |& f9 r( X
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   + N% B% [- B\" l0 f* ]
2. mvar:
   + g) R0 A: f; M& z
3. t=sys::clock(),\" T. A: r\" N3 D$ n& g: o  ~/ C1 _# ^
   
4. s=0,1 m9 N5 ?2 ^2 C) l, l  A( J4 q# z
   
5. x=0, while{x<=1000,
   3 l( U/ q4 \: Q: X+ F) Q0 V6 c
6.   y=0, while{y<=100,& y7 f9 x& O1 M5 ]  [# k
   
7.     z=0, while{z<=100,: Y8 _/ W( |  Y5 W+ e# B9 g\" l8 U
   
8.       s=s+f(x,y,z),9 V: P6 l2 ?4 a3 m
   
9.       z++6 C/ ]/ w* C& c* g0 z\" [
   
10.     },
    . P; T0 ]\" _& H: K
11.     y++4 C8 U; u; ~9 [& ?* O1 g9 b* E
    
12.   },
    + J9 f5 y! b$ s$ l
13.   x++) [. F- U4 E* C7 C5 a( x6 n
    
14. },+ j* ~: v8 Q2 y4 R. k1 i. p
    
15. s;2 j& e# ~! x& d! r
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
( w' r! J8 P% _% K5 T/ o1.109  秒
3 c! b! Q4 L# ^/ c
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
9 z; b0 T8 k0 d0 k5 J

hxftntA

) |& q# ?1 k0 T5 t$ W* p
' B/ T- X: r- `5 @
" W$ t8 [, h5 C! E1 |6 F  x, i
5 H2 S6 l. @2 ]+ ~, T. v

' k+ q, Z( s+ Q
& d) y7 a" U2 q7 \4 G
4 j! ]% X' q0 t, s$ F9 l8 [8 N
. |8 {" f0 }) J; d: T0 Y

8 Z* Z4 L6 U. Z9 O: t" L4 i# A


" |+ b: Y6 Z. V# G. Q
! h1 t: i. i( E2 O! h" A
3 d- V) f; D7 g2 ^- G, U( {: @, H& Y: p
7 o' D; V$ x8 l0 d