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

forcal

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

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

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   5 n7 ?' O: y% e1 N
2. tic;
   ) S; L3 |/ P  B
3. s=0;
   9 f9 r) V% E\" I\" B; k
4. for x=0:1000& a+ ~5 G, t2 S3 G9 u
   
5.   for y=0:100
   9 }7 a* N) l; X# b  ]# A( P( ]
6.     for z=0:100
   ' O- S. Q0 U. G' z/ _1 m8 F
7.       s=s+f(x,y,z);
   ' P- F0 |4 o/ ?\" c
8.     end( Y% f/ W7 c0 }- O1 K
   
9.   end
   8 }. l9 y  [% _7 @( u5 i
10. end6 d! t. r, }: L+ q\" O3 [, Q8 ]
    
11. s- y6 C6 |1 u1 a  N; M/ n$ G* |9 F
    
12. toc
    1 U2 c& @7 }' w% g/ f
13. 7 ~2 ~2 l+ j& z% j
    
14. s =4 {$ j3 h2 V& G# u0 m8 |0 L
    
15.     6.126720600000000e+0095 b( u6 [$ w1 K* f
    
16. 
    5 \, Z3 T- a# V8 \( `3 e5 G$ B
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
+ o! R! q" z3 u: J7 m5 u% |0 Q9 rm函数：

1. %file xyz.m
   + u3 S7 p, D: u8 O/ b
2. function c=xyz(x,y,z)
   ; F. L! V1 R/ f* ~0 T  @
3. c=x+y+z;5 u$ b! q\" K) @, H+ J( S% f* q
   
4. end

复制代码

测试代码：

1. tic;
   3 \- Z& w- Q8 `5 m% z
2. s=0;
   . q3 F5 ]' u7 v- _7 X& r# J
3. for x=0:1000
   : ]% F6 o8 a' t9 m3 W
4.   for y=0:1002 M0 g' f& s7 W\" D\" m
   
5.     for z=0:100
   3 u7 v6 O5 |) A5 Y$ L/ i
6.       s=s+xyz(x,y,z);
   ' \& V6 @$ a( q# i
7.     end9 i\" s* s/ s- l3 p2 n7 p% H. W1 `
   
8.   end6 A7 }4 l\" N- \: H
   
9. end
   # q; p% B+ `& \* m% F' n+ u- b
10. s
    $ x4 Y% G! g3 m$ x2 O+ w- x
11. toc0 q8 {. n1 U0 [
    
12. 
    0 _! S0 k2 ^, G2 `
13. s =& H0 _' {6 D) ~1 @; B  A
    
14.     6.126720600000000e+009
    1 K8 R: ]- _6 f
15. 0 r- p6 F3 k; R+ i2 G
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
' k) S# D/ ^' K8 |1 k" Z
& T  ^+ |0 O( A4 q3 X3 B: LForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   9 L* ~3 b7 f, A+ Q: m! `
2. mvar:
   4 m9 h  b- [: m, D9 t. d
3. t=sys::clock(),
   * {; K0 d) }7 p, Z0 t
4. s=0,
   ( o) P0 j; k8 n\" ~) j2 _( y
5. x=0, while{x<=1000,; T  z: L2 O4 O, b( _\" I  Q# K8 Z
   
6.   y=0, while{y<=100,6 Q8 \/ T! |& ]' M# L0 I
   
7.     z=0, while{z<=100,; z\" Y. G: D6 l2 j! O5 l0 m
   
8.       s=s+f(x,y,z),
   ! ?1 I/ a! |3 F$ x- W
9.       z++/ P* o0 q& h$ b. p
   
10.     },
    # b& @9 U% }! c# F) j) J
11.     y++\" S# z  k, `7 Z\" l' o& n# f
    
12.   },8 J/ W$ Q, h+ }6 w\" C
    
13.   x++1 w# D0 e0 |2 i) o' Y6 D  p
    
14. },* C! e) e  r, \  m( a$ }5 {  a
    
15. s;
    3 l! m4 t1 Z4 @/ f) |
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
2 J# g/ _! n' q: o1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

* q! g) X3 s- l1 U2 M. U: W! s


3 L: y, U' S1 a. R& _8 ?



0 z7 o4 J* `. N" [% u9 s
, p, {  k; n/ U4 f
  o8 F  o% L$ I! m0 y6 G: y





" R* J( Y4 c3 w3 q" B" x$ U