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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
9 e3 p3 I# b) i& C. b9 u* c
0 m9 ?" \1 z, X" l＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   : X  p( k2 w  v# N! R
2. tic;
   2 e' K+ V0 k/ m
3. s=0;( z' t: b\" H' B# J$ J/ n; ?
   
4. for x=0:1000\" e* ?7 q7 |9 r; q' @+ O$ ~  [0 C; w
   
5.   for y=0:1000 P! a9 m7 E- l3 m& ?- A
   
6.     for z=0:1000 [3 q8 F; p6 R8 ~# {1 ~
   
7.       s=s+f(x,y,z);1 s' A: [5 _5 U* v
   
8.     end
   ! |& G# U9 l, s/ U+ t
9.   end
   6 I' y% R: N! o- F+ U* J& A% B' A
10. end
    7 x& I3 P2 z3 v\" m/ E! f, }
11. s
    8 |\" ?: b% U# O; v
12. toc4 Y\" C9 S/ A4 \2 v
    
13. * t. n; |1 J; t0 [6 E4 y
    
14. s =: b6 {- v8 o3 K/ J3 w6 h
    
15.     6.126720600000000e+009% A' b& P( i2 V. h! J
    
16. 
    6 q2 c- a3 [2 E+ w2 M: C% W
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   * B+ K# C\" z  X! P6 X
2. function c=xyz(x,y,z)\" T6 _' {* R' M
   
3. c=x+y+z;- O. Y0 x1 K6 O# r+ L$ I% q
   
4. end

复制代码

测试代码：

1. tic;8 A/ w% S& ^2 F. D# h  ?
   
2. s=0;
   0 N( p/ {9 n) H: l9 h- ], I% O: i
3. for x=0:1000\" }. {0 t7 W+ |1 Z# O
   
4.   for y=0:100
   . f% w7 k6 {: m/ Y/ b6 a
5.     for z=0:100
   2 K\" A: j/ U8 \: P7 c& z! t
6.       s=s+xyz(x,y,z);+ N% x# S\" e) w0 w: t, `8 ^3 L
   
7.     end
   0 e5 x6 d/ {. B
8.   end( D\" M  l! y\" j* }1 d8 c
   
9. end\" Y2 X; X/ j( }9 g
   
10. s7 z+ X0 a& H) H+ `
    
11. toc# d/ `$ o6 C4 k& `
    
12. $ ^0 k0 x- q; I/ r9 r0 {
    
13. s =\" Y& t# F1 i1 |/ v9 V
    
14.     6.126720600000000e+009
      _$ U8 S. h6 B6 z* N) H& s
15. 
    5 c9 G  _3 ?+ U* P7 h. f) q
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;4 N  H, b) v# @' a+ |8 m
   
2. mvar:
   - N. M6 w+ L3 k) W, m8 t
3. t=sys::clock(),
   . R/ Y* ]& Y5 U) c& u# Y1 a
4. s=0,
   ! w* |- R1 m+ W$ t) Y/ F/ m
5. x=0, while{x<=1000,
   , [; m3 t! q' |) _; k
6.   y=0, while{y<=100,% h; F& S' t5 @4 Z' _
   
7.     z=0, while{z<=100,3 f+ \- T4 f) M$ S\" X  |; ]
   
8.       s=s+f(x,y,z),
   8 |' A5 b( s+ b6 f% P
9.       z++& r1 @; J$ m, ]
   
10.     },
    3 K2 c9 L( b. s7 U  Z+ T' i
11.     y++
    6 u0 r# C1 \. I8 J* F6 D* Q
12.   },
    ' ^* [9 Y/ `9 |9 ^, s' Y' t
13.   x++
    & Y. }) l1 E4 E7 ^: o  S
14. },
    0 J/ Z( _& E0 a
15. s;
    # t0 C5 B) j' j* J$ J0 g
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
' f1 ~9 ^9 J. r' c9 o& Q0 P1.109  秒
& {9 V/ c" b' I; }* l
0 G4 [& A' K( R' w: ?4 `二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
  o5 i8 V3 U' J/ m+ B  {) V" l8 H

hxftntA

( L" b) P! F8 ?
# W" `9 |' P* S( \* b# z9 x6 N





) \; g  y4 M" C9 u$ z8 ?
) V, G- C0 D+ }: Z: x; k4 \( n  q


% r2 X6 ^& g( h9 F  c
3 _2 l  p7 P+ k6 W7 O

7 h2 H) W) F7 O
* n5 q: h- C* c0 `% w: T