查看: 5656|回复: 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+ n0 ]. K, ^) ]
tic' w& S* X1 w2 n4 N5 I
s=0;4 ]& B$ I7 M: a. g' G5 Q+ R
for k=1:100004 {4 a7 O0 r$ q, D0 ?3 B. T
for i=1:1000$ l' G- l7 r U) M* E' d
s=s+1;
) V- U5 _+ H4 F, \% i
end
; k0 E' K3 b6 D a+ r
end7 [0 x\" I* j s2 c( z0 W
s
; i) @! k3 F7 |8 ?7 E1 [! i$ k
toc
8 l0 W S0 }) `% q$ ~- F+ E
+ s* k3 P, G9 _1 h5 E
s =
: j$ X; E2 v) R6 K* S
\" d, ?1 n6 K& |1 h. r8 v7 d
100000001 d y8 B6 t1 p+ d
& r\" }( @/ Z4 M+ e/ {2 `1 v4 F6 B8 q
Elapsed time is 0.092808 seconds.

复制代码

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

Forcal没有JIT加速器，代码：

!using["sys"];& }+ x0 `$ _! c. k, V! N
mvar:8 Q d% |7 V1 s/ y: N/ R
t0=clock(),$ `, Z2 D\\" U5 h( X' f; @
s=0,
7 W4 I. D; n% e( c& t. S' C, |0 f
k=0, while{k<10000,
6 v7 ~( X/ j1 i. @, K
i=0, while{i<1000,s=s+1, i++},* i2 ~- b! _6 {/ p4 W
k++$ g\\" M: d- L- Y. i
},
2 C% t: {, } P7 S5 }6 ^
s;# r4 g# U& P) v& l1 l
[clock()-t0]/1000;

结果：
10000000.
0.704 秒

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

2、循环中稍微加点东西

Matlab代码：

clear all
+ {- m+ F' p4 \; x0 ?
tic
* @7 I9 y7 @ r
s=0;
8 b\" K, U* y' x) d/ @- J' {
for k=1:1000; Z' w# F3 ^- z5 G6 A, G$ N1 h
for i=1:1000
6 o3 C! w, N2 B3 q: _) b9 A1 k
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);
7 q* H- q+ ]) r& X+ s8 d
end
1 T9 n! H2 p3 w/ G( d3 H! R
end
$ a: J/ u- U+ E: u* H# |7 s
s
, Z5 |9 M5 t# m0 p R% \
toc3 w: j4 ~5 U/ z, a: m( F( g
: H4 E! D+ t7 R; V- x2 u3 J
s =8 {7 l2 T& s7 i7 a. g( L! Q+ N/ S# F
) |; s! D4 H1 v2 |* ?3 l; q
1.1009e+007/ K2 {# G% Y. I% `. E# x
8 l% S0 Z; r6 F; Y8 H! `' B
Elapsed time is 0.727535 seconds.

复制代码

Forcal代码：

!using["sys"];. t# m2 l0 m0 G7 J1 S
mvar:( P2 I\\" u1 X5 v) u; m; J
t0=clock(),
! k7 ]) s+ t% A/ A( j\\" h8 ?4 j
s=0,# J\\" T2 |4 ~0 Z. J
k=1, while{k<=1000,
! _& }; L$ J- s! u5 U8 N& j
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++},: L' y3 q% ?$ U' v: P\\" o+ [
k++
9 T% c9 Q9 H/ i
},
: G3 A7 F5 I2 M0 `
s;
- H$ M* Z9 C1 t
[clock()-t0]/1000;

结果：
11008953.66602346
1.031

二者效率相差已比较小。

3、再增加一个函数调用

Matlab代码：

clear all/ ^/ K, B& i; p' N. e& u8 ~4 g! m: D
f=@(x)x+1;- J! ?* Z& C/ w5 p
tic
: V. }9 }- d8 H9 Q2 P
s=0;; l' }9 s9 z7 M& j3 }/ Y1 _
for k=1:1000
9 {* z! j/ E+ P: |
for i=1:1000* E% j n( P\" _6 a
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);
$ D* V: O- x, | v
end
m# G: J& Y0 i! w
end9 V8 L. i- S7 o
s
, I4 T5 j1 }( ^+ q1 n
toc
0 \* |/ v: x7 d, n# b# X7 c
7 \ z3 V! F8 {: n
s =7 ~! S\" b6 u6 U1 [8 ~
2 C1 h4 m8 H7 q i. ?) f
5.1251e+008; L* K1 A5 v; T
/ R9 ?# W: O8 q\" c; G! G; t7 o
Elapsed time is 2.173128 seconds.

复制代码

Forcal代码：

!using["sys"];0 K' A/ D# O& V' K& Q: C6 g
f(x)=x+1;
2 `/ z2 u/ ^! w0 b0 A1 e l
mvar:
) x1 i9 S+ u; q, M( D
t0=clock(),7 J1 U# C& z# y* f4 ~- J) x4 l
s=0,
6 @* G; a4 M4 ]. R& b$ o& b
k=1, while{k<=1000,* V& H. ~! l+ |; E8 x9 b
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++},
9 E8 b5 g7 ^, X3 B4 h, V
k++, A% n0 L; v: {, l
},
+ j0 ?7 p4 P3 N) v# ~9 h\\" R
s;1 ]& X/ x) ?1 l+ l+ R
[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
/ |4 t; o. k X j
a = linspace(1,1,1000);; G3 \9 {9 p' T
tic( M7 S# n% N& b) w# O \
for k=1:10000
6 Z7 J0 H7 l% {* d' s6 l
for i=1:1000
3 C6 Y4 @( @8 K
a(i)=2;0 |4 }) a* p: B( [
end
0 V; @% Q2 n5 v+ K
end6 |3 C$ N) w\" x) E
toc8 Z! R; H# b/ Q4 Z- a
6 y: O+ c \& E5 B6 K. T# a) w
tic, N O2 P( [% O9 u0 H
s=0;; J3 D. v8 Q( A
for k=1:100001 a* y/ X9 Q6 _! ^) q T* W# Y
for i=1:1000
, _: ?! T* J# T. l
s=s+a(i); y4 p9 Q4 z- Y M) @' @. p
end
; i& m+ J8 U* F8 h
end
7 o4 S6 t- n3 L; G5 D8 E0 `8 g
s
9 z& H! T3 d/ a( }7 q& K
toc8 M\" u\" @. V t2 T; u& P\" ?' U
Elapsed time is 0.140658 seconds.
7 \5 M1 r& w, Y/ w$ k
' P; z' ~8 p) q* P$ [* H
s =
4 d& |% \* l$ M9 [
+ a. O7 P2 U' M- H/ z
20000000
4 n; Y; F q& L% c) _ U; o
$ X2 ~7 Z v- X6 x
Elapsed time is 0.091794 seconds.

复制代码

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

Forcal代码：

!using["math","sys"];7 _! i8 z2 ^/ e/ ?
mvar:
* I2 @4 B5 C0 y3 i2 F) E
main0(::a)=' c+ ?& G# C\\" S& ]( }
a = linspace(1,1,1000),8 S; X B3 ]& w) e\\" I
t0=clock(),
3 A2 h5 |# `5 \$ r& d, R) ]
k=0, while{k<10000,' @% s; a* Q* N; O\\" S
i=0, while{i<1000, a(i)=2, i++},
* n$ C' C) }% \, U
k++
- o7 j, V, \- j\\" E! ~! d
},
( J* j( ?0 ?$ q
[clock()-t0]/1000;
: x; W/ D6 [( S/ I1 G, Y6 p
1 x5 T) a7 m! n q# q
main1(::a)=
) e) Y/ c9 Q- C/ w
t0=clock(), s=0,
2 R, k- n- p8 ^& a7 A D
k=0, while{k<10000,+ n& P1 ?+ o* n' p& B# ]\\" a2 O
i=0, while{i<1000, s=s+a(i), i++},! _5 N& |# s. c, M3 F/ l4 m( G7 e; n
k++
5 q! y. y- D8 H% G. x8 O7 Z
},
! ~4 P1 G3 _+ i& Z1 n
s;
2 f: J! a, ~1 s2 |
[clock()-t0]/1000;

!using["math","sys"];7 _! i8 z2 ^/ e/ ?

mvar: 
* I2 @4 B5 C0 y3 i2 F) E
main0(::a)=' c+ ?& G# C\\" S& ]( }

a = linspace(1,1,1000),8 S; X  B3 ]& w) e\\" I

