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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
, F; B: R( M5 L8 N4 T# G4 h5 c- B
＝＝＝＝＝＝＝＝＝＝＝＝
. ?6 Y) h1 {: `9 x5 G6 B' v
3 y. M8 x  }/ X$ C# AMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;2 A# ~7 Z2 ?  G9 [% S6 j( s$ n\" f
   
2. tic;
   # `* a% G1 u3 I2 f, I
3. s=0;
   + P5 z  d4 B8 _- u6 L/ s8 M
4. for x=0:1000
   4 `6 K- m3 N7 Z1 `5 y
5.   for y=0:1002 l8 X- @; R9 e
   
6.     for z=0:100
   7 ^8 E5 N5 n/ J$ ?* a/ [+ Y
7.       s=s+f(x,y,z);( w4 y# Z' ~9 k4 I) `
   
8.     end
   ) N\" J- \- L+ y- D
9.   end' [' H1 ?$ O7 L, Y) e9 [% [. F
   
10. end! p5 P6 k. t6 ^1 A# [, _
    
11. s
    7 C1 G- A  N\" p) q5 k\" v$ s- _0 J( ~
12. toc5 i1 v: |; k, E& @( a# b- J
    
13. 5 N2 ?+ l' H+ Q% Y! @) R4 Y1 l/ _
    
14. s =* [) C0 w5 x; i
    
15.     6.126720600000000e+009
    8 S2 g7 N( K3 [  J
16. 
    / g\" W9 |; X) ^# T
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m! N, j) u\" q/ W! I* C+ C
   
2. function c=xyz(x,y,z)6 ]/ u: B( `2 B& \. v6 x
   
3. c=x+y+z;* t* Q+ M& h- A1 b7 m\" t0 e
   
4. end

复制代码

测试代码：

1. tic;4 y7 L# V\" W8 f% g
   
2. s=0;% \1 H% n5 Z+ H2 A! f, I& u
   
3. for x=0:1000% y+ d3 ]/ b6 y* b! E3 `3 |
   
4.   for y=0:100
   \" ]: X2 a. t6 s& ~2 `! T8 h
5.     for z=0:1005 M\" G( u- F; h. R2 c( G  ]
   
6.       s=s+xyz(x,y,z);6 P3 F! ?2 F+ F5 A, ?/ \  U
   
7.     end
   & Q% f3 c) _3 Y/ t
8.   end
   0 \! S+ D. |& E
9. end
   0 _$ C: m1 J' _+ ~
10. s
    \" h. M& f5 W! W9 t0 x7 r  w& {
11. toc7 m# c) Q4 k- C: W
    
12. , c+ |4 p3 b3 j/ H) s0 V# z
    
13. s =
    6 a, D+ N7 f/ F9 H
14.     6.126720600000000e+009
    & A8 @$ v! U; m$ n
15. ( L; {% m9 m% i
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;( T: ^  {# b9 _) C; ^* W
   
2. mvar:
     x- A; A. r  c5 ]) ?; e
3. t=sys::clock(),- Z9 e: P+ [/ S4 z
   
4. s=0,% H6 a! L% G\" b: ~- z' ?9 @' T
   
5. x=0, while{x<=1000,
   . e7 s- G' `5 w6 M6 f1 d
6.   y=0, while{y<=100,
   1 `0 p2 o1 c1 m5 C
7.     z=0, while{z<=100,% W1 A, D4 l! L$ |. K+ n* X
   
8.       s=s+f(x,y,z),# A5 ?9 C' g6 W* j
   
9.       z++7 M  u+ ~  v' `, m
   
10.     },
    * O- O1 {0 h; a9 K& W# B
11.     y++
    . U, I7 ^( N\" z  T2 ~4 _\" Q$ l
12.   },
    ! a# P) O& I\" ^
13.   x++
    5 t, d3 s; z- [\" F. h  l
14. },
    . \( k3 }0 p' d+ l  g# t+ y3 J
15. s;
    ; F0 [2 F1 u1 y0 c' W6 F
16. [sys::clock()-t]/1000;

复制代码

结果：
) ]; E) N% C* @6126720600.
1.109  秒
& I8 {5 O& \) e7 {, k* m; G- Y
% ^3 p9 ?9 ]$ j4 t% l7 h二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
4 i6 f/ X5 ?) l

hxftntA

0 u& W; j' f1 l8 ~+ r" D  I6 d# Y

, A. U% R# `/ F% G' Z
% d3 u/ r% |" P' n5 u




, z, ~6 [' p! s0 J- Y' J+ S
* y& u9 a0 A5 p' p& [* n  _, \
3 Z! r3 ^' ~# s, }! K
2 q; {9 v+ o1 e" s/ w