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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。

( n. G9 r. m: U9 ]% K  d＝＝＝＝＝＝＝＝＝＝＝＝

1 G) V0 ~7 V# b; gMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ; z\" ?% u( k% l' f3 B
2. tic;
     v: v. J  ]' T4 |& h* ^% j0 j1 P
3. s=0;- H3 l; ]* \( X
   
4. for x=0:1000\" ?6 U! e' a' V' ]# ^+ Y
   
5.   for y=0:100
   \" V5 j8 p; h1 R1 X2 C2 e  G\" h
6.     for z=0:100\" k6 y( ]& T7 z
   
7.       s=s+f(x,y,z);
   ! i- V* G8 l& m; V( `
8.     end
     [7 z) z8 @- |8 H; g: [\" m
9.   end; r( P9 ]8 m# A- x0 Q  I0 i
   
10. end
    \" L* F* M3 k/ G
11. s3 ]0 r' k. \+ T% @' X$ N6 V
    
12. toc
    ) K) n\" T$ c. F5 y  N4 T# [
13. 
    ) b\" }8 j. K) k( V
14. s =$ s2 z$ \0 c3 V! l\" v
    
15.     6.126720600000000e+009% u! `+ G9 Y4 A6 }6 A
    
16. , H# [: _/ M1 l8 p, H7 E, [& u
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
, O. U0 ~* X3 e( J# `) x9 @m函数：

1. %file xyz.m  x+ E: _; A' t
   
2. function c=xyz(x,y,z)\" ]5 P/ O0 k* `
   
3. c=x+y+z;
   4 g3 i$ B$ r' E* a! E, K$ _# I
4. end

复制代码

测试代码：

1. tic;
   2 Y\" o0 X4 D% s\" {
2. s=0;( K5 V' d, x' H* p& g! n
   
3. for x=0:10008 }: F/ x4 Q6 _+ {- R9 r4 X* F
   
4.   for y=0:1002 F7 ?& E% }, Z6 S- ]. J$ s
   
5.     for z=0:100
   . ?2 x\" d2 a\" P0 S
6.       s=s+xyz(x,y,z);
   ) x' D: X4 W% ~4 K+ {& P0 m( z3 f
7.     end
   . s! N! s: {1 k\" @6 I
8.   end8 k# G6 O: @& X8 j$ f
   
9. end
   \" H, |5 w, y& E% g6 _# q3 b
10. s3 |1 o+ x' }7 z$ {7 v8 g$ n% L
    
11. toc
    1 w\" @# e8 r7 W2 e; M; _( B' ^
12. : v. |9 R. H, N+ y1 n; ?
    
13. s =
    7 U3 U$ l% j( f7 f: x8 o
14.     6.126720600000000e+009
    5 B7 y1 j& \, W2 W: ]
15. 
    ! c& A% ^; ]4 X1 h\" R$ T
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

% z: Z+ [; Z1 h3 |) \( SForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;  p% H0 t% d# d% \( p
   
2. mvar:
   3 Q' k/ i, z! u7 M. L\" B. p8 Y
3. t=sys::clock(),9 {! F+ q$ j7 @% @  n
   
4. s=0,
   3 y: g0 s% G, b/ ~2 r, R) @6 l2 u* a
5. x=0, while{x<=1000,
   \" r- `2 A$ }- @  K
6.   y=0, while{y<=100,6 \\" S( G/ `  y\" D: R
   
7.     z=0, while{z<=100,0 y+ }$ B4 |; w$ b) F! i3 w2 X
   
8.       s=s+f(x,y,z),' u( n' w; R7 h: E5 m1 ^/ W- e
   
9.       z++
   6 [' ?$ K4 e, j
10.     },
    3 {* z9 G1 h7 I9 B; v3 G: B6 ?' j
11.     y++
    : S% m: [# s; a2 s8 S% M- y7 D2 y) h
12.   },% R6 u& a$ F& g( [& f0 w
    
13.   x++
    5 U\" t3 f, N! w
14. },
    6 u9 X& g: K5 P! u( y4 B
15. s;
    5 `1 J9 m( {9 ]# C
16. [sys::clock()-t]/1000;

复制代码

结果：
  h4 T& v" K2 ~; Q7 G6126720600.
& A( A* \3 ~3 I3 \8 `9 N1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

* o. T. z; \) P* ]5 e% V


9 b. y3 M; L+ N# s6 [5 S
% C1 o4 F. ^" ]9 E8 h& W5 @





( ]0 |% v. g6 E/ A$ F1 g
9 m) y& s0 B0 u& A) ^7 t) S
5 \2 n4 ^7 l% E; y+ y6 m8 j


4 I/ d' R6 Y, [$ ^& |
0 [  ^3 K+ V* P$ m