t0=clock(), 
3 A2 h5 |# `5 \$ r& d, R) ]
k=0, while{k<10000,' @% s; a* Q* N; O\\" S

i=0, while{i<1000, a(i)=2, i++}, 
* n$ C' C) }% \, U
  k++ 
- o7 j, V, \- j\\" E! ~! d
}, 
( J* j( ?0 ?$ q
[clock()-t0]/1000; 
: x; W/ D6 [( S/ I1 G, Y6 p
1 x5 T) a7 m! n  q# q

main1(::a)= 
) e) Y/ c9 Q- C/ w
t0=clock(), s=0, 
2 R, k- n- p8 ^& a7 A  D
k=0, while{k<10000,+ n& P1 ?+ o* n' p& B# ]\\" a2 O

i=0, while{i<1000, s=s+a(i), i++},! _5 N& |# s. c, M3 F/ l4 m( G7 e; n

k++ 
5 q! y. y- D8 H% G. x8 O7 Z
}, 
! ~4 P1 G3 _+ i& Z1 n
s; 
2 f: J! a, ~1 s2 |
[clock()-t0]/1000;

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

－－－－－－－－－－

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

!using["math","sys"];5 p( @# n$ {0 [/ @2 B
mvar:7 b; W' N/ C& p
main0(::a)=: i; ^- F) Z; v* X+ `: T* F
a = linspace(1,1,1000),
: E& H6 X/ J) l( g
t0=clock(),
! P/ l1 B\\" k+ J6 t# N
k=0, while{k<10000,# Q A# D* [6 }
i=0, while{i<1000, A(a,i)=2, i++},/ ?# z5 v* M0 f7 _- c. _
k+++ W3 m, p: Q( w5 Z
},# ^& b! l6 z! K( v5 F! T4 v: ]
[clock()-t0]/1000;
N- ?( j) c0 n4 E! U1 e2 I+ f: H
9 o& \! K4 e: T; [+ c1 W
main1(::a)=2 e; s/ S$ ?9 T
t0=clock(), s=0,
3 h( G+ T+ f! _5 v
k=0, while{k<10000,# {! d7 `, n' O- e+ Q9 X
i=0, while{i<1000, s=s+A(a,i), i++},
2 E, Y0 v$ Q! }+ g0 J& h
k++3 u% D Z6 }! j0 f+ k
},
; @\\" f. ]\\" N\\" F# u8 ^2 C3 U! A% F
s;
& z. ^* K& g4 Y3 f; v% O$ f
[clock()-t0]/1000;

