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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
# @2 j: ?. j9 G1 ]
＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;6 I- T8 Z4 H9 H: @5 c+ a3 B* R/ Z
   
2. tic;! {* R5 n5 Q8 P  A8 b$ s
   
3. s=0;
   + M5 `\" g6 ~1 c8 A
4. for x=0:1000* F# L, N  @  L
   
5.   for y=0:100
   4 j  `' F( u1 h! D! ]
6.     for z=0:1003 C9 {- |9 o# |$ D
   
7.       s=s+f(x,y,z);& o! d  \/ |+ |9 ~, o$ X- O7 X
   
8.     end) G, ~* `7 V7 Z- v: t' A
   
9.   end' L- n8 I# [$ y1 T
   
10. end) |  T3 L$ v\" j. g6 |1 V/ l
    
11. s
    8 s) T7 j8 v. H' \; ]
12. toc
    # s- J  U* H2 v
13. $ r* n% A; k8 `2 R  P/ X- n( q
    
14. s =  p; _! x4 Z1 [
    
15.     6.126720600000000e+009% |% [8 Q5 r# D  R3 Z
    
16. 
    + k( d7 k\" i6 r0 l
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
' A3 r1 y8 T  `1 b' ^5 c& ]m函数：

1. %file xyz.m. n# s, S1 R: u* `+ P: s
   
2. function c=xyz(x,y,z)# Y9 J1 A& A6 b) d
   
3. c=x+y+z;6 K' N\" E2 `2 B& R3 `, G
   
4. end

复制代码

测试代码：

1. tic;1 j; Q, a# Y  Z' E  M' S6 H
   
2. s=0;/ y1 r0 Q0 k) }2 H$ m
   
3. for x=0:1000: \. `% z- x5 o+ S2 N- ^
   
4.   for y=0:100
   1 P3 X( l% H; Y: B) V$ g0 H: }. {
5.     for z=0:100
   3 z7 I\" m; {0 z4 E0 a$ H3 X
6.       s=s+xyz(x,y,z);: Q# g2 e$ c& D; l9 w
   
7.     end
   4 k/ w' Z2 X- b- n9 q' W
8.   end\" E, A\" z8 `# \; I* M( x
   
9. end
   ! i5 O6 P\" l4 u' y. E
10. s$ t3 L4 q, f6 b2 O% h\" [
    
11. toc
    * ]4 [* _. b\" B6 M% z; e: l8 F
12.   y& N9 z( i1 H
    
13. s =# S# \4 i4 r2 ]9 C
    
14.     6.126720600000000e+009* \9 E% z7 i7 o
    
15. ( n% P; J* m! ]# h  b) d
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
" l$ ~* D" [% |+ [& v1 {
8 u; ^& H' Z0 _: fForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;) L* H- m1 y  Q+ n7 j
   
2. mvar:& Q2 x+ H0 F& e6 d
   
3. t=sys::clock(),5 l* o2 X6 s; F7 l8 j
   
4. s=0,% I, p* J2 `$ v  t
   
5. x=0, while{x<=1000,
   ! ^, T7 A- Z  ~% @* Z\" z
6.   y=0, while{y<=100,
     W( t6 i4 S$ Z6 E+ K
7.     z=0, while{z<=100,
   ' R$ N( p# N0 ]2 F2 o! y
8.       s=s+f(x,y,z),
   , \' Z) p: D& d$ v% M
9.       z++: b0 v6 g3 G. P8 B
   
10.     },) Z\" m* B7 n; Z% X  g& z
    
11.     y++* X; E' p3 D# b6 T. R3 W1 Y2 r8 f. R
    
12.   },
    & a/ Q7 O+ t. n0 [7 y) J
13.   x++
    * Z- }$ H6 ^6 U4 I. k# g3 N3 y* m( c  ~1 ^7 v
14. },  [; e! q5 m* i! }0 V
    
15. s;
    9 n# d\" ~  b0 F# b
16. [sys::clock()-t]/1000;

复制代码

结果：
; W9 ~) h' V& ~) B; I* T' U- N7 {6126720600.
1.109  秒
  [' G4 \& l, Q4 ?/ ~+ D6 d0 @
- X3 U& L1 g2 y二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

; a, b& u9 W; Z9 H8 R
2 y, a5 ?2 L. J4 m

  {2 M2 c6 [8 Q# v2 {3 r3 _. T! X
$ u0 Y; G. `0 E$ E



& j9 s! e% u- c, `/ I, `" d9 A
4 r- g# {! Q( L7 S5 i! e) z: R0 V1 E; u

3 T2 |- d8 |! E; K$ o+ V" {
$ x9 Q1 t$ j; ~2 d2 M# F* o  H% _
, M- E" A; y, Q) o  y( a
/ I' `; ^' P! F# N3 E! Y