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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝
! Q/ q, g7 E3 s, D3 h% A8 b
5 L3 y6 [- N/ R, N, Y/ W: {- d" L8 MMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   4 t0 }  H# z+ T$ w
2. tic;
     A3 T3 Y' V8 T
3. s=0;
   + w3 C* X1 `/ e+ \1 _1 B0 s
4. for x=0:10003 N8 G3 e( b/ E  o
   
5.   for y=0:100$ W* b6 o: [# X\" f. `0 Z* e
   
6.     for z=0:1001 Z1 J9 j- d5 }! [
   
7.       s=s+f(x,y,z);3 _: e% z2 j3 E+ [\" _& z8 V4 d2 T
   
8.     end6 N3 l3 y3 N1 f
   
9.   end
   ! e) F+ ^# H: H' [  j4 [8 T
10. end
    3 Q' K5 f# C/ j! g
11. s, I+ D' x1 S6 w
    
12. toc/ Z3 [6 d& x* h0 C4 u5 s
    
13. 
    * t- g, g$ L  o  S7 Z# N
14. s =
    3 u3 v/ h, L+ l4 O1 Z$ A
15.     6.126720600000000e+009: ~' v% h3 k& x4 E9 a6 r3 B; k
    
16. 
    3 `. y1 x* s! z4 Q  D% d1 k
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   + s: l, r( c( c
2. function c=xyz(x,y,z)
   , B. ~5 ~2 k6 y: V* u) E; H
3. c=x+y+z;2 d$ _  k8 w, m% X8 z3 I: v4 d
   
4. end

复制代码

测试代码：

1. tic;
   5 U- E( o) z( u\" |  {
2. s=0;4 o) Z8 ~* |$ J' S\" o8 @& r
   
3. for x=0:10006 X! Z! Y7 W( W5 S4 H2 y& d8 n
   
4.   for y=0:100* t, _; K. Z$ p% l
   
5.     for z=0:100$ Z9 Y# c- |: q$ Z
   
6.       s=s+xyz(x,y,z);
   1 @8 N$ \0 U+ H+ x$ k, e' s& A& l5 I
7.     end
   6 h  u$ C7 S* T. y  q6 A! ]4 U
8.   end  u, N& G2 G' v
   
9. end8 G4 ~# S- _1 L8 F7 I' h0 F- n8 W
   
10. s1 V* L- b6 u* w, T
    
11. toc7 Q% h\" c; @- Y* c8 v9 t$ M9 i
    
12. 
    3 I! v% p+ d( I* T! P
13. s =& @& g: d$ Y, \+ \
    
14.     6.126720600000000e+009' |2 c\" A# p' T1 L
    
15. . T2 F$ P4 a1 H$ i
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
' X9 y# Q  r7 l
  H9 e5 Y. }: x& K* w! e, T+ i7 Q* ^Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   , w2 H8 W+ m: m, `
2. mvar:7 t5 g$ n7 s9 F+ ]. `  _
   
3. t=sys::clock(),* {9 m: ]& {8 {* z
   
4. s=0,
   & S- ^1 ~/ ]) u: h' s
5. x=0, while{x<=1000,
   : {& t# J8 A\" Y! Q2 q. E+ z
6.   y=0, while{y<=100,
   ) I4 O# p. ^; `+ Y- D& a4 s\" K
7.     z=0, while{z<=100,
   . `+ E# O' J) N3 \5 j
8.       s=s+f(x,y,z),
   5 i' A- G8 U! a: n5 T3 `
9.       z++
   . D& t8 G* L  @7 T
10.     },
    5 T0 D( L\" O, e
11.     y++# @$ f6 {  b! Z4 o/ A
    
12.   },
    7 T/ c% \- i% j3 o0 o
13.   x++\" O/ r$ }# R9 s1 k6 \7 j8 k: N
    
14. },
    , [2 d* A: I8 s4 W4 `4 {
15. s;5 s% j0 t3 K  f( \- |3 S
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
7 B4 Q: p+ z6 C) D- N. d2 F1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
$ c, S0 P) l! o2 p' A* S

hxftntA

3 }' j" Q, |( r8 \/ ^
3 j0 @+ t, z! y



9 J( \9 }) B, m8 g) z+ [( ^  P
$ @- W7 ~' @2 K5 S



" g) A6 m8 Y( `& w% S7 v2 }
+ w( C6 d3 V. \

, u- U+ o7 _9 M4 x