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

forcal

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

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

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;) J8 G! \% F' h6 y% f
   
2. tic;* F0 Z0 g# `% r  O* a
   
3. s=0;1 m6 ~) Q) `( r2 n% F
   
4. for x=0:1000
   * O# P; d- A! l% B2 d
5.   for y=0:100% [  t  }  H  m, T3 Y
   
6.     for z=0:100\" y2 @5 L% R/ Q3 P
   
7.       s=s+f(x,y,z);; l9 Q7 @; Y\" e1 H  s' P: A) k
   
8.     end
   9 k\" t! K  H% _
9.   end/ n5 `5 G' X3 C$ }6 {3 h0 H6 _# Z
   
10. end
    , v0 E3 h3 U' h
11. s
    ' G& P# y7 W\" E: Q
12. toc8 W- y' \. n# Z6 J( n
    
13. 5 E1 C$ v* B. H/ B+ g' x2 g
    
14. s =; ^) ?\" Q7 W2 z8 d) T) @' @
    
15.     6.126720600000000e+009
    # J( `\" L7 p0 g& w' N! l
16. 
    : c  {% U4 L& [
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
/ a9 A; k" J% T8 Z# i3 }m函数：

1. %file xyz.m
   8 _! e+ W\" w+ O/ \* K, n
2. function c=xyz(x,y,z)
   4 S1 E: |* }6 p5 d, T
3. c=x+y+z;
   & ~! E8 q3 c8 T0 H9 Y
4. end

复制代码

测试代码：

1. tic;
   ( X7 h9 M  Q1 }& c3 U8 O5 ]
2. s=0;
   # n% E2 k  n- S: i! I/ T- e
3. for x=0:1000
   8 b9 J' m( i0 l/ N5 C3 Q1 \
4.   for y=0:1003 I! V# R) ]5 {, C  N* B) [
   
5.     for z=0:100! s& n# y2 X( g; y
   
6.       s=s+xyz(x,y,z);
   : j& K3 q9 t4 P
7.     end! }. {* |) Z% S4 q3 q3 f
   
8.   end) L, u& C/ d) w$ M  H% s
   
9. end+ Y5 A, t$ _/ u5 G( f
   
10. s# t/ L( R6 Q7 a- Z) N9 m
    
11. toc; i$ g1 W5 {# d  Y: k2 [
    
12. 
    $ c8 b& V; E  \/ E2 Q, S4 I* c+ k; X
13. s =
    $ m! k3 V/ Q# |! `& p4 Z# Q, C
14.     6.126720600000000e+009- j* i& I4 w/ u2 |
    
15. 
    3 |, c* {6 q% w
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
4 a: h) H4 c8 M; d" v; y
$ a/ `  \9 d. V/ [9 f- P9 n  jForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   6 ~1 P8 H9 U' P( x  K! K) m* u
2. mvar:
   $ ~& ^1 h9 R+ B6 P- P
3. t=sys::clock(),
   ( |* v9 a6 a7 s# I
4. s=0,
   ; J' Y! J6 F& C- G2 A7 z; @$ {
5. x=0, while{x<=1000,
   - F\" N2 y' `+ q( q5 P
6.   y=0, while{y<=100,
   : I# H* i  G! ^
7.     z=0, while{z<=100,\" v9 v+ v0 k& y2 I
   
8.       s=s+f(x,y,z),
   9 n* m$ @' ~) d
9.       z++; m9 l# u6 C; f( d! u7 f
   
10.     },6 ^8 o; `# F9 V7 q% ]
    
11.     y++2 v\" j3 M: P9 a; e! Y\" U
    
12.   },. I8 e6 y& h4 F( M3 n, \  ?
    
13.   x++
    0 c! p# s7 ~* R\" l4 z3 w8 H9 X9 r
14. },5 ^. n  j8 \* W
    
15. s;9 L3 Y4 B2 }( D9 ^3 ^3 A; O
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
. U* {+ C' f, a6 ~# c4 G0 D2 z1.109  秒

' X9 N2 R3 R  _9 ^二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

/ C/ L7 `' _. z$ p/ h



3 G# x( M8 g# K. ^  C
  D, F/ C8 }) ^! D: P7 f
' d- C/ q. d9 b
3 @7 @  e) C$ I" G. O
# R7 p7 q  H& S# O4 n: H! K; B$ U
' d3 o9 w4 r6 ~7 T6 ?! ]