!using["math","sys"];5 p( @# n$ {0 [/ @2 B

mvar:7 b; W' N/ C& p

main0(::a)=: i; ^- F) Z; v* X+ `: T* F

a = linspace(1,1,1000), 
: E& H6 X/ J) l( g
t0=clock(), 
! P/ l1 B\\" k+ J6 t# N
k=0, while{k<10000,# Q  A# D* [6 }

i=0, while{i<1000, A(a,i)=2, i++},/ ?# z5 v* M0 f7 _- c. _

k+++ W3 m, p: Q( w5 Z

},# ^& b! l6 z! K( v5 F! T4 v: ]

[clock()-t0]/1000; 
  N- ?( j) c0 n4 E! U1 e2 I+ f: H
 
9 o& \! K4 e: T; [+ c1 W
main1(::a)=2 e; s/ S$ ?9 T

t0=clock(), s=0, 
3 h( G+ T+ f! _5 v
k=0, while{k<10000,# {! d7 `, n' O- e+ Q9 X

i=0, while{i<1000, s=s+A(a,i), i++}, 
2 E, Y0 v$ Q! }+ g0 J& h
  k++3 u% D  Z6 }! j0 f+ k

}, 
; @\\" f. ]\\" N\\" F# u8 ^2 C3 U! A% F
s; 
& z. ^* K& g4 Y3 f; v% O$ f
[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
0 u0 O: S' c\" \
f=@(x)x+1; [! Y4 X( C4 K! U: J
a = linspace(2,2,1000);+ T! e5 _9 i8 Z% j5 c5 [. R
tic
) P ^) k# D: G$ e4 x8 s$ ^
s=0;0 \, F- m: ~9 P/ {
for k=1:1000' Z- h3 c+ u* N1 ^
for i=1:1000
, t1 n% w- a4 v( F S7 |
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));4 ]- M9 Y, R1 _- N& ]
end
8 x! B4 j8 a% u1 B
end
& E- q$ z: Q5 i% H
s
- g! G' A7 F. Z9 C6 W\" s/ K
toc
. o% U1 h8 T. k& @
) m2 |% l8 N, j3 B; x
s =
) A E7 B! P2 D% T. a# R
1 Y/ [; c7 T. e }/ t! c% V: Z# y
2.4002e+007
2 x( S3 ~. l9 x
+ Y; m5 g) a- N8 b5 H- n$ d
Elapsed time is 2.252997 seconds.

