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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;4 j3 o) @2 }+ E& s, l$ G
   
2. tic;, T7 v* |$ H( V. i7 T7 z
   
3. s=0;7 g& C' k. O$ h. ?( R6 Q
   
4. for x=0:1000
   : p- z6 {$ R' m8 b
5.   for y=0:1005 O/ k) Z  {1 `  \% _
   
6.     for z=0:100- h; w\" H* i. u% w: N# ]1 Q7 A7 p
   
7.       s=s+f(x,y,z);
   $ B\" J- }1 {' w  M- O. I
8.     end9 |  X. s: _6 Z/ W) {9 q
   
9.   end
   5 G3 {: G1 D4 L8 y
10. end
    % |: _! b7 E( a+ N
11. s7 D+ t\" P/ Z6 L) ^0 |( Q- m# N
    
12. toc\" P+ \: i0 c5 |! p
    
13. 
    / S8 k3 F5 f) h7 ?) C- X+ W6 U
14. s =
    ; |$ l% u$ Z; e1 h) ]
15.     6.126720600000000e+009. V: s) R0 i' p- S
    
16. 
    . ^6 P0 c1 s\" s4 p! D4 d
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m, g+ E. K\" Y\" d- G
   
2. function c=xyz(x,y,z)
   2 b1 O. j* [2 q( n9 g  h* M7 o0 T! }5 T
3. c=x+y+z;
   * w. t& D& N6 @0 m. `
4. end

复制代码

测试代码：

1. tic;- W4 N' x. O2 b
   
2. s=0;
   ( s* w& ^. `0 c9 C. S  C
3. for x=0:1000; O9 Z- i5 R+ F
   
4.   for y=0:100
   4 F' ]4 m1 M\" j( s3 d1 n
5.     for z=0:100. `0 Z6 O- D) h& f! z
   
6.       s=s+xyz(x,y,z);+ ^  @8 G\" G! B
   
7.     end% k7 A$ w! U6 ]+ {' a$ F
   
8.   end9 v/ u8 H% `8 ?' g5 N, a1 z$ s# D4 C
   
9. end* {, q; U/ o+ s& ~6 A
   
10. s4 [6 G1 \( B% o+ i1 m. Q; h0 q
    
11. toc- j  v/ J- r0 Z\" A
    
12. & s- K1 ^( Z; `6 _& {, L
    
13. s =
    \" b+ F' H0 J  G( i) B( N8 R
14.     6.126720600000000e+009
    5 l  I; t\" |, w2 O( P
15. 
    ' |6 K\" u( F) @* x) P
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
- X5 N8 {0 f/ I
3 C+ u+ m+ R9 zForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   * ?+ Y' M7 j  c
2. mvar:  S( d( T7 o. S% P' r
   
3. t=sys::clock(),$ M' V+ `- j* b/ X
   
4. s=0,6 K1 j8 g: G* J* B! \5 y. r
   
5. x=0, while{x<=1000,4 W* ?; ]/ t' z5 u3 b6 |  l
   
6.   y=0, while{y<=100,
   - |! ?+ W\" ~. z\" K* l
7.     z=0, while{z<=100,
   + }# a7 ^6 o8 E9 O, O3 ~\" D, K
8.       s=s+f(x,y,z),: o9 N: D) h. k  N5 z6 J. B
   
9.       z++
   . ?: x3 W% w6 Z$ |\" W/ c5 O  b
10.     },\" G0 y7 T: Y9 I9 T- q
    
11.     y++
    , J5 m+ y% \( e: e8 n7 j
12.   },
    % `- [7 v4 Y3 C6 N$ P
13.   x++0 u9 b6 k% ^/ A2 b8 K; ^( V\" S7 y; f
    
14. },
    1 o. b. B% C5 U& G1 O4 z5 @
15. s;$ w, p2 L* @: E4 \\" Q2 c: n
    
16. [sys::clock()-t]/1000;

复制代码

结果：
. F. h; e6 }% M6126720600.
1.109  秒

/ [' x. o# N5 ~; @4 a) E二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
# ?; R4 Z- y2 Z" A& `

hxftntA

3 N/ A; U* D6 f' E2 E# o9 G


- }  z4 O* s/ J: D  T- t0 t# m  Y
6 ?& F1 M; A2 f' F& q/ ?
9 O; o3 F4 `0 Z4 P' N

9 V" H4 `& u' W/ Q  x. @3 ~




( n: Z6 V" [3 c, d) o' e
8 N5 ]' `: @# d