查看: 5588|回复: 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 all9 k. ~4 _\" F E\" g+ B- Y0 g
tic7 b/ T7 @* Y+ r
s=0;
! t\" T- A\" h- J0 R9 Z$ K! \7 e
for k=1:10000% N O# c. c+ V3 L8 G
for i=1:1000
3 x\" }\" j) ~1 L9 m
s=s+1;' c\" Z; S2 u& D6 e, M4 E3 E2 e
end5 C# T3 N9 x1 B4 j3 r- ^% c
end
0 Z: |* f- Q1 |( H. \. ~
s
/ t& n- X( T' R7 A, d
toc- W) K/ ^1 M( ~& L+ W' j8 ^4 P
$ {0 t% s5 o3 x, ]! ~0 K2 B
s =
/ G/ c) P- d$ r/ w. p
1 I. m- x' {) c0 n
10000000
' Y0 [( u+ @; ?1 X* d* x2 H6 t
8 r\" p$ s1 c, U6 Y4 l0 ]. z( e
Elapsed time is 0.092808 seconds.

复制代码

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

Forcal没有JIT加速器，代码：

!using["sys"];5 h3 u6 j+ c' v( C( R; M
mvar:
0 Y. g\\" \4 O5 r2 [5 [
t0=clock(),
$ S+ e3 h) w6 P# A4 X5 R
s=0,
1 C' y/ ]$ r6 Z ^( c ~
k=0, while{k<10000,
$ i; ?+ p' C. s
i=0, while{i<1000,s=s+1, i++},1 E2 l! @9 o2 Y$ E6 Z
k++
% U. ~% z) M5 T. m4 ]. X/ f
},
5 u! A) m4 S( c$ w' G8 c' z
s;9 l$ J% x. _% W) H, a9 y
[clock()-t0]/1000;

结果：
10000000.
0.704 秒

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

2、循环中稍微加点东西

Matlab代码：

clear all4 M+ y! c# u& |6 V$ F8 Y2 f
tic
$ Y, i\" i9 V: v& v
s=0;# J8 g, S\" e# \+ W! W- m2 x$ z# X
for k=1:10002 e- r5 r6 A, b/ N/ g% s\" T\" N
for i=1:10001 ?\" J9 P! N5 o5 ?* h
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);
( Q2 ^% P+ F. D- S6 C
end
' x, L/ d2 w ?( S5 P( {
end
s% ~$ V# t$ f& D9 h6 U
toc; [7 n) d' x2 d( w8 Q+ r4 ?2 Q3 {+ ^
3 @$ x9 W! r8 y\" b& {* d
s =7 n9 N0 H% c( Z! O: ]5 S
* q* P\" o6 u) L7 ^4 w
1.1009e+007
' A& t5 S; s' J3 S( h
) L! E) e2 T7 a( M* X) t
Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

