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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
) M3 H! M5 t( l: k, o5 k1 Y9 J4 Y
3 ]% M# n& N2 KMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   ; N\" a9 E+ x8 x* r+ k% h5 ?
2. tic;/ b/ X2 O* z0 c\" [+ a
   
3. s=0;% T) i4 ]9 A$ j
   
4. for x=0:1000
   7 A/ h8 h' A! r8 ]! f8 E+ i! k- {
5.   for y=0:1009 K: q; O3 ?2 q\" [- T
   
6.     for z=0:100
   2 ]! X6 r+ z6 B6 a5 H
7.       s=s+f(x,y,z);
   0 F; a1 X$ J  y- w% a
8.     end- k9 b7 z) q$ a2 F7 l# r7 Z2 R8 v
   
9.   end) K' j! f2 c5 b* z2 ?
   
10. end
    & E7 j! S\" ~& I! ~$ L; P  A
11. s\" ?+ c% n1 i9 H/ b. ]
    
12. toc
    8 X$ l: G2 S  \$ a- r$ `' g% J
13. ) u2 k4 L# `3 E! f! ?5 h
    
14. s =
    / B2 g+ F/ G3 v2 d# B, Z7 O/ D: E# H
15.     6.126720600000000e+009  n* M9 R- P# f1 S1 a\" i# l# k
    
16. 
    $ g+ s: F& d& r9 g& J$ e% [  a
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
, s$ F9 T3 w9 f1 F. lm函数：

1. %file xyz.m
   * Q$ A% N% w, A, y: D
2. function c=xyz(x,y,z)
   4 \4 ~9 r/ K* j9 P7 \! P9 R+ x
3. c=x+y+z;. R3 Q: u- d* F\" [' \& a5 }
   
4. end

复制代码

测试代码：

1. tic;
   + `+ X- w% L+ E: V
2. s=0;/ b. Y- v, y/ c# M2 X8 f
   
3. for x=0:1000
   2 U& ?. s  G4 L; P
4.   for y=0:1000 e8 y\" L3 y! Q! o3 v
   
5.     for z=0:100; M) p( S( D8 g( X& s
   
6.       s=s+xyz(x,y,z);\" p# [! Z: {8 b
   
7.     end1 Z8 R9 J3 Q/ @
   
8.   end3 t) O. T: D. s& d$ m
   
9. end
   4 m$ L. u7 F& O. s6 m) |
10. s% n' x# q/ Y' X; g+ ?& U. ^$ h, T! g
    
11. toc& ?1 d/ k0 R1 i. E  C* P
    
12. . ^! n- j% J% I' p5 `
    
13. s =8 m- \- h0 }# c! |+ g
    
14.     6.126720600000000e+009
    . a% p0 H  [9 N6 Z- R/ q, t
15. 2 s& H& K! O& V4 @9 M; l: ?/ d
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
5 `0 y0 }! X$ E8 L& ]1 T
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;8 t0 v6 m\" o\" Y2 ]
   
2. mvar:
   / z6 Q1 y+ b6 W5 t2 Z6 ]' z
3. t=sys::clock(),8 s$ e# |7 b/ [3 K7 s8 h. U& q
   
4. s=0,# e( a& J5 p2 g% Q2 I1 l3 n
   
5. x=0, while{x<=1000,5 v( Q2 G- J  K# n
   
6.   y=0, while{y<=100,
   0 E# I; B, p' r7 r6 \
7.     z=0, while{z<=100,- p+ [2 A0 C( o3 K$ A# N+ T5 A4 X
   
8.       s=s+f(x,y,z),
   5 f' e& \( H- ^9 ^8 A( V0 B
9.       z++$ S0 n\" T, r& @  t/ \
   
10.     },5 D# F+ H+ K) D- U) [$ f
    
11.     y++1 d7 ?! d& d/ j' f' G2 {2 l
    
12.   },2 l: y8 Y* P9 c$ ^
    
13.   x++
    * \! ?2 M* ]' ~$ r' N5 G. r$ P
14. },
    \" {% C  }! Z3 s% D3 R) x
15. s;: f' v. v, k* Y$ g/ `. V0 P- e
    
16. [sys::clock()-t]/1000;

复制代码

结果：
4 t( i$ o# K, S6126720600.
9 z  G0 x& I6 V. e1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
  i5 C0 ^% O* F: U4 J

hxftntA

2 _& ]' V" k0 I7 i# A

7 B9 R# `# u1 O  c: _. U( l

7 p& j% O. P% g' U  s
" a- o: _! F) c) Q" ?


& Q( N6 u4 G$ W1 e) C( c% o
- z: d; u; r; e. T! i' i* B


( B# K6 O6 ~( n2 b+ Y: q0 k/ u


& b1 I- j: _3 A1 ~* ~8 E