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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
, R/ v  _% m1 b) j1 k$ V5 p
＝＝＝＝＝＝＝＝＝＝＝＝
" I+ l1 n( n2 h4 W9 f; S8 I
; V( `( h4 z9 s4 C* H& @- J" vMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;0 `3 M$ ~$ R8 [; n: f; r# C  q4 `
   
2. tic;
   8 N) y5 U# f5 O, e3 c. \
3. s=0;
   0 c% M, h/ L+ i# ^5 {8 Z
4. for x=0:1000
   % Y- X2 N; T$ i7 v
5.   for y=0:100
   # {7 g( {5 c* h/ S# n0 \
6.     for z=0:100# R# A( C- y5 Z. J\" G( f
   
7.       s=s+f(x,y,z);, z; E7 M! E7 w4 ?3 a( T
   
8.     end
   5 L; S5 m\" D7 M9 B/ k: U/ I  ~7 r0 b9 [. u
9.   end
   # _) i, O( Z$ g
10. end
    / J. W. @6 w/ d7 p2 k& X
11. s
    # d\" X  ?& P* _# p\" R0 G' \% ^& y
12. toc* H+ ^7 e; B) l& O7 ], Y- l
    
13. 
    7 V) A6 {% k. I/ I
14. s =( N0 C5 P/ c! x+ E. c2 a\" _
    
15.     6.126720600000000e+009
    $ @4 s. J5 U2 s6 O! i
16. - Q) n  v- x! b: f/ G
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
' ?( x- u- Q- C# Sm函数：

1. %file xyz.m1 P$ m; ?8 q! O/ T
   
2. function c=xyz(x,y,z)! X0 n9 s0 J/ o
   
3. c=x+y+z;( R' w  _* E, {6 @: \
   
4. end

复制代码

测试代码：

1. tic;
   % o# V3 I) k: R4 u. _0 u
2. s=0;
   ( O% C, T8 J/ H& a( ~: _
3. for x=0:1000
   # \% m0 `/ U7 p3 A; ~7 l3 B
4.   for y=0:1000 u. p' B* E* v& g
   
5.     for z=0:100% k\" K/ f8 H. K4 _$ V5 V# d1 \
   
6.       s=s+xyz(x,y,z);: M' g1 ]' H) W6 D, r$ l1 @
   
7.     end6 O9 u$ S5 Z9 f6 T2 a2 ~) |3 c
   
8.   end
   & r) q4 j; L: _
9. end! U7 |7 p* S8 S/ N4 t# W; ?1 i* o7 i; J6 E
   
10. s, N9 F) p9 \1 Z9 _* N: V5 `
    
11. toc
    ; u! Z# ]  v* O0 Z3 V  A8 F
12. 0 F: ~# T9 ^0 R) {; n$ @+ X) a
    
13. s =
    , Y% Q0 h' X+ u; L0 p2 J
14.     6.126720600000000e+009
    ( X( K. f2 A6 D5 {0 K% _
15. 
    . p; ?; I4 ?5 ~3 Z
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;6 O1 |: c% X5 E( ?! u; U
   
2. mvar:) m. O# L% B0 d! B+ W; S
   
3. t=sys::clock(),5 O# c5 B\" E7 v$ {
   
4. s=0,4 }; L  T+ {$ n( V3 \
   
5. x=0, while{x<=1000,5 r: e: g% W! a( j; w
   
6.   y=0, while{y<=100,1 c4 s# |2 O$ X0 j6 R
   
7.     z=0, while{z<=100,
   # v7 Q, d% X1 X; f! s% {  u\" x: \
8.       s=s+f(x,y,z),
   0 ~4 j; D$ I! f
9.       z++
   ) j2 j- W0 y! Q4 ?. R$ L$ U
10.     },
    - f% |: t& B$ B\" ?$ f8 ?
11.     y++& Y5 I! x! i7 U, U
    
12.   },
    5 m3 |/ _: n, r$ A+ _3 X: A
13.   x++
    3 w! L5 l* L4 Y1 D2 z0 ]
14. },; p$ L4 U4 m5 J4 O3 Z
    
15. s;$ x0 H4 z' u- |8 ?
    
16. [sys::clock()-t]/1000;

复制代码

结果：
( U+ p9 ?9 Z3 H! F1 p: [6126720600.
8 e% D$ A* R( `$ N* Z8 R1.109  秒

4 v9 ?, C# j) N5 O4 O二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
* g2 q4 B$ o* Z9 B: H

hxftntA

6 w" S- S! j( _! h


* y+ f5 I% Z! V" d
# z  [/ K4 ?: k, f/ N6 O
& F1 l1 v4 j/ }1 u3 q" j


" Z0 F$ A' ?  x3 e6 ~

; J, C; _+ K1 f" ?$ u  s

$ G/ H- T9 g( @% P9 E
( S/ @' s1 D  o8 J
  P) X$ z3 w' I4 q! [# n