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

forcal

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

( I6 ?! O$ Y/ n＝＝＝＝＝＝＝＝＝＝＝＝

" h- [, b, c8 J3 R, YMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   \" f: I\" t8 F/ _) l+ E
2. tic;
   & h6 R- N6 |; n2 B
3. s=0;
   # K5 S, v, ^# g. W: ~: ^  D
4. for x=0:10007 O! I9 R* {* K; E, s4 R( X
   
5.   for y=0:1005 A: V% y, Q) ^# D; T: E. }
   
6.     for z=0:100
   1 z+ K& [: o( W5 H* C  ?
7.       s=s+f(x,y,z);
   # |7 F+ O+ a$ ~- c/ X8 T7 @
8.     end
   7 R- ]7 m2 f/ u3 n' O7 `% X2 v
9.   end\" q# w& k$ y\" X) O( s\" v/ ?
   
10. end
    5 X# I4 J9 g4 y
11. s
    # |% C2 q! F& `
12. toc
    & _& U& S  H\" i! j$ M! h) |. B
13. & v3 g9 y( o# Q- P' N
    
14. s =
    . P% F' x) M9 c0 j6 l
15.     6.126720600000000e+009  V% k8 n/ }& N- \( N( c
    
16. 
    . B6 X0 o  ?% q- z
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
m函数：

1. %file xyz.m$ R+ x5 ?2 l\" B# ^
   
2. function c=xyz(x,y,z)
   1 J8 y9 @& }# I
3. c=x+y+z;
   + v3 Q. j: }, s1 ^. C8 _5 k
4. end

复制代码

测试代码：

1. tic;
   1 K& A3 ^( T- v3 U
2. s=0;! C! O% l4 \  h4 b
   
3. for x=0:1000
   ; V4 B7 [3 l: I* v. s9 f- Q
4.   for y=0:100
   % R: u4 D5 ~; x% S) k
5.     for z=0:100
   7 I1 @, @+ V( B, [. L# L
6.       s=s+xyz(x,y,z);
   8 D9 L+ ~( ~& ~/ E5 Q% I
7.     end1 \8 Q) S\" @, v7 F0 w( d' }
   
8.   end- H7 T; s* e; A0 `  c
   
9. end! d# t7 M\" _& X0 P6 K0 h! l
   
10. s1 J; t\" D- Y6 i: I5 Y2 c5 P
    
11. toc
    + ?1 v\" R+ G) @- C, @3 C, a
12. 
    # Q$ L7 K4 a. @. ^2 u7 V- x
13. s =6 V\" ~. S' l' `
    
14.     6.126720600000000e+009  N3 e2 ^. P# s! T
    
15. ' z9 E, s! L( `' g4 _
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

2 e9 B) T" y$ @/ F- ?5 J1 EForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;2 _, M) Q3 P6 P$ n6 V1 O4 ?9 H  B
   
2. mvar:5 s5 H. Y1 h( o
   
3. t=sys::clock(),/ n! p1 z* e/ S/ z) G& D
   
4. s=0,4 ?' t) J- N6 {3 v( s4 `
   
5. x=0, while{x<=1000,\" K9 Q( ~4 w' C7 P
   
6.   y=0, while{y<=100,
   5 o5 U! Z7 e) W+ c7 H
7.     z=0, while{z<=100,' D, C% Y! o$ D; ]$ o7 Z
   
8.       s=s+f(x,y,z),0 M0 ^! ?* x: h( C: m7 u/ e! l
   
9.       z++
   & c7 e9 m: R; ?2 V) I3 n9 @
10.     },
    ( _1 T) w0 K2 e( p% l  K* q
11.     y++
    3 d# s5 s0 S* F9 j, x
12.   },
    7 s: b% k. J$ i' M8 {
13.   x++
    * z6 e, D8 E+ Y  q
14. },9 X( r/ R2 E7 Z9 ]1 e
    
15. s;
    : v* K+ J  I6 ?/ `! Z, `: C
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
9 ?) \# Z8 V' r# k+ ^. e1.109  秒
( J& h4 U' f9 ~9 t3 z
3 D, J" U6 z( I4 R% U" j. L二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

2 g% i. r' i7 h, p3 S% Y! s: Z
) P, R5 j- O& ~  P

# b! U# ?. ]+ s" u7 m
  A/ A  Z' g" H( X% a



5 U' C2 H) i  f8 f; F# G



9 a, x7 I/ h0 R% ?1 t4 ^
; M. {& x! g4 w" K# {+ A% f

* w* ~# q) f# P. T
" u- |/ S! p8 O4 K* v9 E2 Z