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

forcal

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

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

) K# |3 x1 a1 c  g6 C( p7 G: X! P. tMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ) r6 D/ e4 t& K4 n. k, M/ R: F
2. tic;
   3 C& E9 a, [# f+ I. z
3. s=0;# B' u9 h9 X& k0 |2 k# w
   
4. for x=0:1000
   / Z$ W6 x0 |5 {
5.   for y=0:100( \5 C- J* Y' F* |  S8 m+ |
   
6.     for z=0:100+ R  z3 D9 H. \\" u) r8 q
   
7.       s=s+f(x,y,z);
   ; y# C! {8 s  e, |) J\" X
8.     end
   , g5 W% ~0 x3 i9 M7 {
9.   end; T\" _2 A  H. H' N\" F' s
   
10. end2 X* S3 v# K: l- X
    
11. s
    * ]$ H, k\" I! F, `
12. toc+ M+ @+ c# x  L7 f7 v/ l/ Y
    
13. 
    6 A\" t! k2 Z6 W  m
14. s =\" b& ]& p7 ^- r# f
    
15.     6.126720600000000e+009! r9 y# q1 E; h) R+ }$ v
    
16. - V- N& S\" h! k: Q$ @; K( i
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
1 P' a# e2 j0 W$ u# sm函数：

1. %file xyz.m& s  v. d6 i& O' H
   
2. function c=xyz(x,y,z)3 W2 I& S: g7 s5 b- E4 n/ K
   
3. c=x+y+z;6 h5 @' y9 q9 b3 p# G, D0 w
   
4. end

复制代码

测试代码：

1. tic;* ]5 v/ t0 Y: |) g$ r
   
2. s=0;
     Z8 t\" G3 k1 D9 E9 D8 r3 S: s
3. for x=0:1000
   + r3 g7 @5 x5 f6 Y
4.   for y=0:100! _8 L0 m* o2 ^5 h. T/ M9 h
   
5.     for z=0:100
   , _8 F; _6 @$ I0 S$ y, l\" Z% Q
6.       s=s+xyz(x,y,z);% p( U3 e+ m8 r: I0 q% `9 D4 o: J
   
7.     end
   1 A- E6 ^( M+ l
8.   end
   * R8 k& j% G+ ~+ G2 i
9. end8 n3 P+ [8 q- Q8 y
   
10. s
    3 n\" u3 d/ E& c4 Q0 n- {
11. toc
    1 }# ^% t, |$ _+ ~
12. 0 P1 Y1 n; p0 O1 v' d$ P
    
13. s =# V+ [1 I$ D' y' U3 M. O% e- \) T* G
    
14.     6.126720600000000e+009
    ) i6 l5 T$ m# n: m# E; T
15. 
    , p8 s: Z0 y$ X+ g% G
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

% P  _% n2 Q* h8 i/ Y0 IForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   8 E9 r) r  X2 e, Q2 c0 }, l0 C$ ^: L
2. mvar:* s& I0 ]& Q% }) J4 k0 z9 \
   
3. t=sys::clock(),
   # P% W& Z  K8 M2 |
4. s=0,
   1 @9 s: q: i+ ?
5. x=0, while{x<=1000,$ W' q7 l, k1 Z1 l* h, d
   
6.   y=0, while{y<=100,. I! M4 ]0 `$ v
   
7.     z=0, while{z<=100,
   9 G7 e, `: }$ w+ g! N' l
8.       s=s+f(x,y,z),! d% Q8 p' V# a4 p
   
9.       z++0 ~- Z  m# W1 \. P+ Q. _: c  h
   
10.     },. d4 j4 J- o1 C( k. H& k: v
    
11.     y++( j\" |! t5 X( ~$ e6 H# T
    
12.   },
    # i% J- N; k; i/ k
13.   x++1 Y8 \( V& K( O9 \% G# d/ H
    
14. },
    : C: U\" t& p+ K
15. s;7 T1 h7 g$ _1 o. _' P8 c
    
16. [sys::clock()-t]/1000;

复制代码

结果：
" \* w7 R/ d; L7 C6126720600.
1.109  秒
1 `6 N0 \. l3 ]5 ?4 r( _: R
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
/ j( z! G& V/ O$ x6 h' N/ L

hxftntA

* W" c9 V* P/ B+ V

0 z+ j/ H1 M+ Q& E2 f4 c$ ~




/ N: h0 \' W6 [" [. @7 M
" u1 o7 o. Q. Q" t' L1 ~

* M& h# ~+ h' p



; S, w' f* T6 ]6 F! `$ ^' q