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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
  `# P7 u! m( p2 L9 \7 q* m
＝＝＝＝＝＝＝＝＝＝＝＝
. z5 R6 R1 G# l# \! S9 ~0 R7 I
; N  I0 t* _: ?, r; ~! M4 x' vMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;7 u  O' B8 q2 |6 X  e
   
2. tic;
   2 }( _: c  J0 R. I6 F4 X; y  B
3. s=0;
   ( ?( H! F' `2 \( B
4. for x=0:1000
   8 B# t. u* t( y* L- z/ d$ [7 z
5.   for y=0:100' q% o- Z) |, A
   
6.     for z=0:100
   7 B& k; j0 G4 S* \
7.       s=s+f(x,y,z);: G' b# j3 d6 E8 w* N2 Y
   
8.     end
   1 M3 N7 b: X/ A0 k
9.   end$ m& |% z& p- [: D& e' o  T4 p$ q
   
10. end
    5 C\" \; y8 N) z3 \
11. s, @% \7 x) U! @+ ~& Q, B: \* |0 o5 [* c% z
    
12. toc
    8 Y5 J; {( Q0 Q) [
13. 
    5 ?+ d7 `/ o; ^$ f# T3 r
14. s =) N# g8 R6 B\" r, r/ J# ~% v$ r1 E4 {
    
15.     6.126720600000000e+009
    6 m; y& u- T, y) k1 `' ^9 F
16. 
    4 Z. O* C6 H7 b9 ]$ b' U
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
6 ]- r6 k! `: ?* om函数：

1. %file xyz.m
   9 W; D- w, V' a! C; h9 T* ~$ w& t
2. function c=xyz(x,y,z)
   # V# x) u( J0 z6 e# X% T  }\" J- m4 d# ^
3. c=x+y+z;6 A$ I( u2 E2 h  x
   
4. end

复制代码

测试代码：

1. tic;
   ( C0 V- P5 R7 S
2. s=0;
   , h( o* u2 m; @
3. for x=0:1000# H+ e6 s! I: B3 h7 {
   
4.   for y=0:1004 K- i7 E\" U. U: l
   
5.     for z=0:100
   % V* P6 K3 F- h. B: [. y
6.       s=s+xyz(x,y,z);) c4 {2 j- B1 L! g* w4 R
   
7.     end# r, J5 y, h7 o\" u3 h9 b7 H
   
8.   end2 @2 S8 C* c8 M/ Z/ y( X8 m
   
9. end
   . V' m; ^- j4 ]3 z- u
10. s5 ]1 z6 ?4 c5 S0 U8 n
    
11. toc  g$ E3 l7 [- j  n% \% X
    
12. \" V# q! v+ _) Q  L' @, S, U& g
    
13. s =7 l- U* a+ Q/ U, `
    
14.     6.126720600000000e+009' F& Y( u3 c: ]# a
    
15. ! I; S0 G& H: e& g- w& X
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
! B6 I7 h3 F2 S4 z% }
# q; ^: }7 w# L& `Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   & x$ l8 N7 _3 m8 O+ a\" Z  c: ~
2. mvar:
   \" C6 P8 v\" ^6 K3 U\" Q1 r- i! {
3. t=sys::clock(),9 H3 R8 H\" I. o: I
   
4. s=0,3 f+ {. l8 q$ @  \* [5 a, |; n
   
5. x=0, while{x<=1000,
   6 {# v% z% j+ l\" h5 z
6.   y=0, while{y<=100,
   7 g. i: ?* r, J- w; u
7.     z=0, while{z<=100,
   : e* T1 b\" E8 i
8.       s=s+f(x,y,z),3 j\" C4 u( u+ [) k6 n- V
   
9.       z++
   - h+ J  V/ T( }% R
10.     },1 y( b! ?5 K7 E\" n* f
    
11.     y++! d- K- i6 x- w9 V4 w3 m  K
    
12.   },: W  l+ P; j8 @\" z. _- s# x
    
13.   x++; e7 W7 E: e+ |
    
14. },( Q8 a6 l, N  u
    
15. s;$ y7 U; H- o4 I\" _3 [
    
16. [sys::clock()-t]/1000;

复制代码

结果：
+ g9 b. ^  v* L4 o2 L4 e6126720600.
$ N, L" a( h, l" {5 H1.109  秒
+ S" t" {; d8 K, f
) e6 f( h" ~0 v9 u二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

. [4 u8 x6 l3 H
4 t% Q! M; {0 H8 @- v9 q

( Z$ x& C4 N: G9 h( s* C
; i( E' n! v: S  D- b. r

+ c3 I: I" Z% T  u5 n

6 A/ L' I* V9 f& K

/ I" F8 ?* B4 l( d9 N


/ c6 K% O/ ]3 R2 R. A) U

3 k" P5 R+ N# j: Y# Z  G' C