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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
$ O' r. E6 S" m/ N
5 O  A! g, Q& w/ `2 Z+ l8 X＝＝＝＝＝＝＝＝＝＝＝＝
- W" \: i6 n, z% d- [+ A
9 E6 y* H' p/ E  T8 W+ B0 mMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   + t% T2 N+ S9 T  v
2. tic;
   : J; p; p( @8 h8 M/ w
3. s=0;
   2 ~\" ?3 d8 ?) ]' j; s- l( D  _$ _
4. for x=0:10000 b2 v  r) C  }; j: o3 d
   
5.   for y=0:1001 u8 B5 U$ m  F8 s0 V
   
6.     for z=0:100$ Y$ _, f9 j2 h) l  \$ f
   
7.       s=s+f(x,y,z);5 i# z4 c5 l# V& \6 V
   
8.     end7 i: v( i  L2 B. ?
   
9.   end; ?+ p3 {8 A) R) ?+ f8 ^; t
   
10. end
    / r- z& |5 u8 B( k/ f2 l
11. s
    ' E8 W9 j8 ^- }6 Y6 I
12. toc2 V# c4 z# `\" E9 w2 T0 B& o
    
13. 0 c( j$ J  U( i' C& b  u
    
14. s =: Q* i: x- q  ^5 {3 x
    
15.     6.126720600000000e+009
    / Y$ G  Q# o+ W9 o8 s
16. 
    $ [- _$ h) J( L
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   % |\" p% X) M, \) I\" G4 ]' q) x2 o
2. function c=xyz(x,y,z)) |. a- a! A0 m
   
3. c=x+y+z;- I# s* M4 H# x* k\" @! M
   
4. end

复制代码

测试代码：

1. tic;
   , y9 J2 z; \% y6 @+ Q
2. s=0;
   # `8 Q+ L8 \* i( ]) A3 t& |
3. for x=0:1000! z9 K7 r$ }2 g1 w% `  B8 M6 @
   
4.   for y=0:100
   ) l, n# b7 u9 R4 f\" M- l3 e4 `1 }
5.     for z=0:100
   . |( c& i8 b  G0 b3 s
6.       s=s+xyz(x,y,z);# e1 d0 q# b9 U& X; u! r/ Z
   
7.     end$ \+ Y% f( C+ T  a. e2 d
   
8.   end5 K: M' N( {/ f
   
9. end
   ; p! o5 p% D$ V- o8 B  b( t
10. s
    2 d8 r4 c\" ?# I* M2 [
11. toc: ?. E. g% G\" J  V- Y8 {; O
    
12. 
    : G$ Z3 Q1 L* {) ~- K3 r) d\" Y
13. s =# F, }% v' H/ a1 M- C
    
14.     6.126720600000000e+009
    9 S* ^/ S6 @2 V' |
15. 
    3 t$ t; g  p6 q3 h( P
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

* q8 ]: k% Z) e2 a2 O6 `, j9 i: tForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   + C; K: w$ B1 ?* P- }& p( m, F: ]
2. mvar:
   3 C& J  E0 l5 D0 e+ O
3. t=sys::clock(),4 H! _$ X+ z6 ?5 R\" n8 n1 P. n, A5 o( O! q
   
4. s=0,
   $ }1 i# r# {$ l) F\" b
5. x=0, while{x<=1000,3 Z$ s& Q# k  R2 _: u. ^, l
   
6.   y=0, while{y<=100,. Q/ U3 G, t9 w! b$ i6 O+ _
   
7.     z=0, while{z<=100,
   4 P5 u2 C9 H4 }3 N) M- G
8.       s=s+f(x,y,z),
   + x3 S+ _  [1 a
9.       z++6 q, ^( m+ H: l8 S7 l  [5 j0 V) m
   
10.     },
    & E& n) V7 Y8 n. t
11.     y++
    1 K- Z( U. ~: t1 C3 J4 m2 K. O( x
12.   },# a2 F( V1 }5 E; z5 k; ~$ L/ P/ j
    
13.   x+++ [5 D3 Z, C8 r! f# G
    
14. },
    ( q; e! j, r/ k$ s7 h! ^8 `: I- X. [
15. s;
      n6 |: |4 Q/ E% d6 f6 Y
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

6 y* J# K$ o. T. j二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

2 ^) [% L/ |5 p9 p' d3 n: R; U
- y& a3 i4 L3 A
9 I: W+ B# P. Q8 [" ]+ t0 B
* Z* r" ~* i. x' e% @1 s

8 I8 k4 z+ o: @+ _% v* g% o
. C% u4 C; O) a


2 D. Y$ C' Z' F% O& S& s
( c' C1 Q+ d+ K0 [0 `7 N



, B+ ~! t" R7 {/ \  `