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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
  z+ O& x& T9 U) q8 e
4 \- p& M1 ^+ m7 B) h＝＝＝＝＝＝＝＝＝＝＝＝
" J: [5 S. V9 c5 j( s* K0 l
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   1 @& ~& A7 Q: O3 t% z8 |
2. tic;8 o: C! j/ G' g( ~
   
3. s=0;% S2 E+ X# y. Z( O' u2 ^/ @
   
4. for x=0:1000. g: L: Y: q# D/ a, Z
   
5.   for y=0:100! [1 ?; g/ h6 ^% x) U
   
6.     for z=0:100
   & y$ C/ o4 k. J5 U6 f# v$ C1 L. z
7.       s=s+f(x,y,z);
   . O# ?4 n7 n. N  W\" t. u
8.     end
   $ {2 B  h: p1 b' B( c, Z
9.   end
   9 l8 A- `7 ]8 ^2 `( k8 ]9 u! J; Z
10. end2 q/ {6 x4 j\" P4 [/ S1 d% Z
    
11. s5 Z4 z6 m  C, Y( |  i
    
12. toc
    5 B8 L+ h* p  B) ~  _) V7 k4 p
13. # G3 R\" @0 I3 r7 a! f0 R
    
14. s =
    ; N7 s0 I! K* [& W, S  c
15.     6.126720600000000e+009: H3 [& F4 z4 z& C1 B+ P% c
    
16. : S$ T; P/ e/ F8 A) q
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
, {* [" }: n0 F6 Ym函数：

1. %file xyz.m- }\" b* c% W4 c. f
   
2. function c=xyz(x,y,z)
   5 P+ R. u+ `) Q\" W4 G
3. c=x+y+z;0 K7 N% @5 O7 N9 q) _1 P& k$ y+ a
   
4. end

复制代码

测试代码：

1. tic;+ z+ Z6 p3 A\" q
   
2. s=0;
   # C6 K, e4 s0 U, [3 s: N
3. for x=0:1000
   & E* }/ k* V/ ]( v+ I( |! |
4.   for y=0:100
   ; H& O' I2 ~5 V; E; x* u: ]
5.     for z=0:100
   3 _/ Z3 O\" ^- p# O( e/ h0 F
6.       s=s+xyz(x,y,z);4 e0 g3 u6 X& K  q9 u4 U4 u  ^
   
7.     end
   ! D6 k. M0 \1 e
8.   end
   * E5 P! p( I1 t4 K$ i: e( M- Y: F
9. end, t7 h- K( Q* c, {* E
   
10. s9 x8 E8 z( Z% `- X3 s# V( h
    
11. toc8 Q+ D, w, ^6 S& s+ l3 Q, ?+ Z
    
12. 
    * P: a' H5 W* h+ q( l
13. s =2 M7 x& v. ~1 `; [- s$ A9 s0 R\" ]$ {& R
    
14.     6.126720600000000e+009* Z8 S8 M$ o5 r\" ?# r\" \9 Q- |6 C
    
15. 
    1 \% `( J) w. ^! s2 I6 _/ D9 C3 X
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
, m1 \3 O( T, G  F/ s: H! y0 Y
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;% U$ X8 M! Z7 s# m. c
   
2. mvar:( N1 _9 Y0 C2 t% C4 J
   
3. t=sys::clock(),5 K* a& }# h* R8 E7 U2 q
   
4. s=0,
   6 E5 G5 F0 e\" m) i
5. x=0, while{x<=1000,( s0 s0 R+ {% g/ _
   
6.   y=0, while{y<=100,
   , l1 q2 W# a2 v. g
7.     z=0, while{z<=100,) Q* G9 Z+ [* v, g. P/ `
   
8.       s=s+f(x,y,z),1 Y1 e  t- T+ J5 S  E4 K
   
9.       z++6 J2 d0 Z# k# f
   
10.     },
    & n7 [$ t; E2 u' t! O1 w& b
11.     y++, J\" {7 C: M6 A- B\" [+ a( h
    
12.   },\" |9 K. w% t9 L( r+ y
    
13.   x++2 c# B& w) a3 i: }- \1 t# Q# V
    
14. },\" e. F: h+ d4 M# ]
    
15. s;( ~: N! \* U  _, j2 w& a
    
16. [sys::clock()-t]/1000;

复制代码

结果：
% v$ s7 w: a. H. o# B6126720600.
/ w. c  L* Z' P& w4 l! ?1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
) \; a$ N$ O9 h& Q

hxftntA

( U) ~3 X. w4 h) w9 `7 `5 }
8 G' ]; K" R  |8 g$ h






$ M! F/ D; S6 f# V

$ }5 u4 L) T5 l

( R9 d7 o) R# n/ j$ L
8 |+ t- c6 m" _7 m0 N

1 h/ A( M3 G( a! }& C8 B, ?4 o