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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
  C3 @" A! T+ p6 a& w$ E. k3 b" }) \
4 u3 O) y) W" f* \0 r. c＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;' D- L' q* f( p* L; k
   
2. tic;: E( p: S, \5 U- a3 W
   
3. s=0;( u& ?, S# ^  t1 Y& p( z) x$ ]
   
4. for x=0:1000
   1 E( t1 m. B\" ^, N
5.   for y=0:100
   . k7 `9 t0 |! i( E& Y\" ?9 d
6.     for z=0:100
   # }0 y* ?$ Z4 p' `
7.       s=s+f(x,y,z);4 C/ C3 q( d) ?( d2 E5 c1 F% v
   
8.     end\" V* G9 m7 Z8 H% P2 r
   
9.   end
   0 O, a' W3 P$ _, L) a9 b$ e0 S
10. end
    # C+ e. d' E: q& a+ i5 }
11. s1 j  [, i5 z# q- d
    
12. toc
    1 J- o/ n\" _' Q* E0 K
13. ; C: c/ E6 Q& x3 J5 a$ h
    
14. s =
    5 ^$ C) U\" u, m1 J2 F3 ^+ {
15.     6.126720600000000e+009/ N: C* S' z8 V
    
16. 
    3 z& }% N% W9 p* a
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
( }* {1 [$ a$ bm函数：

1. %file xyz.m7 \7 m9 `+ Z3 B  X2 O% {  a
   
2. function c=xyz(x,y,z)
   3 K( C1 h\" U' |, @
3. c=x+y+z;4 T\" U( J6 |( z7 O9 Z
   
4. end

复制代码

测试代码：

1. tic;4 g; z) R9 T+ q6 U' J. e
   
2. s=0;7 J+ H( t9 O* V
   
3. for x=0:10006 V. [% b9 D+ u9 {! p; X9 G. l
   
4.   for y=0:100% U) [0 @% d  r6 N% Y. A; ]9 w
   
5.     for z=0:100
   8 y* o% X- B% Q6 F) @
6.       s=s+xyz(x,y,z);9 B! f: S( C4 C: r6 q& g
   
7.     end$ f$ ^  u3 u\" u9 c. I6 y9 f
   
8.   end\" H) l6 ^5 Z3 Z2 o7 S7 q5 ~' D
   
9. end
   ( f6 P+ J$ g$ w* v( N1 ?5 \: P
10. s
    & c( g0 D, a  A% z6 a# f. g\" v
11. toc
    3 a+ L8 s  k, b1 c3 q1 z8 K
12. \" ?  |, F9 k; N1 n; H# y! A4 {
    
13. s =7 e* x9 B0 C* o! O' P0 |: p
    
14.     6.126720600000000e+0097 y6 k% p( V: V9 E1 G3 H- x
    
15. 
    ) @7 y  C& w* s( F5 P3 B
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
3 ]0 L6 D) Z+ F
0 T5 U3 B0 g# K# u5 ~& F3 c6 t$ tForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   \" j/ P! ~6 y, a* I# c) n: w; o- u
2. mvar:
   8 d6 c# ]\" F$ Y
3. t=sys::clock(),
   9 B8 F0 {  M* {  [/ a6 I
4. s=0,  k' S3 J% h5 z) @
   
5. x=0, while{x<=1000,
   4 M9 }: s# z$ ]\" J\" Z) Q% }1 F
6.   y=0, while{y<=100,
   ! U/ l1 w4 C0 }) V% S
7.     z=0, while{z<=100,/ m5 K$ x- U! l( I& z* V9 l
   
8.       s=s+f(x,y,z),
   5 S6 W! p5 p( e+ ~0 i
9.       z++
   / r8 u+ c  w) ~  Y1 c' ?
10.     },
    ) s/ n% X. j5 t\" x, K: u0 l  {4 y( L
11.     y++
      M0 @4 K8 t# s6 J4 r# E5 j4 \% H/ X# j1 r4 v
12.   },( y+ j2 c/ g\" |6 y- K
    
13.   x++: y; c  m' x' [3 Y
    
14. },2 x( ?/ O5 a  ^
    
15. s;* u- Y- @/ J: e# e
    
16. [sys::clock()-t]/1000;

复制代码

结果：
1 Q2 o" a$ Z9 B# E" ~( A6126720600.
4 z" A1 Z0 z3 m/ W: l0 P1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
: b+ x( w: }% y6 H! n3 E

hxftntA

7 z# j7 f" ?! u$ ^
1 n9 P* C# E4 g- [




4 T  r& H& w) w
5 ]3 F- O1 \. s' {' {1 J" b
% t0 y" O% e4 L" l% Y% D+ T




& `# q5 d# ~6 N6 |2 p