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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
3 j2 [/ [: U. v. f( ]# f' P
＝＝＝＝＝＝＝＝＝＝＝＝
2 L" _4 R1 U: [
: p2 m' i4 K$ `/ f6 p* k. sMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   5 L- D; q5 U5 L4 l* }/ n* c
2. tic;! l, \  d7 W; ]\" F+ m
   
3. s=0;: W  J$ \% Y% j4 C
   
4. for x=0:10003 I( Z# ~; c7 M1 Y2 r$ r
   
5.   for y=0:100% s/ w, j+ X/ L4 y
   
6.     for z=0:1004 \0 s3 q% t  P
   
7.       s=s+f(x,y,z);% ]3 U/ i. \' K( N: n( s
   
8.     end
   1 A3 S9 c! n+ h0 P1 ^
9.   end. V1 S7 R# A5 g8 E
   
10. end\" U% u. T7 \9 x, @7 R
    
11. s
    , ~& r; k' w6 u5 D( m' Y
12. toc1 R( A; e8 T$ M4 T\" r7 U5 R
    
13. 
    ; u\" j; g$ U, i. p
14. s =7 P2 N( P  X4 k\" G8 S% P\" [\" t
    
15.     6.126720600000000e+0097 ~. c( I% ^4 u
    
16. 9 L  @# V\" \! t! \
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m9 J# g, p& o2 \/ b\" |! s# {3 K
   
2. function c=xyz(x,y,z)
   1 o- b3 w& r, n/ R, P
3. c=x+y+z;
   8 ^4 W% i% s3 u$ E) s7 |4 O
4. end

复制代码

测试代码：

1. tic;* R5 t& D, n3 c! f2 Q
   
2. s=0;* O3 A2 j* L8 r5 b. A
   
3. for x=0:1000
   . {, K- [% r' U4 l2 {
4.   for y=0:100
   8 K8 c! Z6 M3 B' |8 v# W
5.     for z=0:100
   4 y5 V; p2 U+ u\" m3 H6 w! L
6.       s=s+xyz(x,y,z);1 D0 T* q6 L3 h3 f7 t4 J, n
   
7.     end
     X* k6 T9 n% i+ r: Q8 y
8.   end
   ! X3 V! j; ?/ E! p7 T
9. end
   0 z7 A1 f/ J0 U* u9 p) C
10. s4 r- q\" A* @1 U- q4 S1 {6 N
    
11. toc
    4 Y- H  z: j( h# `4 A# g* Y( h
12. 
    9 C4 [' a: s1 U0 w4 n
13. s =
    . z' i) k$ b$ K
14.     6.126720600000000e+0096 {% _& i  Q1 `3 V% r+ O
    
15. 
    \" K8 Z\" t) p+ G+ Y6 p
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
9 n$ T  a* d3 H/ ?9 g' G6 G0 V
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;\" K3 o0 R, U2 \
   
2. mvar:$ T! o. [6 d  c\" _
   
3. t=sys::clock(),
   9 M( O. ^; K. G3 @' I( ?( W
4. s=0,
   6 G1 B2 m9 G. T( D# N2 o
5. x=0, while{x<=1000,
   9 y- q2 a  I' r# z7 q+ I  r
6.   y=0, while{y<=100,
   + Q; D4 s2 ]6 b* z, F1 G9 g! b# Q0 w
7.     z=0, while{z<=100,
   # n: v2 @' [4 E. }, y* j$ u
8.       s=s+f(x,y,z),$ Y8 S1 K- `6 u( Q- Z
   
9.       z++
   9 U' t( ^  w- z0 O. n/ I
10.     },* n+ V* V+ Q! q# J
    
11.     y++
    * N( Q# C/ t6 b% S/ P
12.   },
    8 i# V& ^# r& @. d
13.   x++0 ~) {5 D; a# j1 ^7 r, S
    
14. },% o$ y/ h/ B6 }
    
15. s;. i8 U4 f8 s# ~6 u6 ]7 w- {3 z
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

- R1 T7 }1 F6 z# F$ |8 {



- ^: L, @4 B! ?/ [6 W. A

* q# T5 ^' T. i# b0 E

# f; \* k+ z8 |2 U5 |1 I
* r8 }* l2 u( s: v: d


7 a" u0 E6 y1 A1 A: S! F8 C


# i# _, j+ l5 o5 V$ ]* m. H