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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
( x# |" ~2 y. g
8 P7 I1 J6 t' g＝＝＝＝＝＝＝＝＝＝＝＝
* y% O! N' l* J' _5 `5 L+ i
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   7 D( V* _$ D6 z3 c
2. tic;3 @9 y! K) P/ H\" G' m/ h
   
3. s=0;7 s' ^# g, A/ ]# Y8 @
   
4. for x=0:1000& B$ G# n. o+ W9 B8 u: }  ]
   
5.   for y=0:100
     J\" U. s: {: W1 E5 A6 G: O$ G
6.     for z=0:1009 N1 w2 d8 Y# S5 S7 K
   
7.       s=s+f(x,y,z);
   9 C: Z1 K0 [\" l: i9 D
8.     end$ Y3 N& o* C) W4 Q
   
9.   end7 @9 L4 L: v; U7 L
   
10. end9 t/ `. z: d$ t! a; L* M0 @4 E
    
11. s
    \" V8 |0 s- [  e7 l2 i5 U; F
12. toc& f! E) v& D  w& m
    
13. ; l+ k* [* n9 ]5 @+ e  x
    
14. s =
    + E/ R% g4 C- P\" C4 k/ Y
15.     6.126720600000000e+009
    : T0 N' I1 O2 t: ]\" K
16. 
    2 d  P; E& O* X, K% D2 ?
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m4 w+ n+ u7 q0 N& D- z
   
2. function c=xyz(x,y,z)
   # C9 J5 Q4 l\" {' R# V* F9 }
3. c=x+y+z;8 `( [/ r; `) [: M: ^
   
4. end

复制代码

测试代码：

1. tic;
   7 r4 U; a. e# c6 M$ b, O
2. s=0;5 p- z7 L+ n2 V$ q; X\" I% }
   
3. for x=0:1000$ T9 n  g5 Y5 w  P7 |' z! w
   
4.   for y=0:100) ~) Q% ^/ F, N9 X* Y: B. ~
   
5.     for z=0:100+ H1 N0 N! v$ ]' [
   
6.       s=s+xyz(x,y,z);
   $ U3 l+ y% C, h
7.     end
   8 S. l# }' g2 w9 E
8.   end
   ! A  x0 c9 V% E% n0 t
9. end
   0 }' t$ d/ Y! [( M1 ~. m4 o
10. s
    & X# E/ q' H8 H* R( ?
11. toc* U' x6 \. ^& {9 ~; A* l! y+ z8 h
    
12. ) i; A' B0 U+ G7 {3 ]0 S\" a1 ^
    
13. s =
    : E) P) z3 l% P$ w
14.     6.126720600000000e+009* W% ?& C3 y' Y8 ~
    
15. / O- k6 W  `0 o\" Z' A
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
0 S5 U6 \% c! Z
' G* f: o( M8 V, H' L/ hForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;  R9 T( t& i8 _. l( H
   
2. mvar:* K; S' ?4 Q2 Y
   
3. t=sys::clock(),- ~9 B  T# F1 S/ G+ C. f$ O+ |
   
4. s=0,: r' A  `0 m  T3 R$ s3 z6 [. T% ?
   
5. x=0, while{x<=1000,
   6 n' R5 T& I' [+ {: c
6.   y=0, while{y<=100,1 a! t& R# M: t\" R\" r/ T( e+ P
   
7.     z=0, while{z<=100,4 ?6 v' s3 V- |1 ]7 A, I
   
8.       s=s+f(x,y,z),
   ( e. C. b$ F- W% G7 M7 b! U
9.       z++
   2 S# {' \5 K& E' K9 W+ M- {
10.     },8 E) b\" n2 l0 H1 N( a
    
11.     y++( ~. @- E; V! e
    
12.   },# J$ M% v! f% l. h0 F
    
13.   x++
    0 B$ Z! |; B# c# r4 o. _5 f
14. },8 u& S- o' Z9 y% ~
    
15. s;
    * g: g4 S1 i( ~' k, [  T
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒
2 ]: E. p8 Q# I
二者效率有几倍的差距。

hxftntA

! U7 e; E9 k8 n' A" I' z& ^




( V! m* x4 X, C, _4 G$ L

/ j7 L& w& n, z% M



0 M$ o, [3 r( F, N% [( E8 U  |

8 K, C# G' l3 d
: i, I( H5 A0 ^, P, \8 N

ehi28

嗯,不错,支持一下.
+ d& l7 ~  O9 K/ Q* A9 s+ F- j