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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
- i( g  r5 j- c: f. T- S
＝＝＝＝＝＝＝＝＝＝＝＝

Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;
   \" W2 }9 C1 q+ W0 J& H# e2 Y
2. tic;! O9 w) L2 f6 W+ m6 B( c
   
3. s=0;
   4 I; p3 Z- u/ R2 g- }
4. for x=0:1000, b0 U+ R! E! h& t- t: z
   
5.   for y=0:100
   # ~# E4 d+ n8 J) P) ?  y8 a$ Y0 v
6.     for z=0:100
   5 E. R$ d( w\" |, }) X+ f# M
7.       s=s+f(x,y,z);
   5 H$ i- c5 z7 m$ I! ~4 \* ^
8.     end) k( ^8 d\" ~/ x6 A0 Q: w+ I/ Y
   
9.   end( z5 J+ h9 ~1 K, [$ g\" Y; [3 P
   
10. end2 m! {) B: F8 x9 _
    
11. s8 V, b6 G+ L\" T4 _/ ^2 ~! q
    
12. toc2 ^7 w& F- h2 h* \8 ]+ x4 s# F
    
13. 
    * I7 K! H/ D1 k/ Y
14. s =
    ' g6 e0 k+ I7 Q: C. B  `
15.     6.126720600000000e+009
    6 O: a3 V1 r3 O; ?0 `' p) P: _
16. 4 }9 k/ f. k% H' z; T. r. l3 ]
    
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
  g$ M; F7 _4 em函数：

1. %file xyz.m
   : M8 K6 O1 f9 R9 i/ [
2. function c=xyz(x,y,z)
   7 S8 z; t4 X5 s& v( V' [- Z
3. c=x+y+z;) H; F- `) k% X  a5 R8 V4 t
   
4. end

复制代码

测试代码：

1. tic;
   0 K& j- |5 @  R1 _. K4 q
2. s=0;
   ; {! p4 [9 S- N/ G
3. for x=0:1000
   ' ^6 t' M& E7 L7 W
4.   for y=0:1009 B3 v) J  A% L& |\" w* g1 w
   
5.     for z=0:100) `: i0 Q0 \8 ]6 D8 M
   
6.       s=s+xyz(x,y,z);7 S6 D) G- _# M  l  ]- |5 O. R
   
7.     end
   . j5 H# }2 ^5 p* [8 w6 I  v
8.   end
   9 R' r4 C$ v; f
9. end+ N  b7 ~% J& N$ t* K. k
   
10. s3 i\" D; a; Z: B$ v\" v; V* {
    
11. toc
    ( V9 u5 m3 U- e
12. 
    1 V$ D! n0 q2 K$ e/ B
13. s =
    . X; l( k5 Q: r: k: H
14.     6.126720600000000e+009
    * k2 e3 O7 d2 O3 ?0 \5 _
15. 3 y' t. D/ ?/ T% V9 O. ?& ^$ |
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
1 V  R/ D8 f, S) ?8 w
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;
   ( F8 `/ W' Z( F
2. mvar:
   + W( i3 n; Y( d0 [% V4 n
3. t=sys::clock(),
   % r) W0 ~! a8 g
4. s=0,& B6 A9 S/ Q% N& k( ?. p
   
5. x=0, while{x<=1000,
   & C: ~& P& r; @( o, m
6.   y=0, while{y<=100,
   ' H/ n- U\" O: S$ ~5 [& C: T8 }
7.     z=0, while{z<=100,
   ! G2 B\" c2 K! C: s) l
8.       s=s+f(x,y,z),+ \/ J* N' e2 [% E' C$ L* T& f
   
9.       z++
   . t3 I+ h5 q% Z/ I: H# P\" O
10.     },
    % \8 ^8 b2 B( \$ D- Z6 g: N
11.     y++
    + \2 ?, c2 k+ w3 _! C) E0 I1 _. q
12.   },
    # x( l0 u4 Z; o- O
13.   x++
    * a, f) H\" y8 b1 \
14. },
    4 R3 W, z- a) x% q
15. s;$ k% D6 O! A4 p
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
( G0 z5 Z- |+ b0 o+ l6 Z: n1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
6 J0 ]- y: e7 r. W9 B

hxftntA

! J( U4 L9 C. H5 P- ^: a( ~- I


* Y  `) ~8 w' N7 {# K


3 O; p- c, \& F% O



2 I1 Y6 l1 }! r! [+ _/ o0 U, l6 e



: L8 C0 _- M" e$ _0 @& }& [