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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
2 b: V7 y* {  w3 R' m' l1 e3 E
＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ' W; j, J8 r5 j
2. tic;
     C6 i\" E! m0 H0 E1 t% y9 G
3. s=0;: ~% I# Z0 `5 X( W% S
   
4. for x=0:1000
   7 ?# T/ z* F# r7 Y: g
5.   for y=0:100. o+ J5 ]' f7 U, A  e2 N! a5 A
   
6.     for z=0:1008 ?* n) G* k' Y
   
7.       s=s+f(x,y,z);' I$ Q2 Z  V, N& _$ j
   
8.     end3 M0 F+ i8 V. R  k: _
   
9.   end4 \9 `, I2 S3 w' ^
   
10. end! R1 I6 M! U! ?; b\" v& r7 \
    
11. s
    6 }: N% r! Q8 v  x- `
12. toc$ Z7 {, a# S. r0 L( @
    
13. 
    \" E$ {2 C% h6 o
14. s =
    0 \  [2 d\" U0 T0 K2 A
15.     6.126720600000000e+0093 R8 ^- @, m6 J  A4 |
    
16.   Z* \7 A8 s* t& }; e: S# ]1 Z+ i
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
# p; m* x+ E/ _  v2 Gm函数：

1. %file xyz.m
     T# o( ^, Z\" O) X
2. function c=xyz(x,y,z)
   4 D* C% @$ i% b9 A4 A4 p7 S2 i
3. c=x+y+z;2 r% `/ ?6 I& n\" _- F( k$ l5 O
   
4. end

复制代码

测试代码：

1. tic;\" V( H+ K& `8 [& U( e
   
2. s=0;
   + G7 z1 F. d9 `* V1 p5 Z
3. for x=0:10007 V; G3 `/ a; D. C8 [3 g- f1 C
   
4.   for y=0:100
   5 ~  Y) u& T( z  U/ ~
5.     for z=0:100\" O0 V' i1 P; m# P' ]
   
6.       s=s+xyz(x,y,z);
   / Z2 |5 p! v' J: S
7.     end% }7 C, l, v5 @# M
   
8.   end) ~% N3 |( X) O; D/ I1 a
   
9. end
   * h/ x% H  j- ~8 ^& G0 r
10. s+ |  Y$ P1 R. n
    
11. toc7 A* }& a/ e, n9 Z9 S1 o
    
12. 
    & N! Q: {! x' u3 r1 R$ D+ d
13. s =3 A$ V4 }; ]: o7 [4 ]
    
14.     6.126720600000000e+009
    5 Y\" E; P9 U7 U- h$ K4 i
15. + v  H$ l( M+ |, m
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;8 b( j. v7 {: T
   
2. mvar:
   0 [5 y! n% k; }6 @
3. t=sys::clock(),
   / ^. g+ _( L2 m. @
4. s=0,9 [& B% W& L- C$ X5 Z# t& q
   
5. x=0, while{x<=1000,4 f\" W& G. W# u\" Q
   
6.   y=0, while{y<=100,
   ' Q# M\" X5 S, i2 M
7.     z=0, while{z<=100,
   9 F# J; M2 {\" R# o
8.       s=s+f(x,y,z),6 V+ w. Y' x8 O) W) \; W
   
9.       z++; L1 f1 y\" }3 _$ R/ M' f; ~
   
10.     },7 ]) {' P2 ]2 A. Y! N
    
11.     y++
    % q8 d+ E7 e6 g' B! p2 S; R& z
12.   },
    \" c) B8 O) E. F: n9 B
13.   x++
    - J& D  z* D/ [) S% U
14. },
    0 n: g, [$ \, I\" s8 T7 Z
15. s;
    , k% [\" K, F% p- H# U* l5 S
16. [sys::clock()-t]/1000;

复制代码

结果：
9 `! @" p) k2 ?6 B, S) Q6126720600.
1 J! D4 L8 Z4 }7 z) h$ d5 T1.109  秒
" ?0 W6 w: f4 M( X
% @$ @/ z4 u" ~2 c* ?2 b二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
. R8 L! e! S' H( N2 Q( J' X/ C9 z+ w

hxftntA

* ^/ o+ n' N2 ]9 X+ u: n

& x, F5 T% P7 z3 J
+ }7 D) ?; f6 |2 f; S1 u( e$ i
! V3 |* R  ^- {9 u4 Y% W


7 V1 f7 h8 j' L5 o
4 y# |  _! }. u, \% N

$ S1 K* J1 t+ B3 z6 \
* c% G* P) X; a" T- b6 T
  L* _5 p0 s* o% Y* A

, x( }( _4 R0 N* a9 Y; o