查看: 5658|回复: 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* m0 R4 o, J) _\" D
tic$ K0 v, f: s. n
s=0;
\" G- _! b# a9 i3 n& T5 E, [/ a
for k=1:100007 X$ Z8 B# L; b$ P
for i=1:1000
2 ^& q2 G/ @1 z9 H* I
s=s+1;
% L. i, N& h& Z; r* P+ K
end/ C! O: S+ b9 b! K4 v
end% s# |' l' G o# J2 d& H! ^
s
2 C, f, A0 ` G, w- l, J
toc
4 L6 _' ]6 ^# o3 N; ^0 h
3 X1 E8 y. H- ]9 ?5 c# e* }
s =- y. ?5 O( N& Q( r2 s# u. @
& j8 i+ w- }% R5 j5 B2 j$ }
10000000' r. O' Z3 u5 O# f0 A% y
1 \' u0 J g' \4 w. m
Elapsed time is 0.092808 seconds.

复制代码

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

Forcal没有JIT加速器，代码：

!using["sys"];. ?. Q& |\\" c0 d! h h
mvar:
6 x! d8 r/ `% ~9 O3 x+ s
t0=clock(),9 e. C# \' V3 C) X9 x4 K1 \9 c
s=0,
+ Q9 Y* F) s; _! K
k=0, while{k<10000,/ i x\\" ?2 P9 d, A; s) @8 y Q
i=0, while{i<1000,s=s+1, i++},: W9 K5 C% ^( F; `; X2 ]
k++8 k# y& g4 [ R& _( b
},1 B' l\\" ?% w0 ?) Q; ?
s;
% q8 A5 I3 R' @; S* u' q% X# g F
[clock()-t0]/1000;

结果：
10000000.
0.704 秒

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

2、循环中稍微加点东西

Matlab代码：

clear all8 |) h8 G4 M( e' e9 D
tic
0 K2 @' V( v: y) `
s=0;
{ B\" `$ l8 s/ B
for k=1:1000
/ _3 Z2 M$ F: `& C9 m\" \0 z
for i=1:1000
2 W3 ?- L! I4 B6 u& S\" z5 {
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);
& P7 Y3 ~& d/ q5 W
end+ E9 J/ m6 t% }: B' }' G: S
end$ m# ?7 r: c( g. D% [
s$ c1 s3 _3 a! @8 f1 ~& B& P
toc
3 ]/ e. Z) O# R- W* l( Z
. h. S7 z7 n# Z# G1 G. @1 s
s =$ o& z\" F5 k0 t5 ~
$ t j0 D* A N/ A5 P* W1 h
1.1009e+007
9 K6 k\" x9 v( D( f: r
, P$ C l: b8 P {
Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

!using["sys"];
0 ^, } Z3 S5 X; O. }4 e
mvar:
3 D m& N3 C, J\\" p+ |! v% _
t0=clock(),; ], e2 b V2 _6 E5 D& k7 I. _
s=0,
/ ^# j! C6 q- ]; C6 o7 ^% _
k=1, while{k<=1000,
8 b. t/ Z& K- E$ m
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++},, A6 X# {6 X1 @3 T( j, ], q; f\\" G
k++& v, }7 T! R6 b) i; F' g6 I$ J
},; t/ S( u2 Y2 H; T& c2 ^
s;1 F$ E- B2 ?9 ^& s
[clock()-t0]/1000;

结果：
11008953.66602346
1.031

二者效率相差已比较小。

3、再增加一个函数调用

Matlab代码：

clear all
2 s5 I8 L8 O$ G3 ]
f=@(x)x+1;
# A7 Y8 Q3 h F# W7 y& q
tic
+ m& }: D2 H! X; T: @) j
s=0;8 Z1 F* ]. O8 R }' ?\" c
for k=1:1000
5 s9 S' ]9 d- w4 u2 `
for i=1:10007 Z* \) i3 d9 y$ y; G( P# }' }
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);
& l! t, m\" l, _2 I: |4 k* \5 r# `( j
end
. N$ V: d% }4 i( X' ?' q
end
, C7 @7 E6 {2 Y/ R' X ?
s
' R: J% f7 R5 h, e; Q/ B
toc
# @. D1 A9 R) }8 @& l7 H% p7 m
/ A9 l8 d& h( x+ l6 p! o s
s =4 d7 P6 y% T; T$ `- H6 A
2 J( ^5 j2 D% A
5.1251e+0082 U+ R k& n# c5 X1 V D8 s
\" r- ~! Z3 H& ^2 P
Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

!using["sys"];( z) g5 t+ `2 \! a) x0 q6 Z$ |7 M
f(x)=x+1;
: [, U4 B1 N1 ^3 s
mvar:
& k+ a8 `+ t( G
t0=clock(),+ h# L& w. Z+ p4 \ P) G* Z5 Q* {% j
s=0,7 d2 g+ P& T5 p# t9 U( Y
k=1, while{k<=1000,# e5 h- L: k6 H6 H7 _. i8 H6 |
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++},+ \+ X0 A6 K8 P
k++- f/ H# a* z) {2 ^, S$ V5 w1 n
},% U$ g; W/ H* y# ]5 S+ ]+ z0 k; l
s;
# M! \8 R* C7 r( `/ P
[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 all7 H) t3 {+ G+ m
a = linspace(1,1,1000);
5 F* X8 N+ f+ j) T D
tic
4 H) t2 }8 ?8 \- P2 g5 U
for k=1:10000
3 x, [+ ^\" ^8 ^
for i=1:1000 c, ~1 ?' _3 S7 B
a(i)=2;+ S7 j$ n2 R6 n
end* a {\" x3 m: m. s
end/ m1 n/ W6 Y6 W4 g
toc
; Y8 a. T2 C. K0 G
: I3 w# n7 U6 v }) i+ u\" L
tic, Y2 g' F* E' ]. d6 I/ @- g+ W
s=0;* m, A# h X, Y
for k=1:10000
# Y. i2 x: s: M- P% m- f
for i=1:10009 r$ D6 M; u8 I: k `- P
s=s+a(i); ]% J\" L: k! g$ V- B
end; ?8 W2 I\" U; {- W+ p' \
end
3 J\" t+ T6 a' P9 _' I# O4 v0 V
s
U# _3 k, F; }; g
toc V# B% q! g/ O7 Y3 _7 O9 S
Elapsed time is 0.140658 seconds.; v9 o! }0 j# y6 S# I
4 E0 R# q$ o, i\" z$ n
s =. {# {, L W1 p5 x/ P
20000000
7 [4 v$ M- L7 j, T, a
( T$ u8 I% @\" ]' D8 ]; ?4 U: s& V
Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

!using["math","sys"];
* Q5 u+ f, z% A
mvar:
3 W2 t+ w$ g0 S: I
main0(::a)=6 Q. R' S7 N2 U+ ^
a = linspace(1,1,1000),2 K- p4 ]& C1 P8 j
t0=clock(),9 ]+ Q {9 @# J i- W% r
k=0, while{k<10000,! F/ q/ c+ X+ q- n+ X
i=0, while{i<1000, a(i)=2, i++},
( b0 }4 E& _% {4 j- A0 J' ~/ t9 B
k++
1 { K# h1 R' u6 A
},
6 M% A( v\\" Y* Y* }9 M9 s
[clock()-t0]/1000;' W. `) h6 y; y8 E* v- c
! [8 p0 w; c5 i1 a6 j& N$ D, _
main1(::a)=. S6 D3 ^9 X4 n
t0=clock(), s=0,7 a J# Y: K7 D, i* q2 b5 W\\" E2 h2 n
k=0, while{k<10000,\\" f0 v( V) q/ `. ^% ~
i=0, while{i<1000, s=s+a(i), i++},
# z$ x. t1 G) }3 A
k++
% l' A/ v5 I2 D1 V) `4 ]
},. I% }9 f, f6 g* L
s;
2 [4 ^3 ~* m/ A; O1 \% J5 i
[clock()-t0]/1000;

!using["math","sys"]; 
* Q5 u+ f, z% A
mvar: 
3 W2 t+ w$ g0 S: I
main0(::a)=6 Q. R' S7 N2 U+ ^

a = linspace(1,1,1000),2 K- p4 ]& C1 P8 j

t0=clock(),9 ]+ Q  {9 @# J  i- W% r

k=0, while{k<10000,! F/ q/ c+ X+ q- n+ X

i=0, while{i<1000, a(i)=2, i++}, 
( b0 }4 E& _% {4 j- A0 J' ~/ t9 B
  k++ 
1 {  K# h1 R' u6 A
}, 
6 M% A( v\\" Y* Y* }9 M9 s
[clock()-t0]/1000;' W. `) h6 y; y8 E* v- c

! [8 p0 w; c5 i1 a6 j& N$ D, _
main1(::a)=. S6 D3 ^9 X4 n

t0=clock(), s=0,7 a  J# Y: K7 D, i* q2 b5 W\\" E2 h2 n

k=0, while{k<10000,\\" f0 v( V) q/ `. ^% ~

i=0, while{i<1000, s=s+a(i), i++}, 
# z$ x. t1 G) }3 A
  k++ 
% l' A/ v5 I2 D1 V) `4 ]
},. I% }9 f, f6 g* L

s; 
2 [4 ^3 ~* m/ A; O1 \% J5 i
[clock()-t0]/1000;

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

－－－－－－－－－－

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

!using["math","sys"];. X: ~) j* a! v, G4 w' k
mvar:
G) c) s9 x. _
main0(::a)=% M- S4 F& O4 p1 t+ U, M
a = linspace(1,1,1000),
! s$ i9 B/ u& Q2 k* W
t0=clock(),- l) I7 o* O3 @1 e+ |/ j
k=0, while{k<10000,( |5 F8 M+ k* d4 q; A2 k
i=0, while{i<1000, A(a,i)=2, i++},* k' o- X6 [' U. w+ x- H/ y; g
k++
' A1 y( E8 X7 k5 \* b- B
},4 B$ k! n' j) A! F: c\\" P
[clock()-t0]/1000;6 m0 q; M+ \/ C
$ a& ^# B/ B- K2 F, R4 y
main1(::a)=5 B9 T% w. _3 h: V% Z
t0=clock(), s=0,8 L/ o. Y+ Z. ]8 j0 P
k=0, while{k<10000,4 `0 N) }5 q2 g, H9 s1 y
i=0, while{i<1000, s=s+A(a,i), i++},% z) P( _6 ^! t& f$ q% @
k++4 c. R. H1 s8 e' w+ x2 t
},
$ P7 d) D! o8 N# I1 K2 y4 U
s;
( b+ q\\" Y2 Q, D& S% ^4 j' F+ U( e
[clock()-t0]/1000;

!using["math","sys"];. X: ~) j* a! v, G4 w' k

mvar: 
  G) c) s9 x. _
main0(::a)=% M- S4 F& O4 p1 t+ U, M

a = linspace(1,1,1000), 
! s$ i9 B/ u& Q2 k* W
t0=clock(),- l) I7 o* O3 @1 e+ |/ j

k=0, while{k<10000,( |5 F8 M+ k* d4 q; A2 k

i=0, while{i<1000, A(a,i)=2, i++},* k' o- X6 [' U. w+ x- H/ y; g

k++ 
' A1 y( E8 X7 k5 \* b- B
},4 B$ k! n' j) A! F: c\\" P

[clock()-t0]/1000;6 m0 q; M+ \/ C

$ a& ^# B/ B- K2 F, R4 y

main1(::a)=5 B9 T% w. _3 h: V% Z

t0=clock(), s=0,8 L/ o. Y+ Z. ]8 j0 P

k=0, while{k<10000,4 `0 N) }5 q2 g, H9 s1 y

i=0, while{i<1000, s=s+A(a,i), i++},% z) P( _6 ^! t& f$ q% @

k++4 c. R. H1 s8 e' w+ x2 t

}, 
$ P7 d) D! o8 N# I1 K2 y4 U
s; 
( b+ q\\" Y2 Q, D& S% ^4 j' F+ U( e
[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
8 e5 d4 z/ k6 x/ z; W! \
f=@(x)x+1;
$ W( L8 h1 k7 r- r\" R
a = linspace(2,2,1000);& g1 d* R: M$ A' l
tic) H+ H4 }( ]4 O' I! E
s=0;
% L- W, m) E/ I! k4 M+ }; M9 a
for k=1:1000
& ]! O* ~$ {9 P1 I. ]5 N2 Q
for i=1:1000- o! @% e1 B( S0 u7 s( _4 c' a\" V
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));
; o\" L, B0 F) a! K7 @
end
) n\" m; | C) E0 V
end
, N+ X\" n( R+ m5 k: D8 u, q
s
1 c4 H7 I, l7 q( a! \
toc
* J: z; I: G/ Q% _0 A- H7 C
2 ~% ?\" }& G7 ]7 Y$ o1 X, x- I
s =0 r' @; Z3 y- w. A0 G
# ]8 n7 w9 j+ F; y/ L5 j, G1 e
2.4002e+007
* |% Q4 E/ r1 I+ v& e
7 p1 y. n0 H* B& J+ X8 h
Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

!using["math","sys"];* o [; G. [) b
mvar:
/ v* n7 `0 j4 G3 F
f(x)=x+1;
7 c: o6 U' m9 Q& N5 n, X\\" i/ K
main0(::a)=* ~' d# G8 w/ }. ^1 }6 x6 N% H
a = linspace(2,2,1000),+ {1 [; W5 }7 N+ O
t0=clock(), s=0,
6 w4 o3 g- a3 `* }\\" N
k=0, while{k<1000,
+ d P! J3 [/ ~
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++},; S l8 Z9 S6 @1 I- i
k++
4 D# {\\" X\\" d7 R' g2 p
},. w\\" Z, L, |3 g
s;/ @+ M/ }9 m* N+ C% F
[clock()-t0]/1000;

