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

forcal

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

) F5 y, ?, z/ ?4 M+ j2 _) c＝＝＝＝＝＝＝＝＝＝＝＝
+ ~% I! p, p! ^9 b
! e$ S, Y( A8 a* L2 P( B8 b. UMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;' Y7 x5 _: z6 L5 W
   
2. tic;# d\" p) @* Q7 V
   
3. s=0;
   9 Y5 `9 N9 Z4 V$ [
4. for x=0:1000: f; K! ~  h+ E' Z( [% s. J6 a
   
5.   for y=0:1003 H6 r( g4 Y* P
   
6.     for z=0:100
   1 g1 ~; p$ s  D  C2 D; p0 I
7.       s=s+f(x,y,z);. x7 f% ^- z  O6 x% M% V
   
8.     end- q8 f. W3 \& \! D1 x  ?9 O
   
9.   end* u1 Z$ I7 _3 l3 \0 [4 l
   
10. end. C- w8 c# [9 I2 J, p
    
11. s* ~9 t' G1 E* \1 W3 T3 z9 w
    
12. toc
    - C4 D  p\" [\" j+ i2 @
13. 
    $ V3 A2 z; S  m8 D* G\" i' X3 F
14. s =# [4 o4 r$ w* K/ O\" f0 `6 [
    
15.     6.126720600000000e+009
    2 R1 Q9 c/ J0 N- p  x9 t7 v
16. * J+ I$ R+ [& V0 f; @' R) |
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   5 J7 ~; e6 a, h' {; H\" o! L. W# j- O
2. function c=xyz(x,y,z)3 R' l9 r) e9 b6 k7 ]# Y8 _1 i
   
3. c=x+y+z;0 L1 _0 Y/ k% g: b
   
4. end

复制代码

测试代码：

1. tic;
     }9 Q% h$ f. i
2. s=0;( h' N% O5 ?8 K/ Z7 Z, {9 o
   
3. for x=0:1000
   : P( Z  v6 E: i2 ]0 u
4.   for y=0:100
   2 X; F/ P( S1 M' W. S( X! n$ t
5.     for z=0:100
   $ ~& U$ I1 q5 C0 V9 I6 n
6.       s=s+xyz(x,y,z);
   / i\" m' s  H: U5 b
7.     end  M. d\" W' y1 r( g, l
   
8.   end1 t: e( K& }7 Z; j, Q0 s% N
   
9. end
   : P( G6 p+ F4 z: @1 n' L4 [
10. s$ ]) j  ?0 J- k% Y: ]8 v, V
    
11. toc
    ; N# n) V/ q: U6 ^
12. ! e' Z  g, `. @/ a
    
13. s =
    . Q  }3 Q1 x- H2 j
14.     6.126720600000000e+009
    1 B# N\" |% t( v\" Y8 Q2 ~1 Q
15. + e: K# E\" P4 g4 E8 P
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
/ T" D7 k; }. `. Y) m
; C; ]. A( `) W* q8 oForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;4 E9 O7 i/ n3 a  ~1 s0 B
   
2. mvar:2 w8 A7 I% O. q! B- M% @! _
   
3. t=sys::clock(),% C; n! E/ H2 t9 D+ t
   
4. s=0,/ \  o( B, f: }) p
   
5. x=0, while{x<=1000,
   ( [3 t& T1 ~4 N9 F4 t6 s
6.   y=0, while{y<=100,2 t1 V! ^' E4 N. U
   
7.     z=0, while{z<=100,
   1 E; _. ]1 u( M0 H  o
8.       s=s+f(x,y,z),
   # p3 C1 m. x# ?2 y% \1 n& f
9.       z++* f/ `  v3 Z+ m2 D# B/ h4 D
   
10.     },3 t. f1 k5 e! O0 D' s, W
    
11.     y++
    \" r5 {' u% }0 F3 Z
12.   },
    % y4 C; i5 y8 Y8 d; p- y
13.   x++9 e1 n1 `( D+ r
    
14. },; J& S1 Z, v( f' i4 K5 c; k\" x
    
15. s;
    1 v0 @( P: n3 p% v$ }
16. [sys::clock()-t]/1000;

复制代码

结果：
3 ], j4 T7 o% l6126720600.
7 J2 v2 K4 r/ v1.109  秒

0 u2 m' K1 j. k  c4 u) V1 A2 X二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

, \; ~5 G; A5 ~/ s


  X7 K8 J* G; h$ C, r( Y; W) W. r. |


  l. z* X8 N  b: ~/ ?* N  h

7 G; ]; F1 ~( V

! E. A7 b& y' U& J6 Q& l
' m' u" {5 ^' x* l/ T( j

4 r: S* P- y" E5 S+ i  V