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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
$ U6 m: M9 q; q' ?" h5 ^  R$ z
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   % D- g( B! O4 t' [7 Y0 T
2. tic;
   0 J4 `6 c9 t0 K4 ?- M$ x
3. s=0;6 f8 t0 A6 G' i1 z& x
   
4. for x=0:1000; e$ w) O- E& z$ w+ V1 J. {2 ^$ V
   
5.   for y=0:100; u' ^8 f& ~. e7 B, n' Q# R) k
   
6.     for z=0:100
   * o  a) @+ ?! Z1 o6 s$ M
7.       s=s+f(x,y,z);
   ! _+ d. ?8 o+ {. e
8.     end! i9 n2 P; _8 O( ]( R
   
9.   end6 ?; T, k5 `8 Q
   
10. end
    8 u# y* P! m) G. v
11. s4 Q2 j\" O( @7 Z
    
12. toc4 f. T( l: A/ ^\" r8 O\" X
    
13. / E7 n* r6 [! p/ d
    
14. s =
    5 q$ Q1 H9 g4 ]1 n' B. w
15.     6.126720600000000e+0092 K: F# r5 d8 d, J6 c$ B
    
16. 9 L, l4 |9 A8 ~2 `1 r  T
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
( q& S9 K% }: g8 @0 F& V9 y; Bm函数：

1. %file xyz.m
   1 a9 J. @% t5 {8 n+ s
2. function c=xyz(x,y,z)! f% T7 U* k\" }6 O$ x+ z
   
3. c=x+y+z;2 |. d& n8 a; e\" c8 F1 S, |+ ^
   
4. end

复制代码

测试代码：

1. tic;' U& ]6 S* x) }2 ]
   
2. s=0;- f' s. m# O( z1 s4 ]
   
3. for x=0:10009 Z& o& `9 R0 |/ W$ e5 T
   
4.   for y=0:100
   1 L1 [+ {8 W, b
5.     for z=0:1007 k- Z0 t; O; Z; |# H. ^
   
6.       s=s+xyz(x,y,z);- O; Q6 y9 y5 b3 `2 l% `) G$ j
   
7.     end% a\" i1 g- J, B, Q$ Q; E3 n
   
8.   end9 q2 ?, D4 S# C, z* }
   
9. end
   7 X4 c  ~\" g- E9 Z% B! \- r5 h
10. s. h. y, C: H7 z( j% v$ |
    
11. toc
    . x& d% |; ?9 J+ e
12. 
    . n, e1 v\" q$ t+ T
13. s =
    7 n1 a! \  G$ G9 g. p
14.     6.126720600000000e+009$ i2 M+ c# |) |% Y9 C' k5 p
    
15. 
    \" [& f0 C: `  J2 m) A
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
% a; a5 R% ]7 s4 ~( M
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   9 a+ U\" t, u& f  ~) N
2. mvar:* ?/ R6 q( a! f  x8 H( U% K* f: P! p\" s
   
3. t=sys::clock(),9 \$ f' J- J2 k! J
   
4. s=0,
   9 Y/ p, C. x- c* A. |& c
5. x=0, while{x<=1000,
   + A2 V9 A! F- Y& ?# {4 C
6.   y=0, while{y<=100,$ \$ u, K. b+ r0 H9 C, ^! V! o
   
7.     z=0, while{z<=100,
   - ~8 i& \1 y; D- C: G, r; y
8.       s=s+f(x,y,z),
   $ C1 _1 k( _& [
9.       z++3 g' X. H! E3 C* V\" V1 H# ~) B
   
10.     },, P* o3 Y$ E2 g9 t. U6 k
    
11.     y++2 d. f7 t4 r( E8 o' `5 v+ @: `$ P
    
12.   },7 j' }( }% J$ V% M) X2 I
    
13.   x++
    ; W  |\" P/ m  Z5 S6 k
14. },
    7 {. r: j6 M( a& F
15. s;
    ' R' I% A% ]: f
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
0 h: d$ k' W1 N4 L
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
4 S$ h4 E- ~% X: S9 @

hxftntA

3 O3 F3 i* C7 E9 j& I8 x
7 N; q. G$ R6 Y/ `. r2 t5 m' W: G: A



& C/ L* o2 o& D2 }5 x! A/ e  D

4 I0 z. M8 y& v& ?  O# D




% }* e! h# H3 C2 V/ g8 Q* D4 D

+ b0 N; G: L5 H2 ^" }+ t
$ q6 s/ G2 u9 b& l) b2 ^
3 s) G2 |/ f) `8 l