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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

( c8 P+ \4 d# X3 h4 [& eMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;$ z9 Y- G; c' G7 ^
   
2. tic;\" t5 A' N8 T8 e+ s
   
3. s=0;8 ^3 \8 ]2 V' s& K% |4 F0 U
   
4. for x=0:1000  n* u/ s\" c2 w+ H. H: v$ ^0 I
   
5.   for y=0:100
   1 ?+ P: F0 }0 P\" u
6.     for z=0:1009 _# [. d) n; X/ w( H
   
7.       s=s+f(x,y,z);# S& a0 b! o6 k* M- g5 ]
   
8.     end, @( c5 w/ `/ [$ A
   
9.   end
   1 A- c5 N8 \7 z' {- W; }4 x8 @! ]. v
10. end
    9 k% X- E; ~6 U# y( Y& g
11. s9 E  p$ f6 G, v
    
12. toc
    0 @, ^0 C& A3 Y) h0 y( D6 p
13. 
    1 Y* h) T3 {5 Q  {* U
14. s =+ j) m2 a$ s- h# G
    
15.     6.126720600000000e+009; X6 R# \$ R- y& ]
    
16. : w) v9 I  q9 }- g
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
1 P. F2 U5 |5 }+ Qm函数：

1. %file xyz.m
   ! g- n: |: W- k0 H$ [+ A. R! d  _+ t
2. function c=xyz(x,y,z)
   * n) l1 D  k- E
3. c=x+y+z;
   + E+ V3 f9 T- N
4. end

复制代码

测试代码：

1. tic;
   5 A' Y/ H5 ^- R# i9 N/ J. a  h
2. s=0;+ b7 @9 ^+ u8 }% g: o9 t7 q7 x. d5 a
   
3. for x=0:1000
   7 s+ W. \- D' s! k
4.   for y=0:100\" L, f4 R% a: j( ?  [
   
5.     for z=0:100
   - `) [8 j+ I7 @9 [( Y' t
6.       s=s+xyz(x,y,z);
   ; _: V, }( ]/ n: Q* x- z
7.     end, z8 R\" \9 u) i5 E5 ]
   
8.   end3 i7 @& I( c! V
   
9. end7 w, k! a7 _: K+ l: H
   
10. s
    0 s2 ~- t$ c7 m+ T$ u$ j3 p
11. toc
    0 Q\" t8 C/ P3 T
12. 2 Z# U8 T* }5 Q* M4 L2 g# d
    
13. s =
    $ c( U$ k& A* C$ O
14.     6.126720600000000e+009
    , s) }  p* y& q4 p
15. 
    , o: \  y4 l7 e
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

) }# T6 T% C4 nForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ) F\" N+ Q0 b2 G* ^7 I
2. mvar:4 b3 p9 }% j. a
   
3. t=sys::clock(),
   : Y\" S1 c- ?: j4 R- U
4. s=0,
   # n; F2 ?& [- `\" K* g. ]+ |
5. x=0, while{x<=1000,
   6 `+ j0 T$ }1 ?
6.   y=0, while{y<=100,\" U3 Z4 L9 C2 t, C\" R
   
7.     z=0, while{z<=100,
   + b4 w1 S5 q1 X  @
8.       s=s+f(x,y,z),
   : C( K( R& ]) R; ~& `& T
9.       z++, P, [& T& R2 F0 P2 R, W* N
   
10.     },8 [+ u  I: P' P
    
11.     y++- g$ O\" B  X/ e) d: e+ Q) U# c
    
12.   },# f6 D0 u- c( I0 R; q8 N' S
    
13.   x++
      o, I  `3 F% y. X( E7 m
14. },
    $ B! P- H0 X/ s6 Z0 f! v+ l+ p
15. s;4 v  o. m3 G9 ?0 I+ i8 p
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
2 O9 Z7 N: \' E: j9 q& U# `3 g3 l
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

; e) O2 q$ Y% ]9 t* m

4 i% u  T! L5 H3 L6 [


0 ]3 M7 E) e) {% {
( i* P+ J4 _$ O, I
1 U$ |9 i$ M  ]4 E
! G1 a5 }- U  C1 O( T2 t
* P: K2 I& ^6 g  f2 M1 V* |

/ z4 S& Q, ~# Q1 ]  C' c2 F
3 S0 ]' \' A/ _# P# [7 c4 I; K$ B6 y