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

forcal

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

! r4 P% n7 e6 I0 i0 J+ T& G' V＝＝＝＝＝＝＝＝＝＝＝＝

8 Q: d5 ~; F$ dMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   - R# K  u, L8 v\" C! q& e8 }
2. tic;% m3 E$ ]+ ^: V  x0 s4 {3 S
   
3. s=0;- I) d# t0 ?7 S' H. }. ?/ }
   
4. for x=0:10007 j, @: j& O) X6 \: n
   
5.   for y=0:100
   9 F1 E1 R9 z+ [  Q. K- F7 `
6.     for z=0:100# }\" x, Y! N8 ]& k
   
7.       s=s+f(x,y,z);
   ; v; R7 k; y5 X; G- `* k
8.     end
   # c# B4 |# s! L2 E* {, ?3 E
9.   end/ e# q3 `8 P\" `5 x$ ~
   
10. end
    3 h2 S5 b) `& Z! ]  u. W' P
11. s
    9 C. G% p4 x2 R& J5 A1 n  v
12. toc
    : G7 y/ J8 P3 j4 R
13. 
    ( m  s( W6 ^2 w& I: p- Y
14. s =
    4 T) |% Q, `4 `# x
15.     6.126720600000000e+009
    ! `! Q* t' |4 T  K\" R3 B
16. 
    / x& V& R6 B7 n3 _% @
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   $ P8 K% Y+ @& q1 I8 A- G
2. function c=xyz(x,y,z); S  t. K1 J2 S1 i7 S4 K
   
3. c=x+y+z;: K$ W2 x* v% f\" k4 t* L1 d
   
4. end

复制代码

测试代码：

1. tic;
   \" ^& [5 L- T. t+ _4 r
2. s=0;
   5 Q5 G( i$ A- X; }
3. for x=0:1000
   / A# [# P9 T5 ^* w
4.   for y=0:100
   . D* q2 p/ o% E; X# M! ]' s( h
5.     for z=0:100
   + m; g8 u) S' ?% o5 O' f. u
6.       s=s+xyz(x,y,z);/ l* I' K% c) ^8 s. W2 D5 m
   
7.     end( Y3 |, J3 i8 U! J( U7 m9 V8 M
   
8.   end
   7 }2 u. g0 I1 V$ I# O\" b) d
9. end
   6 C- z* o1 U4 d7 V& |* E; j2 s
10. s) r* }2 ?7 R4 \+ o3 O# p5 S+ I( B, @
    
11. toc
    ) M( }2 }# c/ D
12. 
      I) E\" K6 l( }3 F5 r) w2 v
13. s =7 I) u7 D7 E3 a8 ~9 Y. G1 Q
    
14.     6.126720600000000e+009
    ; P3 t\" y# G% s- W
15. 
    # @9 z\" R1 D1 ^/ ^3 m* T. g  j
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ) p/ O+ \8 j$ N1 q
2. mvar:
   : d8 B$ ?& N2 V' M  B  v. _- V
3. t=sys::clock(),
   0 W5 L; J4 D+ i: E& y( y
4. s=0,
   - k2 D) d3 P6 ?/ {& \1 t
5. x=0, while{x<=1000,
   : R$ d/ Z( a7 y! H
6.   y=0, while{y<=100,6 C! t1 u: r* _7 z! }
   
7.     z=0, while{z<=100,8 ^. b\" }6 V. T/ m
   
8.       s=s+f(x,y,z),
   4 t6 N7 E2 u0 Y/ r
9.       z++  {9 {5 X# b8 w3 U! B4 ?! E! X9 t
   
10.     },
    7 i, J1 v  V1 r; Q4 E5 L
11.     y++4 n9 K8 U( Z- N; n) n# g; r! w
    
12.   },2 F$ Y! R* M9 ~5 I\" n5 {: X
    
13.   x++! b. f! q. X; Q8 j+ t' J5 Y! I
    
14. },
    / N! K) }0 {\" t4 T8 D. \
15. s;3 I! \% ]) S\" e  z. f
    
16. [sys::clock()-t]/1000;

复制代码

结果：
  q* m. i7 d2 f8 R8 S8 t6126720600.
! O: w* `& J" O& W, S  R1.109  秒

# I0 R4 C4 Q: \' u% l; y8 C- z  `! |二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
' E8 ~9 v3 T! U1 [7 ^% g# A4 w2 ?

hxftntA

8 b2 c1 W1 ]% c$ `) [9 l

' \* H, q* _# H" B; s3 r  I

. |5 \9 [, m) e
2 b: W4 {8 ?* A4 k5 r1 W
2 _/ |# p( T6 R/ b! I( L


) r) r+ v0 q) R/ C& T, C