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

forcal

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

8 t8 W9 w! e8 o/ s3 G' w/ Z8 e4 b＝＝＝＝＝＝＝＝＝＝＝＝
' ?( F1 J, V( ?* e( R& u% q
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;, Y1 J  \4 b% L/ j) N3 r% ~
   
2. tic;: U3 G1 G/ ]0 Q3 R- j
   
3. s=0;7 G2 p9 U; i! E3 r1 c
   
4. for x=0:1000
   ; e3 s' n2 G# c1 `$ r
5.   for y=0:100* ?9 q& l3 n. ?2 C' ]5 G) [) `
   
6.     for z=0:100\" Y. H9 B* n! _4 f
   
7.       s=s+f(x,y,z);: F! x% a) O  i, p4 W; T
   
8.     end9 d4 I4 C4 ?' w4 |5 t6 y& C
   
9.   end
   ( I2 u% A8 ~# u( c
10. end
    7 f* `# g' W3 D\" T
11. s
    9 I1 ]+ C* h% q
12. toc
    ' o+ X\" E% ~# A' |+ v+ R
13. ! F4 n, d* I3 r8 }. T. T
    
14. s =
    * w+ a7 g6 i: K2 c\" l* I) ?
15.     6.126720600000000e+0097 ]+ v) ~$ O5 V5 B7 {3 c+ i6 q
    
16. * S* z4 F) N: z/ w0 c  g. F% y
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   ) j  m& ]) P( a. }7 K\" h  R
2. function c=xyz(x,y,z)
   8 a: N/ K7 Y% F) l\" W+ H# I
3. c=x+y+z;, ]  q1 o2 z# q/ W% F
   
4. end

复制代码

测试代码：

1. tic;5 b1 M. u0 m6 s  H5 B
   
2. s=0;' ^& l* Q& {$ N* d  W0 }) g
   
3. for x=0:1000
   . h5 P/ o/ f$ P$ s
4.   for y=0:1001 L5 q\" w* ^/ y. {% m1 g
   
5.     for z=0:1003 r' r8 P& b0 m7 b3 x9 l5 H8 R
   
6.       s=s+xyz(x,y,z);
   \" q6 r  n) U7 `- t/ ?! {* k
7.     end, P% p, n+ X, L4 M3 l
   
8.   end  W3 K3 }\" F' j! t6 R8 M' {
   
9. end
   0 ~1 J# k6 Y- }. K- c
10. s9 L* N5 x) J, g& _* p3 }
    
11. toc- I. k$ l& u! Q
    
12. \" \7 p& U$ t! J% E  b
    
13. s =
    ' }  }8 |% A1 N8 Z. @
14.     6.126720600000000e+009; }9 S) L4 ^' P( U
    
15. - H. Q) C: H- U- V6 X  [
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
1 D+ x: n* ?; W6 y4 S, j" y
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;4 ^) V$ u# g8 u# b6 B. c
   
2. mvar:1 v, I4 e+ Y3 U  K- E( v8 ?
   
3. t=sys::clock(),5 {: z$ j. N0 j
   
4. s=0,\" _\" g- f; b; T# x# K
   
5. x=0, while{x<=1000,
   + Y9 g' o( O* E! u. G
6.   y=0, while{y<=100,& y& D/ P\" |& ~/ g( k4 o# J' |, ?\" p
   
7.     z=0, while{z<=100,
     w0 P! G' D\" s; u
8.       s=s+f(x,y,z),
     |2 U) F; N8 M7 u. `2 Y7 L
9.       z++
   # d0 m4 b& g$ s' e6 x% E( L
10.     },9 a. c/ ^& d) I. v5 V
    
11.     y++
    4 J  L0 B8 `2 W, ~
12.   },
    * M3 B/ o\" i+ b- Y
13.   x++
    : T6 y; A2 o6 o3 X
14. },
    + L' V6 q% l\" g7 z
15. s;/ f- }9 V8 I\" q3 m4 P
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
3 B* D+ }9 u* i1.109  秒
, R9 q) ~- a  b- k1 J
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
0 t% k' y+ x1 q7 F% a/ n: p

hxftntA

$ ^/ y0 C3 @) _7 T6 \' o8 k0 |- x3 ^




( l) D! C7 a; f+ M% K

$ s! p- [7 ^; |1 `5 j



4 h: F6 ~& y9 w2 Y1 ^& N
* I  P' g$ |2 u. Y6 s, d9 `
) |- ]- ]: f! u- j, }" p. \