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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
. }, T6 X8 p8 A0 F* S* `. o) N" W2 A
＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;' J/ A0 k: U2 z' E) h3 }4 W
   
2. tic;
   + y% a6 n3 ?: Z: Y, h3 A' e
3. s=0;
   7 A* p4 j: T2 C) ?! ^3 W& I
4. for x=0:1000; K8 f0 X0 t1 d& a\" b
   
5.   for y=0:100
   8 f# Q) Z8 M1 g8 o& g
6.     for z=0:100
   + E7 C\" o+ M% f7 G# [
7.       s=s+f(x,y,z);
   1 q6 m( `; w\" p; m( M' V\" c$ X
8.     end; U- W\" A0 v3 |& O0 a) @6 E3 z/ f6 n
   
9.   end
   ) O6 j: m* J( F+ n9 x
10. end
    6 `4 v* K& d6 G( m
11. s
    . M\" S5 T+ K+ G1 n
12. toc( m: K7 k# m% E$ w4 i% T
    
13. ) `0 ^1 e0 k' o9 r! ?% _/ p
    
14. s =
    1 O7 i! ]% F; M8 g& o. Q: {. \
15.     6.126720600000000e+0098 \6 t/ B6 {: \! R
    
16. ) A; G; i9 J! b0 @9 {7 \
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   ; N( w5 v- D9 V% C0 G% q
2. function c=xyz(x,y,z)
   ! k8 E! ?8 ?( Z% y
3. c=x+y+z;4 C\" J) B4 M. [3 Z/ S5 R
   
4. end

复制代码

测试代码：

1. tic;8 ^) [; Q\" J( p4 ^1 t  }: w2 I7 e
   
2. s=0;: ^) j# z/ v% H' S/ |: _2 _5 B- ?
   
3. for x=0:1000
   / s: K6 g2 L1 p$ R\" H
4.   for y=0:100
   % l' y  z  A7 g4 B( T7 V
5.     for z=0:100
   1 c- U$ W3 S: j2 G2 [
6.       s=s+xyz(x,y,z);
   % ~+ ~* z. A5 @. Q& ?# g
7.     end
   ' |4 _. i- F- g0 _+ n
8.   end( Q! H9 V5 S' l$ y% u+ C\" ?2 F
   
9. end
   ! W: u' P! ]\" C$ ^; _/ Z7 ?$ [
10. s2 r; p( K; v\" M/ ?$ [0 v; ~
    
11. toc( ]7 x2 P+ p! E/ Q
    
12. 
    3 n  ^8 c& s8 o5 M
13. s =
    + Z3 r8 ]# f% C* ?$ l# s  E
14.     6.126720600000000e+009% E1 n% q4 v8 D/ G
    
15. 
    5 x! `  Y8 I+ s/ q; O( z
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ( H8 q' i. r9 d3 f' X
2. mvar:# n( R. n+ n8 o: `9 f
   
3. t=sys::clock(),
   1 R# {( t2 t( J$ W+ ]2 B. Z; z1 M
4. s=0,5 D+ x0 a1 Z) t. C8 D; u3 ?
   
5. x=0, while{x<=1000,/ ~\" S3 A( m' }* D$ J$ w
   
6.   y=0, while{y<=100,
   \" t2 O1 b: M# o& J( ?& E# ^5 {
7.     z=0, while{z<=100,3 O8 G. n0 w; B- p
   
8.       s=s+f(x,y,z),
   ' z% p& u5 k& Y& }
9.       z++6 U5 _; N' V9 x. d. j) G% j
   
10.     },2 M  e9 H\" Z8 V% l$ R. K0 r
    
11.     y++
    4 W, T: p! _! y5 Q: V3 [
12.   },1 m. B0 n6 P+ g- h1 n
    
13.   x++9 }0 n. |! W7 L% F+ M6 ~6 f
    
14. },2 r6 y$ z7 `7 P9 v7 ]
    
15. s;
    ; k2 w/ E% _8 ^; c+ {
16. [sys::clock()-t]/1000;

复制代码

结果：
6 ^( g2 F( u- ~1 q, ?( r6126720600.
, J. N& l! n( r4 O6 f1.109  秒
6 y3 R6 \1 q5 _! `5 a8 {. q
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
3 [  K( a4 R% @6 j6 l+ q- W6 P

hxftntA

* w& g8 d  _2 H; n# Z

; l5 F+ k5 ^/ I# V( K" V



: Z0 S2 p: s, B3 i" c# s
# u$ ~, _3 N( F





0 o. q; \+ f! k- D3 T+ l
$ s9 E1 F  w" y# l

* A2 Q5 j. w# d! m2 V