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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   \" P5 C3 ^; d9 v/ u3 p! ?
2. tic;
   / ?2 @' C1 E& \8 }7 [- c
3. s=0;4 q\" i  h: I, g2 d2 j
   
4. for x=0:1000- |: f9 J5 h( g2 N
   
5.   for y=0:100* Y* G) Z! x  E
   
6.     for z=0:100
   # x9 `; y- S2 I2 c
7.       s=s+f(x,y,z);% z7 ~# N' y- A8 b. {5 i. W
   
8.     end1 |4 N9 \0 v: E
   
9.   end2 v& h6 d( o5 k5 w9 P4 j) G/ X/ u& _7 ]
   
10. end\" d# [# e# T0 o  \
    
11. s2 F, D! F! N  j3 J, p
    
12. toc
    ; A+ Z$ o5 ]4 ^2 j2 v* x% ~9 h
13. 4 ]6 X2 e! H2 O) k+ c
    
14. s =
    ( f5 b0 {1 ]# x. {4 V% J8 l
15.     6.126720600000000e+0091 Y$ h* [; I4 V
    
16. 
    + C! V/ k8 ^# F, }
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
3 W' V' l' t* ^* B6 |, Mm函数：

1. %file xyz.m5 [; L4 }$ ]' ^/ t
   
2. function c=xyz(x,y,z)
   3 Z: t6 ?1 d3 H/ \$ r4 ]% T7 S
3. c=x+y+z;' d/ E: m; J( Y2 @0 |  g% m
   
4. end

复制代码

测试代码：

1. tic;. Q2 T! K* P0 d4 c5 ]: x
   
2. s=0;
   ) H* C2 i( s9 i& W+ v* N' G
3. for x=0:1000
   8 t: G) l4 I  r2 }
4.   for y=0:100
   ( B8 f3 b5 d. A8 x. N- C
5.     for z=0:100
   ! B/ ^8 t' z4 ^) Z2 ~
6.       s=s+xyz(x,y,z);
   : `% x1 L: }/ V$ v* `\" c6 \
7.     end, R\" Z2 e) L8 x5 P\" r, O
   
8.   end
   , E/ a' l3 A1 ?
9. end
   5 K) @5 p& U4 X% l1 l2 `
10. s
    8 W. A3 S5 `5 b+ D& D
11. toc
    0 c\" b- m+ ?7 c: D
12. 
    , n# d7 y. m+ B\" O\" c! k! f
13. s =- b2 y\" {) ]/ s9 y& A  W9 [% \
    
14.     6.126720600000000e+009
    * r* P2 R% G+ g3 g% f$ g8 f3 s
15. 
      i9 z( `, T) q5 \4 B5 z, X
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   6 N\" K3 O3 ^\" y\" d5 V
2. mvar:7 P, v2 ?/ H- T3 D
   
3. t=sys::clock(),( ]. M\" l, [# |
   
4. s=0,
   \" k1 ~: A! S- x( S
5. x=0, while{x<=1000,  m) ?5 S. |1 g/ d5 w\" R* Y
   
6.   y=0, while{y<=100,; w, w5 w/ h0 J, }6 b
   
7.     z=0, while{z<=100,, j8 M5 W# o8 H# ~3 R/ i
   
8.       s=s+f(x,y,z),
   # n: u3 }% f+ @) S2 w5 G
9.       z++' H% x% C7 S8 {# q6 j% G/ Q
   
10.     },
    2 o7 v4 r( Q& G% ?6 ]1 ?
11.     y++
    ' v9 D6 Q: c5 c* q9 l
12.   },% J4 c( b7 O% b; K  `$ j
    
13.   x++1 ]9 t) I. x$ F4 J5 Y
    
14. },
    $ D* d5 R& D, T
15. s;4 f) w* G\" d) ^8 u; [7 E' U% f  ]
    
16. [sys::clock()-t]/1000;

复制代码

结果：
; ?& ]0 G5 W* U3 j" U; A6126720600.
$ e: s* U7 m9 M0 u3 a1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
2 g' _- S8 r! N4 f0 u) ~

hxftntA

0 W% T6 R3 g% N' T0 p7 b( P

$ f5 E. o3 U# J  S  Z' q1 F$ \+ {


/ u  G6 t2 G5 \: D. |
+ h+ [$ L  w$ P! r

5 s$ I( t/ a- p+ U" m
' b, F! C: Z; P( H3 o
; @: ]3 t5 K$ _& V
1 r( y6 M+ S  c, H
$ Y- v1 j- }' B8 g* ?
- [, M1 o8 j, b7 ]+ j7 o! X1 U. L6 x