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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
6 P& n; b! P- E4 j- W4 _1 @1 C
＝＝＝＝＝＝＝＝＝＝＝＝
" O. Z6 ?( M9 C5 O0 l( q* U
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   8 E9 @  j9 p& h\" m4 k
2. tic;+ c\" Z7 a$ P; }\" u! B6 A9 g9 y
   
3. s=0;
   ) D+ r: ^! K2 ~: D
4. for x=0:1000
   9 ?! x* a; x! G( P: y8 P  y9 U
5.   for y=0:100- r# }: Q6 m4 Y; n
   
6.     for z=0:100
   , \2 T9 j\" d% e/ q9 {% z6 q
7.       s=s+f(x,y,z);
   . x\" r$ b% G2 j
8.     end! r$ ~! x) G0 t1 F
   
9.   end' e' ?4 Z! m\" U
   
10. end# W/ P6 j: L) s6 S' ~: s% J# e) b) a
    
11. s
    # X! W  T( g8 P% ]6 Z
12. toc/ v( n& h! _; W, e  K
    
13. / u+ X, Z( O, [3 I
    
14. s =
    ( |6 S. f* c) l& w: G& R
15.     6.126720600000000e+009
    # C; ^5 Z4 s: j! o+ l/ l9 K
16. 0 f0 I$ K4 c+ T
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
8 F, I# v3 l: r' g7 qm函数：

1. %file xyz.m. O/ j1 a) {3 o# N! }
   
2. function c=xyz(x,y,z)
   ' o& Q5 \0 p+ v# F' B5 v4 d8 t
3. c=x+y+z;
   / N1 h7 Z/ b\" |2 o1 c
4. end

复制代码

测试代码：

1. tic;' d& o\" x* h! @
   
2. s=0;7 [  o9 T\" F. N; a2 B5 J
   
3. for x=0:1000
   * v( ^+ W1 ]2 J9 K, Q& r  G
4.   for y=0:1009 ]: X5 R% p  K$ a\" T4 q1 l3 g
   
5.     for z=0:100
   7 U+ i. J+ q3 R& t/ r9 Y
6.       s=s+xyz(x,y,z);+ w& W\" Y3 h! C3 v5 A# X
   
7.     end
   # f. j( H5 ^3 y7 _9 g
8.   end
   & z/ ?- x8 ^8 f# D2 e( x4 M
9. end
   5 Q) [3 U) R9 ~0 q* ?' p
10. s\" [' o$ y  b3 \( @; k; e8 I; Y
    
11. toc
    9 J2 A5 u' {: u1 f
12. : a4 y) @* U- ^# i- y+ u9 n9 O
    
13. s =
    : h6 h- Q$ k+ F7 S
14.     6.126720600000000e+0096 G( ~- S( B9 S5 g/ Y' X
    
15. 
    % [3 n* ^1 x0 b) f& B6 ?* O. b
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
/ s2 d9 Z1 H. ?: v! _
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;! F4 M# ?8 c$ m) h. T6 M! ~
   
2. mvar:9 }6 U, W* O9 J( }& B8 O: R# r5 d
   
3. t=sys::clock(),6 l8 L/ r; w( x0 ]. _
   
4. s=0,
   1 b4 d) l7 @/ w) B9 l' r1 c- O$ r
5. x=0, while{x<=1000,( q1 P. ^9 H' N0 [0 x* f
   
6.   y=0, while{y<=100,
   1 L1 w  W; t; m
7.     z=0, while{z<=100,/ h: e2 @% S2 n\" [2 \3 D$ X
   
8.       s=s+f(x,y,z),' Q- S( r- _+ W( I5 [9 ~
   
9.       z++; Z* e3 j# H4 V7 G
   
10.     },
    - e8 T$ J! I$ t; n( H& N
11.     y++
    2 R- a4 B: x5 Y# |4 x2 D0 K0 P
12.   },6 I- h) Q7 O9 q4 J
    
13.   x++1 o: g; u0 C3 t# {
    
14. },2 K# K/ P& M0 s2 }2 o1 o: |: q: @
    
15. s;8 \\" p' u( v+ F
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

; p% R2 e8 J3 n, c. S1 Q
8 `9 T8 U* J0 Y, L! U# A& E

: D! m) j0 h) K

- u7 `6 k# `4 ~% o
# k& M8 B! Q* j7 k
/ M, B8 y- m6 \& M8 b4 A
! e+ u2 M# I& w

% y4 q2 @0 k0 g- F0 A- f3 x
- }, n$ J" \  d. x0 d! L