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

forcal

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

( Y1 H* l) F; ?. ~＝＝＝＝＝＝＝＝＝＝＝＝
7 w5 U4 M" T" V7 `/ [8 o
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   - D* [1 A; B! m7 a, j
2. tic;
   5 P4 Y- b9 I2 H/ l( _9 U2 I. Y\" [' J
3. s=0;
   * U\" q; Q/ e3 l3 i# f! Q8 M
4. for x=0:1000
   , w- m% @: z4 v% ]5 D
5.   for y=0:100' o5 J' }9 T1 v0 W1 y( Q* K4 S
   
6.     for z=0:100! o2 Z$ h$ E, r/ C& V* Y3 g, r
   
7.       s=s+f(x,y,z);/ ^1 U) B' t8 ?; k5 f
   
8.     end. j2 O- U; e- _( u! B
   
9.   end
   6 l6 M0 a8 B! E. P2 H
10. end
    1 k+ p% w0 n* V) N9 V( Z) w  ]
11. s/ |; w: D' J0 D
    
12. toc
    ; \! V' @% k3 ?6 A/ G0 r4 a
13. 
    ; t2 A+ x\" h( O
14. s =
    ; K6 e4 g$ g# x# P2 z. w
15.     6.126720600000000e+009\" {- @) ?! _8 l6 p5 a. g) D. m
    
16. 
      l5 f2 u' P9 \3 o4 a
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m; i# B3 d) ]- G* P7 W9 u
   
2. function c=xyz(x,y,z)
   ) [0 D9 w  H6 n$ f) \% P
3. c=x+y+z;
   ' C: u# S  |3 F! M0 j* B
4. end

复制代码

测试代码：

1. tic;- E/ \' o* a5 X; B0 P! x' v, @! @
   
2. s=0;  k6 n$ i\" d( L( K/ K! }1 I
   
3. for x=0:1000
   $ X# ^* N* O7 W8 _( e# y
4.   for y=0:100
   2 o: t: C! I: s4 {9 ~\" S
5.     for z=0:100  S0 t1 h\" d$ n
   
6.       s=s+xyz(x,y,z);4 Z& d: ]. Z0 ]5 o9 S! a
   
7.     end8 c' n6 K* U& V+ q, V
   
8.   end. J2 z9 ?. Q  K9 O: S
   
9. end2 [7 c' f' r7 Z; y9 \0 \' ?7 r5 }& N
   
10. s- q  w+ e# r$ H
    
11. toc
    ! c) }( Q2 W& C- l: v, h7 l4 p& J
12. 
    6 i0 a8 X/ _4 _$ Y0 [. Y, s0 t: Q- y
13. s =
    ( G1 Y5 I* k+ M/ Q( {7 h' V
14.     6.126720600000000e+009
    5 H\" C( h: c* l1 C  _8 w0 d
15. 1 r6 p4 o( V5 C  f3 Z
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   3 @3 O( B) f; @% P4 ~
2. mvar:
     [9 V2 U( \% H3 c, o! j
3. t=sys::clock(),% v5 w/ l0 W4 S% u1 T
   
4. s=0,( E: L# K2 O! h
   
5. x=0, while{x<=1000,
   & q! ?# h* a& z! J5 N6 i+ S9 ?) s
6.   y=0, while{y<=100,
   ( y! Q) W+ \* G9 H- G1 P: t
7.     z=0, while{z<=100,' N, u3 N2 V0 v' Y) R/ H
   
8.       s=s+f(x,y,z),
   % [+ h! G; e* i4 ~$ @, Z; t
9.       z++
   % h0 I9 X0 C/ Y& \% [; O# X
10.     },4 ~0 s  ^2 E( Z& S0 d4 P
    
11.     y++% Q) `: Q1 V% D5 l
    
12.   },
    6 k/ P& N! O+ l* Q
13.   x++2 N: {7 z+ I( P% V9 h, x
    
14. },- @  o8 g$ E& H9 K2 H2 Q
    
15. s;
    6 C0 o% \8 W# V# [' V$ v5 r
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
" F5 ^4 i3 w! D# S2 P/ w

hxftntA

& F9 R) g6 f; B8 r
  U( A" \9 q5 d1 J8 R7 W7 [
5 |4 Q0 r* A3 H8 a
+ {0 y% ?' t* p$ u) _8 u8 r

- o; y! v; H1 F: {5 e" Q


" [3 H0 B9 \' W' y: \

3 z' o/ D- O! A# `! r4 d, Q/ N
' {0 h) `9 ]& t3 s  a


: o% S3 F5 D- V4 K5 X( G