查看: 5411|回复: 3

极限测试之Matlab与Forcal有无JIT加速器的比较

[复制链接]

字体大小: 正常放大

forcal

45 主题	3 听众	282 积分

升级 91%

TA的每日心情

	难过 2012-8-27 18:22

签到天数: 1 天

[LV.1]初来乍到

电梯直达

1^#

发表于 2011-8-2 07:31 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

1、简单循环

Matlab 2009a在运行以下代码时使用了JIT加速器，获得了和C/C++近乎相同的运行速度：

clear all
) n1 q, r- V\" Q. l
tic* @0 x1 }+ E h2 e
s=0;3 I- V% c% \* S- Y2 D
for k=1:10000\" A! y. Y- ]( I7 H: y& U; t\" N% M2 K
for i=1:1000
( I! S\" p, }! v( ^8 P5 X. T) i
s=s+1;; Y4 g7 p9 f) a7 I! `& q- ~2 g
end( v; Q. k$ ~' z4 ?
end
0 q7 m H1 }& ~' ^4 ~/ ~# a
s
) f5 s. B; \\" S9 L1 M
toc1 U5 X0 \4 m+ a% C\" S: B, J
9 ~# A# Z7 x+ n0 B/ @
s =
. q\" N& O; l, H4 I3 {
2 m1 _: h+ C) N; D( H- K7 a, B, D
10000000- G- O$ N# q3 J3 r
7 i0 _% I$ q6 C* C
Elapsed time is 0.092808 seconds.

复制代码

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

Forcal没有JIT加速器，代码：

!using["sys"];6 m H, A+ O2 [
mvar:
: R* D$ j0 Y4 p9 { y' F, ~
t0=clock(),& H6 B% I9 p z5 y\\" g. _0 E
s=0,
+ D8 f0 Q+ O$ g! b7 w- [/ ]
k=0, while{k<10000,0 A, H# I0 ]& E2 X$ C9 s L
i=0, while{i<1000,s=s+1, i++},+ }( o( V# w1 j/ V6 w7 z
k++6 l7 ^. [$ j7 P4 f3 ]/ }, r
},& u6 i9 S! e! q: N
s;1 F3 ?8 A4 P7 Q! n- V! G8 J( I/ P8 Q% _
[clock()-t0]/1000;

结果：
10000000.
0.704 秒

二者运行效率相差7倍多。

2、循环中稍微加点东西

Matlab代码：

clear all
; {4 O# c/ f& \6 e# E, v+ \: W
tic
1 Q+ Y\" K\" L& _* v, F0 F- Q0 m
s=0;: C% l$ d6 Z' c. n- _
for k=1:10005 b& p) L3 |; C, [0 F
for i=1:1000; A/ I9 a, {# U+ w2 K; X
s=s+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i);9 G# v$ k% `* A# `/ q
end1 X$ S' U. n3 R, k' [3 ?7 r u
end
' Z: g* J5 u( C$ `) b2 z. ?
s$ r9 `8 [1 E& V% Q; @1 O E
toc
- D. S) m* u6 t, s* [
6 |+ y F) x1 J5 l4 k, V* B. [
s =
; J5 M, e$ t1 w7 Y* o
# c\" X\" I2 v) r$ k4 G) d; c
1.1009e+007
3 @, U( C, F5 U
4 L/ y/ T7 m+ V; o3 Q
Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

!using["sys"];$ E7 l6 x( v7 q
mvar:/ Z( U \- Z+ ?# b; Q
t0=clock(),
$ m7 m; B8 y1 i+ W6 y Q
s=0,
5 p& n% r% e* f+ C8 w
k=1, while{k<=1000,# Y\\" y; E' P. n e' Y
i=1, while{i<=1000, s=s+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i), i++},
5 b\\" M5 c H# T- g( K3 g
k++8 F/ y1 B1 d, O$ L9 s8 y: r
},
\\" z( l/ j# d\\" v) Y
s;$ Z! m2 d3 A4 Z Z& I
[clock()-t0]/1000;

结果：
11008953.66602346
1.031

二者效率相差已比较小。

3、再增加一个函数调用

Matlab代码：

clear all
' o, M! E( ?4 V F% m
f=@(x)x+1;8 X5 ]( M2 k# v9 M9 w$ S
tic
\" U& u* @/ [\" K; z2 m
s=0;\" K. s* a% {+ p' t8 H
for k=1:1000& ?0 r0 ?$ c0 {0 x0 B! {: G
for i=1:10002 H# H* x1 k! @0 P5 {$ `' `
s=s+f(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i);' B% p- l. R; Z# Y0 P+ h
end& g1 h J: D4 I& B* H\" r
end
J$ {' J% _* r6 G; w! U
s\" E6 t V& S+ p @, v* \
toc; @# A1 f& ]# H) S& {. X* ]
+ Q/ _5 `* j9 o3 A9 {* `! h
s =
* e m+ o2 y& h9 |5 f* m
3 k- g* p) b1 _ \: \
5.1251e+0081 A& J* ]! P9 K* \* g- [5 p
: h; K( W; ]* K\" K1 X2 _' d0 c. c
Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

!using["sys"];
2 ]- |% o* [( ^1 b5 S3 v$ V$ M
f(x)=x+1;9 ?/ O2 N' a& }: N
mvar:
1 [8 \; W\\" z! y* l+ y
t0=clock(),
; D# C# ?, Y; [\\" t\\" H/ ?( c
s=0,
1 W4 q/ N6 P) i, ?# ~! m
k=1, while{k<=1000,
m- f1 p* j9 U
i=1, while{i<=1000, s=s+f(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i)+1+sin(i), i++},8 H2 n, T# M! k. O
k++
2 d' h7 w1 W\\" N3 ]# l- a
},8 ~0 n0 n5 D6 h* o4 A+ t' a
s;/ Y0 [7 z4 m0 z& [* V7 Z6 E
[clock()-t0]/1000;

结果：
512508953.6658993
1.078

可以看出，Matlab效率已不及Forcal。
故JIT加速器虽好，但效果似乎并不明显。

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

海水

20 主题	4 听众	494 积分

升级 64.67%

TA的每日心情

	衰 2014-10-24 10:14

签到天数: 104 天

[LV.6]常住居民II

群组: Matlab讨论组

群组: 小草的客厅

群组: 2011建模讨论组

群组: 数学建模

群组: 数学建摸协会

2^#

发表于 2011-8-2 08:12 |只看该作者 |招呼Ta 关注Ta

我也说一句

发表

使用道具举报

forcal

45 主题	3 听众	282 积分

升级 91%

TA的每日心情

	难过 2012-8-27 18:22

签到天数: 1 天

[LV.1]初来乍到

3^#

发表于 2011-8-2 08:45 |只看该作者 |招呼Ta 关注Ta

4、在JIT加速器下，Matlab的数组存取效率非常高

Matlab代码：

>> clear all& L6 N. D: k& G. m' [: b: b' m) W* T
a = linspace(1,1,1000);
& Q3 K9 x$ v% y3 {
tic/ b3 u$ {# ~3 |7 ?
for k=1:10000+ W5 p4 B! w2 L; E/ J
for i=1:1000 Q- K2 G: B4 Q. I8 t/ O
a(i)=2;/ J; @; F* q+ c0 A( [+ F
end% V2 g/ a4 [ h* |' w) N( A& r6 u
end `2 D3 W- s3 ^0 p
toc
6 X' b. E% @2 m% L
5 K1 U0 B5 q3 T' M( S* q, l$ S
tic
8 h9 F( u: N1 L. d: F2 Q
s=0;
2 n1 i/ i: i\" h, X( y
for k=1:10000
+ [7 o5 z( t9 _. C; f3 x2 {& V Q
for i=1:1000; T/ l9 v7 m+ i x
s=s+a(i);
9 T8 n% [9 w b/ w. s. V! B% P0 e\" o
end
! S) l6 q# i3 p
end* ]6 H& _; c& L# b
s
1 O: ^\" I$ B* U: j8 Z8 ?
toc
, T\" I4 k9 t\" r ]+ J
Elapsed time is 0.140658 seconds.
' \+ G+ e1 w5 Y\" p8 j
! K3 A: l$ a# n Y0 c0 o2 c$ s- d+ f
s =8 V0 j& e( ?4 H! z( A! r( V/ |
7 Y9 D9 G- Y! K* ^
20000000
. w2 g6 i\" U) G\" K4 _- ^
7 T/ m0 r7 w# G# u% B
Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

!using["math","sys"];
5 U9 U. a7 O$ G
mvar:
# m7 h9 W5 E4 p5 x! M8 `
main0(::a)=0 d3 Z- K ?) w+ I# R# z
a = linspace(1,1,1000),
9 D\\" N) l& X$ K9 H
t0=clock(),
( A; |- e' }7 w
k=0, while{k<10000,: P+ p3 _1 @, a
i=0, while{i<1000, a(i)=2, i++},
) \\\" [9 H5 D$ R, }
k++
4 Y! G! s, W' |. ^0 v! U; m7 }
},
3 I* ]. v: S t) l% _
[clock()-t0]/1000;
6 G8 R1 k\\" x( h$ m0 d6 t
. G\\" }) _1 S/ |, c; q* X8 L
main1(::a)=- S. m9 p1 ^8 ?
t0=clock(), s=0,
+ O: h6 {# N8 l6 D, C) b
k=0, while{k<10000, @2 y0 y0 U6 I* t
i=0, while{i<1000, s=s+a(i), i++},
; F, b\\" y' S% ^
k++
; ]5 b1 c( X5 a& k, r9 g& v7 N2 h& y
},! @2 P+ x5 U+ B) ?\\" J1 t; E
s;! ?+ g6 ~9 t, c; P5 D4 t; Y3 h- m\\" X
[clock()-t0]/1000;

!using["math","sys"]; 
5 U9 U. a7 O$ G
mvar: 
# m7 h9 W5 E4 p5 x! M8 `
main0(::a)=0 d3 Z- K  ?) w+ I# R# z

a = linspace(1,1,1000), 
9 D\\" N) l& X$ K9 H
t0=clock(), 
( A; |- e' }7 w
k=0, while{k<10000,: P+ p3 _1 @, a

i=0, while{i<1000, a(i)=2, i++}, 
) \\\" [9 H5 D$ R, }
  k++ 
4 Y! G! s, W' |. ^0 v! U; m7 }
}, 
3 I* ]. v: S  t) l% _
[clock()-t0]/1000; 
6 G8 R1 k\\" x( h$ m0 d6 t
. G\\" }) _1 S/ |, c; q* X8 L

main1(::a)=- S. m9 p1 ^8 ?

t0=clock(), s=0, 
+ O: h6 {# N8 l6 D, C) b
k=0, while{k<10000,  @2 y0 y0 U6 I* t

i=0, while{i<1000, s=s+a(i), i++}, 
; F, b\\" y' S% ^
  k++ 
; ]5 b1 c( X5 a& k, r9 g& v7 N2 h& y
},! @2 P+ x5 U+ B) ?\\" J1 t; E

s;! ?+ g6 ~9 t, c; P5 D4 t; Y3 h- m\\" X

[clock()-t0]/1000;

结果：
15.484 //存数据所用时间，秒
20000000.
5.531 //取数据所用时间，秒

－－－－－－－－－－

用函数sys::A()存取数组数据效率稍高些：

!using["math","sys"];
' J4 t! e. X! k' [
mvar:% F; [8 j* \- x! \
main0(::a)=' C& J# Q' V1 S. r
a = linspace(1,1,1000),
9 l/ A3 a6 \! V: ~
t0=clock(),
5 l7 O$ b _1 w7 o; T# i8 L$ ]0 k
k=0, while{k<10000,; E0 Q% j! K# X' A) B: v. e
i=0, while{i<1000, A(a,i)=2, i++},
9 U) ` i\\" G0 U- [! A
k++: ]( B: R1 d0 [1 @$ D- y
},5 U( V; H9 _# Z' n- o; O( E
[clock()-t0]/1000;; W/ U V1 R: E9 w
( d4 r. Q# N% `. O5 f8 A5 ^
main1(::a)=
P- A8 U5 B1 _2 w# m. U
t0=clock(), s=0,5 N& `) S4 _8 B' `; w' H8 {. w
k=0, while{k<10000,
# g# A6 w- t. H, {3 O
i=0, while{i<1000, s=s+A(a,i), i++},
: }4 x5 E0 X7 }\\" s# A N( J
k++4 i$ e! r$ |; O& ^
},1 }+ h% z' W3 w! T! I F
s;
( Q* T) A F8 H\\" s$ v* c1 S3 a/ G
[clock()-t0]/1000;

!using["math","sys"]; 
' J4 t! e. X! k' [
mvar:% F; [8 j* \- x! \

main0(::a)=' C& J# Q' V1 S. r

a = linspace(1,1,1000), 
9 l/ A3 a6 \! V: ~
t0=clock(), 
5 l7 O$ b  _1 w7 o; T# i8 L$ ]0 k
k=0, while{k<10000,; E0 Q% j! K# X' A) B: v. e

i=0, while{i<1000, A(a,i)=2, i++}, 
9 U) `  i\\" G0 U- [! A
  k++: ]( B: R1 d0 [1 @$ D- y

},5 U( V; H9 _# Z' n- o; O( E

[clock()-t0]/1000;; W/ U  V1 R: E9 w

( d4 r. Q# N% `. O5 f8 A5 ^

main1(::a)= 
  P- A8 U5 B1 _2 w# m. U
t0=clock(), s=0,5 N& `) S4 _8 B' `; w' H8 {. w

k=0, while{k<10000, 
# g# A6 w- t. H, {3 O
  i=0, while{i<1000, s=s+A(a,i), i++}, 
: }4 x5 E0 X7 }\\" s# A  N( J
  k++4 i$ e! r$ |; O& ^

},1 }+ h% z' W3 w! T! I  F

