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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
& ~6 y  f/ d; ]3 d
＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;$ E9 u  j& s\" B
   
2. tic;3 G) [9 I. T& E) `
   
3. s=0;
   2 i1 F4 x+ l, a
4. for x=0:10005 `, Z* a, w  \! O
   
5.   for y=0:100
   \" ~. n' ?3 D) j6 ^5 c8 [
6.     for z=0:1007 n  j6 B/ J- m0 C8 D
   
7.       s=s+f(x,y,z);6 ^6 S5 D% g) h2 F7 T# c# N
   
8.     end
   - v! E) y5 g* r  g' b
9.   end( X& z+ M, k, u6 P$ P$ {\" a: v! n
   
10. end
    : R  B1 D1 h* v2 m$ l& o) j
11. s$ ^2 A3 W\" M. l% q
    
12. toc\" Y) m) L' Y) X8 B
    
13. 
    ! _5 ?- ^$ A\" C4 H: k8 L+ X. w
14. s =
    + a8 x2 Z! W- u/ h: {# _
15.     6.126720600000000e+009( i% k; w, M+ A0 s3 `
    
16. # M8 Y3 [8 D. u5 b2 ~0 y3 {1 L
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
4 V( i8 K; W/ i( K) A0 u0 o; cm函数：

1. %file xyz.m8 g0 r. H6 s! [. y) H7 y
   
2. function c=xyz(x,y,z). P; {4 y$ z4 z* \( }
   
3. c=x+y+z;6 v% h4 w( V4 A/ c- H9 F8 H4 `- U1 c
   
4. end

复制代码

测试代码：

1. tic;
   % i( O3 J\" @' P5 ], H2 c
2. s=0;- r+ a3 W# ^5 N\" X! L! ~- n
   
3. for x=0:1000
   5 v+ M4 f' s! @' h. Z4 x: u
4.   for y=0:100
   ' o2 a- w- d5 ]6 _& G) X/ b) A# r
5.     for z=0:100
   1 }/ m+ Z- l6 n9 h
6.       s=s+xyz(x,y,z);
   \" a% J9 T; n( k; _\" v4 }$ l
7.     end$ v# [4 N# s6 f/ t, P, s, r0 r
   
8.   end. G6 P6 W\" z9 w: |5 h4 |5 F
   
9. end6 [- M, ^; {, i/ m8 V
   
10. s9 h) d  y; N+ Q0 U4 W' Y
    
11. toc
    * A$ B- f, ]- s5 K
12. 
    ; h7 B& c5 p8 Z# N
13. s =: X. P/ `+ N7 `6 ?! w( T; u
    
14.     6.126720600000000e+009$ @- F+ h( u* @: J7 j: `
    
15. / V1 J8 p2 O  c9 @/ o& Q' r
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   3 o/ t. \, b) }+ {( [
2. mvar:
     s7 W6 b1 Y6 n$ h$ k( o7 I
3. t=sys::clock(),
   9 Z9 B- ]3 E\" {1 d+ I/ `
4. s=0,
   ) Y3 a$ {4 D2 p( @( `. g: a
5. x=0, while{x<=1000,' Y$ G3 X, N7 a0 u+ C
   
6.   y=0, while{y<=100,
   8 ?4 k' m/ Y/ v# x, q$ v
7.     z=0, while{z<=100,1 B$ t, I# Z5 _* m
   
8.       s=s+f(x,y,z),% u) Z* o& g( U, v) b7 j1 L; r\" d
   
9.       z++3 E7 T5 f  w$ l\" A' k$ X: B
   
10.     },
    2 g- X; k$ S, ]
11.     y++& s* ^6 ~' o: F& V( Z$ o1 d
    
12.   },
    9 e1 Y, M$ e8 p8 K- u4 o6 m
13.   x++
    $ R\" A! _( Y# G
14. },. U# X0 F- F. Z6 E
    
15. s;
    8 k\" L% l# p: O% N$ y5 _\" U
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
+ D# b( O( f! i# d  |1.109  秒
3 L: l& N' D/ c# z( L* J
& P2 S5 ^# ]% t2 P* n2 [二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
0 _, t3 A* R0 A& A

hxftntA

. b& o6 ~9 o+ T% n! ?; i" y# T7 q

/ Y+ @' L/ H3 g
8 Y' P8 u6 S- h$ _: a4 }

9 f& ?. j4 `8 H: E/ F1 |

1 A/ k* p* K2 F! B! N
( B* c) W2 m8 Q' f: g8 A# d% ]6 [
' d6 W7 S& R9 d. j9 Y, s% ~
" F2 ?* m# ?- G7 ^* L$ x
7 p1 G. R, K7 |; F: i! ~



% u. L; u, p5 c  @; O  [