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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
/ p) L; t5 k$ @. o
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   / t, m+ h5 C, c' l
2. tic;
   6 x' e# A. t9 [2 ~8 Y  y$ _6 {
3. s=0;7 [2 t7 r- c' r4 z( j5 I
   
4. for x=0:1000: \4 j' z6 v/ T$ L; Q- w, P4 k# s
   
5.   for y=0:1009 E) }0 b: A% N( ^  M6 h
   
6.     for z=0:100
   ' q5 r. k% [0 H8 w\" [1 Y  e- N! z9 ]
7.       s=s+f(x,y,z);! W7 ?4 ]\" T* R1 H: o2 r- A
   
8.     end
   0 U; c% Q& q6 W/ g* M7 z
9.   end
   / g/ l9 P* c& g3 Y5 B\" A3 W& D
10. end
    ; r* b/ [2 r( z* k; C\" C
11. s
    ' L( k/ h. D& u+ I
12. toc9 G6 M* Z( ~- M6 b# ?& L
    
13. 7 a9 W- q7 c$ f* _' G. A  d7 y
    
14. s =
    * g) F  R- m- B7 _2 P$ z6 g
15.     6.126720600000000e+009& s; E6 E# G) P7 C4 b% l
    
16. + L, @( u2 z3 Y  g
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
  [  f* O2 l/ x) X0 G7 A/ Tm函数：

1. %file xyz.m# L. C& ], ~' Z/ D) t
   
2. function c=xyz(x,y,z)
   7 a5 d4 A; Z) Q. i, a2 `$ p  }
3. c=x+y+z;$ C4 O8 _. ^5 p6 [# K3 C, C& X
   
4. end

复制代码

测试代码：

1. tic;
   ! V8 T& s* v# W# ?\" f# W
2. s=0;  v( M  C, R# j# [# F0 }
   
3. for x=0:1000+ [; A: Y! k& ~% t& A
   
4.   for y=0:1005 [5 [% a& ~% j9 X0 l
   
5.     for z=0:100
   + z5 E8 h; z$ w6 J
6.       s=s+xyz(x,y,z);! W( J3 _8 i/ }
   
7.     end* V5 I6 A  f( \& m* a
   
8.   end& O- N* n9 _. J/ ~' B3 e: d+ s8 M
   
9. end
   : G! A0 L8 d# l
10. s
    - C. c& `+ z8 y/ ]2 k
11. toc: Y: M  f5 x  i  P. X3 |# s
    
12. 5 K) ?+ w3 T) i
    
13. s =2 I; C' f! t0 m
    
14.     6.126720600000000e+0098 Q\" h; Z9 \& Z
    
15. 
    : Z\" g; ]3 E! _9 @
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
% v% |* j! T3 O: ~8 D8 c' I
. ]8 ]& s  a. [: y; S; G; ~Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
     B3 q\" c( j3 q* u4 E& u  V/ W; b8 P
2. mvar:
   8 U/ ~. F/ H9 b. A' D, C  T
3. t=sys::clock(),4 h6 E% K% g) J
   
4. s=0,# x; J7 Q* V  M* A
   
5. x=0, while{x<=1000,4 x6 h\" i0 z4 G8 `  l- S  w
   
6.   y=0, while{y<=100,- X1 Q! Q2 \. f  S\" e
   
7.     z=0, while{z<=100,
   + {7 I0 U5 J$ E' l# L
8.       s=s+f(x,y,z),
   & l6 P( z; h( Y7 f2 ^! ?2 Y- f# @
9.       z++
   3 I  |0 A' w$ I7 y' H
10.     },
    % B: U9 e7 L! k: R& \0 i# ?- T
11.     y++
    $ L6 i9 |' |8 d: X, C, \
12.   },% y' |9 |\" P8 w2 P! h
    
13.   x++. b, Y  P9 }8 n1 D7 E1 C% G/ k
    
14. },  O: |1 f' n, g9 b
    
15. s;8 ?- I' t6 Q4 L) M: a5 z\" a
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

8 x6 P2 ^/ ?9 M. m6 K9 f
8 x% u0 ?/ V( P& d
. ?! ?% G# n  `9 C
& R3 x& H, [1 A5 t
; Y! T& Q( M  H/ m
; g( s5 Q5 D  v4 z9 I1 ]. \
" U5 Z: m0 K( S- N

0 K+ p4 v; b7 X) w; G( |' s1 k
+ V7 T, J. E  E  s& s: q" _