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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

" s% b( P- U% T6 aMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   & U1 n# p8 J( L: m8 X
2. tic;2 m- Y  z6 S5 H8 E6 b( N
   
3. s=0;. a/ R\" y# H3 ~# O
   
4. for x=0:1000) R0 I) j) O$ ^$ p5 d& h7 P
   
5.   for y=0:100
   6 d( b5 t$ k; N& ]6 P
6.     for z=0:100
   : i4 D0 @! g/ z. C( t
7.       s=s+f(x,y,z);
   8 s* p0 j& D) ?- k( M\" L2 a8 o
8.     end
   2 N4 ]\" |; H% f9 m& C. Q
9.   end
   & _7 ?* H/ `* o6 v
10. end7 F$ n3 N; V\" D\" l
    
11. s
    ( D! }8 @  R2 b\" m* F
12. toc$ V9 f3 V2 s; N2 b, h
    
13. 9 ?5 Y* K' c9 {6 D
    
14. s =- S7 @9 E! E; n. l: {9 e  }
    
15.     6.126720600000000e+009- [, V7 ]; p+ ^2 y
    
16. 
    % x3 A' t) o; K7 s
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m3 M. E  A7 o2 k- F$ _
   
2. function c=xyz(x,y,z)
   ( \0 ]8 ]7 @4 T& O
3. c=x+y+z;
   % I9 e5 ~9 u( f; P
4. end

复制代码

测试代码：

1. tic;& k) n) q1 e  `6 l
   
2. s=0;; n: R+ X9 H) L+ d5 w1 N$ N
   
3. for x=0:1000
   & ]3 Z7 _% U0 D! D7 A( w0 J
4.   for y=0:100, _+ x8 V\" a; i+ b' ?% x5 g
   
5.     for z=0:100
   \" Y9 |# E8 u5 X# G* v% v8 }! @
6.       s=s+xyz(x,y,z);
   1 V! e, j# `! `3 V2 H; ~, P. q
7.     end: t4 i8 v5 }. J9 K( p# ~
   
8.   end% a% V9 y+ h\" q! p
   
9. end
   2 R% g  L/ T/ H2 y$ x5 g
10. s
    # |5 A0 V' I% S$ d) V
11. toc
    + `( S0 u. v( H4 B' N! M
12. & u! D2 |& @4 G/ Q# @
    
13. s =
    : x; S! x( z( Q  ~  }( ]
14.     6.126720600000000e+009
    6 g; L2 `) @% W! O+ h' ^1 K* x
15. , B7 S8 P1 v1 C& D( x% i; a
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;+ w$ v' g0 K5 _, `
   
2. mvar:6 @# f2 P* C& N! _# W\" b
   
3. t=sys::clock(),6 V; F* _% W0 p: ~; _$ x0 W* n
   
4. s=0,. @0 X1 ~7 _, a% s
   
5. x=0, while{x<=1000,
   ) U- V$ n! V9 W0 _1 \# B
6.   y=0, while{y<=100,
   ( D8 `/ ~) T. z8 g$ |3 h) _
7.     z=0, while{z<=100,
   ; N1 [! d\" U, G4 f\" p& s
8.       s=s+f(x,y,z),
   5 p7 |9 q# c: }2 x
9.       z++$ {1 X. h! `! w
   
10.     },2 E4 O2 ~; _7 T
    
11.     y++
    ! e6 S% |: [. p
12.   },
    + S# L! `& V6 C+ P7 j% A\" I
13.   x++1 |* ~4 J. |1 B  B
    
14. },
    1 }5 X1 }+ y4 }. _
15. s;
    - w/ ?# M\" _6 ?- S$ r' d
16. [sys::clock()-t]/1000;

复制代码

结果：
. k' n3 v# O2 q( L9 v  y3 ~6126720600.
+ v8 `- X6 |0 E+ n3 d/ j5 c1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

# q1 d: Z2 ?; c0 u6 Q' L

/ ~/ K) H! ]. J8 h( y0 ~. D


5 n7 h& w2 G5 v  v+ Z

: _$ v3 @: N0 |) L0 ~! A' s
! O: Y# s. |" @- ^/ j: f# T
8 @5 c: o# b. m& M) y, ?5 K
' W+ {, S2 {! ]/ _


6 R$ A6 n  g/ T4 b; A1 J1 [

) B* ^# Q% K( |9 T" w& _
' ^1 a5 u: }4 X/ ^& O1 J