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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
# e) k+ x' q( u' F  @6 E
' Y# }! N' z* I5 ~0 R＝＝＝＝＝＝＝＝＝＝＝＝
5 X! U- h, z! S8 J! x: c1 B
9 ~- U. t3 h. GMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   8 f( Z+ D! N5 K! L7 z. ~
2. tic;4 B\" l; @9 T% [5 C$ r. `  J
   
3. s=0;3 t/ |; D! Z2 Z6 C, B- z2 ^6 G
   
4. for x=0:1000$ p/ |3 ^! L7 g
   
5.   for y=0:100- l% I8 P1 j+ I$ r\" K. D
   
6.     for z=0:100
     S0 ~* u' q) ?
7.       s=s+f(x,y,z);3 a; J$ Z  P  ]9 r+ R: L
   
8.     end& R1 c& F- B9 k( U0 a3 Y2 f; u' B: t( y
   
9.   end! E5 m1 w4 ?, `7 H\" Y: h
   
10. end4 |6 s/ p& G0 p9 w! s
    
11. s+ K# l. N5 K2 e( p8 o
    
12. toc
    # h  Y; D  N6 `- ]
13. 
    # s  s( w8 ]7 e
14. s =\" q0 ?$ [\" v4 r8 a: `0 c
    
15.     6.126720600000000e+009
    0 r1 |2 [; k. B$ e& N/ K
16. 
    : v/ E% [. E( {\" s/ K( I3 i
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
& T4 I) |1 l1 R! ~m函数：

1. %file xyz.m\" d# J5 a6 }; [% ]) J
   
2. function c=xyz(x,y,z)
   & o\" v; s! e7 @. p* D
3. c=x+y+z;! y) C3 C; t5 S3 `4 M  D0 K* y* k
   
4. end

复制代码

测试代码：

1. tic;  x- A1 ?% D. M% I; T- j
   
2. s=0;
   - F- P/ z) r* ^+ k# K0 p
3. for x=0:1000\" b& O6 c, }! s. f! _\" L6 m: v
   
4.   for y=0:1009 B6 l# r) k* d0 C' {: E  A
   
5.     for z=0:1000 g- I  M# r7 G, z  W2 B, w6 d
   
6.       s=s+xyz(x,y,z);
   8 W; J8 h, z' O- f( W3 j4 K
7.     end6 ?# j: A5 Y\" j% x2 J1 m1 ^2 ]5 ~
   
8.   end
   % @\" C. l  m+ A' J- Q) K* N' H2 D; B
9. end
   - Q) U8 }$ B) p\" s6 O
10. s  N# A* R7 i/ u\" J; }4 S
    
11. toc9 S  y* W) @3 e
    
12. 
    ( p) M; B7 q% B5 C
13. s =. d5 C: q, o\" y4 y( a# T* j2 ^7 L
    
14.     6.126720600000000e+009
    # F& W+ q5 ?0 G, G
15. . _! W4 A. T7 n\" B6 ~
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
8 w! j7 {( ^) G2 R, Y& b
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   & f0 H: }0 k& f( T( c' i) w: _( z
2. mvar:
   0 Z& E7 E7 f9 }5 Y
3. t=sys::clock(),9 R4 q& _! [6 c0 e- I2 J
   
4. s=0,- j/ Y# q0 ]! q
   
5. x=0, while{x<=1000,
   6 t6 E6 i- ?9 d4 ]9 h
6.   y=0, while{y<=100,
   1 g) R* Y: l! h8 q/ T3 K
7.     z=0, while{z<=100,( X& K0 }  A2 n+ }( l: m
   
8.       s=s+f(x,y,z),
   9 z, v  o% C) w  c4 k0 y2 g
9.       z++
   $ X* ?3 C  j! j1 p4 [
10.     },
    9 t' u1 {/ c7 c% O% P- J+ m2 ]
11.     y++
    - l( Q9 R' e5 |  I/ ^  g8 {' U
12.   },4 r! D- r' G) g, R; S
    
13.   x++1 l  D9 t8 w) j. o: A% X7 L\" e
    
14. },
    ; m* q\" R- m7 P* V6 U& B
15. s;/ O, u+ P( X$ z( @' ~2 D
    
16. [sys::clock()-t]/1000;

复制代码

结果：
, u7 W# ?; j3 k% a* k6126720600.
* E. ~8 Z3 S% s6 q! a1.109  秒

" T) f$ M1 f7 Z& [二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
6 H" J: Y6 V' E/ C8 g( Z

hxftntA

  v% h( f$ q" N' R
- a% n8 R9 u  t  J" L# J" O


" b" K4 E8 b# a0 h5 x' D+ c% ^








  m! o5 v) X) m9 `

" R& m5 W4 K8 Q9 Y