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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
1 i# T+ D0 V* H- E/ M
- Q; P- y. c  N0 NMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;* L9 h2 U\" K4 `1 p7 O
   
2. tic;) t, C5 c; f2 D\" O6 \; o
   
3. s=0;
   ' O, i9 z- f7 t\" |: `' k2 q5 t
4. for x=0:1000
   5 L' K/ F6 C4 u% W9 ~
5.   for y=0:1009 |* r' F1 T% H& d' E$ l
   
6.     for z=0:100
   9 Q0 I  o5 m) h( ^6 @9 o( K
7.       s=s+f(x,y,z);
   : y; P$ f& s4 x6 m# \4 E9 H: y
8.     end
   5 V, A% ?( N% E
9.   end
   ( R7 P: _6 @$ @# _
10. end) Q9 ^4 d: I: J
    
11. s) B/ M. Z; |; w
    
12. toc2 h, j$ T1 b: y' K
    
13. 4 z5 i, i5 B( i4 Y2 ^! ]% v
    
14. s =
    1 F4 D9 a9 g( h# H! W4 p
15.     6.126720600000000e+009
    7 y& {. V3 S/ a$ \2 f( O
16. 1 s3 A3 I/ ^0 Y$ E# b# `
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   2 ^\" Y' \9 D6 w/ I
2. function c=xyz(x,y,z)
   7 g: f7 f\" _: e: k% ~
3. c=x+y+z;. `/ Y& O7 |' r\" b! i2 Q, `
   
4. end

复制代码

测试代码：

1. tic;, i) R# z! S% }
   
2. s=0;# M% J4 b1 U* F( Z, [2 O& i\" w
   
3. for x=0:1000! b* i; `4 F9 t6 L) f) x/ w
   
4.   for y=0:100
   4 F( A8 l3 A8 R3 f
5.     for z=0:100- T4 K! R0 z! X2 [5 N
   
6.       s=s+xyz(x,y,z);
   3 s$ E- B6 \. R; G! J1 y
7.     end
   4 I) P: s) ^* y8 o
8.   end
   $ G& T' S5 F5 t
9. end, p# k6 s8 \0 G# K: y
   
10. s8 o. r7 v& d) m( Q* K, f1 @
    
11. toc: s- u% l# R# |; \+ y* a- `
    
12. : `: m. S% r0 a9 J
    
13. s =
    9 L$ p+ r( h/ @  o
14.     6.126720600000000e+009; v$ T! ~# x1 @; L% ~1 g% p6 q
    
15. 6 Z2 v\" P2 s\" h9 G: ]6 S
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;; P' c+ b3 l6 [* w, R( p
   
2. mvar:\" l+ S# k\" P3 P; `7 l
   
3. t=sys::clock(),
   ) z\" D1 R3 W# I
4. s=0,4 {4 V1 P6 e/ W5 f* C
   
5. x=0, while{x<=1000,
   ' p; ~\" l: ]0 G3 {7 f8 k3 G2 k7 H4 I
6.   y=0, while{y<=100,5 Z  u/ c2 B1 W; K\" b* T
   
7.     z=0, while{z<=100,
   ; i) t7 r2 N* E5 r
8.       s=s+f(x,y,z),
   ; H, b+ W) G# S  |- z2 j
9.       z++
   \" I7 r' a8 [5 y0 j. f: {
10.     },! w) \1 m8 V9 s' j1 v9 c; U/ P
    
11.     y++8 H8 s! Y8 w3 a3 ?0 r
    
12.   },
    % ^8 Y, c& t+ ?& e( z+ y5 |9 I
13.   x++
    % C7 b0 I% y! k  h5 g7 a\" [
14. },5 ?6 M5 U( K/ f! @* c
    
15. s;! b9 w) D2 m) B2 k\" D
    
16. [sys::clock()-t]/1000;

复制代码

结果：
; l" m2 t% n" t. J/ u6126720600.
) Q, w9 d8 X3 r9 V8 U1.109  秒

9 k1 U9 }" h. c& }4 M6 n二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

9 j3 I3 l5 z" h, C  \' P

% _; o# t8 O7 |
6 F* T( r0 R# h% y+ ?




$ L! U; e- C5 E, H) |# _3 q
! y4 T) w0 w2 [
/ I0 I; {9 r1 p( ~$ `  ~& D


3 y! i3 e2 @9 ^) D$ ?

/ ]+ e3 Z2 J& `6 |- f2 r