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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
0 z: x2 |: H$ K1 t
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   \" P+ A- H( L& ^9 p
2. tic;4 f1 d% o# [, U$ S: e9 }1 F
   
3. s=0;
     R8 j, T* m# L  |# Q\" L
4. for x=0:1000
   6 J/ d$ Z\" O7 S. j0 A% E
5.   for y=0:100/ w$ M! L; H! |/ M: U/ @9 C
   
6.     for z=0:100
   ) l& I1 `( I- F$ ^  v
7.       s=s+f(x,y,z);4 q& E; R1 o; ]# |, ]  p
   
8.     end
   $ s( b; l6 A* T$ K; j, z
9.   end8 M0 O7 l- f, f. H9 O, i' g0 F
   
10. end
    , L- j' w% }% o, D; @4 q
11. s  V1 j; y3 z8 `1 W$ D/ K
    
12. toc& C. O$ p7 H) s9 `\" H
    
13. 9 t& G2 r  B+ @, s
    
14. s =
    + _, x! M& _! i5 n4 s
15.     6.126720600000000e+009
    9 `% Z/ e: f. m
16. & T; j4 e- G1 T% P
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   ! }1 Q& ~8 z) U  B9 o* o
2. function c=xyz(x,y,z)! I* ?/ X5 U8 D0 O- o
   
3. c=x+y+z;
   ) [; G5 x/ B6 G; z9 M3 M+ U& X
4. end

复制代码

测试代码：

1. tic;
   3 I& B6 P4 B# Z
2. s=0;
   & U9 n2 d, \. C
3. for x=0:1000\" g6 C9 ?1 M- X: X5 l8 Z/ O7 @
   
4.   for y=0:100* V/ [2 R\" E8 [& p( x' Z
   
5.     for z=0:100
   + A& n+ @! S6 I( I- F1 R; j1 }: N
6.       s=s+xyz(x,y,z);
   7 S\" ?* X, y: |0 M2 P\" z$ o
7.     end
   3 N, y0 b2 ^\" Q3 U8 }0 m
8.   end0 x& Z2 c% q) A
   
9. end
   ! T5 T$ h# q9 ^/ [  X2 c7 N
10. s) n& G( Z' j+ t6 p3 z7 d. W
    
11. toc
    9 }  X2 h0 d8 s( B8 c5 e9 q7 }
12. . N( u( g/ h7 y
    
13. s =. z5 h5 i4 t0 W2 o
    
14.     6.126720600000000e+009
    - `1 [( Y4 g$ D* _3 B9 I. H
15. 0 M  F& D9 f7 L0 J5 p, X- B1 x2 s% G
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
/ L) [+ r9 x7 ^
% B2 t) h3 K! _( s' aForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   $ m  {1 [0 o* o6 f. w9 P( k
2. mvar:
   & B/ q1 v: j* {. D
3. t=sys::clock(),
   ! T1 Z# `/ G; p! @0 b
4. s=0,! D* [3 H1 v6 Z/ y' |: D* f
   
5. x=0, while{x<=1000,* R% L. I\" e! y. L
   
6.   y=0, while{y<=100,* ?: C$ y) m' X+ [& a# z8 O3 a
   
7.     z=0, while{z<=100,( n9 q; ^1 j0 B/ V5 C
   
8.       s=s+f(x,y,z),
   . a7 X( I6 B4 `/ h
9.       z++) k, O. T- T! o# e- d6 H
   
10.     },# l0 b, M. f$ g0 x8 O7 b9 K( j
    
11.     y++, [3 q8 |7 J) A
    
12.   },
    % w' l' _! i, b) D
13.   x++1 g2 [9 r: r/ X; U
    
14. },
    8 U$ f6 C% b* `3 O7 W
15. s;
    ( Y2 ^* F! `2 v- g\" D/ g; G
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

. J% g, T. O2 x' d2 p二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

. h4 W% Y. E! D



9 {# O! i4 t+ N3 o) H1 p$ E

3 z3 t3 ~) a' c$ `  W; ]
0 i/ P0 y: r# u4 {; L7 r! Q

$ o1 u: K( t1 u2 d4 z
0 l/ J9 g# d$ d3 w
' X) Y) r# S" Z

! ?! W7 F4 n" A0 u9 K& e8 w
4 u# d! b" B$ m6 w2 c- ?9 W$ i
0 ^& ?1 ~8 f0 Y6 |: s& J) m