复制代码

Forcal代码：

!using["math","sys"];% x! O$ j: Q- B5 F( N) l7 _
mvar:
. [4 k\\" a0 W% b1 k. p# n
f(x)=x+1;
F5 \- [) k1 W7 I8 I) w
main0(::a)=
+ _* T6 D9 T# P/ i
a = linspace(2,2,1000),
2 D( M& w m4 C& y4 q
t0=clock(), s=0, E Z; H8 X. {0 R: ^
k=0, while{k<1000,
' w; T& ~. d3 \: k' w
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\\" `) g$ u4 T. f
k++
) M8 K+ I z8 r2 I& {- D$ J% }
},
+ z) S\\" ?+ Q( p6 G
s;
5 [' W: X) M\\" I- y$ k9 d5 U; n
[clock()-t0]/1000;

!using["math","sys"];% x! O$ j: Q- B5 F( N) l7 _

mvar: 
. [4 k\\" a0 W% b1 k. p# n
f(x)=x+1; 
  F5 \- [) k1 W7 I8 I) w
main0(::a)= 
+ _* T6 D9 T# P/ i
a = linspace(2,2,1000), 
2 D( M& w  m4 C& y4 q
t0=clock(), s=0,  E  Z; H8 X. {0 R: ^

k=0, while{k<1000, 
' w; T& ~. d3 \: k' w
  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\\" `) g$ u4 T. f
  k++ 
) M8 K+ I  z8 r2 I& {- D$ J% }
}, 
+ z) S\\" ?+ Q( p6 G
s; 
5 [' W: X) M\\" I- y$ k9 d5 U; n
[clock()-t0]/1000;

结果：
24002271.69189384
5.172

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

－－－－－－－－－－

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

matlab程序：

clear all& r3 I: U3 A5 u+ ]. a+ x, O
f=@(x)x+1;
a = linspace(2,2,1000);9 V. l. Q3 F: V4 [# D5 w
tic% |7 G$ U4 C2 J2 A! }7 N6 v
s=0;% B l. K: U8 N* n5 N9 \6 T
for k=1:1000/ |) U) ^9 b& }! g
for i=1:1000
1 s; u0 g' s9 A4 v$ ], F1 h
t=a(i);7 W0 M( R: X' U5 G( l4 u
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);
! C7 ^* H* f4 j4 e$ L
end3 f H$ ~# o% H- Q+ o' _3 x9 _
end
2 D& K9 A. k\" Y6 I0 B8 F9 |\" \& E
s
* Z3 h6 T2 t/ Y( i6 Z
toc
7 I% \9 }/ ?: e\" O/ r+ `
; F6 S9 q9 W* h* E
s =% D% y- F) d5 I* U( o! L* Y
2 m$ _+ o6 r& R
2.4002e+007
, D2 G- `# z1 j' C+ C8 G6 _& l
- {+ H# x, l( X$ R
Elapsed time is 2.043909 seconds.

复制代码

Forcal程序：

!using["math","sys"];( t( t. W/ _( v2 f
mvar:
$ G' F: E) b/ c5 n4 J
f(x)=x+1;
6 q( c. R. m2 g0 C8 D
main0(::a)=
8 t$ Y8 |\\" S' M' [/ t! p) `
a = linspace(2,2,1000),4 r/ f% H; b& E: K
t0=clock(), s=0,
\\" H$ V: ]; D0 U1 t
k=0, while{k<1000,0 |% s: D! u0 I. s, 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++},
9 [1 p8 H) o3 o% J c
k++
, R k! E) I4 d( y( R
},
, G& B& S8 J) f7 c' ]7 M1 v
s;8 S z# U5 q9 R8 I6 t
[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"];! {* D# I0 G7 |\\" w- c
mvar:
. m0 R. N' F2 {+ }\\" k0 T2 i, [ y
a="\&\[1000]"/4, //用字符串模拟静态数组( P$ D: _5 A0 B: o% c
t0=clock(),
) A: W) l/ f5 @+ A7 S. t
k=0, while{k<10000,
. Q$ J/ o/ {* C ^! r
i=0, while{i<1000, setn(a+i,2), i++}, //用函数setn(a+i,2)对静态数组元素赋值
! i) M* f$ |. I! P
k+++ n6 D7 C\\" _: C2 O/ i* f
},7 C) T: b& K4 N9 X) W& o
[clock()-t0]/1000;
\\" ?5 C6 p$ T9 v: K9 L
0 J+ P4 i. ^& T) z2 |
a="\&\[1000]"/4, //用字符串模拟静态数组8 ]8 c2 f3 c* }\\" I V6 d1 x
i=0, while{i<1000, setn(a+i,2), i++}, //用函数setn(a+i,2)对静态数组元素赋值
1 x% j: G$ q% C* v
t0=clock(),
\\" M1 `* c! P4 g% T1 L6 c( }% h
s=0,
! I* }! l5 u+ [) h. K, F
k=0, while{k<10000,6 A3 Q\\" {: Y6 g8 z) f) ~& D
i=0, while{i<1000, s=s+getn(a+i), i++}, //用函数getn(a+i)获得静态数组元素的值
9 I4 d. g( j5 A, ?' @( {+ I- g8 n
k++1 p( e* ?) z8 s- Q* y
},1 e' `1 F* x' R8 ` @
s;7 Y8 p4 b0 B/ c& t5 r/ E
[clock()-t0]/1000;

!using["sys"];! {* D# I0 G7 |\\" w- c

mvar: 
. m0 R. N' F2 {+ }\\" k0 T2 i, [  y
a="\&\[1000]"/4,    //用字符串模拟静态数组( P$ D: _5 A0 B: o% c

t0=clock(), 
) A: W) l/ f5 @+ A7 S. t
k=0, while{k<10000, 
. Q$ J/ o/ {* C  ^! r
  i=0, while{i<1000, setn(a+i,2), i++},   //用函数setn(a+i,2)对静态数组元素赋值 
! i) M* f$ |. I! P
  k+++ n6 D7 C\\" _: C2 O/ i* f

},7 C) T: b& K4 N9 X) W& o

[clock()-t0]/1000; 
\\" ?5 C6 p$ T9 v: K9 L
0 J+ P4 i. ^& T) z2 |

a="\&\[1000]"/4,    //用字符串模拟静态数组8 ]8 c2 f3 c* }\\" I  V6 d1 x

i=0, while{i<1000, setn(a+i,2), i++},      //用函数setn(a+i,2)对静态数组元素赋值 
1 x% j: G$ q% C* v
t0=clock(), 
\\" M1 `* c! P4 g% T1 L6 c( }% h
s=0, 
! I* }! l5 u+ [) h. K, F
k=0, while{k<10000,6 A3 Q\\" {: Y6 g8 z) f) ~& D

i=0, while{i<1000, s=s+getn(a+i), i++},  //用函数getn(a+i)获得静态数组元素的值 
9 I4 d. g( j5 A, ?' @( {+ I- g8 n
  k++1 p( e* ?) z8 s- Q* y

},1 e' `1 F* x' R8 `  @

s;7 Y8 p4 b0 B/ c& t5 r/ E

[clock()-t0]/1000;

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

我也说一句

发表

使用道具举报

返回列表

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

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

浏览过的版块

QQ

电话咨询

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