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

forcal

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

2 f. R/ A  E7 ^6 I1 I8 }＝＝＝＝＝＝＝＝＝＝＝＝
9 J: F8 R7 _4 H$ x& v9 r) Y! V
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;, ~0 \$ e\" r) x) \# J& d2 f, H
   
2. tic;
   & }0 G. j4 V' A  C4 R/ [; `
3. s=0;, O0 e( P6 u& |) [( i$ v% Q
   
4. for x=0:1000. I4 l\" _* G; \% r* \: L
   
5.   for y=0:100
   : H2 s& f\" S5 d% w+ Q  B1 |
6.     for z=0:100, O  u7 ~\" q# B8 V- l$ W& v# }
   
7.       s=s+f(x,y,z);
   1 }' n3 C- c$ \. t& R
8.     end0 Z$ R7 o/ c- G, i1 ]
   
9.   end) R0 T+ c( e' b9 M7 K9 j
   
10. end
    + y% L: t# B- U3 j
11. s
    5 Q: |6 b/ w. T* E! G* t
12. toc
    2 I& W9 j; M* l  s0 C; g2 V8 P9 X( m
13. ; n# F+ T0 d$ V* t1 a, ]
    
14. s =
    4 I  p, ]7 M: u9 Q) r
15.     6.126720600000000e+009
    7 ^( `4 s% `. N) S# x) {9 q# y% P
16. 
    ( v% v% O1 K: }0 i
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m0 l, j\" u. u! E4 x, F6 v! I/ O
   
2. function c=xyz(x,y,z)
   / \- J1 v% u0 E5 [
3. c=x+y+z;/ F# G2 M( e\" _( F) c2 Q8 F$ h
   
4. end

复制代码

测试代码：

1. tic;( S2 x\" u# r7 K7 p7 A6 O
   
2. s=0;, y0 m6 b: F8 K' U& v/ u
   
3. for x=0:10006 T2 ~, D1 H: J* {7 N
   
4.   for y=0:100! |$ J2 K( W7 v  `' e
   
5.     for z=0:100
   ( B1 S: ?7 g4 W. ~) J4 I  w\" i8 }9 v
6.       s=s+xyz(x,y,z);
   4 S: s( g; K. A& P4 d\" ~0 P
7.     end
   \" q8 u/ M% t4 v( s; D/ t2 O: N% t
8.   end! H! H6 |/ J- }4 f
   
9. end- q0 z6 A  x$ |. M
   
10. s
    : I' m/ w* m) J: ~5 b+ Q
11. toc6 V% {9 X  ~+ @
    
12. , D* z  D$ i$ j4 a# ^
    
13. s =
    : w7 s7 W% ^7 C# S$ C
14.     6.126720600000000e+009
    - {( ?5 b# ?, X! J$ J6 z
15. # j/ u8 d\" D$ a
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

2 i. i& T* A3 y" h+ J( }' KForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;, K/ l  @  l! [7 ]0 O
   
2. mvar:
     v/ @6 b! F1 |0 q, c
3. t=sys::clock(),: Z3 f4 D7 n) e# p
   
4. s=0,5 [# H3 I) }. f
   
5. x=0, while{x<=1000,/ P4 C: ?, M& x$ g
   
6.   y=0, while{y<=100,
   . S- G7 X. q2 g7 _+ |% h7 K1 u
7.     z=0, while{z<=100,3 f. _$ e2 ]& k$ F8 @* n8 V; ?
   
8.       s=s+f(x,y,z),% e, A8 c: P. L2 g
   
9.       z++( d( l% o* ?* ^6 D6 L  V8 h3 |
   
10.     },2 s/ @+ T0 F- W% [2 f- R
    
11.     y++
    0 L& G. H+ D- i8 B
12.   },
    \" U2 u6 ^+ I& @5 J8 b6 ?
13.   x++; u) n6 W2 l0 m6 @2 c
    
14. },
      _. J7 N  s5 L2 [; v
15. s;
    ; w- w! C! I& x
16. [sys::clock()-t]/1000;

复制代码

结果：
) M& z' n# S* @3 b6126720600.
$ I' x6 s3 E& h: T4 d! ?  D" F1.109  秒
* Z3 p9 w# c) H$ c' f  E  N- q0 t
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

# F% O9 a. K* N: O  Z

6 R3 T- W8 c- L4 u9 X

7 Y$ p* x! Q; n1 C

8 T9 W9 `2 ?, Z4 q. ~9 q# x: ?

! A6 l' Y! _3 ~5 {7 n
  E" V6 A' |- L  g  O& t( y