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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
! J8 W' G  {& u  G; ~9 b4 q7 ?1 [
9 |4 p" [( F( ~% a* T9 R: e. R＝＝＝＝＝＝＝＝＝＝＝＝

0 B2 t0 [1 X& jMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;  {$ p1 z' i+ ^
   
2. tic;
   / m# j# q5 d9 j* r, H1 C+ Y# z
3. s=0;7 a2 B* _( ]5 d/ D) |
   
4. for x=0:1000
   # F4 f3 D\" V: o# d; j8 z1 o\" O( z) v2 d
5.   for y=0:1000 X% Z$ J5 s7 _. e, |' A( J
   
6.     for z=0:100\" `+ V! C3 ?& B6 Z8 _8 L9 o
   
7.       s=s+f(x,y,z);
   7 T* H9 _1 p' l4 m
8.     end& T4 w0 ^; o: ]) a0 o& r$ L
   
9.   end
   , n0 U8 B1 B2 w# |
10. end5 Y$ a6 c8 K. t\" S$ o
    
11. s
    ) b% B1 ]- m- z
12. toc
    3 ^8 O/ s\" S* I
13. 
    # k  K8 i/ F2 B! ^& Q: H
14. s =, R$ K8 X# y1 L- b
    
15.     6.126720600000000e+009( O) T8 U6 {$ s- z4 n* `
    
16. 
    0 q: R6 K: i/ W. c! F
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m  ^  q7 }5 w7 C# x! V; V+ K2 G8 A
   
2. function c=xyz(x,y,z)
   ) C$ d/ A& R+ `5 @4 R. A' C
3. c=x+y+z;, l( X' q. O3 a, y\" n9 U; D
   
4. end

复制代码

测试代码：

1. tic;: T% z& Z) y! [
   
2. s=0;
   , F, n6 k4 E\" ]) J$ P( A5 x
3. for x=0:1000- w3 F/ X! m' d9 K5 `
   
4.   for y=0:100
   8 C3 ~! }7 v0 ^( ]
5.     for z=0:100' Y\" H8 e. g7 C1 A
   
6.       s=s+xyz(x,y,z);4 \) C8 y  D: @8 R! h
   
7.     end1 e$ {/ d: ~6 e4 L6 U+ X8 v; U
   
8.   end0 _+ F$ D) _+ l! Y8 }- k
   
9. end) e9 n\" k6 ^; x- `
   
10. s
    % E) _, A, t& S: q
11. toc
    # O6 \8 K& A! a+ a, L
12. , @6 z1 d6 c' z& z: c. [
    
13. s =
    + E; {* e  q+ Q3 f. L9 d' o
14.     6.126720600000000e+009
    * E6 C& X/ O7 ?: r& C; b( ^
15. 
    * Y; x4 D( }: S
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   0 f8 I5 n! F+ K# t1 e$ F. L2 c, A7 L\" k
2. mvar:+ F& ]: k) {  M; b- C
   
3. t=sys::clock(),
   2 {, |; r1 f) F5 T; f# i, c+ d
4. s=0,) Q6 x! R$ l6 {. h0 U: S
   
5. x=0, while{x<=1000,
   \" M\" N: B- i5 j1 ?1 R8 w
6.   y=0, while{y<=100,
   + m2 k\" l) ?) X) W& Z- H
7.     z=0, while{z<=100,  O$ H6 q0 o. F& y
   
8.       s=s+f(x,y,z),
   ( J+ G( l9 C9 \\" r0 n\" C$ H
9.       z++2 x\" P( \/ R; K9 J0 w: A8 L8 y' Z
   
10.     },9 t9 V6 X3 z6 N4 E5 t$ a. Z. u8 ^
    
11.     y++1 p& i9 |8 C3 `) C
    
12.   },
    9 t% Y0 @# Q) a! z$ b7 f
13.   x++4 F3 L8 M! i; ~: z
    
14. },
    , f0 t  F( _7 [- N3 {$ _8 }
15. s;
    % O& e$ p7 l. S
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
: S: z+ z5 y: k9 _0 N; B7 X1.109  秒

% b4 e" Q4 Z; D0 _8 \( f二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
3 w) x$ p/ a% v9 x$ C; |$ x. y; ]

hxftntA

  h+ J& a. A. \. h- i1 ?
, ]  n% F* ^# k# A3 V
9 {' h; m, u6 ]: ?6 L8 P

: i% P: D/ G/ ?" T" G
6 s' k3 j2 e4 E: N" g4 j


; g6 E. V- Q+ i3 m


0 J! C, R) A) ]! A