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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
) m3 {" h; T2 ^1 z* M  S( Y$ v" e; G
- K. r7 f5 a% x) ?* o" TMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;7 H1 L4 W  K2 M' S
   
2. tic;
   ' L3 k; i( m. _2 z3 }* i
3. s=0;$ p, ~+ _5 p1 J2 E5 v5 {
   
4. for x=0:1000
   ! |( `\" Z; |\" G( U
5.   for y=0:100
   ' j4 R9 h. Q6 K- \
6.     for z=0:100( X2 u9 C& ^% |6 t3 c9 w: q\" ~
   
7.       s=s+f(x,y,z);+ G; H4 V\" l$ H5 L, J
   
8.     end
   9 S5 W, [2 h! r2 K/ ]# c3 `8 d3 p! ~\" C
9.   end
   3 w6 t: x2 }1 |  O8 Y  P
10. end( |2 Z1 f  n+ p
    
11. s
    7 w' }2 E1 h. F$ R; i5 P2 j3 Z2 G4 e9 U2 p
12. toc$ q9 y; X) p4 ^
    
13. 
    , _( \$ @6 y% h% N9 B' l6 u2 s; |
14. s =( h8 o) F5 o! B
    
15.     6.126720600000000e+009; H, }% t7 B  J9 }\" o
    
16. 
    0 U! P4 Y\" v+ M5 x, m+ z
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m8 Z  A) c  o\" N) w; m- O/ p
   
2. function c=xyz(x,y,z)
   : V4 K2 G4 f1 L
3. c=x+y+z;
   3 y5 E1 Y2 V4 T& R3 Z, v5 Y5 l
4. end

复制代码

测试代码：

1. tic;
   ! q8 W3 U  s3 P% g, E1 B0 [- ^0 j
2. s=0;! {6 Q7 z8 k7 s' g5 h
   
3. for x=0:1000
   / M2 h2 B+ c8 H\" l\" {
4.   for y=0:100
     T, z8 q& B. M
5.     for z=0:100
   ; }+ f0 C# ]& o7 G6 y9 |
6.       s=s+xyz(x,y,z);5 e8 B8 Y# Y# ~- b+ q1 w\" w
   
7.     end, \% s+ ]& L7 w8 J& i
   
8.   end
   ( W# f, D1 k4 V( B2 E8 N% W
9. end6 ?\" `8 |& d- d+ P2 G
   
10. s
    % t/ q) X6 S- x; Z
11. toc7 E+ [! X/ m, Y5 [: k' k
    
12. 6 |& j) W' A! x# Z) \  _
    
13. s =# q0 V0 q3 e5 `1 J, i
    
14.     6.126720600000000e+0099 @( B5 \% \+ O# i  P3 ?
    
15. 
    6 S% y# ]  S- W\" b. x4 O4 g
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   $ m$ O( W) z0 Q2 L8 [0 l: a6 h2 a
2. mvar:
   3 N- F: ^# O9 ]  m* P
3. t=sys::clock(),
   ( K% X1 ?+ N$ T( ^8 E8 N# M
4. s=0,9 C  J7 R8 \+ [- M, o* e
   
5. x=0, while{x<=1000,7 ~\" G: F- ~$ i7 e  ~, j. j8 y
   
6.   y=0, while{y<=100,0 C# L& S+ C3 I' G' m/ D
   
7.     z=0, while{z<=100,
   * Q3 V\" T6 T# Z
8.       s=s+f(x,y,z),1 M& _( S, S5 U3 p+ Z
   
9.       z++5 @) q% K6 e! K. ]( }5 s, N. D: h, {
   
10.     },
    \" A' I9 {. R7 K) b* x. ?
11.     y++8 c- o- ]6 U, u
    
12.   },3 X0 U7 w' f/ k& q. r+ [/ l$ Y
    
13.   x++; ]9 S! V! p7 x6 |# E
    
14. },
    , V0 Y\" `0 Q1 K1 d: Z
15. s;\" W/ Q: C3 r/ i: A
    
16. [sys::clock()-t]/1000;

复制代码

结果：
4 \$ z2 V. \. g  ~, h8 v% {6126720600.
1.109  秒
# O, e9 [3 R" X/ B- }- O  C
2 x& i  M( ~! f; g2 \" c8 h二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

9 ~3 F' B8 G3 P/ F% C$ f
8 j" q$ Y) V0 v4 C3 z; k

7 f. Y# K8 R; o- |* Y$ v6 R
! x; i9 J4 }2 R" D& b5 D



: @* m- |) U( F5 q
0 l3 i5 C8 A* g/ A: ]9 m0 [
( g8 R5 u+ Z3 }7 c) I+ |" E