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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
  p" R/ e) z  O# t
＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   * M' ~; e8 A8 C5 S1 M/ P6 `
2. tic;+ j+ i+ P- m; j) I0 n! F
   
3. s=0;
   ' W7 d2 f+ S- J0 G' w
4. for x=0:1000( ]% l5 G' g+ p) p' J7 C0 z
   
5.   for y=0:100
   % v- G) {6 f5 @! |/ J/ B
6.     for z=0:100, y- l/ F: p# L
   
7.       s=s+f(x,y,z);2 L2 [! P: r& ]; N
   
8.     end) x) U# o6 y! ^9 i3 Y5 _- o
   
9.   end
   , z# p/ y  U( @; v8 o
10. end
    ; Y9 R. R, x- f2 c% }$ e* V\" S4 e
11. s% K& R* U- f$ }8 Y% i, I$ M
    
12. toc
    3 S\" j8 ~) M) @: |\" w3 M
13. 
    ) ?( P! a& Q; r! P% G! x. ^0 u6 j
14. s =
    8 T* d$ n$ o  L
15.     6.126720600000000e+0091 t7 a! R+ N0 l
    
16. \" C8 }. v+ R# w
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   + R3 e. Z5 z0 S% n' U
2. function c=xyz(x,y,z)
   1 p( _/ \$ Q$ g3 t7 `# C! d
3. c=x+y+z;; L2 |& n5 ~' {$ l5 g4 y* ?# h
   
4. end

复制代码

测试代码：

1. tic;
   6 J+ R* j6 y+ T* }0 E
2. s=0;( W$ x. Z2 g. b9 j. X  T3 O; [
   
3. for x=0:10009 @. ?  C7 T+ Q* }' ]2 s  q; G
   
4.   for y=0:100: x% p7 L8 I* l, D% {
   
5.     for z=0:100, i4 _  J\" ]9 j0 i
   
6.       s=s+xyz(x,y,z);
   ( P5 `' u  {% K! n0 d, ^+ F' G1 j
7.     end* V- S\" [, {+ x4 K' l3 o! f; b
   
8.   end2 H7 V( Y/ a/ y\" c7 K; A
   
9. end
   ! d. T3 l5 }* B\" \
10. s- J8 K! m) q0 U  h$ p8 N; G5 x
    
11. toc
    ! `' F. e2 P# u
12. / d6 q5 J) A\" c/ r$ U9 z, ~
    
13. s =
    + t8 l) [( k! D, J( J
14.     6.126720600000000e+0091 I- Q/ M3 L, ^- y2 \/ E6 _
    
15. 7 h$ V% \* \- L* H1 G\" h\" l
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
6 @4 D  [, }1 `2 E4 U: C# ~- q
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   \" b: w8 y0 L. R1 V\" a; k
2. mvar:9 ]9 [2 L$ z: t9 m9 e, w+ B' ]2 N
   
3. t=sys::clock(),2 @+ G# O\" t9 y+ w
   
4. s=0,
   . q\" X1 {! }\" F7 `; ?! h) `
5. x=0, while{x<=1000,# E9 ?- a$ K6 I( y. \: m
   
6.   y=0, while{y<=100,, L7 ]1 e4 S( w* N: Y% j* B% l1 }
   
7.     z=0, while{z<=100,
   7 \# {: }' ]( E: v/ S7 H) F1 T% |& i  Y
8.       s=s+f(x,y,z),\" U) G2 O5 s\" A$ y
   
9.       z++
   6 O+ l) B/ K3 F$ Q3 P* R) K
10.     },+ V0 a. j$ D8 X9 N; x5 @
    
11.     y++
    . m1 O7 X- o( G7 U* E: R! e) y
12.   },
    4 R\" o* H. l# T\" }7 f: D: w5 J% `
13.   x++
    . A: ?6 t  w+ }! `' x
14. },
    & j0 ?\" `, K& Y9 k# m
15. s;; T$ l, B* j4 t. P' k: e
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1 l. Z6 B4 [0 ~4 i4 R( Z( K1 B# i1.109  秒

  m( J3 M: p# A# x二者效率有几倍的差距。

hxftntA

8 N# e2 I* g! R
0 a2 E# z' ~5 D# |& n) Z5 y

. c- n! b+ z( l$ t- j

# k; c1 a8 U5 [: O9 u% O
  z) Z$ Z' Q: u
* \( i$ Y% o" Y6 {3 ^5 h; x
2 U8 y; m  h5 f) U% D' c
& s3 m$ m: m" X* D  D
# \( f5 f% l3 ^0 p
+ S' h6 O! V, J" \9 Q

( p$ b" o" X7 g8 g5 K6 o5 X

ehi28

嗯,不错,支持一下.
: l* y% l; B- e& e9 d