!using["math","sys"];* o  [; G. [) b

mvar: 
/ v* n7 `0 j4 G3 F
f(x)=x+1; 
7 c: o6 U' m9 Q& N5 n, X\\" i/ K
main0(::a)=* ~' d# G8 w/ }. ^1 }6 x6 N% H

a = linspace(2,2,1000),+ {1 [; W5 }7 N+ O

t0=clock(), s=0, 
6 w4 o3 g- a3 `* }\\" N
k=0, while{k<1000, 
+ d  P! J3 [/ ~
  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++},; S  l8 Z9 S6 @1 I- i

k++ 
4 D# {\\" X\\" d7 R' g2 p
},. w\\" Z, L, |3 g

s;/ @+ M/ }9 m* N+ C% F

[clock()-t0]/1000;

结果：
24002271.69189384
5.172

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

－－－－－－－－－－

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

matlab程序：

clear all
2 e) \% d( [ ^. b1 T7 L
f=@(x)x+1;* O9 X7 b8 `* N! `
a = linspace(2,2,1000);
- b, J4 x/ Q/ K( }/ T2 x& z
tic# }% c U, ?. W& ^2 s
s=0;0 d/ g; K$ i! Z. f5 I2 G8 U. O
for k=1:1000! o4 B9 j! v5 z\" O
for i=1:1000; {4 s/ a* w% V& J, I3 O
t=a(i);- M- G0 U7 o: S+ o( |8 Z# x
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);' ?$ y7 k4 P1 K6 y
end. A0 L/ I8 y, x {& }8 g
end
6 M; R' J9 k' x
s
/ N, }9 Z7 A/ Y: D) a+ |6 g$ f
toc; o% {2 V; p# z/ l8 u
' ~6 I7 E4 @- d! r8 O8 K
s =# v4 q& `8 O* w& K% y, ?
+ O3 ~* }& a9 ~, P- H
2.4002e+007 T- Q1 W; _ w+ t
^\" m# W) h: Z7 O
Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

!using["math","sys"];
, e7 p* \8 J' d8 L! b% Y4 m
mvar:
5 d7 ~! F$ s, b2 a& @
f(x)=x+1;
+ u; ]4 m4 P7 l$ g- t
main0(::a)=5 |: ]( x$ K6 ~$ H, k, C6 U
a = linspace(2,2,1000),
$ q) ^: b( w) U7 z& l
t0=clock(), s=0,
2 Y: l* F! X$ G\\" p
k=0, while{k<1000,1 v( ]9 H {# o# o* B7 G) H6 R: E/ o7 A
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++},+ s2 a C( b* k* m
k+++ I7 H% t. T# Y7 _3 m* t7 k2 A
},) y' B2 K' X- F) K* X9 g
s;# z. V5 j2 l3 f1 u
[clock()-t0]/1000;

