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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
% o( Y+ |# w2 M( s) X+ B5 K2 S
＝＝＝＝＝＝＝＝＝＝＝＝
! B8 i7 G( A7 y+ f  M3 Y
  P0 f5 R$ y: \) v& U9 pMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;) V$ i# p# V* U! z, L% P
   
2. tic;4 r% o! O# E% U2 e% X0 W/ k8 c
   
3. s=0;* `1 U& C+ P0 F  Z- d: H6 o
   
4. for x=0:1000
   2 Z& b. W2 o, W7 M
5.   for y=0:100
   0 p1 z. H! Q' T1 H
6.     for z=0:100
   * w+ @  n* Q% p5 z
7.       s=s+f(x,y,z);) K\" v  {5 f. w7 [! b
   
8.     end/ E( X0 W# j7 \/ c  z
   
9.   end
   - I$ `) J0 w\" |% N: d) _, ~
10. end
    7 \  y( l4 T6 ~6 w
11. s
    6 \9 X% T1 T' m: u, s$ N6 E0 N. X
12. toc
    : j1 J5 y( W6 X, \
13. 2 E/ e% O: Y7 s
    
14. s =- A& s7 k! M) [. [
    
15.     6.126720600000000e+009
    2 t. ~4 r2 D\" _# }; o
16. 
    9 |; H7 s& V8 ?* @\" L
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
. v' N! B+ y: ^# c% a& x; m: sm函数：

1. %file xyz.m
   5 D\" @! v  V+ D! v; Z8 D! i
2. function c=xyz(x,y,z)  B% d  ~; A( N0 ?: L, L% j+ I
   
3. c=x+y+z;6 ^, d% \: {* U  Y4 b
   
4. end

复制代码

测试代码：

1. tic;
   ' `0 F3 E& t! H% f' q
2. s=0;
   ; p  k% R. K; E
3. for x=0:1000, T  M  U7 F5 f% Y1 |$ k
   
4.   for y=0:100, i9 ^\" N2 e- L+ h+ T% M
   
5.     for z=0:1005 R0 l& T8 b1 [
   
6.       s=s+xyz(x,y,z);
   % b. v$ F5 G$ T5 ~\" E, ^
7.     end( W7 O1 x) n+ L: X5 _! `4 r
   
8.   end  w! c& q4 b# l* w& U& T4 B
   
9. end4 J0 J7 y1 s5 Q) U6 v8 a5 y
   
10. s
    8 s' K9 y, m% X2 t5 e2 w$ f1 F3 y+ g  m
11. toc
    2 }6 S9 V' k2 R/ G0 v; e: ]
12. % q  {- [; ^; ]
    
13. s =+ w, ]8 h( p; f  K3 y6 n( q# f
    
14.     6.126720600000000e+009  Y0 H- V/ `- m
    
15. 
    \" V. f, ~) D- i& q5 F
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   9 l4 N5 L4 k1 P- m2 U$ V
2. mvar:
   & {+ N6 N* I8 \$ C! T
3. t=sys::clock(),
   9 b, z# R' n; |( F& r5 O, E
4. s=0,
     t; ^0 ^7 k# `# P2 Q( S% T% l1 Z
5. x=0, while{x<=1000,
   \" U4 U# |$ U! a* B
6.   y=0, while{y<=100,1 y1 }0 S8 @7 t# |9 Q4 l
   
7.     z=0, while{z<=100,
   4 a6 Z! V+ Q  ~3 t3 b6 J6 [; V% N4 X) ]
8.       s=s+f(x,y,z),: A\" |0 O/ Q% j- V7 J- y- c. ]( d
   
9.       z++( l! E% b' N8 E2 S* g2 }$ f; I3 t
   
10.     },; B( c8 q- y8 O. V+ D( ]2 N
    
11.     y++8 O1 V% j7 O/ ~+ ~4 u2 X8 l
    
12.   },$ Z) l% R8 K9 Q3 f
    
13.   x++
    , F; F  @\" Z; h3 l) _; C/ m6 U
14. },
    # ]4 W& f$ t( `* _0 D! t
15. s;* r/ v8 ?/ P' v. W3 A9 E
    
16. [sys::clock()-t]/1000;

复制代码

结果：
1 N7 V1 S# U+ o, [# V3 B6126720600.
& C+ t$ E# C2 [* M( r; [1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

* b* Q$ _- k7 N9 F$ t2 ?' i
) |6 H" V. F1 ^" a, N) V


- ~( @: _) E9 i; [




  Y. F9 V/ d3 p% {4 r! q0 {" f3 N
: e0 p% i- W2 z1 G4 M8 ~
$ [6 d6 S' E5 Y% t7 t4 m. F* J
, d- w# C/ k$ C5 [2 f