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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
8 l& W- v2 A3 ^  k/ K8 m) l
, `0 \, h3 _  C" O2 P# g＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;1 j. B* I0 e' O4 F\" b6 k& z2 s
   
2. tic;
   - G5 w1 x$ h  l( \8 k
3. s=0;/ I0 g0 b2 C- Z, U
   
4. for x=0:1000
   7 F; M1 V8 G4 y2 u2 k) C& X
5.   for y=0:100
   ! G& }0 }8 V, Q6 Q$ o
6.     for z=0:100
   9 F' e4 z0 d& ?2 [
7.       s=s+f(x,y,z);
   ' C5 i% ^6 D/ r. Z& S& V& B! Y
8.     end6 i6 B: u: z* B+ j
   
9.   end
   5 d7 ?7 u5 s! w: x6 u
10. end
    5 \; {8 D% N5 i6 J
11. s
    + n, x5 I3 A/ W7 ~8 g
12. toc9 Q  Z/ e# i+ m0 `& c
    
13. 4 P; [5 U! y/ d
    
14. s =
    + a, M3 ~  c5 l+ m- B
15.     6.126720600000000e+009! n6 t/ D9 T  U& A  S* O) K
    
16. 6 w- U' X& _- b) y
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   3 G; G. [/ a; e6 M5 U! {2 \
2. function c=xyz(x,y,z)
   / Y6 F! ^- f  t
3. c=x+y+z;5 H8 ^+ W2 m- i+ V7 U* S
   
4. end

复制代码

测试代码：

1. tic;$ {% P. }0 _+ Y( E
   
2. s=0;- t4 M, Q\" K0 i; j6 y8 R5 U+ v
   
3. for x=0:1000
   3 |6 t! p- J: |9 B! p# _/ c
4.   for y=0:100$ h$ c3 v- c, N+ o  S5 F
   
5.     for z=0:100  g6 H/ a0 S+ V' H1 ]. J- j
   
6.       s=s+xyz(x,y,z);
   3 q' k* b7 f; j# ~
7.     end/ Z0 c\" }6 s% |) T+ W3 y, O3 |0 A
   
8.   end
   6 i8 R  W8 P  c8 @* \) ~, u
9. end$ o8 e7 E; V- z, ?\" c$ c9 {& }5 \
   
10. s8 ]3 T1 l& B5 q: W\" m\" n
    
11. toc( g( O. e) T8 W; f, i6 D+ H
    
12. 
    , `4 [# s' s/ D) w& W0 k% l; S
13. s =6 `. h( j5 P$ i( Y\" w$ m& B0 t* x
    
14.     6.126720600000000e+009
    ( v( ^8 u. h: }8 v
15. 
    6 {7 Z: K+ X% L4 Z6 i\" ~
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

# O+ f! D: Z! z( M2 ^  ^& eForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;) ?; Z8 q2 r( r
   
2. mvar:6 t4 y- [; d& p( c
   
3. t=sys::clock(),
   6 `8 |' l3 a- @- I7 ]
4. s=0,
   $ J9 T. R$ ]: o
5. x=0, while{x<=1000,
   7 T( k% c' a3 ^# V\" Q+ ]
6.   y=0, while{y<=100,
   0 g4 ?2 Y1 e: ]# X
7.     z=0, while{z<=100,4 F& W1 O- Z! d( m3 i\" O' e\" \
   
8.       s=s+f(x,y,z),
   - H2 q' T: R; b: b
9.       z+++ f5 ]! `% _) C, [* ]% M( t
   
10.     },' _1 b1 ]) `9 j4 P) S7 D
    
11.     y++
    % i3 b+ J7 k5 h: V# {2 L- O( d5 G
12.   },5 @; y4 W9 I2 r3 i1 i& w
    
13.   x++) @0 ^8 t\" j& d6 e\" }\" ?$ d
    
14. },
    2 u\" J  q' o+ P' H
15. s;6 C\" `1 v) R  _3 m5 M8 e
    
16. [sys::clock()-t]/1000;

复制代码

结果：
- K$ f' t' w5 F+ H- e6126720600.
1.109  秒
: ^& I* I- E- _  p5 ]2 p* m7 x
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
! T+ u7 W; ~& A. Q7 L

hxftntA

  h1 k& V5 j. k& k1 p# q
1 y; a. p$ h: `
, z7 s. I$ r6 l' R0 [" C2 V
% ~( ^4 N7 d9 U



" P( C7 ^% K: v* u9 `$ R
: y& O) x4 {8 P2 Z7 ^7 g) d3 [: n
3 ^, ~6 j& n( `


$ ?+ T, O% M0 _  @% n' k' ~

7 w3 i( Y0 N0 q9 S" K$ S7 G  Y, a