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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
2 D# N' E( ]0 u" Q9 _7 C
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;/ b: `( H) _& F2 m) M0 |' n/ d
   
2. tic;2 u$ U1 i7 z9 n
   
3. s=0;5 z0 [, K# f0 y- Y% d% S% H
   
4. for x=0:1000
   # z( {4 h( \6 o; L: V( j
5.   for y=0:1008 e: b1 b$ l1 C7 Q7 i- a
   
6.     for z=0:1004 Y. _$ V) ^( F$ y% h1 @' P
   
7.       s=s+f(x,y,z);
   7 p4 o5 f: j, |( ^6 ?7 E
8.     end! p: ?4 ^; x( m# {2 }3 n( @
   
9.   end
   1 i+ o# n2 g8 U  _/ h3 p
10. end
    - t! o: t! p! x8 e
11. s
    * p% A$ ]. W& f9 [9 J+ Z: @
12. toc( N$ v; n4 P, X; U+ {
    
13. - _( N% J/ f8 s) J- s3 D- @. K& @
    
14. s =/ u4 H, t) Q# Z# H0 r
    
15.     6.126720600000000e+009# U( Z! L( V8 z# c0 c, t* R* S
    
16. 
    ) O; u, r4 r4 A( S* U
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
( C- O: g! G8 r9 n- A$ Km函数：

1. %file xyz.m, M, ]  `6 s1 z3 f# h/ m7 a
   
2. function c=xyz(x,y,z)
   8 r) Y$ w- l/ z) L. j+ \
3. c=x+y+z;
   / X% x  r) W$ a+ A/ o
4. end

复制代码

测试代码：

1. tic;3 m6 x7 s! d  @: Z) Z3 D
   
2. s=0;
   \" n9 e. U8 s/ A
3. for x=0:1000
   ; f: G% u: N3 V; M! U
4.   for y=0:100
   & b& @\" Z! ~) T4 F9 D
5.     for z=0:100- ~; k5 q0 r4 }3 O; w
   
6.       s=s+xyz(x,y,z);2 c/ L+ ?8 q% d
   
7.     end
   7 W: w3 O/ T/ u4 \: k+ ~
8.   end
   1 e$ ]% q2 {, T& ^5 y/ e- n7 i3 \
9. end
   - a\" r# M4 ^# @  c3 v; T; Z( p% |
10. s\" B& E- s# C( d( ~9 ^; q( f
    
11. toc
    ' Z% Z4 f  Y4 z6 M! @% s, z  O! ]2 L0 g
12. 
    # v) V8 U/ K% {7 O! d) P
13. s =. ^5 K8 v4 {& C0 @, |7 f\" @) R
    
14.     6.126720600000000e+009
    * o' E$ S' I$ P8 r
15.   _7 Z4 V8 k/ w, O( a0 u0 z5 G
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;1 h4 C; W4 m\" ^( P
   
2. mvar:
   - U: X  J5 K* M) ]% t: X! ]+ K3 K
3. t=sys::clock(),
   - L* l1 {. H\" q5 s
4. s=0,
   . O4 w/ v8 G4 V: h$ c2 Z
5. x=0, while{x<=1000,
   % d* D3 r  p1 X  a) D\" |
6.   y=0, while{y<=100,8 f( \8 i& k3 \& D7 v$ e8 T- c  g+ g
   
7.     z=0, while{z<=100,9 B. L5 [2 P$ }( Z0 u4 ^. e1 P
   
8.       s=s+f(x,y,z),2 `, o, @* O3 ^( l: z. F
   
9.       z++' a  h/ G9 E$ s
   
10.     },5 H2 H9 `& V9 }' G; ?8 [* z
    
11.     y+++ E2 _8 O4 r7 ?
    
12.   },1 p2 @/ A& ^0 I1 @9 c0 c5 S
    
13.   x++
    ! a3 B. J9 z; d' r7 c+ k
14. },% ~( Z: F- V+ ^. B
    
15. s;
    + ?' O6 s0 h; _. ?% z, K
16. [sys::clock()-t]/1000;

复制代码

结果：
/ r8 a9 |3 o5 o: b6126720600.
4 S) u. u$ r6 [$ T1.109  秒

( U! Z" v2 J8 C0 |* q  B二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
- Q. w* q6 {8 V4 X  F3 ^- a( U

hxftntA

! r" J$ L8 X: k! C! a

* [! k+ C1 p; S* Q8 I' F/ a




4 C7 d# F% x7 R  ]( l+ O0 N, W



. s/ y4 {& y7 r  \0 b5 [) K
! d) j- P7 ~2 d
; a- K% A9 Q; T
0 j" ], |+ K; \% X, k! u% X