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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
% K- g& D$ G. o. \+ s7 g1 F: m* m4 i
' l3 D+ z  `% i: Y0 q＝＝＝＝＝＝＝＝＝＝＝＝

/ A, B( a! d7 S  S; ^0 C# tMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;. H\" j5 |  t  A$ S( }
   
2. tic;- a, m* G! p& J1 i. ]
   
3. s=0;( L  K3 o\" {6 t4 Z0 ~! ^
   
4. for x=0:1000$ V, b4 v! L# Y4 J) S
   
5.   for y=0:1004 {4 G5 C( F! I/ E! `$ z4 Q
   
6.     for z=0:100  z5 \, t0 {. E# ?% E
   
7.       s=s+f(x,y,z);! ]& c: m% ~9 U  m4 i  s0 ?* {
   
8.     end/ r- u5 I2 r$ k: R
   
9.   end
   8 ?. V- q# ]8 J  I* Z, {. s
10. end7 l$ {6 r; J5 f# ~
    
11. s% m8 |9 w' `% a
    
12. toc$ x3 b' i: Q% u& \* h
    
13. 2 v! q+ C; k9 `8 P. y( x0 r5 P
    
14. s =
    1 Y  p) [8 P2 @$ V9 W, m. \0 b0 W# a
15.     6.126720600000000e+0097 p! L- b  |, E% h3 u2 f+ x9 }
    
16. 9 U6 E) T0 T4 h
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   ! D) x: ^2 F5 w  F+ I) K4 ^
2. function c=xyz(x,y,z)! c( T/ b# f+ W; B) ^
   
3. c=x+y+z;
   ( N# W3 W0 m  e* V5 H1 n( w
4. end

复制代码

测试代码：

1. tic;3 C  G# T! n( T6 {\" l: p
   
2. s=0;
   4 _1 V! ]+ Z/ n' ]
3. for x=0:1000  K! m8 h, ^0 I! s' M( w
   
4.   for y=0:100: K0 P& Q- n( M4 ~8 O- i
   
5.     for z=0:100
   : S  N- O2 r- b8 O6 b) i
6.       s=s+xyz(x,y,z);. h) g; g8 B6 Y' i) d0 S* D# o
   
7.     end
   6 @( r8 D) g# `4 C4 [1 t/ R5 q* j
8.   end
   ' D3 x( Z8 L: \
9. end
   : A9 j8 ^; u, {9 b+ n5 M/ ^\" f
10. s
    0 Q* x1 s4 O! s5 U\" x( n! i
11. toc1 }( E$ p2 J# K7 h5 p
    
12. 
    $ F2 C# ]( `% a9 W% B1 |
13. s =
    4 b9 G# p- A4 x( X, P0 R, `
14.     6.126720600000000e+009: q- n: _. T- W4 s7 I
    
15. , @; B; ~  [6 x( ?
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
2 L* w$ _! ~+ M% L1 J
  u7 Y' x% e3 w- kForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   & V  P! c! W% r+ v& j0 Y
2. mvar:
   . y; v' T$ }! ?
3. t=sys::clock(),
   4 Q0 Q  A1 V1 z7 C
4. s=0,  _: P* V7 R' O$ z
   
5. x=0, while{x<=1000,6 c6 Z0 I/ O- g( e/ T2 M
   
6.   y=0, while{y<=100,
   , b  `( f% x  H( U1 i6 x
7.     z=0, while{z<=100,
   / w+ d0 B! C6 g  @
8.       s=s+f(x,y,z),\" v5 g  V3 t3 ?5 T# j: J5 a' T4 r
   
9.       z++
   7 ]& c/ F/ {6 c/ d5 {* T, Q6 W9 G) p% r
10.     },
    \" @8 O- g\" t. p, P& K
11.     y++/ g- m! M! U! e5 ^/ y
    
12.   },
    9 `+ ]* T/ ~9 G% N: N
13.   x++
    0 x$ f\" _1 P\" b4 @
14. },' I  i6 d6 L4 _5 K7 s, S
    
15. s;9 c0 |8 l$ m- d* u1 b4 {
    
16. [sys::clock()-t]/1000;

复制代码

结果：
1 c1 Z- M* x; O0 ~" k6126720600.
1.109  秒
* Y5 y# x4 J) Z) l' ], W  l
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
, B5 d4 \& q$ M! \

hxftntA

; e+ \! G; G2 d( A
! U6 O1 y, P( q6 R* |# ~
$ S4 U0 v* A* d  d1 I

2 P: n* `: C0 O9 O4 R7 ^9 T1 V


7 {% }5 N) Z/ F



2 P$ T7 F. Q+ ^! j# ]
3 f) _/ i! K4 ~8 ?0 R8 }

. @$ w& G9 ]( \% q