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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
/ {1 C$ W; v2 @/ l/ ]
% R2 a. v8 [6 \7 g) \) M6 H3 J＝＝＝＝＝＝＝＝＝＝＝＝
( t! \4 c, p2 a" U) d1 R
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   6 X5 \1 z5 Y, e$ K8 ]
2. tic;8 H) f- R& j( n* Y
   
3. s=0;' j  @) [1 |4 ^3 w8 |
   
4. for x=0:1000; v! s8 c\" K9 B+ M5 H2 W
   
5.   for y=0:100& e2 r0 ]& z$ ^9 z\" F
   
6.     for z=0:100
   9 V, L3 ~5 F1 ^) h
7.       s=s+f(x,y,z);
   2 G% l) q: s, `: g: _9 c) y  t& K
8.     end
   \" U# z# D. ~8 u; z1 }$ w
9.   end7 E) _3 s5 A\" a% ^1 B6 J
   
10. end
    % x$ V0 b' X  e
11. s6 u: O) f$ S( v. F
    
12. toc
    \" P- h) H6 ^7 D) o/ }+ w+ W& x
13. 
    8 y! v( A- F3 h/ w
14. s =2 T# U* r4 a, N0 Y$ S. n
    
15.     6.126720600000000e+0099 M. J3 L1 c3 E5 n! P3 V2 a
    
16. + y' S$ q8 [5 t
    
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m4 ~+ G9 U# K0 V# U4 \
   
2. function c=xyz(x,y,z)
   $ W) a2 G9 N( J' k, R
3. c=x+y+z;
   2 M& N) @& _; b, X5 H9 e0 B
4. end

复制代码

测试代码：

1. tic;5 b) A# ]: D2 Y+ |, M
   
2. s=0;
   4 z( L4 t' A/ G. u3 e\" a4 Z  L
3. for x=0:10003 N! I1 z9 @7 f. a- j( ?
   
4.   for y=0:100
   3 {; @2 h7 q: O) c
5.     for z=0:100
   1 d9 ^  @: z$ v2 t3 S
6.       s=s+xyz(x,y,z);
     B7 B5 I/ U4 S
7.     end4 F7 ?# _. t9 A1 r\" d8 ]/ |
   
8.   end
   % B* z8 W! e$ [- R( ^$ ]
9. end& K0 x\" t  m7 d) t9 m
   
10. s+ N; M4 ]+ g2 C( Z3 w
    
11. toc4 `& s9 y% @/ Q; a
    
12. \" B5 j) T) w2 i* x
    
13. s =! |9 e  g' u1 x* R( u4 x
    
14.     6.126720600000000e+0094 r) V. U. m( M# C8 G) m
    
15. ' V5 `5 m1 o. k* q# I
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

9 @! _/ s9 _9 H8 \9 hForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;6 T/ k! \# B$ o1 y
   
2. mvar:4 V; t; z( t# }: I2 O
   
3. t=sys::clock(),
   $ |% o1 D3 A% ~
4. s=0,8 V0 Z! l$ d; ~8 r; |
   
5. x=0, while{x<=1000,
   ) _* l1 }. N$ X\" ~) x. {7 Q/ U4 d  j
6.   y=0, while{y<=100,
   \" q& J, ]* a3 A5 c2 J; @, [
7.     z=0, while{z<=100,
   9 _+ l, B- g1 b: p
8.       s=s+f(x,y,z),: x% }# c4 ?  _/ {  e
   
9.       z++
   0 _  G\" a\" k& O3 L9 R( v
10.     },4 o9 I4 J7 a& R\" j
    
11.     y++
    2 }6 n0 @, ~  b! ^) W' c
12.   },* b1 ~0 E8 f' }: @; \
    
13.   x++
    5 t+ }& ^+ i* f9 s5 c( s+ B
14. },# M' e\" L( P0 N% |
    
15. s;+ _: A\" B9 \+ j! W0 n' v+ O
    
16. [sys::clock()-t]/1000;

复制代码

结果：
9 [0 k- \/ h7 [( {: P0 W6126720600.
" J/ g7 }& j) M1 t: ^+ x4 C- H3 ~1.109  秒
# m2 Y! A, x, @: T
/ n6 a+ @# b' ]; G二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
- r5 ]* l+ J' a9 ^( f; ~

hxftntA

9 O5 P3 a) z9 f( u- v

0 L5 h- l1 `/ A* ~( ]

. d6 Q: U- \. p: w- c
+ @2 Q- d2 i+ y' U  {


9 F. `# R2 e5 \8 `0 C/ c! m

9 m# L# O( i* b7 g1 N  r
0 F2 A! l+ n2 e) j3 p' R% Y

) x6 s, V4 d' E