ә
ғ ү
CRAY C90 ң ( )
ө ң қ ә.
ү
, ғ
- ұ - қққ қ . ғ , ғ, , ү , қ, ө ң қғ ө .
қ өң ү
ң қ , ң қ өң . ұ ә ң қ ғ . ү ә - қ. ғ ү-ұ ң ң .
. ң қ ң ө өң қ; ұ ғ , ғ , қ ң өң ң ұқ қ қ өң .
қ ө ү ғ :
Do i=1,n
A(i) = A(i) + s
EndDo
ғ :
Do i=1,n
A(i) = A(i-1)+s
End Do
ү ()
қ ү ө ү, ң :
ә(ңғ ).
ұқ қ ң қ :
Do i=1,n
ij = FUNC(i)
A(i) = A(i)+B(ij)
End Do
ғ ң ұ ү ү қ қ.
ң қ:
Do i=1,n
CALL SUBR(A,B)
End Do
CRAY C90 ң ( ). ң қ ө ө қ қ ң қ. ұ ң ұ ң ң ө қ. , ң - .
ұ ң ә қ ғ . , ң ө ң ө 1 ә 2 . -, 1 2, . , , , , .
|
|
-, , , , " " 2 1.
, . : l A(i) B. l B(j), B(j+1) ( .1). A .
.1
, A, l , B(j) , , j , B. , , , CRAY Y-MP C90.
, , . , , i0, A B , 63 ( 64, B(j) ). A(ibeg+i*63), i=0... ibeg=1...63 , B, B ! : 63 , , B , A(ibeg+i*63), i=1... ibeg . , 30 (!) . - 2, .
( B):
A(i)
B(j) p
B(j) , B(j+1), j=j+1,
p, i=i+1
|
|
:
(p[ibeg]) 63 A - ibeg=1,63:
A(i*63+ibeg)
B(j) p[ibeg]
i=i+1, j=j+bstep
, ( ) - 128 , . , , , , , .
. , . :
Do i=1,n
a(i)=b(i)*s + c(i)
End Do
.1 CRAY Y-MP C90 , .. n ( - ).
, Mflop/s | , Mflop/s | ||
7.0 | 413.2 | ||
14.0 | 548.0 | ||
27.6 | 491.0 | ||
100.5 | 659.2 | ||
181.9 | 720.4 | ||
301.0 | 768.0 | ||
433.7 | 802.0 | ||
364.3 |
.1 CRAY Y-MP C90 ai=bi*s+ci
CRAY Y-MP . CRAY Y-MP C90 8 , - 8 , 16 . , 8*8=64, . , , . :
Do i=1,n*k,k
a(i)=b(i)*s + c(i)
End Do
k, .. , ai=bi*s+ci n . ( ) 2.
100 | 1000 | 12800 | |
240.3 | 705.2 | 805.1 | |
220.4 | 444.6 | 498.5 | |
172.9 | 274.6 | 280.1 | |
108.1 | 142.8 | 147.7 | |
71.7 | 84.5 | 86.0 | |
41.0 | 44.3 | 38.0 | |
22.1 | 25.7 | 22.3 | |
21.2 | 20.6 | 20.3 |
.2 : CRAY Y-MP C90 .
, . , "" ( , , - ). , . , , :
Do i=1,n
Do j=1,n
Do k=1,n
|
|
X(i,j,k) = X(i,j,k)+P(k,i)*Y(k,j)
End Do
End Do
End Do
, X. , :
DIMENSION X(40,40,100)
" ", . , X(i,j,k) X(i,j,k+1) 40*40=1600=25*64, .. . , , :
DIMENSION X(41,41,100),
. - .
Do j=1,n
Do i=1,n
XYZ(IX(i),j) = XYZ(IX(i),j)+P(i,j)*Y(i,j)
End Do
End Do
, XYZ , (, IX(i) i) (, IX(i) i).
, , , . CRAY Y-MP , , . , , , Ai = Bi*s+Ci.
, Ai = Bi*Ci+Di, . , (memory bottleneck). , , , . 3 , B, C, D, .
, Mflop/s | |
57.0 | |
278.3 | |
435.3 | |
445.0 |
.3 CRAY Y-MP C90 ai=bi*ci+di
, . . , / . , , , . :
Do j=1,120
Do i=1,n
DP(i) = DP(i) + s*P(i,j-1) + t*P(i,j)
End Do
End Do
j DP(i), P(i,j-1), P(i,j) - DP(i), 120*3=360 120 . j :
|
|
Do j=1,120,2
Do i=1,n
DP(i) = DP(i)+s*P(i,j-1)+t*P(i,j)+s*P(i,j)+t*P(i,j+1)
End Do
End Do
60- DP(i), P(i,j-1), P(i,j), P(i,j+1) . , 60*4=240 60 . "" , .
, , . - , , . , CRAY Y-MP C90 : 8 128 . , , . 4 CRAY Y-MP C90 .
64 | 128 | 12800 | |
464.4 | 612.9 | 749.0 | |
591.4 | 731.6 | 730.1 | |
639.3 | 780.7 | 752.5 | |
675.3 | 807.7 | 786.8 | |
.4 CRAY Y-MP C90 .
, . , , . , - 5.
. . , , -, , -, , . 5.
ai=bi+ci | ai=bi*ci | ai=bi/ci | ai=s/bi+t | ai=s/bi*t | |
35.5 | 41.9 | 24.8 | 45.7 | 46.1 | |
202.9 | 198.0 | 88.4 | 197.4 | 166.5 | |
343.8 | 341.2 | 117.2 | 283.8 | 215.9 | |
373.1 | 376.8 | 120.0 | 297.0 | 222.5 |
.5 CRAY Y-MP C90 .
, , , , . , . , . , .. , .
.