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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
5 q: D& j7 r; J* g; |
# j7 {3 r8 B2 X: H$ I& s9 t, @＝＝＝＝＝＝＝＝＝＝＝＝

  p* a- d+ P# M5 nMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   + C: X) F. v+ S
2. tic;
   ! R7 K\" s5 B5 E$ L; U& v( l1 o
3. s=0;: L\" G9 M: e) o
   
4. for x=0:1000+ x+ r/ k\" S0 v\" x+ A& z0 s\" A
   
5.   for y=0:100. C: t\" Y\" R% f. w\" b
   
6.     for z=0:100\" V  w5 v& F/ @* N3 L, H  q
   
7.       s=s+f(x,y,z);- v+ V+ F0 `\" c, @* C7 |
   
8.     end
   % `; D6 z1 I0 }) ~5 \3 s6 ]( g
9.   end
   . S* A% C) s1 [' s
10. end
    ' ^. u- P2 A) g9 U
11. s7 v; L; j9 v' x\" l
    
12. toc
    5 W# C1 h/ D% v, D
13. ) x- t+ H8 o+ v$ U
    
14. s =
    ) c2 V: X7 p5 v* N6 i! o
15.     6.126720600000000e+009) Q  t/ b! a! P1 o1 D# ?3 Q
    
16. 
    ( {. }: b, m+ e, C9 F
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   . @2 u6 |& ^) F: h% G
2. function c=xyz(x,y,z)
   / p' W. X' P& b2 p+ H9 Y
3. c=x+y+z;$ @\" \: P3 G! C* A2 w
   
4. end

复制代码

测试代码：

1. tic;. u* F. ^6 b9 _9 h. C- o  g! |0 I
   
2. s=0;
   3 u: C! A! ^! z+ [7 \, R
3. for x=0:1000  X4 q+ R6 e3 I* {
   
4.   for y=0:100\" [( ?: f  i: E. k
   
5.     for z=0:100- \% }1 S' C$ Y
   
6.       s=s+xyz(x,y,z);
   ; s* S- `! L0 i3 B; p# B
7.     end
   0 o% C4 D7 t1 d. u. R
8.   end1 S& C7 @1 S7 s3 B4 y7 H. L& h: l
   
9. end
   % {4 x  V/ B7 [9 N
10. s; ]% v1 d8 y( O\" @5 ]0 S! F
    
11. toc
    / E\" S4 e0 J) e\" @
12. 
    5 _* K) ?; E# d+ e% l: Z$ y0 a  S
13. s =
    3 T1 x9 M9 H: w% b2 a, T! y
14.     6.126720600000000e+009
    2 ]& @- {% O9 }\" _+ ]9 y
15. 
    $ i9 S$ B+ z9 e& U
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
7 U+ M: i: G& x  e
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   + P, O\" ^8 \* u\" e2 v
2. mvar:
   , c. T4 W; V0 J6 K- R
3. t=sys::clock(),- n. w6 y! L* M7 P& C1 `
   
4. s=0,
   6 Z9 {0 K5 ^2 e8 i4 g. q9 T0 P, q
5. x=0, while{x<=1000,. ^6 e8 e; g4 r: E8 |% o. l) e2 ]
   
6.   y=0, while{y<=100,\" k: E! @7 k, f5 ]) x- a
   
7.     z=0, while{z<=100,
   4 k8 L# |$ ?# ~; ]3 b, {
8.       s=s+f(x,y,z),, y9 P! m+ i8 e+ \) Y. M
   
9.       z+++ X) [( k1 A: U* G* W
   
10.     },
    - F/ V* K! u; w: ?: F
11.     y++
    , e; o/ C: k5 k$ E
12.   },0 f* {1 @  i% [) d\" t\" H/ U
    
13.   x++
    , y/ ^1 f1 |8 y\" L3 v4 l$ {
14. },
    & X\" R' M+ g( L/ f
15. s;+ R  {; ]( [8 A
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
# e' \+ f/ w- ]* n1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
0 A4 r$ B; B2 ?( U# a! i

hxftntA

0 s5 V0 p0 ?! u6 o3 l6 ]2 R6 F
1 y" ]: A! V6 z
9 Q$ V2 k1 ~9 |" `# l
) I) I* k* G, [- w1 {2 ^

9 E. s# l& E/ V: W" W1 a

8 N: O  E: I0 M! H/ P4 ]% G2 o
+ a2 G3 g0 y  r- x
9 A9 c- r3 s5 v# m7 y4 p5 R4 O+ H. {
* e& C! }4 d+ u, B; M
8 D) e+ m% [9 j



$ R: G. z% D& d" ]6 `" A
5 X% x7 `+ ^+ b; E" I, V" {# K