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

forcal

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

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

# Y# o5 s, L* QMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   0 C- T: D( l& z  C% f
2. tic;
   $ M. O& d7 O7 m3 c( r
3. s=0;  m\" V/ R. T& O
   
4. for x=0:1000' I3 b\" o/ o+ F% ?  U/ B
   
5.   for y=0:100
   9 I' X8 R+ H* l  S1 M. ~
6.     for z=0:100( n2 S5 F4 N5 y# {9 H
   
7.       s=s+f(x,y,z);
   0 y( J4 P' d$ F% V* }& Q
8.     end
   ! p: F- F: E% t- w) [* Y
9.   end; C# }1 M0 U  }% a+ d
   
10. end1 A! e! R' t0 u/ n; \
    
11. s) u, j, a$ B9 I+ R* C
    
12. toc
    , @1 G: Y9 Z& T/ l
13. 4 _6 z+ C. A) S1 d9 U! r5 g
    
14. s =
    : w0 d- q5 z) ^
15.     6.126720600000000e+0093 V3 }; x  U( I& B6 F
    
16. 
    3 a3 ?; M\" v8 U
17. Elapsed time is 9.546717 seconds.

复制代码

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

1. %file xyz.m
   5 n5 B: F( C% D* l# c
2. function c=xyz(x,y,z)% o- H$ }$ }5 C+ T5 s
   
3. c=x+y+z;8 e) a: F5 M) i9 }! Z5 i  K$ i  ^
   
4. end

复制代码

测试代码：

1. tic;# z( z# p4 r4 Y7 P: j
   
2. s=0;
   ( X$ K\" ?- c0 ]+ t; v2 [
3. for x=0:10007 h$ q\" M- Z\" I% w4 g
   
4.   for y=0:100) f4 v: a7 o' W; T0 K  }) v
   
5.     for z=0:100
   1 x) _1 F* N( L. ]1 {6 _4 y
6.       s=s+xyz(x,y,z);
   4 w6 u- D) d' c; q4 V5 K7 R% v. M
7.     end
   . K; I, H$ H; h
8.   end
   , I1 w/ K; V9 L6 J& E1 L7 X7 E
9. end
   / g# k' D' a2 \
10. s\" ?. n\" T. H* U4 W( G
    
11. toc
    5 K8 m( X0 t  L2 I
12. 4 a5 \8 p* u- b/ c; R
    
13. s =
    $ ?8 }( `. V/ e: b
14.     6.126720600000000e+0096 U' b# C2 v% `% M\" Y) y# ~, t
    
15. 
    * A# C! A8 S( q( N! }7 [9 h
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

3 Q- z$ K6 e% [+ N! C6 y# ?, l7 j9 JForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ( u4 F( w& z8 c0 I' d3 }' O/ F. H
2. mvar:5 c! P( N\" C* S\" A% H; c) X
   
3. t=sys::clock(),: j* i, F) L: b1 f8 o2 z: g
   
4. s=0,; ~# @' w- v' F: W% P6 t
   
5. x=0, while{x<=1000,8 f( e! A, F1 m, u6 u
   
6.   y=0, while{y<=100,
   5 T0 W( H$ t3 f  |
7.     z=0, while{z<=100,3 R' B+ m. F* U
   
8.       s=s+f(x,y,z),* N/ V% T/ I3 G\" f* Y) ^& |, }
   
9.       z++
   \" j7 y! A% h+ l5 I5 \
10.     },4 r3 _! e& ^& y' R5 H) Z
    
11.     y++2 E4 l8 v# ~  h8 }
    
12.   },2 j: N! [* G5 ?4 k( ?1 \2 b& k
    
13.   x++- C! b( q% N9 @
    
14. },8 N) s1 E; y\" O2 C+ ^/ j3 \4 ?; [
    
15. s;, N# M  }* I8 p+ p% L: C8 w
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

& {# q) \7 U% V8 g& A0 Q* h





( [, _% a5 o; o0 F
; Y& E+ \( N8 w9 U1 T

6 }$ J5 J! l( U; s


2 w5 Y3 ]9 b* O8 e

) P6 S! U$ e; B& r  O