s; 
( Q* T) A  F8 H\\" s$ v* c1 S3 a/ G
[clock()-t0]/1000;

结果：
10.469
20000000.
3.781

－－－－－－－

在不使用函数sys::A()时，Forcal存数组数据耗时与matlab之比约110:1，取数组数据耗时与matlab之比约60:1。
在使用函数sys::A( )时，Forcal存数组数据耗时与matlab之比约75:1，取数组数据耗时与matlab之比约41:1。

matlab的JIT加速器对数组元素的存取效果是明显的，获得了与C/C++近乎相同的存取效率。

Forcal之所以存取比较慢，原因在于动态数组由扩展库提供，数组元素的存取由函数的层层调用来实现。不过，数组元素的存取效率体现了Forcal对扩展库中任意对象的存取效率。

5、在循环中加点东西

matlab代码：

>> clear all5 X# `( V\" {& n9 L- f
f=@(x)x+1;, K* m1 ^' G) w6 u
a = linspace(2,2,1000);5 r, t8 L\" H2 z5 D: C5 Y( c [. I
tic
( X- y& H9 Y3 z6 [) h8 v* W
s=0;
\" _5 W, h# _+ ~0 i6 ?5 A
for k=1:1000
) I% |! V& d' d4 `
for i=1:10002 W$ H& A\" J. J0 S( }! W
s=s+f(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i))+1+sin(a(i));
, @: u& I' s0 ?4 b$ p0 B2 ~7 C1 N w
end4 ~/ [\" Y* l* O1 k5 i- z( c. v
end
2 c5 X: }) Z% C1 [% [9 d: \
s$ s$ e( o( y' m8 w
toc
+ \$ |' j/ g1 m
6 O% z4 }8 Q/ A O' K5 D# `
s =+ h2 |% V1 G) f; n( b
7 S) H7 `$ k, A# r. R! P
2.4002e+007
3 h\" T$ [, l% b: _: |1 @
$ I0 N: q- M( b, o) @; }: N
Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

!using["math","sys"];$ h7 H+ ?; ]8 _) t( l6 a' Y
mvar:9 R; i7 X3 `3 l1 D; V& u J; r/ o
f(x)=x+1; b1 R6 Q( {! }
main0(::a)=
' x* Q$ v\\" o1 o& R$ u) C0 |. A
a = linspace(2,2,1000),
: d4 f5 q+ k0 X; ^& x; X
t0=clock(), s=0,3 K4 v3 \7 O\\" l% I
k=0, while{k<1000,# @6 i! @! ?4 ^+ _+ P& t [
i=0, while{i<1000, s=s+f(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i)), i++},
: x' q$ o- N% y4 c
k++
$ k! v1 v0 z3 }9 V9 P$ Y1 d8 m
},
( W; z4 D* {2 Q7 L# O w
s;6 o1 s4 E) v\\" Y9 C* E# S
[clock()-t0]/1000;

!using["math","sys"];$ h7 H+ ?; ]8 _) t( l6 a' Y

mvar:9 R; i7 X3 `3 l1 D; V& u  J; r/ o

f(x)=x+1;  b1 R6 Q( {! }

main0(::a)= 
' x* Q$ v\\" o1 o& R$ u) C0 |. A
a = linspace(2,2,1000), 
: d4 f5 q+ k0 X; ^& x; X
t0=clock(), s=0,3 K4 v3 \7 O\\" l% I

k=0, while{k<1000,# @6 i! @! ?4 ^+ _+ P& t  [

i=0, while{i<1000, s=s+f(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i))+1+sin(A(a,i)), i++}, 
: x' q$ o- N% y4 c
  k++ 
$ k! v1 v0 z3 }9 V9 P$ Y1 d8 m
}, 
( W; z4 D* {2 Q7 L# O  w
s;6 o1 s4 E) v\\" Y9 C* E# S

[clock()-t0]/1000;

结果：
24002271.69189384
5.172

二者效率现在只有一倍差距。

－－－－－－－－－－

不过，显然matlab和Forcal的程序应该这样写：

matlab程序：

clear all$ f6 y4 N9 T\" E9 p
f=@(x)x+1;
! g8 Z6 L* e6 d- M
a = linspace(2,2,1000);$ N& F. Z! p6 @; n1 y9 V: g6 f
tic
' i, ~2 l5 X% D# d
s=0;
4 E7 K. R6 j( i; t
for k=1:1000
- A' E8 X2 C' u/ o, Y. @
for i=1:1000
4 i# t- Y: |6 d# N1 Q3 R: ]3 D
t=a(i);; k0 \- W/ R4 A) ~0 ~
s=s+f(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t);% S* }! `0 b4 t1 P
end
; M& T! k5 T0 n2 {
end
4 L5 p0 K, a: D: E\" @8 b
s
/ x5 T, \1 n* Y
toc# k, X3 ]9 b( J/ ?1 ]) C
\" n. b( X! h; g( s `
s =, S8 k- q6 D& W7 m6 ` s6 \* e
* ~! g5 {; h# L\" I/ U( ?6 F
2.4002e+007# L+ ?) R9 k* J8 m
: z& U7 L8 X5 g) e4 J/ ?# y
Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

!using["math","sys"];
$ L5 j) \\\" \) _ S: y `: z: F0 d
mvar:' W. p6 D2 v4 G0 a, y, i
f(x)=x+1;
5 y9 N7 b; Q* K
main0(::a)=5 j\\" s1 Q# b9 @7 ^! n3 f
a = linspace(2,2,1000),8 v* f% c, S$ Z& E5 w; J* z
t0=clock(), s=0,7 g' z) ]8 `- |5 }% A1 t# H, y8 @
k=0, while{k<1000,- J7 q2 T4 ?' m
i=0, while{i<1000, t=a(i), s=s+f(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t)+1+sin(t), i++},
) r J5 r- x' n8 V
k++
. f0 @( ~0 q+ a
},
6 t) N$ w3 h2 ?/ d7 V
s;
2 w+ Y. L. b- A4 @ p
[clock()-t0]/1000;

结果：
24002271.69189384
1.625

matlab效率已不及Forcal。

提示：目前Forcal对数组元素存取效率较低，要想提高效率：（1）用临时变量暂存数组元素；（2）用专用函数sys::GA( )和sys::SA( )一次存取大量连续元素；（3）对数组整体操作。

我也说一句

发表

使用道具举报

forcal

45 主题	3 听众	282 积分

升级 91%

TA的每日心情

	难过 2012-8-27 18:22

签到天数: 1 天

[LV.1]初来乍到

4^#

发表于 2011-8-2 09:27 |只看该作者 |招呼Ta 关注Ta

在这里顺便说一下Forcal静态数组，静态数组是用字符串模拟的，但数组元素存取效率非常高，约是C/C++静态数组元素存取效率的1/5左右。

Forcal代码：

!using["sys"];; x9 N& O. d, h# b3 X: Q
mvar:( P6 F8 ?# n9 ?& ~
a="\&\[1000]"/4, //用字符串模拟静态数组
# s* }* |. r; R8 B\\" _9 f
t0=clock(),2 z( n; ]9 o; V; h8 j4 N/ y
k=0, while{k<10000,
% q; S7 o& p8 w6 t
i=0, while{i<1000, setn(a+i,2), i++}, //用函数setn(a+i,2)对静态数组元素赋值5 L: G0 ?\\" c9 n# Q, u* ^
k++- J9 `6 v, S7 n+ T; T$ m6 }( b
},5 Y2 m# w; B7 H' T
[clock()-t0]/1000;
, a- n1 v+ L1 t( j8 @6 Y6 e; R6 o
j: O! w& G5 [\\" Z F
a="\&\[1000]"/4, //用字符串模拟静态数组2 u8 d! z3 L2 E) \2 u
i=0, while{i<1000, setn(a+i,2), i++}, //用函数setn(a+i,2)对静态数组元素赋值1 ^# ^* d; G3 X: `) T$ Z
t0=clock(),& { {6 q6 {' s7 s6 g0 j
s=0,
6 V! m4 a0 Z- F. _\\" y
k=0, while{k<10000,
$ `# U/ q1 R. ^5 O
i=0, while{i<1000, s=s+getn(a+i), i++}, //用函数getn(a+i)获得静态数组元素的值
$ K$ S. d! M( W8 m5 `
k++
- t0 _7 m! l% u3 O0 S2 ?
},& T# K/ l$ o6 j( G% c
s;! n- ]6 `6 w* h2 g7 w3 E1 @& \4 ^
[clock()-t0]/1000;

!using["sys"];; x9 N& O. d, h# b3 X: Q

mvar:( P6 F8 ?# n9 ?& ~

a="\&\[1000]"/4,    //用字符串模拟静态数组 
# s* }* |. r; R8 B\\" _9 f
t0=clock(),2 z( n; ]9 o; V; h8 j4 N/ y

k=0, while{k<10000, 
% q; S7 o& p8 w6 t
  i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值5 L: G0 ?\\" c9 n# Q, u* ^

k++- J9 `6 v, S7 n+ T; T$ m6 }( b

},5 Y2 m# w; B7 H' T

[clock()-t0]/1000; 
, a- n1 v+ L1 t( j8 @6 Y6 e; R6 o
  j: O! w& G5 [\\" Z  F

a="\&\[1000]"/4,    //用字符串模拟静态数组2 u8 d! z3 L2 E) \2 u

i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值1 ^# ^* d; G3 X: `) T$ Z

t0=clock(),& {  {6 q6 {' s7 s6 g0 j

s=0, 
6 V! m4 a0 Z- F. _\\" y
k=0, while{k<10000, 
$ `# U/ q1 R. ^5 O
  i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值 
$ K$ S. d! M( W8 m5 `
  k++ 
- t0 _7 m! l% u3 O0 S2 ?
},& T# K/ l$ o6 j( G% c

s;! n- ]6 `6 w* h2 g7 w3 E1 @& \4 ^

[clock()-t0]/1000;

结果：
0.5620000000000001 //赋值用时间，秒
20000000.
0.579 //取值用时间，秒

我也说一句

发表

使用道具举报

返回列表

帐号		密码		只需要一步，快速开始		注册地址	找回密码

极限测试之Matlab与Forcal有无JIT加速器的比较

QQ

电话咨询

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务|