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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
" W& X" ~! N, g) |8 L) S
1 o& f; M$ |3 F( _+ o( O/ T& l＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   # |6 K) j' O! H/ P, M8 z
2. tic;
   $ s2 r& U! K8 y1 Q- {$ |
3. s=0;
   ( A1 w6 c; ?4 [' f  d4 U
4. for x=0:1000: E1 j+ e: d! h/ I2 p  i. j7 f
   
5.   for y=0:1004 b3 e. E& v; _1 ^& C( E
   
6.     for z=0:100; o% T2 ]6 Z5 q# {( P. B! @
   
7.       s=s+f(x,y,z);9 D5 n3 h1 U+ a5 I/ K& `1 R
   
8.     end
   ! ]' E& d' s$ O+ w\" h0 r4 x' ~. g& @
9.   end- F, p3 E% Z% f8 \
   
10. end
    1 a9 g4 ~1 Q+ P+ |& O4 T
11. s9 m; k* [; I; m: v+ V4 x+ S% I
    
12. toc3 P9 w/ ^# _9 V& S' \
    
13. * c$ G* w3 ~\" [: J6 j
    
14. s =
    8 `2 Q% }/ ]# ?\" K* `# B  w
15.     6.126720600000000e+009
    3 B, e6 S\" H7 y; z  P
16. ! w  G$ s  H- S
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   9 i& z$ w! K# t4 Q
2. function c=xyz(x,y,z)
   5 a% X* ^\" _/ h0 [2 b
3. c=x+y+z;8 X. B% [0 M+ k* B! W4 x2 _1 A
   
4. end

复制代码

测试代码：

1. tic;8 D. b2 P/ M6 L
   
2. s=0;2 T! L# H\" Z6 F\" ~1 y/ x
   
3. for x=0:1000* L# K/ k1 O/ {\" I9 \8 z
   
4.   for y=0:100( Y% R+ Y/ O: O+ G) E1 n' s\" H
   
5.     for z=0:100
   7 G  `+ l9 b% d* d
6.       s=s+xyz(x,y,z);$ M' O- f  i. t5 j; \
   
7.     end
   9 K! o2 M1 Z\" o* M, n; v
8.   end
   ) y8 }+ c0 ?$ k% D& Y8 i
9. end. ~) P1 R9 Z; g! Q
   
10. s( M# }' N, \* \$ [
    
11. toc- m# d) ]) y6 k5 [
    
12. 
    3 Y' c# M3 k! Q8 Y' A; T2 B9 d- V( L
13. s =
    8 |) F' Y& X2 h  Z
14.     6.126720600000000e+009
    3 K, c; D( t, @# h
15. 
    / D, I  T. _/ [3 P- o) y# P\" R% i% @
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
! `2 U5 g7 R8 o2 V* U
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   - S\" ?) R( I* I5 c' d: \
2. mvar:
   2 n0 e) A7 ~& I: I5 \/ E! _$ E1 l( c
3. t=sys::clock(),
   ' L) J$ `7 Q4 p% j; Q
4. s=0,  W/ m4 t! l- ]0 R
   
5. x=0, while{x<=1000,
   6 w* [4 M% ]. N6 a0 p3 N: x, I
6.   y=0, while{y<=100,. n% p, N& b; P0 `2 t5 B
   
7.     z=0, while{z<=100,
   % y0 D6 J! u( e8 Y* y\" w\" X2 n5 H
8.       s=s+f(x,y,z),
   7 H9 Q; ]7 Z2 ?! J( m( e6 @$ ~
9.       z++9 B  D$ I% E8 n/ E; U. a
   
10.     },, e9 F; d* V1 D/ e3 [\" [6 s
    
11.     y++& @# `+ _\" V/ L& w- s8 t6 c* A
    
12.   },2 N5 r# c\" Q3 e
    
13.   x++
    + e& i6 c& s+ j. d$ T, Z# o% r
14. },5 K' E4 r  i# ~2 |7 n; E( Y
    
15. s;
    / w4 p; B) K4 {% x
16. [sys::clock()-t]/1000;

复制代码

结果：
% [' B% a% R: n1 b0 g, n6126720600.
; n% _, U; ?1 j1 J! ?) o1.109  秒
! G2 F6 r+ x* L
5 P7 R8 F' T0 o二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

9 D/ d2 f& z' z, R/ {( ^
0 N& n( k% v! g- q





. ^  [% ?+ U% u8 L) u% r( {
& z1 f' @* ]2 @2 p
$ q/ ?1 z4 X' U2 X- q0 N
$ [/ g% `! h/ G& q  J; K) o
3 L$ q: H* u) U) [, Y, z. N

  C; S3 v# f! d5 ~( b
! s5 `2 B# i- F( g
( |! v/ D- y4 K6 Y2 q6 k+ R" \- `
0 C# i, U) R6 M) `