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

forcal

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

, ^! I0 ^+ U/ H＝＝＝＝＝＝＝＝＝＝＝＝

: ]# p/ ~+ K& R+ X) @* V3 q& ?Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;, W) j! u6 c. _9 F5 b
   
2. tic;\" u4 e0 o9 m/ P3 m1 f* X* r9 F
   
3. s=0;
   + ?; ?5 J- R3 z& P. r
4. for x=0:1000# ]) N' M% \0 Y; |# N
   
5.   for y=0:100
   . C) t* }) `- M4 P
6.     for z=0:100) x1 L) \8 t4 z# H, K* H\" \
   
7.       s=s+f(x,y,z);( J) v4 e3 }! b\" e9 C; Z
   
8.     end
   ( X  u6 J  U2 |) C- w
9.   end+ H6 Z& i5 q0 \$ M- i
   
10. end& w9 [7 i; w- y9 {' n
    
11. s\" t2 u0 p\" h9 p* a
    
12. toc
    9 v( O4 M# W4 ?8 S! }5 f- M9 @
13. 
    7 `( A2 c' {, U7 \; P0 s/ f
14. s =2 l  U! H0 ]2 P
    
15.     6.126720600000000e+009
    5 |1 }5 r) c6 t3 X* N
16. 
    * B, H- B; J3 `; g0 I- w3 ]
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m% X1 t# Y/ X) [; B) }. ]% A6 }
   
2. function c=xyz(x,y,z)
   ; }2 a  s  S( i8 ]/ u
3. c=x+y+z;
   0 ^! D; ?2 L2 Q\" l3 Q
4. end

复制代码

测试代码：

1. tic;+ P! @# f% J: J7 W
   
2. s=0;
   / u- f3 X- R  l  e' b& D) \
3. for x=0:1000
   + \4 T9 u- M' y) r$ E# w; m
4.   for y=0:100( c+ {, E3 @3 y' o
   
5.     for z=0:100
   , H2 h. h6 j/ d5 O, ^( `2 v% c( ^
6.       s=s+xyz(x,y,z);3 n! [\" [7 D6 }8 e4 Q
   
7.     end
   5 b1 M\" D! I2 d; `2 \\" C
8.   end9 g\" Y7 ?! L7 g7 n
   
9. end
   8 T8 P$ U  C7 t\" s6 G3 N: u
10. s# Y! G& y% D  O2 ~' E' q
    
11. toc9 l) U! I9 ~7 }. U& K6 V1 s
    
12. 
    ' _  }9 F- q; O: @% J; O9 W. o. s
13. s =2 |' c8 I1 S  y0 a4 n% u/ V0 q' ^
    
14.     6.126720600000000e+0093 U/ k8 K* J' p4 L
    
15. 
    6 u3 e& V\" Q- K8 {$ w
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

$ Q+ K  O0 A  dForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;* ~. e4 i8 b3 a2 W3 Z- @) s/ @- a& Y
   
2. mvar:) o' n3 g2 i9 |
   
3. t=sys::clock(),5 q$ o* E2 o: H( x7 G1 g7 X
   
4. s=0,) u8 F1 v* S' f
   
5. x=0, while{x<=1000,
   5 d1 r\" ~! K$ L6 T! X/ B
6.   y=0, while{y<=100,
   5 z/ w' `; n. ?! \! g! s
7.     z=0, while{z<=100,  A: `* K3 w/ q, w) }; t
   
8.       s=s+f(x,y,z),9 c1 O/ X7 x: s1 f7 w6 Z
   
9.       z++
   $ L6 f3 K  T- {' C/ T+ ^7 m4 S2 a) h
10.     },\" X6 I- |% @: c3 W( ?/ w. R! m
    
11.     y++
    \" z; ]- T\" F7 U
12.   },- n4 C2 n  f! D# t- U! f3 I$ y\" F
    
13.   x++8 T( p7 Z, I) [
    
14. },* D, ~, A+ Y/ c# s2 D3 b  [# G. x) g
    
15. s;9 w3 y% |/ g2 h! H
    
16. [sys::clock()-t]/1000;

复制代码

结果：
, Y5 A0 {7 X0 d2 _1 ^6126720600.
1.109  秒
- T4 p. T6 W9 r- X
* ]; z+ ]1 L9 A0 K# ]8 X二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
6 F8 Q7 ^& l4 [$ B& z0 V# @

hxftntA

, Y6 N2 ]' w, K) z! `3 f
' F% Z3 ~- V7 N: [* n! |
' z4 U) n3 h$ p


- V- U- R4 n( {4 V+ I- q


3 z2 D, B7 m1 F0 Y% a
$ |' }! o8 a- N# m' [& ^


2 a9 U# F6 ^9 H( I7 S
9 `/ e3 r+ s% `+ o

  n. @( u. d3 ?+ N