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

forcal

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

/ s& T0 M, c; j9 I8 U7 [＝＝＝＝＝＝＝＝＝＝＝＝
8 L( K( Q1 X; m2 ?! \( n
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;$ e- m8 A6 K  p1 ]
   
2. tic;& i. @! r% E$ r% `
   
3. s=0;% o- I* P, Y9 K3 s3 o
   
4. for x=0:10006 }, j6 V. Y' s1 y
   
5.   for y=0:100, Q1 m- B- j. Z3 |
   
6.     for z=0:1005 j- X5 C% a4 B9 ?. {4 y& k
   
7.       s=s+f(x,y,z);
   6 e  @* i. r6 M3 }7 d
8.     end, M- c) V- X5 @9 [
   
9.   end! I* B- `% z+ {! _2 G& i
   
10. end3 r8 h$ |* ]3 Z\" W# ]- u+ y
    
11. s
    , L  x0 b1 k8 w, b# r0 g
12. toc
    7 L) h6 ]8 g( t' E$ o# Y3 H& r
13. 
    5 E3 W, P\" d# K; f$ K
14. s =
    9 Q7 N\" z, `$ j0 `0 X0 U9 f
15.     6.126720600000000e+009
    4 I1 G1 M% a) @/ \6 P
16. 
    ' J! J5 J, J6 K
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
4 x, V- n& t3 P& S- N, mm函数：

1. %file xyz.m
   9 v# w( a& W4 I9 k9 N% h! f' x
2. function c=xyz(x,y,z)
   - C$ @. a# f# V: u0 `- L4 D\" B! P1 b
3. c=x+y+z;+ v) t+ s0 T. o# w
   
4. end

复制代码

测试代码：

1. tic;# ^7 f% W0 a+ O% `' O4 c
   
2. s=0;2 H: ^( C, i$ u# \. N, U
   
3. for x=0:1000
   9 p) I$ l! Z! R4 y+ a; y; Q% t: j
4.   for y=0:100% c% v% u\" T8 s4 x2 f
   
5.     for z=0:1007 i; }8 l  l: D: O4 {+ A9 {, w. w/ Z
   
6.       s=s+xyz(x,y,z);1 j' N5 Z. i3 a2 t# c
   
7.     end
   & H8 x* v5 u# P7 g1 a+ |
8.   end3 c' Z\" T: C7 W$ p, d
   
9. end: ]; Z% b: ]% n( p8 o' Q
   
10. s
    : K9 D* M\" X3 C' _. \+ P
11. toc
    + I! _$ w( Z, D# i6 J% z+ u1 x- o; i( [
12. & _9 v! O4 V; \% y$ G4 F% v; _
    
13. s =+ o& b2 o/ e& Y  ~# A/ a
    
14.     6.126720600000000e+009: l6 c% @; Y5 y) E' w
    
15. 
    * p9 N0 }. _  s
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

  t. F" ]! B5 \2 t# YForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;3 ?' F& r2 t- M, Y
   
2. mvar:' R# v+ q$ K- @
   
3. t=sys::clock(),
   1 `- G  |( ^' [. c. d; ?
4. s=0,
     Q2 ^- Q2 }. }2 Q: K
5. x=0, while{x<=1000,
   * ~* V- S0 R! F# |) a
6.   y=0, while{y<=100,
   % D8 T' ~* k; o: Y' q; f) r, `
7.     z=0, while{z<=100,; ~4 p9 l+ \% g! C
   
8.       s=s+f(x,y,z),/ D0 b6 P/ `# A3 l  c' w: Q
   
9.       z++
   4 y; D$ H7 j6 i% M8 X$ Y\" E. k# T
10.     },+ r4 K( J0 L, C# m8 ?' n2 `
    
11.     y++
    ( Y% `; I3 U4 K% h9 P! u) ^0 y' m
12.   },7 l* I0 e2 U  C* b2 W
    
13.   x++
    ( a# ]2 E* W7 W9 T  _5 q
14. },, [6 H  O, `' ]) Q& }
    
15. s;
    ) f, |$ I3 O) S5 A) W
16. [sys::clock()-t]/1000;

复制代码

结果：
% y" E$ z! Q7 O; G6126720600.
1.109  秒
$ x1 x& Y; _/ D  e
9 k! l7 }* x# t; ?8 P二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

- T" C) n2 J/ b2 R1 l8 z3 @
+ G9 i. H$ Y& Y- }1 m





( u! K! g' N0 f5 i, w

  [5 e* G. [/ |* k
  o+ r: n- B. h& f6 X7 l8 o; \




1 A2 h' V3 x9 E1 \# {5 k% Y