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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
( w" L* I6 U4 y5 a; r/ L6 i, T
5 D9 k4 R) H/ Q' g* J+ ~' Z＝＝＝＝＝＝＝＝＝＝＝＝
# @8 S: ~* x* R+ v6 n, T9 n
( S9 I9 j) C  l& a6 z; E4 l0 t2 yMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ) u% d- K6 j+ w' n
2. tic;
   - X+ E1 r2 F2 g% O8 L
3. s=0;+ B' r7 ?& _# ^9 ~9 r: Y6 b
   
4. for x=0:1000: k2 e* k7 F; \+ d5 M
   
5.   for y=0:100$ ?1 N+ \* z* n' W, _
   
6.     for z=0:100
   % l% v! q/ C8 `3 u
7.       s=s+f(x,y,z);# m$ w/ Q9 Z( w7 Z3 [* d3 i9 {
   
8.     end
   / Z7 A; C  o/ h3 Q
9.   end# g, ^* A; r  G0 D9 d6 X
   
10. end
    # [% y2 V3 q6 o4 n) s5 c0 c, q
11. s& }: C, @$ B\" T
    
12. toc! m' _/ B& w8 T
    
13. 9 p( ?. b\" g; D2 X
    
14. s =
    # U' F% V\" h3 G( X; G% ^- y  b
15.     6.126720600000000e+0096 R6 c1 |. u( j$ l- w9 D
    
16. 
    $ ^\" N. O& i3 O# A1 }5 R, n
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m0 n  u* d& ?& A8 K
   
2. function c=xyz(x,y,z)
   4 S6 ?% E' d0 }6 n; L) h
3. c=x+y+z;
   7 N! [7 x1 @0 B: z
4. end

复制代码

测试代码：

1. tic;
   : N7 h( l% ]& w$ a' d
2. s=0;
   6 c5 N  M0 V) d7 N( c) ^
3. for x=0:1000
   : v7 T! z4 N3 G% h2 S, l* {
4.   for y=0:100
   . L& W7 j* J3 ]1 H  t; V
5.     for z=0:100
   \" A7 X& V' J( k( K) v# s
6.       s=s+xyz(x,y,z);) H! Q1 o: n- C: k
   
7.     end
   , F7 H/ Y5 x( W
8.   end
   6 V! Z1 T9 A! t5 A# P( l8 y* G
9. end. E7 l+ I$ v1 U1 V$ j
   
10. s
    ! Z' ?; z& H4 \. b0 v4 o
11. toc! d- e' e8 J6 p
    
12. 7 G$ F0 ^2 C) n) c+ n
    
13. s =1 O, ?* o/ r\" h( `( f/ e
    
14.     6.126720600000000e+009- U2 g4 Y  A* L
    
15. 
    7 s* N$ |% S2 A, J, q8 u; c
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
( `$ K. M8 ?" v- r5 N+ w+ S! K2 g
7 M, Q# e/ K) [2 `Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;5 [5 y4 F7 K) f9 v' E* ~. y
   
2. mvar:# }\" y. h# c& |1 _; \
   
3. t=sys::clock(),
   , s% ^- g5 `  N3 ]( H
4. s=0,\" P- x\" G\" q6 x- V5 D
   
5. x=0, while{x<=1000,
   3 d\" S6 M: T3 J- @7 q0 M5 }
6.   y=0, while{y<=100,
   9 b; k\" `/ v3 p9 T# |( f
7.     z=0, while{z<=100,
   % W: \( x* Z' V) F( E% Q
8.       s=s+f(x,y,z),
   1 N8 N2 c' i3 K6 N* O3 w
9.       z++
   6 L1 Z  U* J! t
10.     },
    3 n7 m! F; x! u6 b
11.     y++
    4 P, {- K; H\" X6 n\" n* e2 ?# ?
12.   },
    ) o! a, l# M6 v/ c
13.   x++
    $ s\" L9 p/ z# Q) [+ z2 i! y
14. },
    9 G1 H$ E1 |+ m7 o
15. s;. p/ c% a0 s7 t) A: a! U
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
' E  _! L2 |! H1.109  秒
5 @' E8 Y( X/ k/ c7 F# V
4 w" x+ X- C# h$ o, u# i二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
+ I0 K7 h' v( w

hxftntA

1 n; A, R7 z$ c


: A/ ^4 W* N& y5 N

( d0 H, v6 D3 |+ ?" }9 H& ~

. F1 r( m% `! q

4 c% b  v! Z5 W2 d) R  B$ a0 M

  c- z) v: U9 w7 `0 k. l6 n

. b3 E0 t- ~, W; U# c" M
+ Y! q8 p" }' w5 e