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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
7 L) K0 b, Y7 x' z- t
0 B  Z+ B1 K* p4 R7 S1 `9 f＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;% A  p; ~: ]+ L% `; N
   
2. tic;
   7 ^- G5 P( r5 r# C* ?
3. s=0;\" u0 S' Z1 u7 O4 J; u, ^
   
4. for x=0:1000% y1 G3 L$ E0 u% H! p% i6 \) X$ r; N4 u
   
5.   for y=0:100
   5 W9 M3 _6 I1 E$ S% b
6.     for z=0:100. p5 q! h9 S2 ]% A5 j
   
7.       s=s+f(x,y,z);
   \" b7 x+ e# v) x$ j4 p, \$ D
8.     end7 V+ {$ o. y3 w  h4 S' ?$ ?\" k7 h0 _
   
9.   end
   * @: E( V4 D  W, ?; q5 V% g
10. end/ a! H' s0 \+ o
    
11. s
    ( J9 c0 h* b6 d7 V
12. toc, k- @5 i- {  g4 v, X) T
    
13. 8 D# B. V  Z, n; g
    
14. s =, l* ^8 U; S6 H4 T) k
    
15.     6.126720600000000e+009- @8 J& j8 F1 T: A/ o/ z/ L8 ~
    
16. 
    * X  @3 R5 ?) c
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   $ ?1 [% C2 [/ `+ a; |
2. function c=xyz(x,y,z)7 q1 Q/ U1 b) W% H' R  L2 C+ d
   
3. c=x+y+z;8 S$ _5 R; C7 i2 F9 S' g
   
4. end

复制代码

测试代码：

1. tic;0 V- C1 o/ _, [# d) _( Z
   
2. s=0;
   - f, W% J5 l# u- W! V
3. for x=0:1000
   & f+ ~. y; Q5 a, e% K6 c
4.   for y=0:1007 s. ?' K* ?( {% [
   
5.     for z=0:1005 ?! \) @# m9 N; E' D2 o
   
6.       s=s+xyz(x,y,z);
   * b$ D% [6 v6 W6 |' G) _& U( c+ ]
7.     end
   4 P! V. F+ {2 ?7 U3 K* ~\" [4 R
8.   end
   ( L2 ?4 M3 k! @
9. end- L! q0 J4 Z% S4 d( {3 a  k
   
10. s
    \" r, R; X8 c\" o7 n# o+ h! t9 j\" t1 ~
11. toc, I4 }( v2 |6 v: A2 _) V& w! e
    
12. # @4 m- O4 ~* ?# J+ a) E- B
    
13. s =\" Z2 g$ b% u: M# n% P3 e' P
    
14.     6.126720600000000e+009
    ( c' ?# Q( o% y% l
15. 
    7 r9 i3 O; ^/ F- s. c1 `
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
  c* D; Q+ k5 b# [0 L$ b
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   6 l% T( G$ M0 }+ V
2. mvar:6 c) ]; {. R/ v3 n
   
3. t=sys::clock(),( J3 D7 [) @2 B! x  h2 c
   
4. s=0,* L\" I8 m  G: h\" n5 N) G  c; C
   
5. x=0, while{x<=1000,
   1 }, u4 L& i1 r' @
6.   y=0, while{y<=100,
     y\" J, ~. @* z. j0 i7 d* E
7.     z=0, while{z<=100,) R* y( @2 U4 i+ W  K8 X% b
   
8.       s=s+f(x,y,z),' x, R* O- `! V) Q6 Z
   
9.       z++
   $ k# P1 f; @4 }2 I) V  W' F5 s
10.     },
    , P: U% f0 G2 [
11.     y++& b2 `' @4 ^# L  T7 b' p/ L  c% `
    
12.   },% N- k$ _, d8 `- E2 m8 V8 a
    
13.   x++# c\" g; G1 J- X
    
14. },8 S* d* e' |\" k! t/ Z
    
15. s;
    - ]6 Y, y2 k5 n  c) I3 R
16. [sys::clock()-t]/1000;

复制代码

结果：
3 m* {( D# d: q% b  O6126720600.
2 y  b: s4 Q* \- G1.109  秒
) O/ L: O+ L' Y0 Z8 j
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

* k- q; O% L- o& z5 x0 q7 t: V






) |1 H8 _& L: {' l. Z( Q/ o$ \( q7 n
3 l" Y' W0 B! X  [% D( a) r+ h
2 I. z" [6 B7 Y. f# s4 U5 ?

$ `$ C7 R; |2 F+ X2 I/ Y
, d! n- q7 D3 c  V. W' m: r7 |

. Q0 l* G' J. Q) Q1 ~) e