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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
  E7 E  e6 y% B' G2 h" b
＝＝＝＝＝＝＝＝＝＝＝＝

) N" z$ y* X' F4 ^* HMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ' ?% I4 }, v9 [% `\" W2 `
2. tic;- \4 P5 L3 a; \0 T6 r
   
3. s=0;) g* p3 t- z+ k5 a% T- m
   
4. for x=0:1000* @9 e2 j, `8 |' S1 y
   
5.   for y=0:100; L$ W% G7 }\" t' Y
   
6.     for z=0:100
   7 D' T4 J2 S3 b% ^- J+ Q$ W
7.       s=s+f(x,y,z);
   # y9 l( r) a3 u  q' `+ R3 f2 c
8.     end
   0 p# ^- y; P5 I- x( E
9.   end\" V# i2 F) G! n- L( @  V
   
10. end
    ) g- x& B  U# L4 E( `
11. s
    9 n0 w! U9 k9 Z5 ^
12. toc\" f/ [0 U1 g! b, X- u$ [: {
    
13. % N, U1 b. t- X: E5 d9 T
    
14. s =* F4 U6 Z( ]; h\" \( U, \2 n6 F0 N
    
15.     6.126720600000000e+009
    , e% w+ N8 c3 J) I% W, Y$ B4 \
16. 
    2 P0 s9 j5 q+ p$ v& M/ T
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   : ]\" b+ s3 n; \% q% D3 ?) g
2. function c=xyz(x,y,z)
   7 \/ s! Y5 v# K/ {4 b- a$ t
3. c=x+y+z;! s; d5 y3 m2 k7 Q7 O
   
4. end

复制代码

测试代码：

1. tic;
   ' g; c+ t& b' ~* t$ ~9 H2 q
2. s=0;$ B* h. z; I/ o) {- w
   
3. for x=0:1000( ]+ {; E- L+ v; I4 A8 \
   
4.   for y=0:100$ E. b  H2 X2 D* }) B  |
   
5.     for z=0:100
   . u6 R! b  D# x\" \
6.       s=s+xyz(x,y,z);
   & _# h\" x7 H7 k8 \+ y/ v
7.     end0 c) b. ]% S+ b! ?0 k3 f
   
8.   end  Q% z  [3 j6 p7 K1 g
   
9. end
     R+ y\" X& o; R( \5 ?& [
10. s+ R, [* g3 l( k; W7 @
    
11. toc+ w6 D% I' x8 ]$ G( l1 C5 {
    
12. 
    ) \5 `, i/ L' Y' ?
13. s =
    : n: P\" W& s. c. ^5 R9 G
14.     6.126720600000000e+009
    3 T+ Y* e5 C5 }8 f& R0 u0 P$ Q\" b
15.   e7 {2 g, b% A# q
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

$ i- J# P5 e9 K  _2 T' |Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;8 Z8 H$ f: c3 G\" G3 Y$ i
   
2. mvar:
   - a, s8 y8 w& t: U
3. t=sys::clock(),
   9 k7 t6 W! \  \; l8 `
4. s=0,! M) k, S- y7 N
   
5. x=0, while{x<=1000,* u/ M) T5 i) N! C0 E% ]6 D
   
6.   y=0, while{y<=100,\" ]! M1 c( K$ c) |: `9 }
   
7.     z=0, while{z<=100,3 ^5 W# x\" {/ L0 x, }( p: ]* h  m$ J
   
8.       s=s+f(x,y,z),
   3 E7 f! d& \+ J! d+ O$ d& u3 Q\" E
9.       z++\" J0 i5 ^8 f1 E$ s# h7 Y  F* c+ Y0 B
   
10.     },
    8 z. \9 r# ~& h' q
11.     y++) _6 J, U& `8 b4 ~! W2 M
    
12.   },6 M$ q- |0 y  ?' A; A0 u, j; V
    
13.   x++
    7 u/ v0 _4 J) v+ ]
14. },
    \" S; f# h( X3 v1 v$ A% v3 h- F
15. s;
    , D4 u8 @1 @1 m- V7 X: E8 K
16. [sys::clock()-t]/1000;

复制代码

结果：
; n+ a6 o% n4 r. N* [6126720600.
: \: J  r0 Q6 l* @) j  Y9 [1.109  秒
4 A1 M9 I/ }! |" j: D* Q. l
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

) U2 f' Y, d! B7 t
$ B( o. g: W# j8 f( y: {
8 d8 j1 Y) B9 `
- f$ }* p1 O) H% b9 F

5 a* c" {+ G# @! d9 D& `7 m
" b) ~9 q) l% t4 U3 V4 ]# g

3 p) b: S, o% t2 t5 m6 S
! n( e2 I+ `% t: ^  N# V! J( C
' Q4 O- |  `3 p5 n# B
' {: f' s1 s' ~& ^  Y8 Q" }
9 Z. _$ `8 m9 J( n5 E( z9 D4 r
5 q+ I- v. n* T8 `$ F
1 m" d5 b2 V7 G2 ]