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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
9 S" ]% R: t! C8 Y8 g3 \
' b  q% M* ]5 L1 x  `＝＝＝＝＝＝＝＝＝＝＝＝

: Y- q- `5 S3 C: x+ {Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;\" ], S1 p9 J$ Y- ^\" A1 }
   
2. tic;
   6 m4 I% Y0 B6 q, T
3. s=0;% P1 r% Q6 O2 [
   
4. for x=0:1000
   1 e2 s  Y6 o$ ?/ y  i9 b3 u) m( j
5.   for y=0:100
   # \8 ]3 {, l7 t1 E* `; m: }+ g
6.     for z=0:100, }( d2 C3 u6 O1 p\" y0 F
   
7.       s=s+f(x,y,z);9 p, s, @: x  |+ ]1 s\" k
   
8.     end' ~' p$ P0 o2 z  a9 h3 I2 `
   
9.   end. j- l2 @7 f+ Y  n6 n
   
10. end7 X0 w8 c) W, N\" i\" H4 ~' Q
    
11. s
    4 O0 o3 [) J6 Q9 g& l
12. toc
    . A( g2 P3 S( A\" o8 {
13. 
    & o1 z0 c% r- _8 z$ w1 s0 [. }
14. s =
    . u+ A- ?- w/ v& n3 W% w. t
15.     6.126720600000000e+009! |8 _\" W1 J+ k( O) [\" p. ~
    
16. 
    9 k5 j7 n2 y2 @/ z
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
  [! m4 O" i5 x% N+ c. Gm函数：

1. %file xyz.m
   ! e- k8 n2 l1 G6 A
2. function c=xyz(x,y,z)9 U2 o/ z; v; U9 M( X
   
3. c=x+y+z;
   ( K# ]% ~7 x! C  \7 G- V; ]9 U\" J. V
4. end

复制代码

测试代码：

1. tic;9 ]: Z; }! K) u9 ]) }
   
2. s=0;
   6 Q. }5 Z# h! D6 b& h
3. for x=0:1000+ o\" ~\" a- b% E% i- m4 N) C! N
   
4.   for y=0:100
   ; y( F1 q  a: ?' g4 E- A4 Z) z3 S
5.     for z=0:1009 D. D7 K\" M! F
   
6.       s=s+xyz(x,y,z);4 g0 M! Q. f- Q- _
   
7.     end
   6 J3 w+ }0 j2 r: W# w
8.   end
     f+ x# H9 t) K9 D# B4 H
9. end- _  r6 P# r% ~, r$ \' A\" q( ~$ z
   
10. s: d3 `9 N& X- ?! z7 N
    
11. toc
    & ?* X  K+ \4 ?/ I
12. 
    $ e, k4 C. ~  @4 @% j: |+ b
13. s =
    ( B, U* |' N( G9 s
14.     6.126720600000000e+009
    9 |$ \' h\" @6 v: z# I/ s  i  |! b
15. 
    . Z! w2 ]9 R7 u5 j
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
3 E8 k3 X. I9 d) Y5 E/ \
6 X4 Z" w1 _7 I) F  o* C7 M; hForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   8 M\" B5 ]1 M# M5 g& G$ t. A
2. mvar:
   + z! h0 q* C3 F4 _0 I- y* Y: H0 W
3. t=sys::clock(),
   2 T, ~& ^- b5 m# d
4. s=0,' [& W3 I# J/ D3 Q9 w
   
5. x=0, while{x<=1000,
   8 ~\" ?( J+ y# u6 v; g
6.   y=0, while{y<=100,! l$ o6 u/ \  J* [9 p0 c
   
7.     z=0, while{z<=100,6 u9 ]2 O4 v7 t
   
8.       s=s+f(x,y,z),
   # k6 K# P: T, T6 L0 q
9.       z++
   ( M, B5 W2 G\" V  t  F  \, d+ L
10.     },7 _6 k2 y9 S& e\" c$ l/ S  g
    
11.     y++
    - Q7 C0 T* p) e# z
12.   },# T' S5 a- I, g$ ~, @' _
    
13.   x++6 c+ _& D/ O/ F0 {' c
    
14. },* L7 [) z- ]) n. o
    
15. s;
    4 w% N5 E9 r1 C
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1 N. z$ A' L# U& e- w; _' I1.109  秒

9 P* O. F& N3 p5 M' W8 D% l* Q二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

: q" x" Z5 L8 w0 s* b+ ?1 N4 @
' h5 }! ~0 C3 X  z  o9 [) m9 |: b


7 C, F* t( e# c; l% I: A+ e9 _

. {* e2 h2 ]6 h" t. P5 H, J
+ w; C, L% s. I8 q8 u! h

: ]% k9 {3 c2 V# k0 _( ]* `



6 l5 v) i% c: X/ @( j) q+ O
$ G/ x! B4 j. k1 s+ w) k
0 f3 _- m4 F4 E& c+ t
' b4 L; Y1 `9 {( ~