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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
7 S# s3 o1 z0 T' F
; `8 `$ z  ^/ O; d- f* Y＝＝＝＝＝＝＝＝＝＝＝＝
8 f5 R* O* {* }3 E$ s, j. X
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   : b6 l1 V& @1 ?. E3 n3 I) m
2. tic;% W6 C* N, l! t/ S: \
   
3. s=0;) J/ d4 x% u. g5 N\" _: |
   
4. for x=0:10005 u- o% y8 a3 P* O
   
5.   for y=0:100
   # A2 L& Z5 h/ h( t
6.     for z=0:100
   7 I- ]( L3 X  `) \8 |
7.       s=s+f(x,y,z);6 O. ^7 G. W! d; I( H* ?
   
8.     end( P+ e8 s( i0 o  \/ p
   
9.   end
   2 ?( x0 T) L* L% m
10. end7 l) i( Q7 R- B, H  C
    
11. s
    0 A# T4 w) A6 Q6 I3 I# \0 J. m
12. toc( O$ F5 Q: ]& o- x+ U) M2 x
    
13. 8 ^' N4 Q- R- N: y* ~0 V
    
14. s =
    0 m& ]  n4 d0 e* F! w/ L
15.     6.126720600000000e+0094 s: X8 }' v' I0 O2 C
    
16. * ?* B4 H0 ?6 H# R2 p
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   ; `% |2 {, s3 v% n, |
2. function c=xyz(x,y,z)$ _( t, W8 |2 `8 a. ?\" S\" C
   
3. c=x+y+z;
   5 r2 Z' ?3 j5 t: D! d
4. end

复制代码

测试代码：

1. tic;
   ) z$ c' z, F0 ?1 [1 G
2. s=0;$ [1 l5 M* ]+ R. ^( s/ \7 r: q5 q6 M
   
3. for x=0:1000, _2 S) c\" o) d2 I! C  |) Z
   
4.   for y=0:100
   # K% ^$ }9 Y+ b- Q  c6 }  ^0 n
5.     for z=0:100* w5 ?7 j3 ~' j, R
   
6.       s=s+xyz(x,y,z);
   0 K8 F8 E2 ?' r% W* O
7.     end0 ]2 _9 N4 T7 ~
   
8.   end( u6 B  k  l( K  l
   
9. end# J\" E1 P* H! u  G6 U3 a
   
10. s- Q5 l1 ]4 F0 _5 M- L( J3 x9 J
    
11. toc) M2 e/ b* L, T# C6 z* y
    
12. ) I! Y* T3 V* m( X& {- y
    
13. s =
    ; o! m0 o: _& U8 `
14.     6.126720600000000e+009
    ! V: B! i$ U; ?1 w4 I5 ~2 [& B) U
15. 
    0 @  {  P7 b9 ^
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;; F8 |7 y' F' z
   
2. mvar:2 l$ M% Q/ e) f' ^* H% S; a/ i
   
3. t=sys::clock(),
   \" X$ {2 h+ m: ]* ?2 P) {
4. s=0,' @) h7 \# _  S/ j: K) M8 H
   
5. x=0, while{x<=1000,
   3 V( |! a( n) p. f) x) L
6.   y=0, while{y<=100,
   % Q2 Y1 J3 _4 h( q! z
7.     z=0, while{z<=100,
   8 r1 s2 d8 @! b# G. v+ ?8 V- {
8.       s=s+f(x,y,z),
   5 h7 Q  ~) q( o% C  Q% d
9.       z++0 P5 B) D8 H' l9 {$ m2 g7 U
   
10.     },; Z7 j3 R. L6 A2 g6 Q
    
11.     y++
    # N; n( U9 L! B7 l4 w3 {
12.   },! Q, `0 e, x0 W7 k2 z0 I
    
13.   x++
    ' P% e0 _: v! L; s5 Y  A) N, ?
14. },  e  n% D: f4 T0 q( S
    
15. s;
    . d% X) ^, `% K; F; M
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

6 l1 p, D9 }: j8 y4 d. |. g
$ u2 x, g* J; L; S+ n; R5 e& I% u- T
3 M8 Y, O$ h: ]7 t; ?; c



: X2 q+ n" z' C
/ l; [, e/ a3 B2 s2 @
. l! O# `- P3 A% g  \, f3 m% ?
. ]% ~$ j/ S' k; p. Y


+ p+ W; {) ^- v" T
0 i  g  A( l$ f- [1 {: k
* j) n1 ?: q  J* i
/ l# Q1 P# m6 H. g. ?