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

forcal

数值计算，例如微积分、微分方程求解等存在大量的函数调用，故高效的函数调用对数值计算软件来讲是至关重要的。
* A  C9 F8 E; A6 d8 g
＝＝＝＝＝＝＝＝＝＝＝＝
2 k1 ?( P4 R8 B- ^- B' R4 J1 A
Matlab 2009a的测试代码：

1. f=@(x,y,z,t)x+y+z;8 ?# Q2 t4 L4 K' E
   
2. tic;/ T+ A) J8 e; E+ o5 x
   
3. s=0;
   5 A& S0 W2 Q3 _4 p' @
4. for x=0:1000
   8 P: `8 Y1 b9 |9 w; y, x
5.   for y=0:100* x2 j5 ]! j4 `% x+ C' o
   
6.     for z=0:100
   3 @& b2 U' K7 o( C) R
7.       s=s+f(x,y,z);$ V, N$ R5 n\" ~
   
8.     end
   \" M0 r+ A7 h; {: r4 B
9.   end
   \" B. u# v- h! @( R/ @% H, b$ n
10. end
    $ x' e6 ?6 z7 r3 W8 D. y
11. s  S: ?0 X- y7 l' W* o
    
12. toc, c4 t1 Y' c% t; O/ T  `
    
13.   X9 D, U: [/ l( h1 B, B: H
    
14. s =
    0 }* Y1 t) v5 w+ Q3 P
15.     6.126720600000000e+009
      o: z/ d* w4 k7 x/ i
16. 
    : E7 b# k8 m* Z$ h2 i\" n0 K
17. Elapsed time is 9.546717 seconds.

复制代码

发现将函数写成m文件后效率会提高，如下例：
# j9 A* g) J& ~m函数：

1. %file xyz.m2 \2 r7 K8 s\" G
   
2. function c=xyz(x,y,z)
   3 z& ]3 o6 V- C( a) Q; p+ E! P
3. c=x+y+z;
   2 z, W2 f1 Z! P5 ~7 H
4. end

复制代码

测试代码：

1. tic;
   , x1 t: H\" t  b' y( j
2. s=0;
   7 X- ?' f5 ?& Z\" P
3. for x=0:1000
   0 M( X% ^, d  k7 a, H) T
4.   for y=0:1009 a# z4 a+ L\" ^% m+ |
   
5.     for z=0:1006 W& U& J4 ~3 r+ F1 c) `! M; p' T( n! }
   
6.       s=s+xyz(x,y,z);
   . N$ X- M\" v( K4 a2 a+ d% P
7.     end
   \" }/ V0 `' i. W, h! I* j6 H
8.   end
     [, H9 g& Y3 f
9. end
   $ h# W8 x& Q# o8 g\" Y
10. s
    % a7 }6 d6 V8 m/ R0 K$ X
11. toc7 o' M2 h- d( @4 ~- S' i6 N6 W
    
12. 3 A5 H, |8 K\" T8 B* b
    
13. s =9 I( w9 u( Y: `9 o
    
14.     6.126720600000000e+0090 ]- u5 E% }. e) v' h: M
    
15. * K0 J4 g5 W8 l/ ~
    
16. Elapsed time is 4.724592 seconds.

复制代码

＝＝＝＝＝＝＝＝＝＝
  Q8 {6 c( F) f5 M: o' C, E) F
Forcal 9（OpenFC演示）代码：

1. f(x,y,z)=x+y+z;\" ^% f& y, J; Q; N8 z  T2 [
   
2. mvar:
   - ^! y& _\" I1 ~- i4 V
3. t=sys::clock(),
   + i. m1 c- e6 C+ I8 q1 }7 t
4. s=0,
   9 [6 W- y) a, T
5. x=0, while{x<=1000,
   9 ~6 C0 z& Q6 G' }
6.   y=0, while{y<=100,
   * O7 v- Q1 j8 i. X' o
7.     z=0, while{z<=100,; I5 }1 o! @, Y( }) I
   
8.       s=s+f(x,y,z),
   7 W1 R8 }\" D+ j# Q0 V4 F  v
9.       z++  q+ s# W/ \( _; P
   
10.     },8 k. O+ L3 n* a9 S8 B$ B+ r) D
    
11.     y++
      U/ r& D% s4 f+ A
12.   },
    ( C- e6 x. B, ]1 O$ ?
13.   x++6 ^- F6 Z) a- v- j, s
    
14. },7 B0 G9 J\" R: v2 L6 R+ T1 M8 l
    
15. s;* L$ n  W5 }2 [! c
    
16. [sys::clock()-t]/1000;

复制代码

结果：
6126720600.
: i; [0 Y" g$ a& q+ Q3 Q0 J6 D6 O$ j& V1.109  秒

二者效率有几倍的差距。

ehi28

嗯,不错,支持一下.
+ q0 N" y+ p. p$ `

hxftntA

, G- @$ c6 _0 K6 }: m





9 R/ @, }7 b" x; G+ Q$ }7 W! P
" n  y# S* H: A
. ?4 K8 e! a* q! j( A





( S: {+ T3 |+ d7 F

; K/ Z1 C) i) k2 f