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

forcal

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

: \2 d+ L3 {0 t" i7 p; m4 q7 s＝＝＝＝＝＝＝＝＝＝＝＝
; G& P8 a/ ?! `+ N  U9 n% N
# P, y( `+ K8 n! V4 G, PMatlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;/ P# G+ b* G6 ^
   
2. tic;
   % ?/ u: }4 @1 F% k2 P7 ~: [
3. s=0;
   ! t( ]2 T, X( g. e7 g' q2 v
4. for x=0:1000
   9 {$ ?* f3 O: u9 a, \) i
5.   for y=0:100, S( u& p) g% @9 o5 X
   
6.     for z=0:100  J1 `$ u/ F+ C  ~
   
7.       s=s+f(x,y,z);
   1 f- J# Z4 i1 B( o! q; \( o2 C
8.     end
   + ]! B5 P. ~* N' M( N% |
9.   end- ?! a; R3 Y6 s2 `  s
   
10. end
    5 k: `* C# \0 L% }! N
11. s- l# H. h. F; \5 n9 Y
    
12. toc1 M4 G  ]( h\" C% ~
    
13. 4 b: ~! B+ n3 T) Z\" \: L
    
14. s =
    , ~1 e5 h3 q. F9 `$ G* E  e
15.     6.126720600000000e+009$ f* j* }8 |2 J
    
16. ) ]9 f/ M9 R; z& T/ r# R! _. g
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
0 q# N( h; \) f: A4 G+ {3 um函数：

1. %file xyz.m
   - x( E$ q; o+ s3 @( a4 |+ W
2. function c=xyz(x,y,z)7 N5 _. d; \% t
   
3. c=x+y+z;; n: J* q5 R' i8 \
   
4. end

复制代码

测试代码：

1. tic;! ?5 k7 L1 G0 \
   
2. s=0;\" y! ]% ]; g& U* B8 S6 Y% E
   
3. for x=0:10008 s/ O' X6 b& r( s, L8 x
   
4.   for y=0:100
   $ L2 j2 Y  u; `. U4 A, q4 F
5.     for z=0:100
   8 `3 S/ v7 P% |- ]# Y# X\" q
6.       s=s+xyz(x,y,z);
   \" x( Q5 I6 C) M: Q5 C) }
7.     end  s; x6 y) d& |2 M8 n
   
8.   end% G\" }- q0 `1 E3 X; v( C) |/ F1 r* P
   
9. end- A\" g( t- M/ h8 c. y6 x
   
10. s
    ( _- f; F6 \% f3 P% i% I
11. toc- Y2 a8 q, ^- t1 ?( O
    
12. 
    * g4 o+ W& @\" t1 @  B
13. s =
    7 i' T& {* e. u5 H% v) B- b6 K% T
14.     6.126720600000000e+009
    2 k5 [5 U8 J/ {6 ?0 \( O. @
15. 
    & n! A# p. d' }  N& p
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝

; q9 z# S/ I# q$ j$ aForcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;. K- I( l3 I, I4 J; e* I! H
   
2. mvar:
   2 ^7 W, T. @0 v) n* B
3. t=sys::clock(),
   8 [9 O  w, i$ j
4. s=0,9 e' N2 Q& R+ `4 a
   
5. x=0, while{x<=1000,
   2 O, \! ^; U+ S* S2 |- w
6.   y=0, while{y<=100,
     c2 m  E# O2 l$ y2 V7 U
7.     z=0, while{z<=100,
   0 @. i5 F- b: M  i8 |
8.       s=s+f(x,y,z),
   ( a* V4 M\" U* Q, }6 O! a7 i- H
9.       z++
   + A) M( Q8 F\" ?0 L8 n5 g! l2 G
10.     },2 z/ g\" y8 n. T1 n/ s
    
11.     y++
    5 Q& V; E) ^1 s9 U$ u3 h
12.   },3 C; w3 F2 d7 u; l% p' E) F3 R3 @
    
13.   x++
    ! ~, E' @' t4 s\" E
14. },
    # R2 \4 d3 y: o: I9 @
15. s;
    5 k, {( i. `\" T9 B: z' M+ y% J
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
1.109  秒

: q% s# \. c8 J% [4 W二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.

hxftntA

; k/ {- U! d# U4 ^+ @
4 Q9 v, L6 a$ u- W( l5 l& q0 v$ m
& I" l) {5 K/ |" ~* g( j
9 v% _. Q: g% t+ ]6 c4 y- E
- C7 y( z; V* d0 N2 v- a' S& g2 e
+ l/ x' F; {! P" d! M/ _- D% i  K


( P( H1 _! _7 j, l: u4 E* l
% w9 C* ]+ e5 o5 F
9 _( m8 b5 U. X4 m5 n* A' _& l
4 s. K+ [$ @7 A7 {' K  I1 A4 [
7 x2 V$ V1 J9 ^3 u9 d' w


9 J) C% s) Z1 e) D  T