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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
- N' p( I- b- W- J! O
$ R' f9 u; ^' Q! \, `0 `& Z7 o; h＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ! [+ E; S4 P6 a5 {2 D! b4 u
2. tic;
   6 G\" ?6 S, L3 H4 {+ D! p; R/ z
3. s=0;) o, }0 }  c1 ~& t. ^$ W  t
   
4. for x=0:1000
   4 o: E8 D, D$ v& F4 X
5.   for y=0:1009 F1 ]  j' A# U: T
   
6.     for z=0:100& w9 H% d- V2 ~4 [* I1 X1 w
   
7.       s=s+f(x,y,z);
   0 N+ k3 Q: u' e2 j4 k, |3 R* T
8.     end5 B- D! @2 `; T2 M( t; V
   
9.   end+ `3 L9 o& g6 v) ?
   
10. end: ]! m% }# `- U. h! F
    
11. s
    * S* r% X* ~5 d1 C9 Q! `# w
12. toc
    5 U5 J) ^2 A+ D. G9 O0 h
13. 
    5 S+ ]$ {$ S& O6 t' S# q
14. s =
      z* P2 b! ~6 k# c\" ~5 U# d' I
15.     6.126720600000000e+009) X- s2 R  }9 G# {4 b+ {' R\" F
    
16. 5 `7 Z6 S8 _3 L9 v\" E
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m/ \# H+ D( a$ j& c0 z
   
2. function c=xyz(x,y,z)) I, r- g8 X, h
   
3. c=x+y+z;% O- F\" D6 ]5 s# A3 F
   
4. end

复制代码

测试代码：

1. tic;
   6 `8 V' W3 W\" e) m' t, O\" X
2. s=0;
   % O9 U8 [- F2 R$ h+ u
3. for x=0:1000
   3 G. q, M- z+ \- x
4.   for y=0:100
   ! S9 T4 E  E+ P# ^
5.     for z=0:100
   - }2 W$ ]! c( A8 S- g7 y
6.       s=s+xyz(x,y,z);
   6 e& J7 m9 G1 z- q5 w: o
7.     end: b* m6 l( I9 E% v) A, H0 b
   
8.   end
   ' e\" j  l& @0 d
9. end- P$ W5 e. C0 F% p( h8 o
   
10. s\" ]2 [+ X$ q) z' k9 Z9 }; |
    
11. toc
    % R+ M9 ?0 N6 O5 `
12. 1 J3 J& [; \0 x9 f: w
    
13. s =
    . I5 H$ J0 s  J+ |
14.     6.126720600000000e+009  E, F: e; `8 K  J$ }( M
    
15. % r4 u! P8 g; \! X& @! i8 J
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;! O, a$ t4 ~7 [2 a6 m' D, U9 l/ C
   
2. mvar:9 ]! [1 @6 }: i% Z. b' Y  i\" u! N
   
3. t=sys::clock(),
   ; v7 p! J' Y* d* @% g8 c+ O: @1 `3 f
4. s=0,- r8 I- q: r\" u9 _4 b  V
   
5. x=0, while{x<=1000,. l: K1 X\" k% M( J
   
6.   y=0, while{y<=100,4 `) g& B) j: S3 j- x
   
7.     z=0, while{z<=100,
   & g' `9 C, }9 ]2 l8 e- N8 v$ K
8.       s=s+f(x,y,z),
   \" Z+ N; R* R! D6 p/ z8 o7 V
9.       z++. ]6 D* o6 k  z. O! Z3 V, {+ B, J. O3 x
   
10.     },) J( \& Y# W; L; N2 G
    
11.     y++2 d& r. b3 K/ Q# I! K0 Q! u+ W' n$ Y
    
12.   },
    0 ?6 q0 f5 m  a
13.   x++6 t$ v7 f% ?$ F3 |
    
14. },
    - Z  D& G8 n7 d+ |
15. s;- g; l! r6 ?7 O
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
6 F# r) R% c+ F6 n- ]; C! u# |; ^3 M
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

% U! L! J" B* }9 k7 p
! d& n7 ^. ~& r


: w* v4 P5 K8 u2 }3 r* u



2 h; O# m0 z+ m, ?- w

$ Z1 T7 o3 \" v6 L! u& B5 p
) {5 \( f5 s8 h
2 N' T4 O; i0 |: h, x; P8 e
. P9 m; V) S$ y

1 m) w  B3 Y. t' z' N' ~
/ u0 c: _8 {, L: h  o