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

forcal

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

, U9 y  Z  H# I＝＝＝＝＝＝＝＝＝＝＝＝
- s1 Y4 _& ]8 }
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;$ Z\" i6 {# G; ^\" ~8 c& R  L0 t
   
2. tic;
   3 M! _0 i8 e) @/ n7 k
3. s=0;
   & ]7 G7 f9 \8 m
4. for x=0:1000
   5 E$ V$ O( D8 u, X! P# \
5.   for y=0:100
   $ }. {\" s# }0 ^! n
6.     for z=0:100
   1 N* @0 X$ k: p  U& f\" K9 r# |0 a
7.       s=s+f(x,y,z);
   + L- B+ ?\" M7 F: f4 f; L& x9 n
8.     end
   - B1 w+ X- Z6 ~& X+ R* x# X- A$ B
9.   end
   5 l/ n6 q$ A. S% }/ n
10. end/ B1 k. p\" A7 P6 W# p( ]\" X
    
11. s
    1 A) t% }# j4 x
12. toc
      |0 r4 G9 q! n7 f: r0 ?. W. \
13. ; F0 L; R, _+ H* v! I6 ~- \3 W
    
14. s =4 L0 `\" k6 ?; K8 w1 s
    
15.     6.126720600000000e+009& w1 C! A' M\" I
    
16. 
    $ y) A' r; R8 h8 D+ H1 _& d
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   4 ~5 a; w9 Z& \
2. function c=xyz(x,y,z)
   ) P8 \3 H; I% G$ ]4 v$ q! w* g
3. c=x+y+z;/ r6 Z\" x0 O5 Y% B7 `& w/ c
   
4. end

复制代码

测试代码：

1. tic;
   ) z. s1 n* ^5 |0 m# u& A
2. s=0;
   ' `& J6 |$ l, T! n$ j: ~) Q
3. for x=0:10001 E4 P: A1 C/ {
   
4.   for y=0:100% I2 T9 {\" H2 L4 L0 D0 U# I( w6 m
   
5.     for z=0:100
   - e: @( y9 T; M
6.       s=s+xyz(x,y,z);
   7 q1 R5 x+ b5 a3 S
7.     end
   % D' j+ D! t# h4 Q# G
8.   end3 ^1 A7 n7 x1 n% s( A( r
   
9. end3 l$ s! n* q8 L! x6 F3 k7 a2 c' q
   
10. s- k5 @8 u\" r1 y$ W* m) w
    
11. toc+ b* W) U2 |) G
    
12. 
    0 x# T' U: m- P: o  k3 g: u; b
13. s =) p+ i; {5 U2 G! b( y
    
14.     6.126720600000000e+009
    2 L0 H3 H2 B! r; S( E: q) Z+ P6 n
15. 0 n  i6 \' q; s. t3 H$ |/ g
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
3 ?6 e4 P3 H3 c" L# H/ h  ~6 l
) [2 E7 U* W6 l$ a+ U2 }- `Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;\" }. [% \\" g- g( Y  W
   
2. mvar:
   5 Q6 w2 z; l% O
3. t=sys::clock(),
   + O! D1 E* x# v\" m\" m3 w
4. s=0,
   7 u! Z8 |5 c1 |4 H0 B2 R8 l
5. x=0, while{x<=1000,\" D- @$ _5 r7 ^
   
6.   y=0, while{y<=100,: s\" k( h0 I) ~
   
7.     z=0, while{z<=100,
   ; Z- g, X% f& F\" K# Y$ x4 n! J
8.       s=s+f(x,y,z),; E* f& U+ V# X9 F7 b3 ?0 b
   
9.       z++
   ; g. u2 _& U- F
10.     },/ l. U( u: I/ ?8 B$ a2 k! Q
    
11.     y++( B# J. m/ l' ?7 i$ w& s  d
    
12.   },$ T\" @3 K: n( Z& }* e9 }
    
13.   x++
    # [' w3 L\" b9 b: L$ K, Y\" q9 p6 E
14. },
    / }: g0 i; r' v
15. s;9 E+ _5 P0 r6 H' ^
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
0 }& F2 @! d' Y% z1 ^* \1.109  秒
& i+ x' g$ {4 a
- w9 H1 _& B- {6 D二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
4 [2 j: q# X3 b$ E. F" _

hxftntA

. |$ [6 ]8 p9 I$ g9 k2 r& U. A
7 N) H# D5 g! O
( `4 k5 C5 N# W( r" P  F
" e9 g1 C  F$ @9 \


. f3 H* D4 s1 ?7 O, w% \
  E" Y$ P( Z3 S0 k8 w. A0 E
( l$ u1 l/ e: e9 d/ }


% a# C3 e/ X. _7 O


! M7 o: ~6 o  V; t5 b
. b# n; e3 t$ P* Z) B