Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 1
EECS 252 Graduate Computer Architecture
Lec. 12: Vector Computers
Krste Asanovic
(krste@[Link])
Computer Science and Artificial Intelligence Laboratory
Massachusetts Institute of Technology
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 2
Supercomputers
Definition of a supercomputer:
• Fastest machine in world at given task
• A device to turn a compute-bound problem into an
I/O bound problem
• Any machine costing $30M+
• Any machine designed by Seymour Cray
CDC6600 (Cray, 1964) regarded as first supercomputer
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 3
Supercomputer Applications
Typical application areas
• Military research (nuclear weapons, cryptography)
• Scientific research
• Weather forecasting
• Oil exploration
• Industrial design (car crash simulation)
All involve huge computations on large data sets
In 70s-80s, Supercomputer Vector Machine
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 4
Vector Supercomputers
Epitomized by Cray-1, 1976:
Scalar Unit + Vector Extensions
• Load/Store Architecture
• Vector Registers
• Vector Instructions
• Hardwired Control
• Highly Pipelined Functional Units
• Interleaved Memory System
• No Data Caches
• No Virtual Memory
Krste
CS 252 Feb. 27, 2006
Cray-1 (1976) Lecture 12, Slide 5
Krste
CS 252 Feb. 27, 2006
Cray-1 (1976) Lecture 12, Slide 6
V0 Vi V. Mask
V1
V2 Vj
64 Element V3 V. Length
Vector Registers V4 Vk
Single Port V5
V6
Memory V7
FP Add
S0 Sj FP Mul
16 banks of ( (Ah) + j k m ) S1
S2 Sk FP Recip
64-bit words Si S3
(A0) 64 S4 Si Int Add
+ Tjk S5
T Regs S6
8-bit SECDED S7
Int Logic
Int Shift
A0
80MW/sec data ( (Ah) + j k m ) A1 Pop Cnt
A2
load/store Ai A3
Aj
(A0) 64 A4 Ak Addr Add
Bjk A5
Ai
320MW/sec B Regs A6 Addr Mul
A7
instruction
buffer refill NIP CIP
64-bitx16
LIP
4 Instruction Buffers
memory bank cycle 50 ns processor cycle 12.5 ns (80MHz)
Krste
CS 252 Feb. 27, 2006
Vector Programming Model Lecture 12, Slide 7
Scalar Registers Vector Registers
r15 v15
r0 v0
[0] [1] [2] [VLRMAX-1]
Vector Length Register VLR
v1
Vector Arithmetic v2
Instructions + + + + + +
ADDV v3, v1, v2 v3
[0] [1] [VLR-1]
Vector Load and Vector Register
Store Instructions v1
LV v1, r1, r2
Memory
Base, r1 Stride, r2
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 8
Vector Code Example
# C code # Scalar Code # Vector Code
for (i=0; i<64; i++) LI R4, 64 LI VLR, 64
C[i] = A[i] + B[i]; loop: LV V1, R1
L.D F0, 0(R1) LV V2, R2
L.D F2, 0(R2) ADDV.D V3, V1, V2
ADD.D F4, F2, F0 SV V3, R3
S.D F4, 0(R3)
DADDIU R1, 8
DADDIU R2, 8
DADDIU R3, 8
DSUBIU R4, 1
BNEZ R4, loop
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 9
Vector Instruction Set Advantages
• Compact
– one short instruction encodes N operations
• Expressive, tells hardware that these N operations:
– are independent
– use the same functional unit
– access disjoint registers
– access registers in the same pattern as previous instructions
– access a contiguous block of memory (unit-stride load/store)
– access memory in a known pattern (strided load/store)
• Scalable
– can run same object code on more parallel pipelines or lanes
Krste
CS 252 Feb. 27, 2006
Vector Arithmetic Execution Lecture 12, Slide 10
• Use deep pipeline (=> fast clock) V V V
to execute element operations 1 2 3
• Simplifies control of deep pipeline
because elements in vector are
independent (=> no hazards!)
Six stage multiply pipeline
V3 <- v1 * v2
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 11
Vector Memory System
Cray-1, 16 banks, 4 cycle bank busy time, 12 cycle latency
• Bank busy time: Cycles between accesses to same bank
Base Stride
Vector Registers
Address
Generator +
0 1 2 3 4 5 6 7 8 9 A B C D E F
Memory Banks
Krste
CS 252 Feb. 27, 2006
Vector Instruction Execution Lecture 12, Slide 12
ADDV C,A,B
Execution using Execution using
one pipelined four pipelined
functional unit functional units
A[6] B[6] A[24] B[24] A[25] B[25] A[26] B[26] A[27] B[27]
A[5] B[5] A[20] B[20] A[21] B[21] A[22] B[22] A[23] B[23]
A[4] B[4] A[16] B[16] A[17] B[17] A[18] B[18] A[19] B[19]
A[3] B[3] A[12] B[12] A[13] B[13] A[14] B[14] A[15] B[15]
C[2] C[8] C[9] C[10] C[11]
C[1] C[4] C[5] C[6] C[7]
C[0] C[0] C[1] C[2] C[3]
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 13
Vector Unit Structure
Functional Unit
Vector
Registers
Elements Elements Elements Elements
0, 4, 8, … 1, 5, 9, … 2, 6, 10, … 3, 7, 11, …
Lane
Memory Subsystem
Krste
CS 252 Feb. 27, 2006
T0 Vector Microprocessor (1995) Lecture 12, Slide 14
Vector register Lane
elements striped
over lanes
[24][25] [26] [27] [28] [29] [30] [31]
[16][17] [18] [19] [20] [21] [22] [23]
[8] [9] [10] [11] [12] [13] [14] [15]
[0] [1] [2] [3] [4] [5] [6] [7]
Krste
Vector Memory-Memory versus CS 252 Feb. 27, 2006
Lecture 12, Slide 15
Vector Register Machines
• Vector memory-memory instructions hold all vector operands
in main memory
• The first vector machines, CDC Star-100 (‘73) and TI ASC (‘71),
were memory-memory machines
• Cray-1 (’76) was first vector register machine
Vector Memory-Memory Code
Example Source Code ADDV C, A, B
for (i=0; i<N; i++) SUBV D, A, B
{
C[i] = A[i] + B[i]; Vector Register Code
D[i] = A[i] - B[i]; LV V1, A
} LV V2, B
ADDV V3, V1, V2
SV V3, C
SUBV V4, V1, V2
SV V4, D
Krste
Vector Memory-Memory vs. CS 252 Feb. 27, 2006
Lecture 12, Slide 16
Vector Register Machines
• Vector memory-memory architectures (VMMA) require
greater main memory bandwidth, why?
– All operands must be read in and out of memory
• VMMAs make if difficult to overlap execution of
multiple vector operations, why?
– Must check dependencies on memory addresses
• VMMAs incur greater startup latency
– Scalar code was faster on CDC Star-100 for vectors < 100 elements
– For Cray-1, vector/scalar breakeven point was around 2 elements
Apart from CDC follow-ons (Cyber-205, ETA-10) all
major vector machines since Cray-1 have had vector
register architectures
(we ignore vector memory-memory from now on)
Krste
Automatic Code Vectorization CS 252 Feb. 27, 2006
Lecture 12, Slide 17
for (i=0; i < N; i++)
C[i] = A[i] + B[i];
Scalar Sequential Code Vectorized Code
load load load
Iter. 1 load load load
add Time add add
store store store
load
Iter. Iter.
load 1 2 Vector Instruction
Iter. 2
add
Vectorization is a massive compile-time
reordering of operation sequencing
store requires extensive loop dependence
analysis
Krste
Vector Stripmining CS 252 Feb. 27, 2006
Lecture 12, Slide 18
Problem: Vector registers have finite length
Solution: Break loops into pieces that fit into vector
registers, “Stripmining”
ANDI R1, N, 63 # N mod 64
for (i=0; i<N; i++) MTC1 VLR, R1 # Do remainder
C[i] = A[i]+B[i]; loop:
LV V1, RA
A B C
DSLL R2, R1, 3 # Multiply by 8
+ Remainder DADDU RA, RA, R2 # Bump pointer
LV V2, RB
DADDU RB, RB, R2
+ 64 elements ADDV.D V3, V1, V2
SV V3, RC
DADDU RC, RC, R2
DSUBU N, N, R1 # Subtract elements
+ LI R1, 64
MTC1 VLR, R1 # Reset full length
BGTZ N, loop # Any more to do?
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 19
Vector Instruction Parallelism
Can overlap execution of multiple vector instructions
– example machine has 32 elements per vector register and 8 lanes
Load Unit Multiply Unit Add Unit
load
mul
add
time
load
mul
add
Instruction
issue
Complete 24 operations/cycle while issuing 1 short instruction/cycle
Krste
Vector Chaining CS 252 Feb. 27, 2006
Lecture 12, Slide 20
• Vector version of register bypassing
– introduced with Cray-1
V V V V V
LV v1
1 2 3 4 5
MULV v3,v1,v2
ADDV v5, v3, v4
Chain Chain
Load
Unit
Mult. Add
Memory
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 21
Vector Chaining Advantage
• Without chaining, must wait for last element of result to
be written before starting dependent instruction
Load
Mul
Time Add
• With chaining, can start dependent instruction as soon
as first result appears
Load
Mul
Add
Krste
Vector Startup CS 252 Feb. 27, 2006
Lecture 12, Slide 22
Two components of vector startup penalty
– functional unit latency (time through pipeline)
– dead time or recovery time (time before another vector
instruction can start down pipeline)
Functional Unit Latency
R X X X W
R X X X W First Vector Instruction
R X X X W
R X X X W
R X X X W
Dead Time
R X X X W
R X X X W
R X X X W
Dead Time R X X X W Second Vector Instruction
R X X X W
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 23
Dead Time and Short Vectors
No dead time
4 cycles dead time T0, Eight lanes
No dead time
100% efficiency with 8 element
vectors
64 cycles active
Cray C90, Two lanes
4 cycle dead time
Maximum efficiency 94%
with 128 element vectors
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 24
Vector Scatter/Gather
Want to vectorize loops with indirect accesses:
for (i=0; i<N; i++)
A[i] = B[i] + C[D[i]]
Indexed load instruction (Gather)
LV vD, rD # Load indices in D vector
LVI vC, rC, vD # Load indirect from rC base
LV vB, rB # Load B vector
ADDV.D vA, vB, vC # Do add
SV vA, rA # Store result
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 25
Vector Scatter/Gather
Scatter example:
for (i=0; i<N; i++)
A[B[i]]++;
Is following a correct translation?
LV vB, rB # Load indices in B vector
LVI vA, rA, vB # Gather initial A values
ADDV vA, vA, 1 # Increment
SVI vA, rA, vB # Scatter incremented values
Krste
CS 252 Feb. 27, 2006
Vector Conditional Execution Lecture 12, Slide 26
Problem: Want to vectorize loops with conditional code:
for (i=0; i<N; i++)
if (A[i]>0) then
A[i] = B[i];
Solution: Add vector mask (or flag) registers
– vector version of predicate registers, 1 bit per element
…and maskable vector instructions
– vector operation becomes NOP at elements where mask bit is clear
Code example:
CVM # Turn on all elements
LV vA, rA # Load entire A vector
SGTVS.D vA, F0 # Set bits in mask register where A>0
LV vA, rB # Load B vector into A under mask
SV vA, rA # Store A back to memory under mask
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 27
Masked Vector Instructions
Simple Implementation Density-Time Implementation
– execute all N operations, turn off – scan mask vector and only execute
result writeback according to mask elements with non-zero masks
M[7]=1 A[7] B[7] M[7]=1
M[6]=0 A[6] B[6] M[6]=0 A[7] B[7]
M[5]=1 A[5] B[5] M[5]=1
M[4]=1 A[4] B[4] M[4]=1
M[3]=0 A[3] B[3] M[3]=0 C[5]
M[2]=0 C[4]
M[1]=1
M[2]=0 C[2]
M[0]=0
M[1]=1 C[1] C[1]
Write data port
M[0]=0 C[0]
Write Enable Write data port
Krste
Compress/Expand Operations CS 252 Feb. 27, 2006
Lecture 12, Slide 28
• Compress packs non-masked elements from one
vector register contiguously at start of destination
vector register
– population count of mask vector gives packed vector length
• Expand performs inverse operation
M[7]=1 A[7] A[7] A[7] M[7]=1
M[6]=0 A[6] A[5] B[6] M[6]=0
M[5]=1 A[5] A[4] A[5] M[5]=1
M[4]=1 A[4] A[1] A[4] M[4]=1
M[3]=0 A[3] A[7] B[3] M[3]=0
M[2]=0 A[2] A[5] B[2] M[2]=0
M[1]=1 A[1] A[4] A[1] M[1]=1
M[0]=0 A[0] A[1] B[0] M[0]=0
Compress Expand
Used for density-time conditionals and also for general
selection operations
Krste
CS 252 Feb. 27, 2006
Vector Reductions Lecture 12, Slide 29
Problem: Loop-carried dependence on reduction variables
sum = 0;
for (i=0; i<N; i++)
sum += A[i]; # Loop-carried dependence on sum
Solution: Re-associate operations if possible, use binary
tree to perform reduction
# Rearrange as:
sum[0:VL-1] = 0 # Vector of VL partial sums
for(i=0; i<N; i+=VL) # Stripmine VL-sized chunks
sum[0:VL-1] += A[i:i+VL-1]; # Vector sum
# Now have VL partial sums in one vector register
do {
VL = VL/2; # Halve vector length
sum[0:VL-1] += sum[VL:2*VL-1] # Halve no. of partials
} while (VL>1)
Krste
A Modern Vector Super: NEC SX-6 (2003) CS 252 Feb. 27, 2006
Lecture 12, Slide 30
• CMOS Technology
– 500 MHz CPU, fits on single chip
– SDRAM main memory (up to 64GB)
• Scalar unit
– 4-way superscalar with out-of-order and speculative
execution
– 64KB I-cache and 64KB data cache
• Vector unit
– 8 foreground VRegs + 64 background VRegs (256x64-bit
elements/VReg)
– 1 multiply unit, 1 divide unit, 1 add/shift unit, 1 logical unit,
1 mask unit
– 8 lanes (8 GFLOPS peak, 16 FLOPS/cycle)
– 1 load & store unit (32x8 byte accesses/cycle)
– 32 GB/s memory bandwidth per processor
• SMP structure
– 8 CPUs connected to memory through crossbar
– 256 GB/s shared memory bandwidth (4096 interleaved
banks)
Krste
CS 252 Feb. 27, 2006
Lecture 12, Slide 31
Multimedia Extensions
• Very short vectors added to existing ISAs for micros
• Usually 64-bit registers split into 2x32b or 4x16b or 8x8b
• Newer designs have 128-bit registers (Altivec, SSE2)
• Limited instruction set:
– no vector length control
– no strided load/store or scatter/gather
– unit-stride loads must be aligned to 64/128-bit boundary
• Limited vector register length:
– requires superscalar dispatch to keep multiply/add/load units busy
– loop unrolling to hide latencies increases register pressure
• Trend towards fuller vector support in microprocessors