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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
9 T1 F$ ]8 n3 O% [& L
+ B% y/ d! H6 V; w, ]0 ^: g＝＝＝＝＝＝＝＝＝＝＝＝
  W* }2 q8 Y2 d  R3 o3 b& K" F& S& m
- ^/ \) ?6 Q5 R6 [7 tMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;9 }\" G! B9 Q) p/ k- P/ R
   
2. tic;
   $ \% d/ i1 K0 O! O; S
3. s=0;4 B! Q( F+ t) {; V( P, m
   
4. for x=0:1000
   ( `( l1 D\" O, g2 P# F
5.   for y=0:1000 F( K' e\" ^( x) ^( X! B
   
6.     for z=0:100\" ^: r- K* w' U. ]4 c
   
7.       s=s+f(x,y,z);
   3 V- k  V+ }- D8 s. d6 V0 p
8.     end; K3 l4 u1 s/ N+ R; p
   
9.   end
   0 B5 a5 y  Q! @2 J
10. end
    ( U2 Z- ]. l8 {  h* ?
11. s
    - [3 X# b# {& t9 G5 g
12. toc( f3 V( r0 e5 I0 Q$ s# T- w
    
13. 
    * I2 m( P) j7 G1 b2 O2 F- J
14. s =
    . w% m4 o+ ?: G6 \
15.     6.126720600000000e+009
    3 P2 K, a* Q  g6 |# c( h
16. # ~# ?% J\" R7 _2 S0 j
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   . z2 {2 ^  e; B& T: r' l. @* W5 d6 W
2. function c=xyz(x,y,z)( w7 W. B1 M8 p! q: g1 n9 s
   
3. c=x+y+z;
   $ i/ M& p, s\" V+ X; e: o, |* N+ F
4. end

复制代码

测试代码：

1. tic;- G8 N- r( s  T5 B9 z' l
   
2. s=0;0 C* [: P. ]* j; j, L: j0 v( i
   
3. for x=0:1000
   5 r; ]! v9 T! R& E5 ]3 c
4.   for y=0:100$ r7 i' U2 A, F  J
   
5.     for z=0:100
   ( d- o4 K9 v+ s/ u: e2 C0 w& s
6.       s=s+xyz(x,y,z);\" W6 [/ Z# |% E% J) ?7 X6 B# i( n
   
7.     end
   1 u+ u. H( s5 @7 j
8.   end
   \" V) u% C4 h& |5 }! ^
9. end
   5 {0 Z0 k) ^) e8 p6 e( L
10. s& c# \* X6 L' F) t. V. p
    
11. toc
    9 W! G7 Q. j- A( h+ ~9 c
12. 
    * l& J7 h+ u$ v
13. s =: b* c- ]* n% T- u1 W- g3 `
    
14.     6.126720600000000e+009
    : K  t8 X( Z1 I, G: Y# X/ i
15. 
    2 X5 z+ U) a) z\" I* s( C
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
$ R  c* N9 y. M  [' @
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   % X% H' O\" J( i
2. mvar:( p3 Z0 h' T% |' t! T
   
3. t=sys::clock(),
   - ]+ a) D  `  C* ~/ d
4. s=0,3 F\" I' |0 {* O
   
5. x=0, while{x<=1000,4 ]* T8 n! Q' a  ~# ?- r
   
6.   y=0, while{y<=100,
   1 L' U- L3 M. ^, S4 R! P
7.     z=0, while{z<=100,0 J$ o$ m2 C3 T7 `. @( q4 o
   
8.       s=s+f(x,y,z),: }1 m  q$ ~3 g9 d! y
   
9.       z++
   ; G+ a) ~+ }# }5 Y1 ?1 B5 I8 S
10.     },
    ) e7 e6 M+ m$ F* J$ s! d: @
11.     y++' y9 k3 M$ L\" l1 t
    
12.   },
    + e* q- g) t5 X\" w
13.   x++( z5 J2 ~# b  O
    
14. },! I+ B* D5 F$ Q
    
15. s;9 }! D. H4 V9 Y2 s+ _& S% J
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
1 K1 U: `6 s3 |' v) ~
1 \4 q6 t5 u8 \3 p; j5 I2 l二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

7 y5 t+ }* v" Q9 E


& R' q" \/ D6 y, S






6 [- a4 p. z9 \% ^( ~- H5 I
" j8 B. G4 {% g1 k  }

, c# @! p8 s+ m3 ~& A4 g; t6 |