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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
% R# n0 a+ n0 y$ w
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   / t7 Y+ n; Z9 H. H
2. tic;
   6 ^# H. L/ ?0 e) r$ u' M8 q/ @
3. s=0;' v% I$ S8 k9 x, {7 ]  {
   
4. for x=0:1000* H* ?5 ?* X8 e0 W6 f$ g
   
5.   for y=0:100( @% q! V1 v( Z
   
6.     for z=0:100
   9 y/ X. E; ^8 ~+ B' E
7.       s=s+f(x,y,z);
   ' [8 e+ c: c5 _. B1 n  U* m2 [* ?
8.     end
   * F1 {4 y4 d5 Z; |\" b
9.   end
   ! j8 }3 \3 Z4 u1 T\" L
10. end
    3 [/ j, |3 g) p) q) z
11. s
    % R7 u( f4 m! z/ z/ C
12. toc
    8 c2 |4 Q( I& u& M/ d$ p. g
13. 
    $ J) t9 `* P( R- g+ y' T4 {# d
14. s =$ _/ Z& C3 Y5 F
    
15.     6.126720600000000e+009
    , u9 s' v6 |- a* L( p
16. . b4 ~' _1 K+ B' X1 M4 k
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m0 {6 i, ]6 o6 k# Z- _6 T  q( l
   
2. function c=xyz(x,y,z)
   4 J! K; u  o9 u! o0 U
3. c=x+y+z;
   , h; D! P6 I! M
4. end

复制代码

测试代码：

1. tic;$ m( e) @( c  F* o
   
2. s=0;
   0 `9 b( V+ V' j0 |% F( s5 a5 d# O
3. for x=0:1000- h# @9 ~8 M; N/ F( Y. f
   
4.   for y=0:100- p/ }# v, S8 F( L
   
5.     for z=0:100, A* X/ n& l9 t, }! ?) r
   
6.       s=s+xyz(x,y,z);
   - b3 _% p8 n; D/ t
7.     end
   ; Z- S1 F, X: {5 Y- ^
8.   end
   0 v0 h\" l6 q1 G. ]1 ]' G
9. end+ _8 O/ _- H4 v  }- \/ a
   
10. s
    1 a. z3 X, _  n! n9 \! J
11. toc! i5 X  s8 o- W, {! m5 E
    
12. 
    6 r/ K  a' S* x8 y6 c: }
13. s =
    , I- }7 p) N3 c4 p6 e2 D# Q6 ^& I
14.     6.126720600000000e+009
    6 d8 i# x4 `! _0 p/ Y9 C. i
15. + X/ l; {4 b0 x
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
) S( Q) }1 s' ^- D+ C* j' U
9 x1 K! h& @% V, dForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;7 q, m/ y3 N7 ~
   
2. mvar:
   ) R0 p& D* y0 V\" i\" o/ b& h7 x( Q
3. t=sys::clock(),6 b. c* ~3 V1 S\" F
   
4. s=0,
   1 v) }1 C* `' \3 k/ d; C
5. x=0, while{x<=1000,
   / O! v3 u( }\" a, B7 h
6.   y=0, while{y<=100,
     q6 i' c+ }' ^
7.     z=0, while{z<=100,8 W4 d8 W; c+ y8 a
   
8.       s=s+f(x,y,z),
   # @* ^' N9 X5 z* k5 q: I: f7 Y9 R
9.       z++
     w) l  L+ g4 t! s+ a& P
10.     },' \' S& ^/ D\" G' v& m1 ]) ^
    
11.     y++\" @) ], O\" p- A' Y* d1 |& T# D
    
12.   },
    1 O3 z, @) Y! f- [+ \
13.   x++4 H5 L1 ^6 a) U+ A% W2 s
    
14. },
    ( h0 U+ L3 e  L' n# c\" B
15. s;( V$ x3 p5 F! ^
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
4 e  R! n+ ^# w* R4 u0 W1.109  秒

9 b& Z3 _+ ~- r& g! R二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

2 o0 X, [1 S! x( E1 I) B% H

  H2 \1 |# h1 M2 n* D& t( V8 v7 X* r; _
& N! d. l2 S5 w* k; t6 x
) P& T2 F$ D% z/ `


% K2 b/ \% H' b( ~. D! o3 B
1 ~# r6 D3 p+ P0 V6 m8 e% q& E' u

9 o0 J7 l+ ]- ?. H# u7 }

3 h0 p( S& F) K+ c
7 d* `6 T3 v" J4 r8 e7 \