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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

! d: q! S  d7 N7 V( [: R) KMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;! \) e6 g! V. A: e
   
2. tic;$ v* |; L\" U, r( C5 R1 a$ j
   
3. s=0;
   8 r1 T- e; {9 ^( F! E  q9 l
4. for x=0:1000
   , y: v3 V! z5 K7 L
5.   for y=0:100
   ' Y' o  |  }- a$ E3 }$ J
6.     for z=0:100
   / l! C0 P9 L; w7 U9 w, d
7.       s=s+f(x,y,z);
   # t# X  R\" }+ I/ S2 m
8.     end
   $ L% [$ n* C3 Q7 M
9.   end
   4 a& M% y- F% l
10. end  X\" c8 L7 u: ?3 L3 x
    
11. s; L& W4 x+ |' V- v+ }, e
    
12. toc
    2 P: M\" L2 L) p+ ~6 ]/ {5 _. P3 i8 W
13. 
    6 \$ U3 L. ]* B2 R9 J. ]! q
14. s =, L% z; e$ y7 a  Q( k2 R, P
    
15.     6.126720600000000e+009
    7 Z8 `6 ?& _# t( {) V/ J# j
16. 8 C6 n8 K  s% u# A
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m# u' O; @: i! k! B
   
2. function c=xyz(x,y,z)0 L9 Y$ y  E8 j$ O: |
   
3. c=x+y+z;
   : u; }0 s% E, M, N
4. end

复制代码

测试代码：

1. tic;7 @4 Z/ v; W3 ]6 W  M
   
2. s=0;
   \" p6 F0 g$ J4 E( u% a1 z7 m, H
3. for x=0:1000
   \" M3 e$ t- m3 K& m! k' c' j7 _
4.   for y=0:100
   # r/ |* A9 l( D# ^
5.     for z=0:100$ ]\" v; x( |2 P; s- Z7 F. |
   
6.       s=s+xyz(x,y,z);2 N. g  w4 X# g4 T\" O# R
   
7.     end, T6 a3 i2 l2 v7 K
   
8.   end
   8 f/ ~: M- T& h( t2 W
9. end
   # U) @  M, B2 I* m8 J! _0 ^. r
10. s) K2 X  x4 T& e/ B
    
11. toc
    ; U8 Q+ r; @; r* g6 M& E\" b
12. & k8 i4 L  ^, ^& h
    
13. s =
    2 x+ X- N0 O# x3 I& N  Z
14.     6.126720600000000e+009; _7 Z: D4 A* ~% w\" }/ F. V) I
    
15. 5 N1 n' c9 d) q, ^3 ]! }
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
- G2 w* v! j# q/ w# u$ T
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;8 g) d) O7 B/ h
   
2. mvar:
   2 Z1 z2 e+ f, @
3. t=sys::clock(),
   4 w/ T7 c4 Z' E8 t/ q# i$ K0 F# O\" z
4. s=0,
   ! o2 V& J) U& J, ?, O
5. x=0, while{x<=1000,
   # ~* S) _* j: |, }
6.   y=0, while{y<=100,
   / [2 C+ `/ \/ l
7.     z=0, while{z<=100,% G- k8 M+ b5 h3 D( e/ Y2 d' d. b$ b
   
8.       s=s+f(x,y,z),3 f' k  a. U% O\" K3 @) [. y$ \8 p
   
9.       z++
   2 p. V+ X$ E! ^3 [$ ?8 S+ Q
10.     },
    ( r\" y8 y- @1 q
11.     y++: ^. C\" B3 Q: \/ D1 ]0 {
    
12.   },
    ! o0 _$ N! [8 w% l+ D
13.   x++) b\" u6 V\" q$ E0 a5 \+ e) s
    
14. },! H) N3 b; y' H. w9 o* X7 m# p
    
15. s;
    ! I4 T3 q# M' ^6 q3 ]2 h5 A; \' Z\" ^
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
2 k. T" G: i3 u2 }# P& I1.109  秒
# Y  b7 b+ F* i. p  S) R3 C
二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

+ X1 M  ?1 o7 _$ i

8 @/ H+ d4 w* Q, |: D1 W& J



; \  u- B7 A9 Y


4 w7 B- {8 m) m( s- @/ a
. B  j$ `' `# R) \2 c3 p9 e  Y, t

( b: l$ W/ P5 B' b0 G
' x" w- p0 F0 y% C1 {& n
7 D* p7 F8 \' i) P2 K3 e$ m
. R3 k" K, K! t2 o, p$ e
4 F$ A. q6 T$ ^! f7 D