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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
( L. f5 ]3 [% D6 w6 e
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   3 Z$ I9 Q* j+ V0 |\" o+ i6 E
2. tic;& J, l6 Y# n0 B$ x
   
3. s=0;
   * q! i0 j! U7 E! w, v! k
4. for x=0:1000
   / J# a2 A\" y8 L+ J
5.   for y=0:100
   3 H$ ?9 t! q5 C4 v% ^/ @2 w
6.     for z=0:100/ K) L8 h$ w4 t7 [
   
7.       s=s+f(x,y,z);
   \" d3 L7 C9 }$ X
8.     end
   ; u6 C* D7 S$ J* G
9.   end/ c& \5 {8 l# u/ X8 U
   
10. end
    . T6 m& d  m5 l. t! m3 J% \( ~
11. s' g) t3 w4 v\" Q0 w) Q
    
12. toc, H6 P( _2 N! @9 c' ]
    
13. 0 K4 F6 E- y0 S' o
    
14. s =
    4 ?) D  B& C3 m5 D
15.     6.126720600000000e+009
    $ i4 u# j) N9 E! A/ f
16. 
    ! D- B; f& ]9 u2 f. t, D
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   \" y' @0 S5 |/ T
2. function c=xyz(x,y,z)
   3 S$ g# D. Z$ A
3. c=x+y+z;
   8 l/ V; T# W9 a
4. end

复制代码

测试代码：

1. tic;( j\" B2 B, I1 z. j
   
2. s=0;( A0 w+ x; N: f, z
   
3. for x=0:10001 c. H# A. l+ ^( J( e. m/ f
   
4.   for y=0:1008 E4 V3 b4 Z\" r3 @0 g8 ^- E
   
5.     for z=0:100
   2 p  A5 U4 z2 W/ I  L, K% @% w. E
6.       s=s+xyz(x,y,z);
   \" m3 ?/ Q7 A; C
7.     end
   ( S0 n5 E/ E. x. G
8.   end: l+ s0 _% N; l( e) m
   
9. end, Q* F; E  a7 O1 V
   
10. s
    % ~8 n% j  z# {  i4 `' \  P
11. toc
    ' U) m$ ?- B2 M; S
12. 
    ) k: e4 h% D- t' F8 @8 Z) p3 n' C+ F+ h
13. s =$ M: {6 I1 w2 {. M9 V
    
14.     6.126720600000000e+009
    : p, `3 r4 E8 ~
15. 
    $ X! C+ d* E6 o9 Q9 {
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
/ C% f0 M) w, P
4 W* J' M7 g9 B/ w% aForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;/ T9 @: y6 n0 h
   
2. mvar:
   : b1 x2 p) G0 h2 K- W
3. t=sys::clock(),
   ' m& K' X$ M  f* r
4. s=0,
     Z& b' a\" ^1 {+ _. t5 H
5. x=0, while{x<=1000,3 x, d; R7 M2 [$ `\" W$ D
   
6.   y=0, while{y<=100,4 r  y* N0 @$ p3 J
   
7.     z=0, while{z<=100,
   ' w# |' ^6 g& v6 b# n& Z
8.       s=s+f(x,y,z),5 `( l: F  s& |\" h* b( o' f
   
9.       z+++ `: p, J\" o8 o! |; N
   
10.     },1 Y# l2 c3 C! r. U) g4 s; x
    
11.     y++
    4 _) z! `\" Z8 \
12.   },7 X\" Y; b, O2 g6 N& Y% D1 X% M
    
13.   x++9 B+ |/ L' [; V  V
    
14. },  V8 t1 X7 l7 g\" N( @! f/ ~
    
15. s;
    . f: l/ u; G1 c. H
16. [sys::clock()-t]/1000;

复制代码

结果：
" y$ F2 r5 m# D$ U6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
6 r" w) O; H  O& S* a; `6 a

hxftntA

* _% z) \; s) v; a) u/ b
, n8 V, C1 o; |  J

* p* W( j  C0 v




6 A# {0 ~) Z8 _


- H! Z( C; C; T. }6 ~/ u; R( ~6 `+ D- |6 D


! }1 _; B1 G6 Q+ n