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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
3 @9 y# i( Z/ b* z
1 x  J5 q( m0 X0 @＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;5 V' n8 t/ d) P! E7 e
   
2. tic;1 Z) E  `' W) q5 o; Z. V' T9 L( Y
   
3. s=0;& U2 J' q1 x3 s8 R
   
4. for x=0:1000
   3 t5 S) P7 w9 I( o6 p2 m/ u
5.   for y=0:100
   ; _/ W/ I. c& r, L: t& _
6.     for z=0:100/ N5 K- y% o$ ~
   
7.       s=s+f(x,y,z);
   3 n- [: F$ J# Q
8.     end7 m  x: U8 l5 j8 a) Y) [
   
9.   end
   * \2 U( P5 a% y$ v. M. z
10. end
    % d+ T6 M, C& }! f+ z; q8 N
11. s! u5 a' n8 d/ O& k
    
12. toc
    5 K( i# w4 L9 |3 i# g; I. W
13. 0 y% S5 z\" L* e, T  R
    
14. s =1 K( I# ]% Q5 {6 {
    
15.     6.126720600000000e+009
    ) w% m) o7 M7 X+ m
16. 
    + r0 C5 @7 C1 @# C  E! g' `  l% S
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m0 w: T( O) v; |3 g& C! U+ q
   
2. function c=xyz(x,y,z)
   1 i  [0 g# C& v: Y/ u  N- ]
3. c=x+y+z;
   * x: ], e! I  c  z/ e8 G3 z  U8 S. X
4. end

复制代码

测试代码：

1. tic;
   1 D7 |, E: a% i# F
2. s=0;* b' f7 M6 [$ V, {# ], ?, T1 h
   
3. for x=0:1000
   2 |; a\" }# |* k- n9 S2 I
4.   for y=0:100% i+ R6 ^\" i& j6 x
   
5.     for z=0:100
   5 Y+ Z) x2 w% V\" d- W/ @
6.       s=s+xyz(x,y,z);: S9 H$ j. y5 s# p1 t$ f5 {) Y
   
7.     end
   / i- _0 O+ ~. W, S1 j, c( B. O
8.   end& @6 @. |6 j9 k9 b4 m6 ?
   
9. end
   ; A3 R. [, b3 R; s: t) \
10. s
    9 o% W. G; e\" r* A9 }
11. toc8 R. w7 W' Q+ }0 a# v7 n\" z
    
12. 8 {! }\" V9 R3 ^
    
13. s =! I3 v4 M' \; M) O+ o  E: q
    
14.     6.126720600000000e+009
    2 Z* m6 H, [! a* t# _: w$ O
15. 0 R5 K1 O% _% C& d- Y
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;% M; p, b# P: k1 I' }
   
2. mvar:
   ) f2 h3 r& P9 j: t& A
3. t=sys::clock(),
   + g- @: m+ B3 M% h4 O2 \- t
4. s=0,
   6 T* L! M\" S+ V* J9 v
5. x=0, while{x<=1000,' J, B  F3 }\" R; W/ [
   
6.   y=0, while{y<=100,+ _' G- ]( {4 t. X7 y/ z
   
7.     z=0, while{z<=100,2 s6 {7 l4 S, v
   
8.       s=s+f(x,y,z),
   / L4 A5 @0 M* t( d: n; p5 I
9.       z++8 x7 A8 A% \0 M
   
10.     },
    ; P- D2 l\" H% q  b3 ]; Z4 b, q5 Z
11.     y+++ O* I& X! w' s/ N
    
12.   },
    : _# d& m) D4 R: \) N. w5 L
13.   x++# ]0 x2 F9 a  w2 j2 Q: e& ~, a. Z
    
14. },
    ) X& e4 P# g6 b2 z/ `
15. s;
    - h\" N$ f3 W! `& {; f# i( _9 n3 `
16. [sys::clock()-t]/1000;

复制代码

结果：
; l7 y8 y( C; S  f' u& g+ b6 I6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
1 z* P* F6 v+ V& y7 l/ O. G

hxftntA

% X- ]+ B! l3 p2 O

6 |. C3 H( f7 L% _' S% i+ Z
, h& |# l; r! A. O/ G  K3 B
7 \/ E, r# t7 |2 O
/ o/ D) U; F# u3 k3 D* N
* t0 e$ W( T5 T3 c  C3 W# U
* W2 F0 z' D! I: a
  C, Z9 P# T* F! e8 I1 e# B
+ f) c2 h$ J( y. l
: z4 B% D$ b" J/ q3 ?) J! c. R
  _! v" p, r  K
8 g" _+ D4 D0 ~, E5 {0 l
) k2 g! V; o# q- w9 u

. ]! n1 Q( X, U2 E: R: |9 a! ~