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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
* ^/ S' _6 b0 S" a
2 Y7 k6 R$ B- r: ]4 ~Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;6 h, x+ b# [# r# e
   
2. tic;
   ; F8 n+ Y/ c4 ^  T1 Q
3. s=0;
   ! U3 _7 V0 e% a% q; H$ Q+ V. T- Z
4. for x=0:1000
   # K/ T7 X  V, d2 Y
5.   for y=0:100
   : `/ Z4 S( U, Z4 B6 ?% x
6.     for z=0:100
   , ]+ X% x9 j8 M
7.       s=s+f(x,y,z);& `. d; f  W, B# W# ~' V4 a5 H, H/ A
   
8.     end* j9 u9 v- O3 i; ?4 @- q# [
   
9.   end
   : z! a# p: _5 N- \1 G\" v
10. end
    * @7 X# {& l. x& i
11. s' g2 _  k9 n+ U- q. _3 ^
    
12. toc: Y1 B4 K1 c1 e0 e
    
13. 
    * G8 k7 }' j2 y% g* w
14. s =
    - u( Y' n% i6 j% E, z8 }
15.     6.126720600000000e+009
    / m/ N  q7 \3 J* Y' ^1 I; a1 I& N
16. ) q% J+ y) G9 T+ K% `
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
* i3 B5 _; x" [8 j7 u3 y! wm函数：

1. %file xyz.m4 v( I! }  ^+ t$ ~2 y$ _
   
2. function c=xyz(x,y,z)
   ' n7 F: `; h/ L$ v
3. c=x+y+z;: E0 D; g. H! n1 S) l$ e
   
4. end

复制代码

测试代码：

1. tic;
   + j$ h% _\" s' K6 r6 z4 o( V, x% p3 t( t# j
2. s=0;
   : M) z$ K0 h! B  o2 W$ Z2 Z
3. for x=0:1000
   4 X/ d' G) X, v, g5 j  R
4.   for y=0:100
   0 m$ q5 M, c, b8 ~4 `# _& r! k
5.     for z=0:100* e4 K% h9 Z* W& X, i
   
6.       s=s+xyz(x,y,z);
   2 I2 L5 D& C, |5 H% F
7.     end
   4 b4 e- j3 {- W- ~* j4 Q. h6 C
8.   end4 [7 ~$ [8 i0 O; |
   
9. end' k$ f! U5 Q) x. h$ U: |! i- u
   
10. s\" ?3 S9 g% l( w\" W6 K2 ^6 a. [1 V% o
    
11. toc
    8 n1 t/ Q) `5 x  N% V1 x\" d
12. 
    * F, K2 c. E! h6 z0 w, O6 ]& L1 v6 b
13. s =; D/ h& S4 k8 d  t( V' Y. ~
    
14.     6.126720600000000e+009
    9 b/ W' C) G# ~$ _
15. . _# p8 V# i7 u! c/ \- o& ~
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

$ t- J, I$ l" @) [6 C/ k- ^Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;; P, d$ W( y: e
   
2. mvar:, r* U. G  }7 J. ]9 B9 G0 g  }
   
3. t=sys::clock(),! S: Q! C4 @) y\" O
   
4. s=0,& g5 p- O8 ^( t# ^! o
   
5. x=0, while{x<=1000,7 \8 R# l' p5 P* k  x
   
6.   y=0, while{y<=100,6 Q9 w- p+ n; ?2 n5 t/ T6 s
   
7.     z=0, while{z<=100,
   : h+ K# t; f, N% M7 c! I# G\" A/ L
8.       s=s+f(x,y,z),
   7 l3 r: Z( f) e$ g; I- m1 W\" s
9.       z++9 Q3 s\" S) ~8 v\" |: c
   
10.     },5 g* |; i. H9 f$ E\" W% ~
    
11.     y++
    & u+ @, }  ^. U4 A
12.   },\" k' c3 m3 q\" |* W! m7 f# c
    
13.   x++! P/ z4 [1 I. a) |: R# }8 F  g
    
14. },# J& W( p* E. y* x- b' a5 L
    
15. s;
    0 E+ q& v: x$ M9 l4 g4 |# I
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
* a4 Z# L" n- t3 C/ e2 J9 t: A0 c- E1.109  秒
6 `% t; L0 H, c$ C7 I  \: J
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
$ Y+ ^( X) d& ~* G4 c8 j

hxftntA

$ p0 P; G1 U% _




) o3 p  `8 z4 S- v
, S' G! H3 n$ X/ m7 }
2 I! o; h- `5 U" k1 y

. E3 p& z3 N1 }/ e3 \' Y0 K6 o$ y  D

" k+ \3 L! t& O, i5 ~7 T
0 u; U4 M) j/ ?$ d. R# @