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

forcal

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

＝＝＝＝＝＝＝＝＝＝＝＝

' u; ?& A# J/ ZMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;4 ~7 E; l6 c& S+ I5 I7 a6 m5 Q1 E
   
2. tic;
   : n. u' C7 `, S0 O9 g% |' G# s
3. s=0;3 P, r! M6 }# H3 C, F/ F% p
   
4. for x=0:10005 ], E* @/ f/ X9 b, r
   
5.   for y=0:100% ?) s6 Y\" t- K/ V, L
   
6.     for z=0:100
   2 ?& z4 S) u1 r- c* S
7.       s=s+f(x,y,z);\" y\" X3 N: M5 g) B
   
8.     end
   + a9 A8 Z3 n9 D5 d3 x$ S
9.   end
   2 _0 H7 p2 W8 e3 p+ L
10. end
    \" ?( q4 m. C0 K  u
11. s3 m& b% F3 n0 G1 y# @/ q
    
12. toc
    \" H3 t$ E7 K, r+ n
13. 
    # E. G\" B  y- Q$ n
14. s =, g8 @' P4 O4 W
    
15.     6.126720600000000e+009
    2 R, U. c2 p5 r* j. V; v9 I# d
16. 0 i; i! {6 \3 x* ]7 |* V9 u$ N
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   \" q7 P# i2 t% i9 D, [
2. function c=xyz(x,y,z): J/ b  t$ x, C0 k8 n3 E: b
   
3. c=x+y+z;. s$ _  M. U; z
   
4. end

复制代码

测试代码：

1. tic;6 d2 N4 A& h4 R
   
2. s=0;; i! x2 N) |4 p  t
   
3. for x=0:1000
   8 n7 i; b( t. W- n* d
4.   for y=0:1002 r( t/ @& Z# n7 _- X, O: ?
   
5.     for z=0:100- X' ?6 y) i* L% d
   
6.       s=s+xyz(x,y,z);
   $ j8 Y/ N0 T. Z6 \) a
7.     end
   2 v' U! L5 c  e
8.   end6 }4 C7 U6 ?% M& V( @& j1 F
   
9. end
   / M0 a3 n6 h8 i! q
10. s\" ]4 t6 Z: v) e
    
11. toc
    + @: q% J$ `1 ^% H. G: P) Z
12. 
    4 m\" [$ E& @$ R, C
13. s =7 p) B; ~. `3 [' ?/ t; i/ s2 n5 K
    
14.     6.126720600000000e+009, r( V$ _% M( X* {/ O. e! f+ U/ d
    
15. : G- i9 Q& Z+ O1 W8 z: U8 t
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
. T" e. t. Y4 L/ |; d- L
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ! N1 S. W9 L7 u2 X\" i9 {! X; v
2. mvar:/ l0 R! u- K/ N
   
3. t=sys::clock(),. q; J2 s( }; d& [& l3 N
   
4. s=0,. I4 [; V4 C* t+ k1 ~
   
5. x=0, while{x<=1000,
   3 F- |( K1 D, C# H2 Q
6.   y=0, while{y<=100,
   2 \* c  l% v5 x, i, Q& N  b
7.     z=0, while{z<=100,
   3 V9 p\" k$ }1 G6 P
8.       s=s+f(x,y,z),
     C8 q% G4 |1 P
9.       z++
   $ j) g0 R) K\" A# b1 c8 F
10.     },
    ' e& Z. V0 `3 q% a6 [
11.     y++& \0 Z% G* m0 e( i( D& J& j2 Z
    
12.   },- k# l5 A. ]4 l( G
    
13.   x++
    4 g; ~, z  s+ I
14. },- R3 U+ Q* L% a/ i2 c
    
15. s;
    $ O; O; }, x) R8 Y5 n2 ]4 F9 A
16. [sys::clock()-t]/1000;

复制代码

结果：
7 B# H+ G  `# H' J6 x$ H! _& S6126720600.
; q2 q- ?7 d. j; \2 t2 I. b* `! M1.109  秒
5 b* h) F$ {( b( U$ j
8 c/ W5 R/ v- U% |6 ?二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

/ A2 a: G+ v1 L4 k5 X! E! G

. K" T: U( b; u3 w) s& Z

# L2 J# b4 H! c! Z

/ M* U5 [# @# r- b# t- p
/ |/ i2 t3 Q3 A7 @8 R8 l, a


1 z2 q  p7 t) g* ~8 v5 E9 b# E
" r6 q  M/ T1 C