!using["sys"];8 I! r/ U; p$ `9 @0 R. @8 W9 Y
mvar:; g4 L1 w# ]8 H |8 A/ k7 U3 M! f% X
t0=clock(),
3 g* F3 ?) W4 ?% K
s=0,
. ^9 A# f7 }: ]# f. n
k=1, while{k<=1000,' Y7 s2 [\\" {/ y& W# t L8 [9 T
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++},
; R/ p; H* h3 w3 S' s: G5 O5 x( a
k++: D' u+ Z% g! e; |' r
}, S6 C7 l9 ?% R; |
s;' k2 r, o7 M9 S6 b
[clock()-t0]/1000;

结果：
11008953.66602346
1.031

二者效率相差已比较小。

3、再增加一个函数调用

Matlab代码：

clear all, p9 J/ j: u6 a& I2 `
f=@(x)x+1;# r! a- ^: i1 [: I
tic
- K$ z. g' a# f) P
s=0;+ f. ?# u3 Z) ~\" a) {: R' `
for k=1:1000
* h4 U e: Q1 r# r5 c7 W
for i=1:1000
3 g( F+ C2 w4 v! j3 k
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);) V! H l! T ]0 k; \
end
$ L6 B+ I\" Z+ v1 l) ~
end
2 x7 e2 @9 p) {
s
9 i( C7 S( g7 ]' z
toc
) A2 _2 C( `8 f! h4 U
+ O( r3 z8 c) S* w4 l: _
s =2 ]& P! ]+ Q/ \\" u- G1 {8 Q
( N/ w5 B, @! a
5.1251e+008* q! q, X8 `- M* u6 q9 ^
`9 ]! l3 h8 L
Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

!using["sys"];) ^6 t\\" g\\" x. U- r
f(x)=x+1;
) {# s& [. i1 Z/ T! R! D0 J
mvar:# R% u4 \; l( ]4 J( k: {
t0=clock(),3 @3 H) {+ r u+ X/ K. }+ N' C
s=0,: R; h' @' p' X( ? d4 `
k=1, while{k<=1000,0 E3 N. o* |\\" N% `$ I
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++},
( T3 G5 g7 T, r2 b9 o# ? j
k++
$ J! i5 [/ T! y$ |1 x
},
% }- H5 q9 g: o! O5 R
s;
; \; Y, t0 N+ u\\" d0 L- p+ v+ K
[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 all6 D& Q% A5 f- K0 f# z6 D
a = linspace(1,1,1000);
, }3 \3 d+ Z/ v
tic, S0 }5 r& S$ j6 c5 c5 k) p
for k=1:10000
& ~) |) x& i% l8 ]- B
for i=1:1000( C- j4 H+ @: [0 q2 Q
a(i)=2;
9 w2 f\" o5 p' l) q( } M
end E: F! ^% [6 }0 A9 [) Z, |6 Y+ M9 F/ J
end
1 Z- E# ?# H. J7 u
toc
$ Q% P w3 A' {* z
% D, U; b$ P- w8 m/ A9 s# u
tic2 I* O% p% B- k! G% y
s=0;1 ^+ ]1 S' a4 {0 K# @7 [- u$ V; e
for k=1:100000 r7 i2 V2 F* I: L7 `0 T5 S
for i=1:1000
' u+ n1 M) C; K
s=s+a(i);
4 @# |\" @7 X* A( Y
end
1 ]' i2 u0 B8 n7 e1 L o
end
! \( ]+ Y( `3 D! G6 g: p
s
7 D( f) E _6 S- p/ u
toc
: V' C5 P/ e6 K0 P5 I! @2 O3 _: o
Elapsed time is 0.140658 seconds.
; \' @8 I4 [1 A/ q! P1 a- _% s
1 O' f# ^$ f. a. v8 @- e! `+ Q; R
s =8 i\" I/ x8 U# N( V& ]+ m
4 c8 v1 K1 n1 l; L& u( h
200000004 }$ d' ~( \* d$ |/ p
: M4 c3 u0 Q, ~
Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

!using["math","sys"];: _& B. G. ^8 O) Z1 J
mvar:
; H& \& G. ?7 n8 w\\" q0 [
main0(::a)=
Y! i/ `5 z% q- P
a = linspace(1,1,1000), N- @9 ?3 {0 @
t0=clock(),
) z# A+ f# U/ o+ R' \
k=0, while{k<10000,( h: t3 ?( Z& A
i=0, while{i<1000, a(i)=2, i++},* O' q% N; P, j* |. B- l
k++
$ q2 q. ~+ P& B' }9 Z
},
$ v. q! Q, a# U# t+ I4 z0 c, Q
[clock()-t0]/1000;/ X$ l _# H1 K( W/ s2 M
6 E; F: |' ?) s\\" ] x
main1(::a)=
& C6 @1 d' X9 y Q2 C
t0=clock(), s=0,5 \' a* k2 n/ k' h/ [
k=0, while{k<10000,
! w# s) I- h1 x. F- v+ c0 P
i=0, while{i<1000, s=s+a(i), i++},3 r4 ?% c\\" }: ?2 ~+ \
k++7 U3 d4 u; S7 L4 _\\" h' q
},, M! S; O0 u x' Y0 a' f [
s;
! \0 ^5 |! v& j- w$ p3 N' ?9 a$ Q
[clock()-t0]/1000;

!using["math","sys"];: _& B. G. ^8 O) Z1 J

mvar: 
; H& \& G. ?7 n8 w\\" q0 [
main0(::a)= 
  Y! i/ `5 z% q- P
