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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
! g  f# {4 L) T! Y/ c+ o1 W% {) L
＝＝＝＝＝＝＝＝＝＝＝＝

/ n% W2 h3 g( G5 G/ K) KMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   \" v\" w3 S; }' H
2. tic;
   ! q: K/ w8 G& ]. N5 U+ j8 v
3. s=0;
   # J/ F) N\" o1 i' R' q; m4 d
4. for x=0:1000! w: o+ Q# c( _7 H
   
5.   for y=0:100
   7 l# o8 c\" h7 z( s, H9 Q7 U4 G3 B
6.     for z=0:1005 _5 Y$ v+ V  y3 x- @! w
   
7.       s=s+f(x,y,z);\" k3 w  S5 `# v' c
   
8.     end
   \" F' R! B  C9 |
9.   end
   6 u, n- @/ o& T# e! ^\" X/ J, }( L
10. end
    5 A6 ~1 Z+ ?( p; ^  W
11. s+ }2 W6 y/ x! s% S; ]+ l, b
    
12. toc: }% [* s' @0 R) P. P; l3 `
    
13.   Q. S9 W7 @' m6 @
    
14. s =
    2 ?( h1 ~' `' V\" v& V! N
15.     6.126720600000000e+009
    ( b8 j\" c' {0 f. y: l; A$ _
16. 
    . S( N. H& i; E, s! G  o
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
( i8 }9 i9 Y  Im函数：

1. %file xyz.m6 p  v# C- I; u/ y; L9 h1 G
   
2. function c=xyz(x,y,z)
   4 @8 P, M: o3 }5 W$ y( [% O
3. c=x+y+z;) P- k% ^0 \& E6 M& d7 I
   
4. end

复制代码

测试代码：

1. tic;3 C6 c, n0 J( g9 y
   
2. s=0;
   9 Y+ x' Y+ x+ D8 y\" y4 `- P+ T* r( F
3. for x=0:1000/ }- D& i- a- q' K6 m7 X/ d
   
4.   for y=0:100
   5 `; @' Q\" D! |$ S
5.     for z=0:100
   ; f7 W. q% R4 `1 ^1 H
6.       s=s+xyz(x,y,z);* U) c6 k! ^. h3 E6 K8 H* k
   
7.     end
   & \% ?. c% E2 e+ r$ }/ T
8.   end  K$ v; g4 O) d; g
   
9. end  U: i4 o2 S5 j0 S* w
   
10. s
    ' b8 g+ i% F' }5 Z% a3 U
11. toc5 }. [. r; ?\" H5 k% s; K- N
    
12. 3 U$ S# u, b( ]/ P\" n/ M3 B
    
13. s =
    2 O( x9 F' I/ P( ]$ C; j& k
14.     6.126720600000000e+009% d+ V1 D1 ], Y
    
15. . M6 i9 m- |1 U4 {0 O
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
# w! b8 g6 a  S) b  L$ {
- X9 `$ Y- b, p5 V8 T. YForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   . Q. n; ^* S' h6 W6 o
2. mvar:
   2 Z( f2 C6 z# l! B5 F) i
3. t=sys::clock(),
   , ?1 l! Y9 K9 ]3 s4 C
4. s=0,
   7 c# t8 D; `\" ?9 D7 b/ i; B
5. x=0, while{x<=1000,
   $ _4 ^; `; u8 i1 K2 E8 k- G- `0 X$ `( x
6.   y=0, while{y<=100,
   & p8 Y1 d6 p+ H$ L# I
7.     z=0, while{z<=100,. `4 I0 u5 F7 g* v# z# K
   
8.       s=s+f(x,y,z),, `' V2 @6 g8 T6 F$ k# N* u$ C
   
9.       z++
   8 \' i7 k# f8 v3 ^- K
10.     },
    \" U  W' y2 M1 F& K# Q8 G& v; J+ ]
11.     y++4 i( r* X$ ~# ]( {7 @  C
    
12.   },
    & L5 [: Z2 K. ^8 L
13.   x++
    $ u* z! ~; z; A* }! c; V# |
14. },
    ! t) b8 Y/ H$ ]
15. s;
    + ^1 W% {' e. r# E
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
$ r7 O7 a: ^/ \" F1.109  秒
' s- N# |" o! O0 K& e" j3 N+ W
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
! H- T0 _6 u  b1 Z; [0 |

hxftntA

2 p, B9 m; s- P* K( G
; r4 `0 Y2 }' M: ~3 p7 e- J6 {




% a0 _# t9 R; p8 j4 N5 b; Y- K3 X
  ~# \3 S# j: Z0 a1 L% g4 z2 u/ Z9 m* f

: Y6 g' E* b7 A' z9 d
. g* e  Y2 K1 X' ]
7 }$ v! w+ K: c* G" R



' V( c, A) f3 c2 i) e) I! w
* O7 }5 r8 i/ V" {6 d