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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;0 V, T1 _5 p  w7 @; W& Y
   
2. tic;
   0 i+ k+ ~$ c3 a8 R$ y
3. s=0;2 h: \, x) x6 [# B! ?\" K3 X5 X
   
4. for x=0:1000: _: Y; Z6 K/ p\" p4 G- V/ z0 B\" G& h
   
5.   for y=0:100' t! O  p- p8 B- g* G, ?
   
6.     for z=0:1004 O9 ^; R& A3 l$ i. g
   
7.       s=s+f(x,y,z);
   5 D: d; W1 Q( H; r1 j
8.     end6 _0 n  z+ n& p
   
9.   end( D: M! i% b9 ^* L/ Y: D
   
10. end1 a\" e5 a  Q\" n& K, d; A. t
    
11. s' O3 u1 \9 W/ ^. K) U
    
12. toc
    ( j, Z; y6 D$ G
13. ; L/ L2 d\" F- x7 k
    
14. s =\" U/ H1 U: {- n) @; s, B
    
15.     6.126720600000000e+009
      Y1 B4 W6 D5 Q& J  N
16. 
    ! K5 X. b9 _* M
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
+ a! r# m3 O  q% ?m函数：

1. %file xyz.m! K+ U/ F\" ^+ T% L0 x$ v- P1 f
   
2. function c=xyz(x,y,z)7 S8 w/ T# m7 n, E9 T# U! G
   
3. c=x+y+z;& ~/ L6 j6 t2 X
   
4. end

复制代码

测试代码：

1. tic;
   \" a3 {) R5 p\" ?5 X
2. s=0;# R\" O2 ~2 `( N6 \7 |, A
   
3. for x=0:1000
   2 W8 i1 g3 j# k% m% m7 ~; b) x# ]
4.   for y=0:100; o3 }4 Q1 J& x1 [* s% C) |  D
   
5.     for z=0:100
   , r' ~- v. s9 {( l' k3 g
6.       s=s+xyz(x,y,z);
   # v/ E  o5 x. F% w7 Z  b
7.     end
   # b1 o: A+ @- J4 U; A1 q4 e
8.   end
   0 P/ a. P7 Z1 p0 [
9. end
   ( q) o, |5 f8 Z& a
10. s2 k& Z1 h. n7 e
    
11. toc. X  A  S+ P4 B4 _6 ?
    
12. 
    $ E1 s: y% h0 W- K$ }7 c
13. s =7 [* E- o# K) S1 L: B
    
14.     6.126720600000000e+009) V, ]* J& g% E  e6 D
    
15. 
    , G# Q( ?\" m# d5 }* i! V
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
& J8 H# d( P" e, x+ K( z8 f+ e
: E' ~' J- \+ v* o) gForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;. H, d, z; M; n. E
   
2. mvar:
   + m5 c6 B2 e: K8 G. Y
3. t=sys::clock(),8 b. o* z3 q- c/ D% y& ]
   
4. s=0,, \1 }0 X& ]- f: V& Z: r\" S( |
   
5. x=0, while{x<=1000,
   1 M, x/ U  H3 W( [4 T
6.   y=0, while{y<=100,& r7 G4 c, f. K, \, l7 o6 [% q
   
7.     z=0, while{z<=100,
   : U7 V6 a# H/ o+ ~
8.       s=s+f(x,y,z),
   ( ^/ |/ W/ l  [* r( h7 }
9.       z++) L' B5 j* U' [3 v+ M
   
10.     },
    $ p$ C# S+ b6 h: Q5 u/ G3 d
11.     y++* T! ~) o8 o5 D/ t( R% D
    
12.   },
    # F\" n5 W: X, W- D2 r# ^
13.   x++
    & K3 Y& a$ Y( l( Q% F
14. },; v& p8 ]* O! M\" Z0 o/ ^6 R
    
15. s;
    : B. I& z  d\" v
16. [sys::clock()-t]/1000;

复制代码

结果：
& L( c* K0 ~) t- s9 A6126720600.
- r/ U4 Y+ P+ D& w2 D1.109  秒

: b2 l$ M  L8 q: m. q/ }3 F1 U二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

& C7 m; Q7 R- L' s

! M3 S: w. `' f! V+ ~3 q

8 k; p. a" j+ k8 X


2 S& t" I  {: A6 T; ]
" Z% M% n' h! L0 x
* @/ |( N2 R( u% M% l
3 q4 E* g" I* t6 t" q8 l8 l