a = linspace(1,1,1000),  N- @9 ?3 {0 @

t0=clock(), 
) z# A+ f# U/ o+ R' \
k=0, while{k<10000,( h: t3 ?( Z& A

i=0, while{i<1000, a(i)=2, i++},* O' q% N; P, j* |. B- l

k++ 
$ q2 q. ~+ P& B' }9 Z
}, 
$ v. q! Q, a# U# t+ I4 z0 c, Q
[clock()-t0]/1000;/ X$ l  _# H1 K( W/ s2 M

6 E; F: |' ?) s\\" ]  x
main1(::a)= 
& C6 @1 d' X9 y  Q2 C
t0=clock(), s=0,5 \' a* k2 n/ k' h/ [

k=0, while{k<10000, 
! w# s) I- h1 x. F- v+ c0 P
  i=0, while{i<1000, s=s+a(i), i++},3 r4 ?% c\\" }: ?2 ~+ \

k++7 U3 d4 u; S7 L4 _\\" h' q

},, M! S; O0 u  x' Y0 a' f  [

s; 
! \0 ^5 |! v& j- w$ p3 N' ?9 a$ Q
[clock()-t0]/1000;

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

－－－－－－－－－－

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

!using["math","sys"];
; t' |4 r+ x9 u' f' ]
mvar:' S) i! r. I# _- O* w
main0(::a)=
\\" P& o# x- I& y/ S; D: j% `* u
a = linspace(1,1,1000),/ d. ]0 f5 K! o, `6 q& D
t0=clock(),\\" ~, `+ Q7 e! }* x% ? q$ `+ _
k=0, while{k<10000,4 z: ?5 X. G+ D7 j' x* E% L0 i5 a
i=0, while{i<1000, A(a,i)=2, i++},6 r* x7 S8 f# k U, R! r1 l
k++& G. M$ ]7 F$ P% W8 V
},$ b! ]! Q7 O9 q1 d9 n
[clock()-t0]/1000;& [5 r' Z D) k
, Z' C% V+ T2 g/ O
main1(::a)=
# f\\" D2 B1 h7 m. z
t0=clock(), s=0,/ v& w* U4 a+ |, F# L
k=0, while{k<10000,\\" O% D6 D5 J* ]7 l- r
i=0, while{i<1000, s=s+A(a,i), i++},
% B! k' S% B/ o
k++
2 w8 o: D$ r# Y) q+ K, C2 n% H
},. u5 w8 s8 m. S7 _
s;
! N* I# f! i C* |2 J
[clock()-t0]/1000;

!using["math","sys"]; 
; t' |4 r+ x9 u' f' ]
mvar:' S) i! r. I# _- O* w

main0(::a)= 
\\" P& o# x- I& y/ S; D: j% `* u
a = linspace(1,1,1000),/ d. ]0 f5 K! o, `6 q& D

t0=clock(),\\" ~, `+ Q7 e! }* x% ?  q$ `+ _

k=0, while{k<10000,4 z: ?5 X. G+ D7 j' x* E% L0 i5 a

i=0, while{i<1000, A(a,i)=2, i++},6 r* x7 S8 f# k  U, R! r1 l

k++& G. M$ ]7 F$ P% W8 V

},$ b! ]! Q7 O9 q1 d9 n

[clock()-t0]/1000;& [5 r' Z  D) k

, Z' C% V+ T2 g/ O

main1(::a)= 
# f\\" D2 B1 h7 m. z
t0=clock(), s=0,/ v& w* U4 a+ |, F# L

k=0, while{k<10000,\\" O% D6 D5 J* ]7 l- r

i=0, while{i<1000, s=s+A(a,i), i++}, 
% B! k' S% B/ o
  k++ 
2 w8 o: D$ r# Y) q+ K, C2 n% H
},. u5 w8 s8 m. S7 _

s; 
! N* I# f! i  C* |2 J
[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 all& v( B. k) F' r5 {% ]: q
f=@(x)x+1;\" ]7 d1 h% \1 m3 ~/ x
a = linspace(2,2,1000);. T) T8 k2 ?5 e( ` m) Q
tic\" c c3 Q6 N! z. d* q8 W
s=0;
\" y% ]/ c/ |9 w. B( D: `\" u) {
for k=1:1000
\" s& L$ U\" k U2 B# b! ]1 x
for i=1:1000; c/ D* e+ g# `/ d* P
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));
/ e6 c( m; f/ g. _
end; C: \. U\" R6 L
end
9 v* }; d$ A! a8 v( n
s+ ?0 {9 |, T- Y7 q: y% x. J* b
toc0 Q6 G% G- W0 n. W- a0 E
: i6 \5 Z* U- R- u
s =1 N% }5 r O$ A. [) f
* P: G7 M- V. z0 h, {
2.4002e+007
$ t7 j7 b7 A5 [' W
: L9 f7 J0 r. M0 v- l\" ]7 Q
Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

!using["math","sys"];) M6 M, C# ]7 n; M5 X3 M
mvar:
& ?1 y* I2 a6 w' f( j3 s
f(x)=x+1;% l, F) E& r/ o/ Z# Z/ u* u9 ]/ o7 ]# D
main0(::a)=6 A' h( F* U0 A2 K
a = linspace(2,2,1000),
; g f( q% F. t
t0=clock(), s=0,
1 K+ a+ q# T% H0 u J
k=0, while{k<1000,* {) k/ e# M0 ]5 @) f
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++},4 o! Z4 U. Y8 R7 H0 T% d2 [+ W
k++) o0 G& V K2 C
},& s+ X/ q7 \7 n: d
s;
7 ]8 o\\" z\\" o$ C: U
[clock()-t0]/1000;

