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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
8 w# j# x0 f( d; {/ q# A
＝＝＝＝＝＝＝＝＝＝＝＝
7 F' Q  y% t. L, O) m1 ]8 Q( }
9 `4 A2 K1 T0 _Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   3 U7 L* p& R* s
2. tic;  K6 `/ S( m) n8 D; N2 h' j
   
3. s=0;
   : @# o+ l$ A# ]6 z1 i$ a
4. for x=0:1000) L4 i* u, k' i# p
   
5.   for y=0:100) \/ w\" ^3 _! \* H8 [
   
6.     for z=0:100
   / r5 j8 @1 M5 j+ p! v2 B& D. H( t
7.       s=s+f(x,y,z);
   + O# m* N( _+ o8 [- K! h3 D( e
8.     end
   3 w5 C* o9 L: `# `' U\" {# o0 X
9.   end
   $ E! r1 }, a2 B6 I5 t  C
10. end9 z9 Q) b- N, s  w
    
11. s
    0 F. q( V* N& C) `0 ]
12. toc; }; \. _$ G\" W% b8 x
    
13. \" W/ t9 S/ d( o8 o( i, _9 y
    
14. s =1 L. p/ _5 T/ F4 g' G
    
15.     6.126720600000000e+009+ ~+ i5 F3 r4 ~0 Z
    
16. # W# `9 p' `! v+ Y6 D- w( G- c
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   \" b! z6 j8 i+ M: o0 \3 X2 [
2. function c=xyz(x,y,z)  C: [% z) |9 ]
   
3. c=x+y+z;# k2 D2 m\" ?. L( ~\" @- s( O
   
4. end

复制代码

测试代码：

1. tic;
   5 s. m\" e3 K& ^# Z\" T, V
2. s=0;4 E- N. P/ c) z# o
   
3. for x=0:1000
   5 O4 g  n, Y# _; w/ u
4.   for y=0:100) N, U+ s; Y* ?
   
5.     for z=0:1003 p3 ~/ c# l7 _' p7 F/ v& u
   
6.       s=s+xyz(x,y,z);: z& @. m, s+ x1 X; R* C
   
7.     end
   7 ^1 w  r( O. s! F0 w1 y3 W; o
8.   end
   : E( |\" J; R1 E5 v* h
9. end
   & l8 S+ `4 d0 x7 @* Q) a6 |) E
10. s
    8 F# `1 }, d6 @2 q
11. toc& r7 e3 j1 e\" D7 B6 Q
    
12. 
    ' `. ^, V9 e3 ^
13. s =. l8 W1 v, b' i6 [1 S, H1 b/ I+ B
    
14.     6.126720600000000e+009' d# ]4 E( U\" c: N, V+ j
    
15. ' o# K1 l. H  y4 ~\" ?4 p
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
# e$ G& i: K" l, D
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;  W7 c( ^7 s. ?7 c  N5 L5 i
   
2. mvar:
   * e7 l\" z: ?' k1 H3 A
3. t=sys::clock(),
   8 E/ T3 |- p5 ?4 Z
4. s=0,
     x  g! \\" [/ H4 T( w
5. x=0, while{x<=1000,
   - \0 R: x* E; K0 i, e1 M, G
6.   y=0, while{y<=100,8 G) u  ]3 T2 F1 ^
   
7.     z=0, while{z<=100,) p+ l, a4 y- o! v7 d
   
8.       s=s+f(x,y,z),
   , S* G0 m; f) C% T, R9 P3 X
9.       z++
   ' _9 h  w2 V. _. f& ~+ N
10.     },
    4 a8 q) B$ m) f7 v* `
11.     y++
    ) m. f# @+ H6 N  l; m
12.   },4 I) j\" F$ N/ @1 t2 T$ t
    
13.   x++* `6 c\" c- f- m  Z/ N  x$ N% K
    
14. },6 \) [; d& O- l2 d, E# k
    
15. s;: h' L  v, [0 |) h7 T
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
/ p5 A# O$ K5 k8 }7 L7 ^+ P
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

% O" ?1 y- C0 a
: {; a9 B! j8 o6 M. A& V
5 g  M* W( q( X; X3 @7 \
2 a. m0 X! ]4 u* Z8 z. l' F+ @1 W

5 h$ B+ V$ l, K& y

; K3 ?6 `5 A# a, m/ ~5 L


2 m2 E4 |2 e+ V" B9 {$ q* k
' @8 d+ {" y) R; L! Y' P2 _! i

# ?/ D8 G$ a& W0 U/ B, w9 j

, [# b4 W/ h7 c' A  U7 U7 b
! K$ i+ s" r4 l2 P. [* h