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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
' H# E+ P- K* k& p) o
9 d# z7 f' s; g＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   $ h& `! d6 k* Q- a6 C8 W1 D/ X9 {
2. tic;
   ! ^) z' N: F' g. K8 e% m8 ]& x8 P0 V
3. s=0;+ s4 u7 t4 F( L) U
   
4. for x=0:10006 f8 B& C3 V2 y+ V% m: Z
   
5.   for y=0:100
   - `5 u) `( P\" e/ u
6.     for z=0:100# s- U  t\" z  P  y
   
7.       s=s+f(x,y,z);2 B7 w6 c! n; ]* K' O
   
8.     end
   ; Q8 y\" v# ?2 m* ~5 Y
9.   end
   7 s9 c. Z- z9 _5 Y% \
10. end
    9 |1 P: g8 z( c1 b4 y
11. s
    ) R9 M1 L6 q  j2 \9 a
12. toc, ?% Z) l- {* }( {, U- P\" v
    
13. 
    & k# t% W! |9 r9 |; s' J) a/ Q
14. s =
    / \; i: K5 R7 v; p
15.     6.126720600000000e+009\" u2 y' B8 B3 _0 h% w
    
16. : D) g1 s; W4 c* Y  V5 j
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   & z9 n6 o- J3 _
2. function c=xyz(x,y,z)7 p1 q3 c7 A7 U5 y9 t
   
3. c=x+y+z;
   ) T8 M- {; R; B
4. end

复制代码

测试代码：

1. tic;
   1 i) t' Z4 h3 K. v
2. s=0;& o1 p  S' R) H- D2 y5 H
   
3. for x=0:1000
   ; J: h$ U1 k3 r/ A, L: A
4.   for y=0:100
   . c2 }% ]4 a: q4 \* Y
5.     for z=0:100- M/ j4 e\" v5 Z4 M4 n5 C
   
6.       s=s+xyz(x,y,z);\" G- Z- x  r5 a/ u6 h7 y1 {
   
7.     end
   # v$ y# m9 k5 @% T0 i
8.   end
   - V' F: m8 _- a: O9 a6 X
9. end% q9 G6 q1 P1 r$ |3 X
   
10. s, c: m* Q* e( `2 d4 ~: ?4 u1 D& t
    
11. toc+ j) J( ]3 @* y+ c
    
12. 
    3 i5 C# N4 O6 K$ V, [- c% x
13. s =
    3 N- S! x9 D: Y8 O7 f6 n* E2 x\" v
14.     6.126720600000000e+009
    5 M* q+ A: w' Z1 p  m
15. 
    4 M) E* L& G0 S& t8 G) H
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

0 Z0 W% m; W" ?1 bForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   & A. I7 U$ i2 u* j$ @7 z
2. mvar:8 ]; h) _, S& R0 O. n
   
3. t=sys::clock(),
   / v9 f! l+ b) O8 s5 F
4. s=0,
   ( B- c% I' P( L$ i4 B
5. x=0, while{x<=1000,
   5 \) h$ X9 Q% b0 R' A
6.   y=0, while{y<=100,; i: U  t4 s* q% U' C% k. y
   
7.     z=0, while{z<=100,# B2 A8 O3 g\" Q- l1 J
   
8.       s=s+f(x,y,z),
   9 `, x3 ^( K\" X0 I% }
9.       z++\" Q7 n' ?$ j5 a
   
10.     },* }! }5 R% K1 W4 ]% P  q
    
11.     y++, d6 D' n) f4 r2 w
    
12.   },% O% b2 L! C; \1 N2 t
    
13.   x++! {1 q% w3 I% N7 E) q
    
14. },* t; i6 R: H! D- B# `
    
15. s;
    ' N* N1 ?0 t. D8 w, e- A: ?
16. [sys::clock()-t]/1000;

复制代码

结果：
' L4 Q# B/ s! e+ S$ f% U6126720600.
1.109  秒
* T, n- t+ l6 s7 R
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

$ S* D2 V6 u" _6 l


6 L3 L- ^  Q6 @$ g/ Q5 B& x3 B0 z










0 n5 R: c) y3 {  p( F

) g/ ]& r6 V5 w7 l4 W