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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
6 _2 q7 s& d) H2 Q2 Q
& O- B/ K* V6 h$ t＝＝＝＝＝＝＝＝＝＝＝＝
: J9 Q! v" L) m+ x
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;- L2 `: x+ q4 r$ A+ w
   
2. tic;
   * e, Y5 w\" j- O, Q. {4 U\" {7 M  b, _
3. s=0;* ]8 J( h+ t: P* ?  Y' `
   
4. for x=0:1000
   ! Y( @2 |. u7 p. l; p4 X
5.   for y=0:100
   \" F! s' a! s7 D  b' r8 s
6.     for z=0:100
   8 |, ]' d: M& O. q8 Q  }
7.       s=s+f(x,y,z);
   $ ~' V1 o& j9 G: }- v
8.     end) _2 D* m) a- G5 }
   
9.   end
   5 R* r# V$ v$ B! O- o0 [( k6 ^
10. end$ I. `, Q% a: T
    
11. s
    6 _0 ?$ j( ^9 _+ A; p2 {
12. toc2 T; _0 @* B4 J# n2 m/ P& z
    
13. / ^8 @7 Z4 ^& m: G: Z
    
14. s =
    ! C# _, A! A) ]5 ]+ q
15.     6.126720600000000e+009\" w3 H) N- H+ }( [
    
16. 
    - Q5 W+ }7 o2 q& p4 S
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m6 ^3 m5 D' n! l* j  q8 u- h
   
2. function c=xyz(x,y,z)+ {; `( P\" X9 i: h% @; U( `2 F' I
   
3. c=x+y+z;9 [5 n. g  e0 a\" b
   
4. end

复制代码

测试代码：

1. tic;3 n2 n0 K1 f6 {6 M3 S  ~3 d6 O
   
2. s=0;' s) G- ~: ^+ V7 w  S
   
3. for x=0:10009 `/ m. o2 |3 J2 w$ d7 h
   
4.   for y=0:100
   $ T' p2 G1 f9 H
5.     for z=0:100
   % v8 V/ o) q) h4 E
6.       s=s+xyz(x,y,z);
   # ~0 k/ V\" y; G# k7 k) n
7.     end) r, P  }3 }4 A) E
   
8.   end* ^' j/ @3 d- V/ U% t6 D  M3 |- V
   
9. end! F2 z. {0 [- P
   
10. s
    ) b, B. m& i; v+ @8 i0 |  @
11. toc; ?; o+ [, {) P
    
12. 0 D* k9 Z* C6 P. H0 x4 i
    
13. s =
    ! Y, @2 I) j& i
14.     6.126720600000000e+009
    2 w( T6 F; U3 F& J$ b
15. 9 ?% Z8 l* J* X$ f2 Z( x
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
6 t6 x, d0 e/ Z' {. x. {) F
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;2 {& v! ?; `* z3 M
   
2. mvar:
   : n- y; @\" c8 O
3. t=sys::clock(),
   0 }: u; U9 u8 O) ^$ i% k7 J
4. s=0,* ~; ?1 A- n& g$ _; }0 e
   
5. x=0, while{x<=1000,  x& a3 L: F6 Q: W8 c& [; E
   
6.   y=0, while{y<=100,
   & c, z' ^9 \# O- k
7.     z=0, while{z<=100,
   # Y/ H! A; [4 u+ Z9 N4 i
8.       s=s+f(x,y,z),. R\" {9 W6 `( E/ `# k
   
9.       z++9 h+ }0 e5 _3 s& i! z& Y
   
10.     },
    , d& h5 N! o7 n, x
11.     y++6 M\" J! U9 O% V- m) E: M3 W
    
12.   },; u. d, W4 N6 o) h+ r& t
    
13.   x++1 X1 N' ~& E\" h\" \+ V) M! g8 [) P
    
14. },
    + J* H. F* s5 Y7 F
15. s;  d& R. ]# ~5 l% d: R' J
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

" Y7 @' g  A' K. ^二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
; N# L' \2 A' J- z% A- y7 W

hxftntA

; u8 c' \1 a, p2 U# O9 H



# B. I, c+ w; n! S  V: w





( z" N# v& ]  B. y! e5 H

! O* J" v3 u! r  u% B- y

5 A  Q' u6 `9 `" H7 K