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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
& ]. I( w4 C) h) l% b5 e8 F3 k
7 ~5 R5 G) L8 l. o& R3 D" F, |/ Q0 C& g＝＝＝＝＝＝＝＝＝＝＝＝
8 E6 x) ^% o3 ^* x
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;$ `3 A$ {7 t% @\" u% P5 B
   
2. tic;
   & ?2 M# |+ g% L4 Q1 f6 C/ q
3. s=0;
   & c! H. y) x. n+ g
4. for x=0:1000, m' C8 O4 z2 z. t; C$ H& E* s
   
5.   for y=0:100
   6 t3 ?( A- w! {. }+ u4 l
6.     for z=0:100
     k: `! G. ^3 F, w* K5 f- }
7.       s=s+f(x,y,z);
   0 b3 g% E/ d\" i0 X! e$ ~+ s* a* S
8.     end) w  e, M: Z) n& t0 z
   
9.   end% O+ @+ N$ }# K/ k1 y/ k* N
   
10. end* [% |  r& m: O9 q3 h\" z' ~, ~
    
11. s/ T, ]  X0 Z% D) D
    
12. toc- q' f0 X* \/ y9 `0 @) g
    
13. 6 [; e7 [: n+ Y6 ~+ Y; S  k
    
14. s =
    ) ^' X\" E* [- U
15.     6.126720600000000e+009% R\" B% Y% E& N# J, y& Y3 V9 U
    
16. 5 i2 W4 ?/ d4 ~- i% J9 _
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   + {' `5 p# h6 Y
2. function c=xyz(x,y,z)# O4 i0 Z! J6 J1 V9 l& h; E% f& ?
   
3. c=x+y+z;0 _; G6 u: e# Z6 ^$ }
   
4. end

复制代码

测试代码：

1. tic;
   ! N\" i+ t' [/ A, S# B9 R! w' C0 Z
2. s=0;- _9 r# O; n  }) {& L: ]% \0 i
   
3. for x=0:1000# R9 D$ `6 ?' Z7 |7 ^! u
   
4.   for y=0:100
   \" c0 R; T+ [& Z
5.     for z=0:1001 u\" I& I9 H) c& b+ E; d2 v
   
6.       s=s+xyz(x,y,z);
   7 a% q- f! k& W5 p
7.     end
   ; Z$ d5 l+ M4 _( o: l\" e
8.   end$ }\" C4 L$ Y5 B& @: e  D
   
9. end. T, a- c; H2 [9 g
   
10. s
    / A( y5 h8 s7 r+ S2 K
11. toc( C# n. s# C2 L4 J
    
12. 
    - z; u, w' A5 @# v: t- r
13. s =
    3 l2 W  `! O\" I$ K4 {' X# M
14.     6.126720600000000e+009; o; J# {$ h3 z  s+ \. e* L
    
15. 
    $ {9 i/ K: A0 q  j1 i. R4 g8 p
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
2 V" l; U$ M9 o* h
3 y6 E2 `# X# I" {Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;8 P$ J  R6 D( G4 B) @! M
   
2. mvar:. g& b: e5 K* Z5 T
   
3. t=sys::clock(),
   . c! d* m6 T0 w3 \- U
4. s=0,0 J; b) s1 E# h5 H* ^8 U
   
5. x=0, while{x<=1000,
   6 \7 v+ {  {9 w8 C
6.   y=0, while{y<=100,
   % ~9 P+ B( m8 \& c! `/ Z! C2 M0 q
7.     z=0, while{z<=100,
   9 ]6 |4 o- a4 L* a
8.       s=s+f(x,y,z),
   , j! P) b4 M) d2 Q7 \/ s3 g5 V
9.       z++2 j4 t* N& ~# t0 \$ ?# z; R0 `8 `3 h
   
10.     },/ n) ^7 T1 v9 r7 A
    
11.     y++3 g/ ]% p+ m1 _2 ]
    
12.   },9 \2 T5 C+ j. B; V
    
13.   x++
    : R  c' m' z+ k( M) i
14. },. J) Y+ _. q8 B# K0 r3 A
    
15. s;
    6 W+ S\" c. {8 F( @1 s9 t, A
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
6 _3 W% E" D* G; G4 x+ |
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

- ^, A% J8 b/ w( d7 S0 k( A
! F: {$ i. j& B( ~: C






; m% m2 B# N6 e8 j& t8 q3 [

1 X5 _, m9 x. X. K/ J6 K
4 g. U- ^: J9 w, M) R+ A1 {
( V$ U0 Z2 }" [
9 J& Q$ P$ W3 l5 g4 `" Z, ]

) N+ L( Q1 C' q- y, b