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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
+ k6 U* K7 e/ C9 h4 s
, a5 d5 E7 D4 Z; k- H＝＝＝＝＝＝＝＝＝＝＝＝
9 O, A8 {3 E% x6 d/ z% T
" j! C# [* ^/ T* O9 V" mMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;\" P! H) a8 R* L/ d$ Y0 E$ D
   
2. tic;* N( B& R$ J6 ^# @% R: F  [( q
   
3. s=0;
   ! H4 J0 L- w; a# S
4. for x=0:1000
   3 f7 l1 p& k* K6 ~* K1 u1 z
5.   for y=0:100. A( N$ [5 h; i( V, S  G7 K
   
6.     for z=0:1002 U0 }+ T4 f1 B- x+ t
   
7.       s=s+f(x,y,z);9 I2 i. K& j) h2 Y; H
   
8.     end
   8 A2 G; B$ z7 O) {
9.   end$ d: c1 `' i2 q( Y* z5 T
   
10. end
    & b) K* v2 E# A) \/ r9 @5 K6 `; y
11. s
    , H7 C2 r( b# i1 J. l6 o- s
12. toc
    . _5 D$ f% b. [$ O5 m
13. 
    8 D3 k  p\" {5 O
14. s =
    ; w! k& a. G, Z$ c( f* ?
15.     6.126720600000000e+009
    5 v. C0 p/ h% T6 l4 f- l5 W4 P) p
16. & z. R  ~& [3 j) V7 S$ \5 ]' M
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
7 Z' D$ n8 d4 U$ A# ]$ }, z- Zm函数：

1. %file xyz.m) m& w$ T- L' V$ A! `
   
2. function c=xyz(x,y,z)1 L% z+ [, N3 T3 Q* N5 p& C- S! u4 }
   
3. c=x+y+z;
   1 a! H  N, C( s' m$ }
4. end

复制代码

测试代码：

1. tic;
     P4 ?  B; D+ s/ m: X: P0 a2 D  q3 M
2. s=0;
   \" ~\" D5 q6 W  R
3. for x=0:1000
   % N2 O: p% O\" @) H6 E
4.   for y=0:100
   + j1 l) f# r  p\" }3 t6 G
5.     for z=0:100) b( x( u$ P' P\" a2 q0 K
   
6.       s=s+xyz(x,y,z);) F' V! L6 G2 [2 X$ H0 B
   
7.     end
   ) ^1 u' W8 @6 o- N, W. u
8.   end
   \" @1 G+ u/ ~. |& Y0 C7 u, A
9. end
   6 H* n' Z. k2 b4 Y/ t
10. s! ~6 B# T! c; _2 `: `7 t
    
11. toc# g+ m! ~5 g) p3 t
    
12. 
    6 }$ ]& g' H1 V1 [  N\" e
13. s =3 f6 s) H2 B/ x) i+ N
    
14.     6.126720600000000e+009
    ; u8 J* f% M2 E3 h0 n# h\" z( w
15. 
    ( Z9 K0 l7 q1 d3 T: b5 H6 \0 [
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
5 z& j0 d0 Y3 i" x% P% W& U
  \2 Q. c& ]1 _! Z$ r  r3 U; o; eForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   \" _+ t# o. x& ~4 n
2. mvar:
   ' m6 T) ]\" u' D
3. t=sys::clock(),
   3 {/ X4 ?# y  I2 V6 Y( G
4. s=0,
   : C+ [\" `) I) t8 O; G0 @! ]8 Z  S8 o
5. x=0, while{x<=1000,
   $ N( R0 h+ l1 q) y( |1 \
6.   y=0, while{y<=100,( }/ ^# m+ G7 |0 J. n! U1 q0 D' |
   
7.     z=0, while{z<=100,
   . x% _9 d5 g, o0 `5 A3 n5 [% ?& p
8.       s=s+f(x,y,z),
   $ y\" M( Q1 j+ V9 r. k: m1 I
9.       z++) F! G5 h- v\" |1 B# t5 ]
   
10.     },5 G+ f3 L8 u2 U
    
11.     y++& w' p4 U' h' v1 g8 l
    
12.   },
    , Q; D4 ]3 H( q
13.   x++
    + Z& f5 X\" n; y% M* C' v( r
14. },2 B4 k! [- Q2 `\" }
    
15. s;( l, Y0 W& ]- ^6 L9 F: y, f8 M( v
    
16. [sys::clock()-t]/1000;

复制代码

结果：
5 O' Q7 M, w# ]8 B% {+ v6126720600.
1.109  秒

! g0 o- N5 {9 M二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

* K( C4 u8 Q) ~8 l( G1 L

4 ]* X5 U( u3 Y. `; ?
; I9 E1 b; R8 `" {; [5 Q; g* ^6 }
5 z- Y; |9 t! z) d. H) Z


. u( k+ M1 v* ~* C( l
5 f1 G! I" y% A$ p" k; S) y
% \  K: Z4 ?/ A8 o4 z
3 P9 ?8 E1 M1 n& m$ X
# K. q6 l6 w1 s& X8 L



8 R, N. m1 {- Q# m: O! x! b
, ^  L$ S2 i' r  k4 ~