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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
4 k2 s, G  u* [7 B7 b: u% d
/ {1 Z0 L9 l+ g+ i- x- `＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;\" i# J- t+ g5 F* K, R. C! Y
   
2. tic;
   \" |) u4 T+ P/ m( e
3. s=0;
   . `6 ]  L\" u! ~: @: z% j
4. for x=0:1000
   2 v0 ?7 M% `# U, S% @6 q) O
5.   for y=0:100- Z: T. n\" m- z. Y
   
6.     for z=0:100
   % g6 {$ e: t; G- h/ ~1 i0 D2 Y; n
7.       s=s+f(x,y,z);
   ( x- n5 q8 I0 i  e# ]! }
8.     end
   / h' q6 b  _5 M- R
9.   end( j, t0 m5 \3 {& B$ S
   
10. end2 h0 ^: v8 \( O& A1 @
    
11. s) s/ M0 _8 r' _  a3 l\" W
    
12. toc3 Z; Z& [& J$ B: F
    
13. 5 z' o; O% F  Y, S' I
    
14. s =
    - V1 Y% v1 q% Q0 `
15.     6.126720600000000e+009
    ( ]' ^& M) E$ t6 C2 J
16. $ a7 @8 S3 v3 ]! @' {4 t
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m% `4 n% e\" p& Y6 ?
   
2. function c=xyz(x,y,z)
   2 }% E+ A; u7 N# i/ y& l4 {
3. c=x+y+z;' n4 o/ X9 R2 N& L
   
4. end

复制代码

测试代码：

1. tic;
   0 X4 w- k2 r7 {! b! c; p0 r
2. s=0;\" {7 W0 v\" R$ {( ^& S! A6 r
   
3. for x=0:10003 b* N4 l, m: x( _. l# g\" \
   
4.   for y=0:1008 `  k. `# |* k; J5 O  M: Q* r4 g
   
5.     for z=0:100$ s1 \% [1 s( b4 J\" j# G7 E
   
6.       s=s+xyz(x,y,z);
   8 n. a0 b: k# M+ ]/ s! X
7.     end# n# S4 \' N2 `6 p
   
8.   end+ e9 J  \- }1 B: k4 {) v
   
9. end: `1 H# E6 s) j' \- E' F/ v
   
10. s) U3 [$ F! \! l
    
11. toc) R: \+ o6 \( M% C. O! m- G
    
12. 
    ; }* O  q2 M* M% D& G: z8 y
13. s =* k1 v0 X. [\" s
    
14.     6.126720600000000e+009& q/ b- j+ u+ y% F& J
    
15. 
    9 T0 e. m; E- a; M
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
, G) U' B6 A5 F/ W& w% I' U1 @
0 c4 y  q3 P% cForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;0 z7 P& K# Z- a\" j1 O8 v  Q
   
2. mvar:
   * m% J* ~7 _4 N6 n5 ^0 U$ t8 D
3. t=sys::clock(),1 p; r+ y! n3 ^+ \& {, k. o# M
   
4. s=0,
   6 q9 L3 S\" g! ]3 G
5. x=0, while{x<=1000,
   5 a1 s* B+ \; z5 M
6.   y=0, while{y<=100,& {( G  r, Z8 r, d# M2 P% A
   
7.     z=0, while{z<=100,
   & `- {* y4 Z  N9 c/ G9 V
8.       s=s+f(x,y,z),
   1 h# D2 l8 h0 O5 ]
9.       z++% J$ h0 a3 {6 I$ x2 z) G0 x1 s
   
10.     },. U9 V0 \9 A8 l* ~3 Q9 q
    
11.     y++
    : q# k/ }; u7 ?% b( P
12.   },2 g; x& q& H3 h  X7 G7 c5 z1 A+ b
    
13.   x++
    ; M+ l# a. X6 a* }8 L9 ^
14. },
      N( ^7 s6 S9 `  p6 O
15. s;
    - X) R# r5 |$ C: B# i; R$ r
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
8 [/ F% q- B  b4 {' T5 L$ x" M1.109  秒

% @3 `, b; s( \& [二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
: R1 Q  H3 D" c. q3 o; e- B& t* U

hxftntA

5 F+ F- v, `. Q: n  b

$ I- A! u5 I) Q" ^7 r/ Z& T' G
; y% Q/ x. }" p
! g& j: A; H4 F; O' \; i

' ^. y" I8 T6 }1 u! F5 Y6 `
" r/ u+ j; ^$ X$ ~; B






2 u3 [% _% f4 [' G" U