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

forcal

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

/ r, N! Y3 ^- r4 f＝＝＝＝＝＝＝＝＝＝＝＝
6 |. ^. i0 {9 X* P
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   % l. C- o8 y% W9 X! s0 I
2. tic;: h* O7 x1 A& M, K* c! w4 }0 s/ v; K
   
3. s=0;- h8 h) n( F, }; h; ^: w
   
4. for x=0:1000: V: t* `+ ]# b8 ?. s
   
5.   for y=0:100' X8 s( Z3 W  |- n) Z: U. l& d/ i* A
   
6.     for z=0:100, l, K; ~/ x* ?9 V2 i
   
7.       s=s+f(x,y,z);
   ! k2 b! i+ k! F+ c+ ]
8.     end' k* G) M. e  o
   
9.   end* `/ e$ g' V' Q+ T3 r\" x
   
10. end4 V% d) ^+ S3 L
    
11. s! o( c4 d3 b9 H4 J/ y1 g
    
12. toc- M- q6 D! d: U' W  E/ u; i  S1 B
    
13. 
    3 ?/ J9 j' M1 M2 m0 P/ m% _$ J
14. s =
    ) C' z) @& H: p: g
15.     6.126720600000000e+0090 U2 \( O) }8 U7 ]
    
16. ; V, B2 [6 R) I- @9 v1 g+ Z
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
0 X- J+ s9 y9 f6 D1 T1 G4 |% wm函数：

1. %file xyz.m0 [$ A/ _+ ^8 O5 P6 n1 N
   
2. function c=xyz(x,y,z)8 ^. M* O& D: [$ ?6 P1 r
   
3. c=x+y+z;7 O. C7 x& |# D' f\" X4 n* W( [
   
4. end

复制代码

测试代码：

1. tic;- ?* H7 N- m; h
   
2. s=0;
     o8 v% s. t8 j
3. for x=0:1000
   ; k' [; ~: ^8 m1 P3 }( K$ ]
4.   for y=0:1000 u2 M+ }5 C5 a- E. W7 x
   
5.     for z=0:100
   4 \( K! J& ~) x
6.       s=s+xyz(x,y,z);1 X+ C% M5 j) m0 @  `$ ]\" W2 Y
   
7.     end
   7 f) c* T\" u# j$ M
8.   end. _8 F6 V6 Q/ h- I! \* Z
   
9. end
   ! d  D\" |/ a$ a
10. s$ |: ^2 f; T  N/ s7 G/ K
    
11. toc  e( o, V% Q1 |' R  t$ f0 b
    
12. # f! U/ w8 u0 |( {/ }
    
13. s =% o) y5 B. ?- b1 k+ l2 p& N
    
14.     6.126720600000000e+0091 l; }- [$ Z2 H' I% g5 g
    
15. # N\" F, H+ Y/ {
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
0 v$ `4 C/ {# T
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   0 o* b. w8 k: U1 A7 O# R; h\" _
2. mvar:; U& X. b& F7 c% O5 F
   
3. t=sys::clock(),
     H& z( v8 t2 j- [
4. s=0,
   / ^! K8 B- I% b9 t
5. x=0, while{x<=1000,
   : G' c8 s\" A8 I1 X& y) J
6.   y=0, while{y<=100,2 i. n* b) r/ B\" M4 z  L\" ~$ v
   
7.     z=0, while{z<=100,; y) N' k$ S& |7 Q$ S/ U0 G$ ?
   
8.       s=s+f(x,y,z),
   3 B- m! P' S0 \, Y3 a! E
9.       z++
   & M& w( K\" ^9 v: R: {- r
10.     },
    3 h$ `2 }1 F  ?4 \1 Z8 a% N2 X
11.     y++; ~) p: I+ E* Y
    
12.   },5 b! H  C6 {+ y; U
    
13.   x++8 h7 ?! F6 J- R  t
    
14. },5 m+ Y# s5 b  Y
    
15. s;9 E3 f* x/ x: U) g+ a
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
6 `5 w, G9 _: \  t4 P
0 Y' i  H: t/ z0 J, {: D二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
* f) D# \! W. V1 e

hxftntA

; v& O2 Z" J) d* z5 ?2 C


) o/ g5 t& K( U( X( D+ [3 x2 _# x' a

+ {! u9 h* a  k& X. h/ V( w* f% k

$ D9 [% t# ^9 p8 K3 v5 K

) W0 ~: P4 }, q- `6 M8 G
, l$ M+ p: @  Z' K
. _( ?9 ~0 b$ i( O" v+ P- \