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

forcal

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

6 q$ v  ?. j9 c; F& \/ B＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;* `- q3 n7 o. r! |! B% k1 a
   
2. tic;
   ( e6 i/ Q4 V\" T% c3 r) Q7 p; a
3. s=0;# @6 w% G% N! _6 t
   
4. for x=0:1000
   * X0 c% L. `7 L0 \8 `: R
5.   for y=0:1003 f7 G/ W0 I0 |8 J$ K2 f  G$ U7 k& Z; z
   
6.     for z=0:1000 [7 a9 s: s( f! ]# q
   
7.       s=s+f(x,y,z);
   / t\" f3 i. S\" \  X1 z! f
8.     end$ n\" s( Q\" f, l6 N- Q
   
9.   end
   ) ^) K) F9 y. @7 |: I- @
10. end
      O$ G\" Q% T( ~. g4 [5 y) a
11. s
      ]6 X) p7 o) z0 c3 H
12. toc
    9 V' |9 e) s* ~. I
13.   H8 `; `6 b' n: P. j2 ?  S
    
14. s =  s+ x\" V+ Z; c) M2 r
    
15.     6.126720600000000e+009
    + L) p  ?\" a  s( k/ \- M
16. 
    9 [# z- S! W* b
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
+ T. y( H% F5 e$ K4 G. @m函数：

1. %file xyz.m- F' ]: X9 Z( |8 J. I0 h
   
2. function c=xyz(x,y,z)
     }+ u. t# Q$ M# @
3. c=x+y+z;
   : j. Y* z; {$ }- [/ ]1 m
4. end

复制代码

测试代码：

1. tic;
   * u* U, j6 K) K% Z
2. s=0;7 [  w( K4 O( t  I0 y8 w0 L
   
3. for x=0:1000
     O+ H* Q+ A% I8 V, a: H
4.   for y=0:100
   % a* D) W2 M2 F* q( \
5.     for z=0:100
   # o) d/ o% p4 D; o( j9 C) k8 B
6.       s=s+xyz(x,y,z);
   ' \# j, N. T5 z# H
7.     end
   9 h% \! U2 t/ y
8.   end
   ! X/ B# F7 a9 @5 y6 \/ {
9. end
   ! ]# y4 P8 ?0 M3 `9 ], \8 d  z1 y
10. s
    + y5 b# Q  q3 X* }6 l
11. toc1 `5 D9 d- |0 x. J- L. a\" z
    
12. 
    3 Q5 _\" a$ Q\" Q+ w0 d6 X
13. s =3 k% k, x4 j! ~8 Z
    
14.     6.126720600000000e+009
    \" S+ a5 r0 ?& s9 F$ T3 h' l
15. 
    1 P% A2 }9 U( J( k  p
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
0 ?+ z0 u3 }# t2 j# }# ~' g7 \
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   / \  |- r) K2 }
2. mvar:8 S! s; a) f+ @. A; X3 ?
   
3. t=sys::clock(),$ x$ x0 H7 s& x1 {0 n# _
   
4. s=0,
   , O! u$ g5 c8 [3 {4 f
5. x=0, while{x<=1000,
   - {. W: Z. T& ~3 ?  d
6.   y=0, while{y<=100,# e5 J\" N. o' ^: _, O\" Y3 N
   
7.     z=0, while{z<=100,8 i; g\" i% z! c% D& b\" v  L\" d. w
   
8.       s=s+f(x,y,z),
   0 Z$ J: Z7 R2 Z  l/ W* c* W; }. T
9.       z++
   / `- y; W7 l- v9 L
10.     },' G& t0 l6 e. s  E! D
    
11.     y++
    1 U; q; ?& I# U; N+ r/ L
12.   },# V: }& t8 d7 @8 z8 v; P( e
    
13.   x++
      ~6 B\" x- J- Q# ^' f' U3 I
14. },# R* S2 k; a1 @3 ^& f3 w\" `* S
    
15. s;
    $ H( ^& b) L1 l/ p( T
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
: h. S1 `6 x" g5 [; |
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

- E( p2 N+ [" W2 m& A' ^
' k& y# w2 h& Q, n+ ]3 K5 i$ r
" s1 _! P/ K1 I: b% B3 S1 V

0 f! p) O" }: F) z

  r6 ]3 U/ P& E+ H; `6 |2 U
: Y, r9 n. u; }$ L
+ B" P. J" T7 T2 ~; `8 y; C

6 z; ?! n0 a$ \! n

8 U: _( ]; p9 s