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

forcal

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

1 s" T8 m6 Q3 _6 `＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   5 \2 j/ Y) ~* t+ m' [
2. tic;
   , ?\" e( z( }0 e9 }( ]
3. s=0;; W! P0 _8 ]7 U
   
4. for x=0:1000: p# c9 |4 U8 S1 N2 H
   
5.   for y=0:100( X2 M& X3 Q7 u0 `) ?& @0 o# ^
   
6.     for z=0:1006 n; c) [( C3 S- A2 M: Q# @3 ?
   
7.       s=s+f(x,y,z);8 O; s9 S- t( e5 N
   
8.     end
   * D. r5 S% F\" d/ Y
9.   end
   ! t. h5 M, W0 N! M; P
10. end
    \" M2 n6 w7 L8 T% ?: x; ?
11. s0 u3 X$ ?+ r) l, K4 @
    
12. toc1 _; K2 x: R. E! [3 g
    
13. 
    0 C/ |) c2 z% k
14. s =% B0 Z! Q. x$ V& {9 I3 e0 [4 p
    
15.     6.126720600000000e+009
    ' b- G: M! O& w# U- h! D4 ^3 ]
16. 3 P, G9 u- y9 A6 f1 h
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
. G2 u9 L$ Y2 \& c( l: @' A7 H, `m函数：

1. %file xyz.m2 M6 {) v( Q6 W. _
   
2. function c=xyz(x,y,z)
   & r  [0 b- K3 o: c6 J
3. c=x+y+z;& `( T5 Z$ `: E+ i
   
4. end

复制代码

测试代码：

1. tic;4 t7 U1 I: \  j0 i0 D+ }2 ^2 ~
   
2. s=0;
   ! \/ m+ ~! M0 J9 }* {8 M* u
3. for x=0:1000. }* \3 O* E9 R3 n
   
4.   for y=0:100\" ?3 v2 y% s/ f' [& e2 U
   
5.     for z=0:100
   ! D- z6 Z& G: F3 r: t2 g1 D
6.       s=s+xyz(x,y,z);' ?) x  |0 F; O, x
   
7.     end
   \" T0 z- P, {5 s6 {3 x; c
8.   end
   1 w5 Y: C4 w1 m+ h$ l
9. end1 O: u7 `6 [$ a. q' U
   
10. s9 _. g& m& X. }' p& O
    
11. toc
    6 R2 o% \/ N& {+ Z6 j6 M0 k
12. 
    7 w5 N6 H+ b5 W( U
13. s =3 e  T$ A\" E$ g
    
14.     6.126720600000000e+009
    & M6 H2 L0 h  l; q+ Y6 k/ Q
15. 2 \6 \& Q- q0 F% ?8 M! s: J
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;- S1 v6 i8 R  {! j6 P1 j
   
2. mvar:
   8 e- C; k5 g$ P4 I. S! I
3. t=sys::clock(),
   + W: S6 g2 e; D$ g! f
4. s=0,
   , G7 M* E. W' n, e
5. x=0, while{x<=1000,% U\" G% f7 l/ g\" ?: h
   
6.   y=0, while{y<=100,
   4 _2 w. s. I/ A9 ?  ]
7.     z=0, while{z<=100,) I( t' r% h3 Z  N
   
8.       s=s+f(x,y,z),) t7 h) c9 I- A
   
9.       z++& n2 N/ V8 u& Z
   
10.     },
    # _% w$ T# J* M) ~5 q
11.     y++( V. h( m) `7 n9 G- s3 T/ p  n
    
12.   },2 A0 V) W( ]+ I- {
    
13.   x++
    : y! k/ h4 _/ K* Z
14. },
    2 ~& m\" i! }+ ?7 t- ~  L
15. s;  d3 @/ W& [' a! ~
    
16. [sys::clock()-t]/1000;

复制代码

结果：
5 P# |% S  z" h! G# R+ B6126720600.
% m4 p' N0 K! q1.109  秒
! R( U" V) Q% L9 I& s$ W- T
3 A+ x* |! {% Q8 U7 v二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
% L# r4 q& q, l4 ~0 s( e' m

hxftntA

, a, t+ B4 U3 ~: T; G/ j) {. t
( j4 q1 }+ ?( @9 v/ T

7 S# ^0 C. ^6 K. c* Y+ ]

  t4 A! ~3 r9 [





0 h/ E' m- d! O+ e" L