!using["math","sys"]; 
, e7 p* \8 J' d8 L! b% Y4 m
mvar: 
5 d7 ~! F$ s, b2 a& @
f(x)=x+1; 
+ u; ]4 m4 P7 l$ g- t
main0(::a)=5 |: ]( x$ K6 ~$ H, k, C6 U

a = linspace(2,2,1000), 
$ q) ^: b( w) U7 z& l
t0=clock(), s=0, 
2 Y: l* F! X$ G\\" p
k=0, while{k<1000,1 v( ]9 H  {# o# o* B7 G) H6 R: E/ o7 A

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++},+ s2 a  C( b* k* m

k+++ I7 H% t. T# Y7 _3 m* t7 k2 A

},) y' B2 K' X- F) K* X9 g

s;# z. V5 j2 l3 f1 u

[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"];
2 j\\" k- I+ M5 ?5 D+ n; I
mvar:
9 J- e, s- D6 O
a="\&\[1000]"/4, //用字符串模拟静态数组
% }) N+ L. i% Q$ i( Y3 W& Q# f4 j0 C
t0=clock(),3 ~6 A& T- k: p0 C\\" ^% n+ I& H, L
k=0, while{k<10000,
' `1 F7 N, D) W2 D
i=0, while{i<1000, setn(a+i,2), i++}, //用函数setn(a+i,2)对静态数组元素赋值; j( _, ], x( x
k++1 S9 p\\" G7 F' P8 V) P, a
},7 w4 [; s+ B- p: y y+ A
[clock()-t0]/1000;
+ ^7 T; X8 m- a: S a
. K' w7 D( J' f) M+ ^
a="\&\[1000]"/4, //用字符串模拟静态数组
1 v3 {- X! j& l* W* `
i=0, while{i<1000, setn(a+i,2), i++}, //用函数setn(a+i,2)对静态数组元素赋值- m4 q2 J, \1 V' w9 y9 Y8 s
t0=clock(),4 Q' E$ s' }/ d9 ^7 g. \\\" F3 M
s=0,
\\" d+ g\\" t5 O2 l7 |- b. D
k=0, while{k<10000,% s! F% k4 t9 A3 E
i=0, while{i<1000, s=s+getn(a+i), i++}, //用函数getn(a+i)获得静态数组元素的值
3 P! Y8 H$ Z! n, J! I/ s: G1 I: M
k++; \8 U Y% F& W
},! @0 E7 M2 F* @3 e7 W
s;
: y: I8 H, x% x; f- ]1 l
[clock()-t0]/1000;

!using["sys"]; 
2 j\\" k- I+ M5 ?5 D+ n; I
mvar: 
9 J- e, s- D6 O
a="\&\[1000]"/4,    //用字符串模拟静态数组 
% }) N+ L. i% Q$ i( Y3 W& Q# f4 j0 C
t0=clock(),3 ~6 A& T- k: p0 C\\" ^% n+ I& H, L

k=0, while{k<10000, 
' `1 F7 N, D) W2 D
  i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值; j( _, ], x( x

k++1 S9 p\\" G7 F' P8 V) P, a

},7 w4 [; s+ B- p: y  y+ A

[clock()-t0]/1000; 
+ ^7 T; X8 m- a: S  a
. K' w7 D( J' f) M+ ^

a="\&\[1000]"/4,    //用字符串模拟静态数组 
1 v3 {- X! j& l* W* `
i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值- m4 q2 J, \1 V' w9 y9 Y8 s

t0=clock(),4 Q' E$ s' }/ d9 ^7 g. \\\" F3 M

s=0, 
\\" d+ g\\" t5 O2 l7 |- b. D
k=0, while{k<10000,% s! F% k4 t9 A3 E

i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值 
3 P! Y8 H$ Z! n, J! I/ s: G1 I: M
  k++; \8 U  Y% F& W

},! @0 E7 M2 F* @3 e7 W

s; 
: y: I8 H, x% x; f- ]1 l
[clock()-t0]/1000;

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

我也说一句

发表

使用道具举报

返回列表

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

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

浏览过的版块

QQ

电话咨询

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