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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
6 E' u5 G) ?) M- @
＝＝＝＝＝＝＝＝＝＝＝＝

% @* V) a1 ?2 e# n( W; M* y  qMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   - L1 l' I5 \. @0 F\" p. D& U
2. tic;. I- o, b8 x& {
   
3. s=0;7 E, }' Z. z. k
   
4. for x=0:1000
   ) v0 N2 R7 @/ Q. Q3 j$ k
5.   for y=0:1008 b+ ^$ P, T9 c6 f5 A! f
   
6.     for z=0:1006 A+ I- K1 V8 N
   
7.       s=s+f(x,y,z);
   6 m+ l8 \7 e& Y0 ^* R4 A* N, I
8.     end$ q! o5 n  n$ }5 ]
   
9.   end
   $ s' Y4 y; c6 g7 k
10. end
    2 f\" M0 p! U' O
11. s
    ! J- ^) u$ K+ U( m! ~
12. toc0 u% a6 T9 f! w) g
    
13. 
    $ [. }4 _) g\" J
14. s =
    & i. v# W- d1 ?; ~1 w9 X7 {  m9 p4 }
15.     6.126720600000000e+009
    2 P4 ^# {7 {/ `( U8 R\" x$ t5 U
16. % a. {$ q* y, T) H2 g\" Y. a
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m2 M/ z. e3 S6 `# H2 m8 U8 x4 c
   
2. function c=xyz(x,y,z)8 T- e\" ], t6 L4 O1 C\" u
   
3. c=x+y+z;; F% N+ C. Q\" X( B. i
   
4. end

复制代码

测试代码：

1. tic;1 c; L) f. \9 {& M: v\" q
   
2. s=0;$ U$ {% t. _$ n3 a: a
   
3. for x=0:1000+ v. j\" d: |\" a# M' R/ I
   
4.   for y=0:100
   $ P8 s7 Y5 V5 y# `% p3 B
5.     for z=0:100\" `5 p$ a+ p& T- J9 }- f4 p\" c
   
6.       s=s+xyz(x,y,z);/ v# j5 O: ~; @6 X, x
   
7.     end# ?\" H1 e4 i& @
   
8.   end
   . [2 X2 j: m( ^3 l, ^4 Z; C
9. end
   ; ~& C. Y8 A8 B8 n% c, t, f- R
10. s
    ! j7 r  Z3 m5 ^% S& x+ h. }
11. toc
    ; W6 r# G( b8 t7 {% N' J
12. 
    4 h* j, Z- h$ _0 M/ W1 b9 g\" V
13. s =
    6 Q\" w, }( ?( Q\" {' k6 _
14.     6.126720600000000e+009
    1 h6 ^: X' T. I# L
15. 
    . w' G& D4 A3 \
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

, Z3 g9 J  B+ |/ `Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;5 O, W# n8 ^# ^1 N  w
   
2. mvar:( y3 h: w\" \) N9 p, f
   
3. t=sys::clock(),
   ( K: n; X9 J4 U2 l\" A
4. s=0,; F1 V4 }2 X0 {  s) u6 z8 i2 \; |$ t
   
5. x=0, while{x<=1000,* U5 h$ C% ?6 d6 t+ T% U
   
6.   y=0, while{y<=100,0 T6 P4 i  \; {) a! c1 Z
   
7.     z=0, while{z<=100,- d; j( E2 F\" u- H\" |8 R9 U1 H6 F
   
8.       s=s+f(x,y,z),
   ) S+ _( p* s* E  f\" ~5 d3 o3 g. Q
9.       z++
   ' _6 }8 U4 b+ k; @
10.     },; d4 b2 z3 G& `* h4 k) }
    
11.     y++8 H2 }  @1 L# u1 p# n7 n
    
12.   },. Z4 y* M\" X8 ~
    
13.   x++6 ?6 A; s' M8 S\" J3 e
    
14. },
    $ B' D& F: S. Q# o0 H7 T, F% F$ U
15. s;$ i! u7 q% n0 I. o/ c9 C2 V
    
16. [sys::clock()-t]/1000;

复制代码

结果：
7 C" Y; o, h6 |8 W* v. |, h, b4 K9 B, U, O6126720600.
1.109  秒

: Q( M0 w3 r/ m& K二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

- a& e9 e- O8 S- G! j3 O6 r
+ |& W  j6 I1 O; R; ^, k8 g
2 r# H7 J6 P- Z+ y, |

. }' l- s6 W6 _. B
# s& B  m& I- q0 V


0 ~5 n2 @* L! O# U1 b

. f1 Z% `8 C/ n% z
( X; b: z3 X! w' j* p" x/ U
5 u  Q  V0 I. t' r$ H4 }