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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
2 p: i% O( e+ }) `3 y
＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;. G% E8 z* p0 b0 Z
   
2. tic;
   \" R. j) K! a/ ?  q9 K. u& p
3. s=0;
   4 V% Q. B# g5 S5 R
4. for x=0:1000
   ( r, k+ T! }\" G: b$ i1 B9 p
5.   for y=0:1001 \- }5 ]8 f, d- N
   
6.     for z=0:100, G( ~; m/ N+ {- `
   
7.       s=s+f(x,y,z);' G9 u# N9 d& B
   
8.     end
   2 k\" I3 J5 G1 d- {$ L3 t
9.   end. W  _* h. q2 S6 ~5 ^, T5 m
   
10. end
    3 w. q/ T# s% d. D
11. s
    - l' R) B\" D# s# F& n3 z
12. toc) X7 ]5 j) b# k' {
    
13. 9 T  Y& |3 S* Q8 w: ]2 e7 ?1 r
    
14. s =
    ; V9 B7 |* X1 h; I- a# w- e& p
15.     6.126720600000000e+0090 y7 M1 I* U0 }& k& N, ~: i2 }
    
16. \" w; Q( \! j+ h7 i1 I& d; F9 k6 s
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   0 l* v; x% [7 p+ C( i% P3 Q* ?
2. function c=xyz(x,y,z)
   & O& D$ D\" u( P- f( J
3. c=x+y+z;, G9 K* `4 w$ P+ h: j7 [
   
4. end

复制代码

测试代码：

1. tic;
   2 q; j7 q5 o+ y
2. s=0;1 J' i1 E' T$ a+ T
   
3. for x=0:10001 Y4 \8 T% O( v/ ~2 P; N) p
   
4.   for y=0:100
   \" G4 W0 e* a5 u2 R' c) `. |
5.     for z=0:100
   - F\" ^. D* O0 ^# r7 ~
6.       s=s+xyz(x,y,z);
   3 f' M, R' e5 B4 ]' \0 m
7.     end; @& _- _  T4 m7 u$ |7 ^
   
8.   end
   2 r: {' g8 k& C. U7 N/ ~5 z0 V# @
9. end
   ! Z4 J, P% z2 V\" s1 b
10. s8 T\" v3 ?8 J5 A
    
11. toc
    + F) A' c3 o; d+ N& U
12. 
    : O\" \+ t\" ~1 w6 M& Q# q+ d; P
13. s =
    $ Q7 }! @8 E  t2 T: G\" O8 S
14.     6.126720600000000e+009
    \" |& d4 C$ u! l' ~& T3 e! C
15. 
    - ?7 b/ S8 b1 Z( `8 [
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

# ]% g* j+ \1 C# a2 VForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   + Y: {; G0 }! X\" v- }4 g
2. mvar:
   9 P- ]! N8 k; D& O3 J
3. t=sys::clock(),! _6 r+ y8 z% Y$ A; Y5 u2 ^
   
4. s=0,
   6 |& P+ l8 r: |5 z
5. x=0, while{x<=1000,
   0 [* I- r. I' k
6.   y=0, while{y<=100,
   5 p* {+ g! z% D$ S
7.     z=0, while{z<=100,
   2 X  v- R# {2 Q
8.       s=s+f(x,y,z),\" {5 n6 O: `$ c; j) i
   
9.       z++
   % ]3 ]\" E0 D( W4 g+ s
10.     },
    1 V8 I- r5 _+ i
11.     y++2 b) U4 X) U# u! [- k; e
    
12.   },
    8 `3 ]! M$ f: M. }; l6 Y
13.   x++
    \" c/ |- g: m' b  z+ a) c5 k& n+ x
14. },
    1 M2 ?2 N$ Q* r
15. s;\" R- L0 |\" \\" v! M9 \
    
16. [sys::clock()-t]/1000;

复制代码

结果：
7 `# ?* ]& |% p! y/ ^$ m+ R6126720600.
1.109  秒
' X( V  V6 x1 e# R
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
2 j  Y, F8 D3 H9 e+ p- j) J5 b8 e

hxftntA

( W, G" x% T; ~" n  J

* A1 v8 g1 ]# ~" y2 ^


& z3 }7 o# g7 S3 X





% m1 Y' @8 i! a4 G7 s/ L$ b# G: A
; B$ q1 g( W- |: u; [, y1 k3 O0 B
' ^5 M1 z5 v( `+ b