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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
2 T1 U: O8 g1 S. I
8 D; s4 J; z1 [) k- {＝＝＝＝＝＝＝＝＝＝＝＝
" s2 r: F; f0 \& M
5 n& u& b2 w# ZMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;# A2 z/ t+ i\" M
   
2. tic;4 \: i7 G9 ?0 g' L; r
   
3. s=0;5 ]; `' ]; w% z9 L# g5 X
   
4. for x=0:10001 Q7 T; C; o0 p  x/ k
   
5.   for y=0:100/ v4 ^$ j; O' _
   
6.     for z=0:100/ C' X% q7 Y# Y+ s7 F
   
7.       s=s+f(x,y,z);$ v$ g. i& \! n3 F  h8 D. z8 }
   
8.     end4 b3 {6 X# R. V- L
   
9.   end
   ( D- B\" r7 o4 [7 O$ a
10. end  h& J$ ^4 j& J; U  a
    
11. s
    ; L% ]: }7 _: W8 z( r. H0 N0 m
12. toc. ?( e1 J& z2 J2 D/ I9 q' I
    
13. 
    & d; j( O4 W% d% i$ z
14. s =6 d9 A: H# H' h& y
    
15.     6.126720600000000e+009
    ! u\" G# \\" S% B3 n# g# F
16. 
    / T# B4 {4 a8 j* T/ B
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m0 q( m- z0 Z6 g. m0 g8 M. N3 F
   
2. function c=xyz(x,y,z)
   : Q& ^0 u5 l/ O+ M$ N; o  _8 d
3. c=x+y+z;
   ) A! f* {+ A  W8 X0 W. t' ~
4. end

复制代码

测试代码：

1. tic;
   * `! _' G5 p, e6 ~8 P
2. s=0;1 [2 ^/ O& i8 p3 F6 D6 Z6 e
   
3. for x=0:1000/ X: e/ O2 E# m2 S* z2 Q
   
4.   for y=0:100' k! S* N3 @( g2 Q! Z
   
5.     for z=0:100( i1 r# M; ~4 O5 |9 a5 i
   
6.       s=s+xyz(x,y,z);: a* U; b$ S, }) p
   
7.     end$ l\" W2 P9 }5 R, P\" ~3 g& e8 L/ k
   
8.   end
   \" W, b6 [  {\" n% m\" {
9. end
   ' d/ u3 q# Y+ [
10. s' a/ S+ t8 u. s7 X) _* ^
    
11. toc
    # S9 F$ e8 {* w+ Z  q
12. ! N5 [/ N4 z9 n& p/ C/ A& r9 J
    
13. s =4 ~9 }# I5 M1 l+ X
    
14.     6.126720600000000e+009$ V- X) f3 ]: S& F
    
15. 5 N; q, I9 U( U1 n
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
. F3 B1 W( R3 P7 V* X
+ Q' C  J$ y/ SForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;2 R% k: m$ T\" O\" K6 A
   
2. mvar:
   8 |, V( _2 P  G8 j) j8 \# q+ ?
3. t=sys::clock(),
   - M) I5 u! d5 D, ]/ `, O
4. s=0,4 t& j/ n6 D. X7 S, b# F
   
5. x=0, while{x<=1000,
   \" E- T+ M9 Z3 E3 R% m! D
6.   y=0, while{y<=100,: W* _5 I' }- N, o) ?, @
   
7.     z=0, while{z<=100,; O- K( Y; S  c4 a; E2 y
   
8.       s=s+f(x,y,z),0 j, A7 e( P0 w0 q7 t
   
9.       z++5 I& P) _4 W% o9 D
   
10.     },
    2 ~3 @# d( g- K% }; g
11.     y++
    ( J3 H: z8 ~2 S: W+ P
12.   },
    * O- t* n\" @% c2 c% ~
13.   x++
    ; C- k  ~8 |; ^% O9 S) ]/ W6 |
14. },
    ; H: u0 X: l, I: x& g4 R
15. s;; w; M$ d9 W  D+ u) i+ s: V
    
16. [sys::clock()-t]/1000;

复制代码

结果：
+ r9 w5 A3 K8 _. F) Q6126720600.
1.109  秒
, X* D# h5 y9 B6 w' L
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
7 j2 R7 R. y' a& t

hxftntA

: I. i8 a, W& L

6 ?  Q9 Y* X& g( t. h! J
( r0 y, C# Z5 h  A- l0 R4 R

8 E: e2 i4 c2 D6 p, h. @


% a2 |* k; B& [( [: U


. E) i# Q9 f/ o, z+ ~: \7 a# e+ o- r
3 y. ~: ?3 s/ y4 [0 g' a0 I

5 q& J: M+ Q$ o
7 v* j- n- Q/ C- ^/ B8 }