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

forcal

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

, V: z/ t9 ^. a; f' B＝＝＝＝＝＝＝＝＝＝＝＝
  f2 V( H/ D2 B5 k2 |1 J- E6 w: ?7 f
. Q8 A: `. T' g' T9 N% m4 oMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;9 I7 {' y$ r5 u' S% B( M7 O
   
2. tic;' n/ j) Z6 D( s# y, K8 p) K& z& Z
   
3. s=0;& S& I# }% k; |- S
   
4. for x=0:1000
   / b, ]: z8 }0 n7 i- w\" {4 i6 X
5.   for y=0:100. i+ J, ]: v0 J) O; ?, ?
   
6.     for z=0:100
   / y. H( U- O, \/ Q- j1 o: v
7.       s=s+f(x,y,z);
   ( Y# M& R  Y- k
8.     end
     u8 F1 f4 B- V4 u) W* Z
9.   end
   3 f5 H6 W$ {, Z
10. end9 y- W! f1 t2 m* n: x) L4 r# P% f
    
11. s: A\" w6 u) _% H
    
12. toc
    , R# N) q: ^# P
13. 4 p. y  C3 l8 I/ g9 J# U2 h
    
14. s =# b6 V4 X  s. n, n
    
15.     6.126720600000000e+009
    * I- d1 g) a. \' A: m2 x
16. 
    # O' {  c; L1 z, r0 m, m7 V, k, Y
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
$ `! h% R0 ^1 u' c5 q) x3 r9 u' B9 b6 pm函数：

1. %file xyz.m3 _% e) s. f% [
   
2. function c=xyz(x,y,z)
   $ n6 W/ R* _1 a* N\" D' F; C
3. c=x+y+z;
   3 n' @\" p8 N  Z
4. end

复制代码

测试代码：

1. tic;3 n+ H& t. I$ U8 |, C. W% I
   
2. s=0;
   : S3 K3 c  @) R/ ^5 x
3. for x=0:1000
   ( [4 n% C! d' \$ r' G8 a
4.   for y=0:100* e\" y# A$ ^8 a2 A  a' V* {
   
5.     for z=0:100
   ! K5 \1 _6 B( D. ^4 P6 x
6.       s=s+xyz(x,y,z);
   # f2 u& r: F\" V9 O0 q
7.     end
   $ o- n\" ^$ j1 C: H, [& K+ K
8.   end+ ?2 ?4 c* c: k: @, ?
   
9. end
   ( Q0 J8 J; k8 k\" c* `, C' }
10. s9 K# u/ e1 ~7 i% }
    
11. toc; k) i# A9 \3 y$ w' E* O! m
    
12. 2 ?4 Y/ }, P% N( n0 \5 p9 X
    
13. s =
    8 V6 L# T' G+ ]( P$ N0 w5 @5 r
14.     6.126720600000000e+009
    1 T( R- d' X2 T
15. 3 Z- h& H\" [$ s& D1 r: B
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
- F: K* O8 v! T2 t* m8 w
5 b. o% v: ~" t7 G2 lForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;, R/ r+ r, M0 l( h
   
2. mvar:! H, i: g5 Z9 ^: `  l) O
   
3. t=sys::clock(),9 m3 R5 j! T+ ~3 L\" F
   
4. s=0,
   2 u% z: P1 S( D8 D. G! |, \0 ~
5. x=0, while{x<=1000,1 @' P: r' ]) m5 _. Y: r% _
   
6.   y=0, while{y<=100,
   2 ]1 d' E! q! z* E. }
7.     z=0, while{z<=100,  V3 d8 _) O7 W5 ?/ M& q. Y& @
   
8.       s=s+f(x,y,z),
   . R: N! I6 T, \% g
9.       z++
   9 a( O& d+ N. z5 R$ T* d
10.     },: Y4 I1 Z8 k4 i/ g- k+ Y
    
11.     y++4 Y& `2 N; Q( N/ P; D\" R  K6 B
    
12.   },0 u2 q; K9 O3 B% H
    
13.   x++
    ; w+ s2 q) |- _/ q5 {% R7 z
14. },
    : K4 |) w$ D* D
15. s;+ Y' a1 E$ E\" L$ v- ]) W1 t/ X
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
9 K5 R4 X! t1 c' A' n, |; m& w1.109  秒

二者效率有几倍的差距。

hxftntA

* d9 J9 I* O  C
9 u5 r0 t& n0 ?7 w' j

; \, ], p( i; E% X7 t5 [5 W

) f2 K8 [9 y1 i$ I0 }7 h/ g% d
. ?- K  p* I/ J3 `8 y
6 H  w2 b4 K+ E6 q# e+ F: f( P

" c/ u8 e+ b8 l& w8 s2 @( }7 O, k- [
8 r0 A) b+ ~7 H. C

  S! @$ J) S, k0 Z( \3 s1 J' r
  ^6 L  V. I( ]( v3 K
4 F2 C8 \1 v) {9 E
! }+ _2 b, G  o+ _* Q* _7 a
. N0 [( m3 D( R* J3 P5 T

ehi28

嗯,不错,支持一下.