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

forcal

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

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

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   & J6 e4 G8 t\" }+ e8 I5 x9 J3 s* [8 h* p
2. tic;6 o\" D7 o; y! X) Y# X
   
3. s=0;6 _( c) I7 w* l4 E8 d: y- ]
   
4. for x=0:1000! Q0 X) r6 |7 r
   
5.   for y=0:100
   1 X0 \5 D7 c; h2 E+ j. n
6.     for z=0:100
   ! N! E8 I. Y6 I1 ], Z
7.       s=s+f(x,y,z);( T: n' D$ m( J/ ], b0 N\" X
   
8.     end4 O1 L8 p7 y1 r' D) @) _
   
9.   end( b, y. x9 ?9 B6 G\" s, ~- `
   
10. end
    9 o# l8 \+ r/ m, W* S
11. s
    5 O7 }  I$ I' w) S
12. toc
    $ J' \/ n* @/ N5 |3 e4 c0 K
13. 
      Z4 k3 N1 X\" k3 E  ^7 f3 O
14. s =( _) g. p  N3 S/ L9 \
    
15.     6.126720600000000e+0098 w) k0 f. [9 w2 A( K3 }) I
    
16. 
    2 b  s' v; P, S4 |, H
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   ' ^3 }+ }) k' i% E7 g5 d( r5 L
2. function c=xyz(x,y,z)
   & j4 D2 E; K% W7 E0 t8 W
3. c=x+y+z;
   9 ?7 ?; \0 r2 F, w% Z4 ^3 G
4. end

复制代码

测试代码：

1. tic;
   , n% K. g+ \) M2 B. Q7 D
2. s=0;
   ( |3 |* x# |; w. q9 V: j
3. for x=0:1000
     p2 h8 z, {6 S
4.   for y=0:1004 x\" j+ p1 A! F
   
5.     for z=0:1008 ^, ?7 Y4 g& K5 g4 Y2 K* t& b
   
6.       s=s+xyz(x,y,z);
   8 N  F' v# c5 J! w& p; R' W+ C
7.     end
     Z$ W8 r: z4 k5 K/ _7 N; g8 S3 }5 y
8.   end9 t$ {( u1 Y1 \7 |& S
   
9. end: L' ^/ A6 V: u( G9 z6 t' g
   
10. s' e+ R! H8 m' \4 K) u) ?5 j\" G
    
11. toc
    1 G* W! E: N5 c8 c. {5 w
12. 
      ^+ _$ g+ K# n7 Z\" r
13. s =5 r8 N8 [/ T: K, O
    
14.     6.126720600000000e+009
    # ?6 ^$ ~. w' j  m& I
15. ! g6 `; U  B7 X: s
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
) o( k: Z. p- O. J: G/ f6 k
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;) a# v9 X, A8 n) v: T
   
2. mvar:, N( H7 N+ R, F5 }2 w
   
3. t=sys::clock(),
   8 ]0 ^8 Q: ^, m
4. s=0,
   ( U3 D2 \/ f% \- n7 Q
5. x=0, while{x<=1000,
   \" [# E. K1 v0 [( k\" u
6.   y=0, while{y<=100,) [4 C/ C+ J4 k! m0 `3 K* @* S
   
7.     z=0, while{z<=100,
   7 J5 ~: v* }8 y. @
8.       s=s+f(x,y,z),
     F& p3 ]6 M- q\" b4 Z7 l3 R
9.       z++
   ( S; r8 ]' f! l, ~  l
10.     },
    2 ~9 n9 S  d) |& C. w
11.     y++' q1 a- _2 ?/ \
    
12.   },
    ; e/ G' g4 [- `$ ~% h
13.   x++
    7 \8 F. P7 Y3 w. S/ d1 w/ p' ~
14. },
    ! a$ f9 K) S$ Q. `4 w6 d. T3 K0 Y
15. s;
    , _$ O) C# X\" X: f6 g' X
16. [sys::clock()-t]/1000;

复制代码

结果：
5 U8 g9 h1 P+ Q7 w% `! k5 o6126720600.
+ t* [3 u& q% J5 ^4 O6 z1.109  秒
7 y& i8 [" J8 h+ f# y% _
/ c, T& N* F  p9 p二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
* R: I. }8 E& D' e4 z  t- U1 N" d, z

hxftntA

1 \7 P# u1 p4 s& U! ~  }9 z& t+ V7 N
1 q. q  r9 F! Y' u4 o! Y
+ [& A' }' E6 I. ?: k
7 W1 {6 N$ z  W& f
' y- ~/ H/ |9 ^8 ^; D; M5 a( T2 I% A

& J3 ^/ ?8 j9 K: u+ b
8 @& [+ H) ~5 T* ]* ?4 V! o

( i& K) w4 }$ v$ f" L' Q+ @, W

6 h5 O. ]) M( _7 B9 j
+ j+ k* n* J# S5 E8 A7 y
7 D: Q' K+ |( F0 _" `0 H8 v

8 c) c! z3 n$ _' J
1 X$ I# o* f' m1 N3 Y% a: U