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

forcal

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

7 _% E  j8 r# N4 f. v* N＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   4 f, v: v( V9 s% M. r1 s1 H# r
2. tic;8 q+ F% P; X/ @- v
   
3. s=0;3 K+ ^, g- d- K& a1 Y
   
4. for x=0:1000
   # ?/ r) {4 J: C& Z* b+ B
5.   for y=0:100. F, _7 I. x, ~8 ^\" @, ~1 ?
   
6.     for z=0:1001 v! k' Z* t4 R) ~2 B+ M: F
   
7.       s=s+f(x,y,z);( y6 n8 G/ F! o  c\" Q6 b, }! {2 y
   
8.     end' s* b; W5 n% a  c3 K& j# J
   
9.   end
   , X& V) Y, R& V- Z0 G: F\" z: f9 R
10. end
    ! n( x' U\" e7 e$ h4 d
11. s. [, I. [  u4 i\" ]9 h. @
    
12. toc
    ) l2 }, G0 ?6 \' p  g
13. 5 ~  b! F' `/ s& ^6 t. `/ Y
    
14. s =. Z. o6 f1 l% P! C3 X
    
15.     6.126720600000000e+009
    & o3 e6 ^& f- }+ L( p: @1 j
16. & a8 P' ?& H; X# J' Z
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m* l3 i\" n8 K# I' _4 X6 k- o, s
   
2. function c=xyz(x,y,z)7 t0 [' U$ h\" M! C+ X
   
3. c=x+y+z;+ V) l2 }  U- L
   
4. end

复制代码

测试代码：

1. tic;# s) g( T( m7 ~3 `7 F8 }* ]
   
2. s=0;
   , k, ~0 U8 Q6 b7 t& t, d  u, S9 k
3. for x=0:1000
   & Z\" W: m' k& }' k, i  M: {
4.   for y=0:100
   & q; i& ^9 m8 ?3 Y8 b8 ?  O
5.     for z=0:100' |% |9 |* ^$ _: |2 Y  f& J9 ^% D
   
6.       s=s+xyz(x,y,z);: j, f% v* O* J, k$ n
   
7.     end
   2 H\" [' [+ p& D7 y& I- w
8.   end
   % G7 V0 L\" S# [
9. end
   , g5 k+ I6 z; P
10. s
    ' y7 r# D8 k: U2 W/ ?; x
11. toc  n0 T6 n7 B: l1 i& e' P( g& }* E
    
12. 
    ! q: o. E% c( c# G0 Y
13. s =% i2 Y0 t% N+ B! s4 Q; `
    
14.     6.126720600000000e+009( {8 l2 x  \1 _+ i- Y
    
15. 3 S5 e* m& H9 f
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
" {! C/ m; q5 A4 p
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   * W; B2 E) G; ?- Y\" E4 q+ z/ A
2. mvar:
   - \7 Y# x% b* [
3. t=sys::clock(),
   0 i: h) `' d1 {% c% e. P- Q
4. s=0,1 N4 }8 y! c) a$ Y+ }3 L8 O7 ^
   
5. x=0, while{x<=1000,5 A4 b, m/ }# V& f9 n6 y
   
6.   y=0, while{y<=100,
   - ?& [8 N/ C  s2 l% [4 L
7.     z=0, while{z<=100,
   ) U$ l! m. y3 i5 q
8.       s=s+f(x,y,z),
   * K% i0 x. D# U, S. H! _\" [5 k2 @1 q0 K
9.       z++6 ^7 `: x$ B\" ^
   
10.     },
    ) j1 J3 p2 ^' z
11.     y++
    : k6 c# W8 {& @% N6 y2 q
12.   },, I* |, t# N) ], h8 C6 Q% J0 @( o8 L9 @
    
13.   x++
    & d: i9 a6 N0 J! U\" Q
14. },
    ; [5 ~3 j8 ^: F\" D+ \  O9 {9 [
15. s;
    - B, ^7 B9 u$ ?& s9 j
16. [sys::clock()-t]/1000;

复制代码

结果：
9 W# l% k4 D! s1 V( ^. Q9 T6126720600.
8 r2 B" U8 w' B( ~1.109  秒

) n: a; j7 Y9 V/ {) C, D7 B; g, D二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
( Z) b5 P( v! t  m6 v- l

hxftntA

% c+ J8 b$ S8 S+ m



+ I+ D% Q7 K; t% A. y) h! T: I: I
% i# P/ h9 b: |, B4 I7 X8 M

8 o0 Q$ B6 A5 K: G8 F( e& F




( h9 H5 ^7 ?: t9 X) Y3 f
6 t/ Q' x( b' T6 A4 \
2 M$ D! L, D8 O) a5 f: Q) {( r
+ [0 _$ o6 U, X7 m1 f6 V* L