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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
: Y5 a, Z& ~. {$ ^
- ^* S5 g9 q/ [- B/ P7 zMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;+ A; e! v% Z% t
   
2. tic;
   $ y( V5 G+ s3 l0 j\" h
3. s=0;
   3 m% D- S\" e% B) T6 w$ ^4 i
4. for x=0:1000% I/ ~; O1 p  j$ a
   
5.   for y=0:100
   ; w, d8 d5 K; U  J  f/ s- Z  L3 i! W
6.     for z=0:100
   : K! F5 o/ W9 o9 F1 D8 J4 b5 U* _
7.       s=s+f(x,y,z);+ ]& @' _( {9 k: }1 n9 O- D! \1 X
   
8.     end$ X! H- Z1 X' M( x4 p6 S7 t& K
   
9.   end( x* O3 ~! J/ d' L1 R
   
10. end
    & r( F- |7 z6 x% k\" D8 [5 l
11. s; `. l9 F9 }# t- {$ O+ {7 g% a
    
12. toc
    % n- s. J9 W% Q. M2 q  b8 `
13. 
    6 `4 N5 E, s' U# Z# m- n4 I
14. s =4 w( N* ^1 F# T  v
    
15.     6.126720600000000e+0092 Y\" l  F* n0 j0 q% n  F
    
16. 
    2 T  _  {6 x0 t* `6 J, u! K# x
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   2 X/ Y/ ~7 C) B) c3 @' S
2. function c=xyz(x,y,z). P2 R6 R# b7 a! S
   
3. c=x+y+z;- e1 c5 Z( H* ]0 Q
   
4. end

复制代码

测试代码：

1. tic;& c\" r2 F/ R% l# z1 C8 K, x. k
   
2. s=0;. y1 o, V7 I. S( T* P( f7 U
   
3. for x=0:10000 z# s3 X' x# s4 x/ C9 d
   
4.   for y=0:100% L, g. a! r; Z; T
   
5.     for z=0:100
   9 ~& i, c; n/ U. K
6.       s=s+xyz(x,y,z);
   3 ]) D4 p, i( W- q2 V+ @0 O; z
7.     end
   \" R, a; {0 i* e1 h, [
8.   end
   7 e8 D6 t, G( E\" I( l8 U/ Q/ f
9. end
   + S6 z, W0 W# _1 b! M
10. s
    3 \) A7 t8 g5 \  J& p) e. ]
11. toc
    5 y* j1 o+ E# }; x' u1 ~
12. 
    , o\" k2 ]8 c; Y7 N- `& u$ d/ Z
13. s =
    : q5 o  f( S  F) t7 Y* u8 d
14.     6.126720600000000e+0094 U8 d; H6 l, g- L
    
15. # S$ U3 W# k/ M8 W
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
0 @7 G# ]6 z. T
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   & S' j, \/ W2 x\" a# @% i+ ^: `' @
2. mvar:
   , W- O9 m) @. V0 g/ v5 P! x
3. t=sys::clock(),! X  [/ E8 ?; Z6 \\" x\" s
   
4. s=0,. f- c+ y5 U# C( W5 V2 B
   
5. x=0, while{x<=1000,6 d' q\" {1 M; F0 A, P
   
6.   y=0, while{y<=100,7 s# Q9 N2 O, U7 S9 K$ k  b+ p
   
7.     z=0, while{z<=100,
   # D8 [( `2 T/ b% a) w
8.       s=s+f(x,y,z),
   ) B9 N' K6 A8 x( J* K/ m
9.       z++
     ]8 f& i; P' e' j- H\" J! c
10.     },
    * u2 E( z5 E8 z0 E3 o/ s
11.     y++
    & r1 ^2 J3 a( p) {
12.   },
    7 w4 ~& h; j. X; ~6 @8 r+ n
13.   x++* \$ `! b1 I! x: ^9 |* e, q
    
14. },' o& S& M6 T1 _+ O, T  x. M
    
15. s;
    * P' ~, }' v' u8 p5 m& C
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

3 Y% T* t5 c/ ~$ f* D. ~
7 @" w1 F  Z/ y* S



' A5 ^( S7 L8 `5 Q9 }! v$ l2 x
9 m/ X1 S& f! {0 |. K' G
/ B8 Q1 D1 \$ m4 s* r: k
3 s4 K( j; Q! W: r
4 O( p! G5 a2 R! v6 @



& V1 x  L9 v4 V+ r- y" d, c