!using["math","sys"];) M6 M, C# ]7 n; M5 X3 M

mvar: 
& ?1 y* I2 a6 w' f( j3 s
f(x)=x+1;% l, F) E& r/ o/ Z# Z/ u* u9 ]/ o7 ]# D

main0(::a)=6 A' h( F* U0 A2 K

a = linspace(2,2,1000), 
; g  f( q% F. t
t0=clock(), s=0, 
1 K+ a+ q# T% H0 u  J
k=0, while{k<1000,* {) k/ e# M0 ]5 @) f

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++},4 o! Z4 U. Y8 R7 H0 T% d2 [+ W

k++) o0 G& V  K2 C

},& s+ X/ q7 \7 n: d

s; 
7 ]8 o\\" z\\" o$ C: U
[clock()-t0]/1000;

结果：
24002271.69189384
5.172

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

－－－－－－－－－－

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

matlab程序：

clear all
! C. T* e H/ T3 X8 ~1 G
f=@(x)x+1;
- {! {! T& m3 I9 ?; X9 h
a = linspace(2,2,1000);
3 T: P: l) Q4 |8 a# R
tic
\" p j% ^; K' X3 U/ E/ g% ~8 D
s=0;* r& I\" q8 A: G$ t0 h& J2 T- ^* ^
for k=1:1000/ p8 Z# V) h( @\" P1 M. h
for i=1:10004 E1 I* r! s; G3 `- _% w9 W
t=a(i);
2 w, l: ^ @* w. _* ^
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);
+ w; \\" J% ~9 ~, e\" L8 i/ Z0 j* o' C
end5 b5 c; J* P) e$ x# j; q' S
end
6 t6 _. y1 g( L
s
# a. c+ W' P( n
toc
! Q5 d1 w. Z6 }1 i8 f. m8 W$ z
3 K3 f, m3 c+ X
s =
& T- Z9 G$ _& X
& S! r: R& W+ l7 O3 P. `; ^) h1 y
2.4002e+007
1 ]8 t( o3 c. m- k3 W' S) C
5 {: M; l- N+ @* D( r
Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

!using["math","sys"]; Y/ `. L9 [) M0 n n7 W0 e' S
mvar:
/ k# ^0 j9 W9 m! x' ~% p9 f
f(x)=x+1;
4 C `1 i9 }( c4 ^6 M
main0(::a)=
% Q$ ?. g6 \6 ?. z\\" @
a = linspace(2,2,1000),; u8 {2 ^* n. s( f `+ \# q
t0=clock(), s=0,
9 s5 y; a: \2 l8 L9 u$ k8 M. r\\" `
k=0, while{k<1000,: _( |- A% L8 K$ k5 Q; Y
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++},
/ x$ ~0 z9 ]0 z- [3 U: J x4 y! O
k++1 x& _3 H+ d9 I( m' l6 T
},+ ~; i! @2 r; Y# P
s;8 o! |: d2 w7 }: Z ]
[clock()-t0]/1000;

!using["math","sys"];  Y/ `. L9 [) M0 n  n7 W0 e' S

mvar: 
/ k# ^0 j9 W9 m! x' ~% p9 f
f(x)=x+1; 
4 C  `1 i9 }( c4 ^6 M
main0(::a)= 
% Q$ ?. g6 \6 ?. z\\" @
a = linspace(2,2,1000),; u8 {2 ^* n. s( f  `+ \# q

t0=clock(), s=0, 
9 s5 y; a: \2 l8 L9 u$ k8 M. r\\" `
k=0, while{k<1000,: _( |- A% L8 K$ k5 Q; Y

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++}, 
/ x$ ~0 z9 ]0 z- [3 U: J  x4 y! O
  k++1 x& _3 H+ d9 I( m' l6 T

},+ ~; i! @2 r; Y# P

s;8 o! |: d2 w7 }: Z  ]

[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"];
3 x, j& `$ d4 p
mvar:
6 l: j; O3 q0 ~4 a7 C8 J0 u1 i
a="\&\[1000]"/4, //用字符串模拟静态数组\\" x4 Z1 d' i; V1 L\\" m
t0=clock(),6 [8 K5 a* T: z8 y4 W( D
k=0, while{k<10000,
, w0 I\\" G/ J+ r\\" {
i=0, while{i<1000, setn(a+i,2), i++}, //用函数setn(a+i,2)对静态数组元素赋值
5 A3 T. g1 ^1 b2 Z4 x9 y
k++! ~$ [\\" Z: }, ?! m
},
: i* R4 O' }4 H! D4 a! U! ~
[clock()-t0]/1000;# h6 o, R1 @8 O! B- Y( ^+ l5 O
9 |% s4 A# A2 [
a="\&\[1000]"/4, //用字符串模拟静态数组: _6 P6 H$ j% F1 L7 w6 P
i=0, while{i<1000, setn(a+i,2), i++}, //用函数setn(a+i,2)对静态数组元素赋值
t0=clock(),8 b- j& U3 y t# X* Z* L. M/ F9 |
s=0,
7 o- o) V0 H7 ~0 u/ p/ E/ F- Z
k=0, while{k<10000,
$ z, r! }$ M, l3 I8 v: E1 G\\" D
i=0, while{i<1000, s=s+getn(a+i), i++}, //用函数getn(a+i)获得静态数组元素的值
% y! L- _/ m4 c2 ^ o% g$ t: w
k++: c1 S, _4 Z6 o4 [8 A
},
\\" V; i5 j, f; W. q$ m8 `5 T
s;! c4 A% ?7 Q9 R( `5 Y+ O
[clock()-t0]/1000;

!using["sys"]; 
3 x, j& `$ d4 p
mvar: 
6 l: j; O3 q0 ~4 a7 C8 J0 u1 i
a="\&\[1000]"/4,    //用字符串模拟静态数组\\" x4 Z1 d' i; V1 L\\" m

t0=clock(),6 [8 K5 a* T: z8 y4 W( D

k=0, while{k<10000, 
, w0 I\\" G/ J+ r\\" {
  i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值 
5 A3 T. g1 ^1 b2 Z4 x9 y
  k++! ~$ [\\" Z: }, ?! m

}, 
: i* R4 O' }4 H! D4 a! U! ~
[clock()-t0]/1000;# h6 o, R1 @8 O! B- Y( ^+ l5 O

9 |% s4 A# A2 [

a="\&\[1000]"/4,    //用字符串模拟静态数组: _6 P6 H$ j% F1 L7 w6 P

i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值 
' Z+ a/ f3 O+ d\\" r5 i\\" I1 k$ ?  E( t
t0=clock(),8 b- j& U3 y  t# X* Z* L. M/ F9 |

s=0, 
7 o- o) V0 H7 ~0 u/ p/ E/ F- Z
k=0, while{k<10000, 
$ z, r! }$ M, l3 I8 v: E1 G\\" D
  i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值 
% y! L- _/ m4 c2 ^  o% g$ t: w
  k++: c1 S, _4 Z6 o4 [8 A

}, 
\\" V; i5 j, f; W. q$ m8 `5 T
s;! c4 A% ?7 Q9 R( `5 Y+ O

[clock()-t0]/1000;

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

我也说一句

发表

使用道具举报

返回列表

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

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

浏览过的版块

QQ

电话咨询

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