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

forcal

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

. v" T: C9 i! ^& q" c＝＝＝＝＝＝＝＝＝＝＝＝
& w1 z1 `  H; Z9 o# M
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   8 V7 Q' H! u' e* y: K( n2 C
2. tic;
   8 i& K5 j1 u1 a+ {1 r
3. s=0;+ X9 w+ a- B* V. y+ l, _
   
4. for x=0:1000
   & Z! m+ `7 `4 s( j9 f
5.   for y=0:100' N' e% o, O6 ?
   
6.     for z=0:100
   7 W% s& p# O3 j9 y\" D. z
7.       s=s+f(x,y,z);
   ( O1 z, n) R0 P# u* Q2 y
8.     end5 j1 P0 Y: u3 m% x! Y( ]+ ~
   
9.   end2 y* v+ Z4 f! u
   
10. end$ s  n; {8 a2 R* ?- Z, k5 j
    
11. s
    ) m\" S/ b' Y2 Y# \$ y9 j
12. toc1 b( v$ w; X( R5 G& A
    
13. 
    / u) h% R# m  b6 ?/ z
14. s =
    ! P/ R8 G; x) P' J
15.     6.126720600000000e+009
    7 {2 Z5 Q- ~% `/ r
16. & D1 A4 p! o5 w% _6 c
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
0 g7 R* y) i) Z4 o% R  {$ Y2 k, }m函数：

1. %file xyz.m$ }1 n; y; x: |' M; W& z3 M4 m
   
2. function c=xyz(x,y,z)6 m/ e. \3 T, e* O+ T; B4 w! x
   
3. c=x+y+z;# e% d; Q8 u3 j3 R4 a6 W1 C
   
4. end

复制代码

测试代码：

1. tic;/ U5 \! F) d2 O$ {2 T' F% R8 r
   
2. s=0;1 |: I* S; V6 j4 M$ L
   
3. for x=0:10008 t1 D0 @; q1 ^: G4 U2 h; I
   
4.   for y=0:100/ I& C' W8 C+ E% ]) f, x8 l
   
5.     for z=0:100, W/ ?5 |! w# S( C$ q. y0 r6 b( ^. H
   
6.       s=s+xyz(x,y,z);6 s2 s. l% [& s. W
   
7.     end
   ; T! c9 t3 i6 i! p+ O# l
8.   end
   & I6 y. t$ @7 O2 Z. D2 G2 I' Y* d
9. end
   , \7 u8 S2 w9 I: l7 ]
10. s
    ( h1 H: a, q% r; F9 h+ e( A) @( F1 K
11. toc$ \! w) C8 a1 ?' @, A4 T
    
12. & F/ H! C, T- s2 ]\" J% U
    
13. s =
    ! E$ F& e; l/ o. ]0 a4 C- v
14.     6.126720600000000e+009- H7 T0 s+ R' Z
    
15. 8 H, C: U( ~5 q  `
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

( }: d9 T4 _# m* H# e# Q! G. BForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;4 j( g; P, K* l: L: |$ D! ^
   
2. mvar:3 V7 s: ?7 H# a2 i2 o( {9 T
   
3. t=sys::clock(),' N8 _7 F. g6 s$ t0 v
   
4. s=0,6 y5 j1 U( H# I& E
   
5. x=0, while{x<=1000,
   2 q- U1 d; V1 B) x9 R
6.   y=0, while{y<=100,- Y\" c! [2 z8 N- o* b' `4 h/ v
   
7.     z=0, while{z<=100,
   ; }+ l2 `. L( [8 j1 o
8.       s=s+f(x,y,z),3 J6 J0 V7 J/ \
   
9.       z++
   ! c$ M5 ?  Q5 k! D
10.     },
    8 M$ e7 n2 F( ]- O/ A
11.     y++0 C( i( q6 O; L! t3 \; ^
    
12.   },
    \" i: f% T& F7 \9 o) \7 {( X
13.   x++  w+ _. V: U4 r2 X* t  {+ J
    
14. },7 c$ c9 d) w7 [5 x5 l
    
15. s;
    , ]% [  K; K2 D$ A8 k
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
+ P$ R9 H8 y6 g& B  Z1.109  秒
0 j& T7 p( V" z# P1 B  e
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

2 }! p5 x, E& [8 ^1 H- t8 `





3 X: I4 u( K; }8 l8 b
" J4 a7 P3 d3 @- W, }; t



( j. I1 A3 D: G# p

" \$ Z5 \3 A$ X+ Y6 k2 B% K

: B  {, ~4 u$ T" ^- A! L8 Z" L
. j$ b% U( ~" ^" u