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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
  {7 r- e/ _' X; p5 e
＝＝＝＝＝＝＝＝＝＝＝＝
* A, W9 {; j3 N0 [6 f/ M
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;# ^8 d( U9 a, e; O+ @+ U7 \, Q
   
2. tic;
   \" ~2 e# u5 J, L$ ^- Z4 a
3. s=0;
   + f' Q( u3 e; N8 D  n* m
4. for x=0:1000
   3 U: b' {6 P) `9 |2 i1 ~
5.   for y=0:100
   9 B- j/ q7 n: g1 m2 u. p3 O2 e3 U
6.     for z=0:100# f# i4 d% A6 ?
   
7.       s=s+f(x,y,z);
   3 T0 [% a, {: ?! [7 G
8.     end
   , B8 r3 M6 P2 Y, O, P0 h\" [
9.   end
   7 S* E. _/ `* d, i1 f0 `; A
10. end
    ; D; U: n\" ^8 \  g
11. s8 R# v  A0 X; d, r3 E
    
12. toc, r# E4 \* X\" A  v3 ?
    
13. 0 M: D  r\" \3 X# n\" G: s, y
    
14. s =
    $ v! B5 g# z& i' |( L2 V
15.     6.126720600000000e+0092 u2 p; S6 [! e\" X
    
16. # U* a# p; U9 W
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   * x' j\" ]/ E\" P6 N% u
2. function c=xyz(x,y,z)& {0 u) U\" P0 K# X# Z
   
3. c=x+y+z;
   + h: R% U2 E4 M  a
4. end

复制代码

测试代码：

1. tic;
   \" k+ p+ I6 Z; R# j% G\" k
2. s=0;
   7 T) G3 h& Y& L- n3 A
3. for x=0:1000
   3 r+ t/ T2 S( |/ u
4.   for y=0:100
   ! \8 ]4 H2 X+ n- ?
5.     for z=0:100\" ?& g8 N# B- J; g( M
   
6.       s=s+xyz(x,y,z);
     t\" D- K6 H% n3 \$ k1 G! ]
7.     end
   . k3 t7 E' P+ g' N/ m( G5 I# ?! a
8.   end9 W8 C- \0 U\" g) u+ m) T1 R
   
9. end
   $ N+ |3 {' q# s* l8 [9 c
10. s
    ' C/ X6 _+ V8 t+ k+ l: s8 H2 h- K7 ]
11. toc$ ?* [& K$ b- G9 Q
    
12. ; m% S; G) O$ g+ k, ~- k; b
    
13. s =( y. m- {+ u, @- a& L- l
    
14.     6.126720600000000e+0098 z2 J4 z) u\" i8 ?, H$ J
    
15. 
    + r% s. p9 V% S3 |. k( S- ]
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

3 ?7 s4 {7 X3 l' X& s  C* KForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   * L* \$ P3 d3 {0 Q& L% K
2. mvar:
   ( l1 R9 I7 {: g7 P! U& z
3. t=sys::clock(),# `: w, m% x2 h/ P, }2 G' g7 y: F
   
4. s=0,
   ' x1 q4 K- P0 s7 \' ^9 N# w
5. x=0, while{x<=1000,* O8 v  e* C4 [2 z0 b3 _, j
   
6.   y=0, while{y<=100,6 A4 q( y# x  g' n\" o2 `
   
7.     z=0, while{z<=100,
   . ~1 c& Y! E# J8 h5 H2 j% N\" `! A
8.       s=s+f(x,y,z),
   $ `; |, {8 f1 a7 R3 p) i$ s* }- b
9.       z++
   * \! q6 F\" N0 \! F3 N+ e
10.     },, l8 d& R, y# g
    
11.     y++1 M! C& j% Z( ^: Z- m
    
12.   },
      |7 X+ ?+ O- x1 g4 H4 `
13.   x++5 P! D3 h$ u- ~3 y3 F* q* l
    
14. },
      E4 e; S! N! g  m
15. s;+ s9 [) b2 [, M/ e- C0 K/ d
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1 l  ?0 X- Q0 Y# V' f! V+ L3 a1.109  秒
' _9 p6 ~* R6 B9 B' ?$ `
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

! e" u5 X+ e& P' d/ @9 A5 M



8 W; n1 t5 |) u+ Y0 Y
) f0 _% m/ e8 r( @# h  j7 ?% J


" M: J5 M6 k! p8 O, l! l
" b/ e2 J/ S2 t

  b' @+ Z8 q; G' _8 ~