Vector
Vector
5381/7381
Review
Speculation: Out-of-order execution, In-order commit (reorder
•
buffer+rename registers)=>precise exceptions
• Branch Prediction
– Branch History Table: 2 bits for loop accuracy
– Recently executed branches correlated with next branch?
– Branch Target Buffer: include branch address & prediction
– Predicated Execution can reduce number of branches, number of mispredicted
branches
• Software Pipelining
– Symbolic loop unrolling (instructions from different iterations) to optimize pipeline with
little code expansion, little overhead
• Superscalar and VLIW(“EPIC”): CPI < 1 (IPC > 1)
– Dynamic issue vs. Static issue
– More instructions issue at same time => larger hazard penalty
– # independent instructions = # functional units X latency
5381/7381
Review: Theoretical Limits to ILP?
(Figure 4.48, Page 332)
60
56
Perfect disambiguation (HW), 1K Selective
Prediction, 16 entry return, 64 registers, 52
50 issue as many as window
FP: 8 - 45
47
45
Instruction issues per cycle
40
35
34
30
IPC
22 22
20
Integer: 6 - 12
15 15
14
17 16
15 14
13
12 12 11 11 12
10 10 10 10
9 8 9 8 9 9
10 8
7
6 6 6 6
4 5
4 4 4
3 2 3 3 3 3
Program
5381/7381
Infinite 256 128 64 32 16 8 4
Review: Instructon Level Parallelism
5381/7381
Problems with conventional approach
SCALAR VECTOR
(1 operation) (N operations)
r1 r2 v1 v2
+ +
r3 v3 vector
length
5381/7381
Operation & Instruction Count:
RISC v. Vector Processor
(from F. Quintana, U. Barcelona.)
Spec92fp Operations (Millions) Instructions (M)
Program RISC Vector R / V RISC Vector R / V
swim256 11595 1.1x115 0.8142x
hydro2d 5840 1.4x 58 0.8 71x
nasa7 6941 1.7x 69 2.2 31x
su2cor 5135 1.4x 51 1.8 29x
tomcatv 1510 1.4x 15 1.3 11x
wave5 2725 1.1x 27 7.2 4x
mdljdp2 3252 0.6x 32 15.8 2x
5381/7381
Components of Vector Processor
• Vector Register: fixed length bank holding a single vector
– has at least 2 read and 1 write ports
– typically 8-32 vector registers, each holding 64-128 64-bit
elements
• Vector Functional Units (FUs): fully pipelined, start new
operation every clock
– typically 4 to 8 FUs: FP add, FP mult, FP reciprocal (1/X),
integer add, logical, shift; may have multiple of same unit
• Vector Load-Store Units (LSUs): fully pipelined unit to
load or store a vector; may have multiple LSUs
• Scalar registers: single element for FP scalar or address
• Cross-bar to connect FUs , LSUs, registers
5381/7381
“DLXV” Vector Instructions
Instr. Operands Operation Comment
• ADDV V1,V2,V3 V1=V2+V3 vector + vector
• ADDSV V1,F0,V2 V1=F0+V2 scalar + vector
• MULTV V1,V2,V3 V1=V2xV3 vector x vector
• MULSV V1,F0,V2 V1=F0xV2 scalar x vector
• LV V1,R1 V1=M[R1..R1+63] load, stride=1
• LVWS V1,R1,R2 V1=M[R1..R1+63*R2] load, stride=R2
• LVI V1,R1,V2 V1=M[R1+V2i,i=0..63] indir.("gather")
• CeqV VM,V1,V2 VMASKi = (V1i=V2i)? comp. setmask
• MOV VLR,R1 Vec. Len. Reg. = R1 set vector length
• MOV VM,R1 Vec. Mask = R1 set vector mask
5381/7381
Memory operations
• Load/store operations move groups of data
between registers and memory
• Three types of addressing
– Unit stride
» Fastest
– Non-unit (constant) stride
– Indexed (gather-scatter)
» Vector equivalent of register indirect
» Good for sparse arrays of data
» Increases number of programs that vectorize
5381/7381
32
DAXPY (Y = a * X + Y)
Assuming vectors X, Y LD F0,a ;load scalar a
are length 64 LV V1,Rx ;load vector X
Scalar vs. Vector MULTS V2,F0,V1 ;vector-scalar mult.
LV V3,Ry ;load vector Y
ADDV V4,V2,V3 ;add
LD F0,a SV Ry,V4 ;store the result
ADDI R4,Rx,#512 ;last address to load 578 (2+9*64) vs.
loop: LD F2, 0(Rx) ;load X(i) 321 (1+5*64) ops (1.8X)
MULTD F2,F0,F2 ;a*X(i)
578 (2+9*64) vs.
LD F4, 0(Ry) ;load Y(i) 6 instructions (96X)
ADDD F4,F2, F4 ;a*X(i) + Y(i)
SD F4 ,0(Ry) ;store into Y(i) 64 operation vectors +
ADDI Rx,Rx,#8 ;increment index to X no loop overhead
ADDI Ry,Ry,#8 ;increment index to Y also 64X fewer pipeline
SUB R20,R4,Rx ;compute bound hazards
5381/7381
BNZ R20,loop ;check if done
Example Vector Machines
• Machine Year ClockRegs Elements FUs LSUs
• Cray 1 1976 80 MHz8 64 6 1
• Cray XMP 1983 120 MHz8 64 8 2 L, 1 S
• Cray YMP 1988 166 MHz8 64 8 2 L, 1 S
• Cray C-90 1991 240 MHz8 128 8 4
• Cray T-90 1996 455 MHz8 128 8 4
• Conv. C-1 1984 10 MHz8 128 4 1
• Conv. C-4 1994 133 MHz16 128 3 1
• Fuj. VP200 1982 133 MHz8-256 32-1024 3 2
• Fuj. VP300 1996 100 MHz8-256 32-1024 3 2
• NEC SX/2 1984 160 MHz8+8K 256+var 16 8
• NEC SX/3 1995 400 MHz8+8K 256+var 16 85381/7381
Vector Linpack Performance
(MFLOPS)
Machine Year Clock 100x100 1kx1k Peak(Procs)
• Cray 1 1976 80 MHz 12 110 160(1)
• Cray XMP 1983 120 MHz 121 218 940(4)
• Cray YMP 1988 166 MHz 150 307 2,667(8)
• Cray C-90 1991 240 MHz 387 902 15,238(16)
• Cray T-90 1996 455 MHz 705 1603 57,600(32)
• Conv. C-1 1984 10 MHz 3 -- 20(1)
• Conv. C-4 1994 135 MHz 160 2531 3240(4)
• Fuj. VP200 1982 133 MHz 18 422 533(1)
• NEC SX/2 1984 166 MHz 43 885 1300(1)
• NEC SX/3 1995 400 MHz 368 2757 25,600(4)
5381/7381
Vector Surprise
• Use vectors for inner loop parallelism (no surprise)
– One dimension of array: A[0, 0], A[0, 1], A[0, 2], ...
– think of machine as, say, 32 vector regs each with 64 elements
– 1 instruction updates 64 elements of 1 vector register
• and for outer loop parallelism!
– 1 element from each column: A[0,0], A[1,0], A[2,0], ...
– think of machine as 64 “virtual processors” (VPs)
each with 32 scalar registers! (- multithreaded processor)
– 1 instruction updates 1 scalar register in 64 VPs
• Hardware identical, just 2 compiler perspectives
5381/7381
Virtual Processor Vector Model
• Vector operations are SIMD
(single instruction multiple data)operations
• Each element is computed by a virtual
processor (VP)
• Number of VPs given by vector length
– vector control register
5381/7381
Vector Architectural State
Virtual Processors ($vlr)
VP0 VP1 VP$vlr-1
General vr0
vr1 Control
Purpose Registers
Registers
vr31 vcr0
$vdw bits vcr1
vf0
Flag vf1
Registers vcr31
(32) vf31
32 bits
1 bit
5381/7381
Vector Implementation
• Vector register file
– Each register is an array of elements
– Size of each register determines maximum
vector length
– Vector length register determines vector length
for a particular operation
• Multiple parallel execution units = “lanes”
(sometimes called “pipelines” or “pipes”)
5381/7381
33
Vector Terminology:
4 lanes, 2 vector functional units
(Vector
Functional
Unit)
5381/7381
34
Vector Execution Time
• Time = f(vector length, data dependicies, struct. hazards)
• Initiation rate: rate that FU consumes vector elements
(= number of lanes; usually 1 or 2 on Cray T-90)
• Convoy: set of vector instructions that can begin
execution in same clock (no struct. or data hazards)
• Chime: approx. time for a vector operation
• m convoys take m chimes; if each vector length is n,
then they take approx. m x n clock cycles (ignores
overhead; good approximization for long vectors)
1: LV V1,Rx ;load vector X
4 conveys, 1 lane, VL=64
2: MULV V2,F0,V1 ;vector-scalar mult.
=> 4 x 64 - 256 clocks
LV V3,Ry ;load vector Y (or 4 clocks per result)
3: ADDV V4,V2,V3 ;add
5381/7381
4: SV Ry,V4 ;store the result
DLXV Start-up Time
• Start-up time: pipeline latency time (depth of FU
pipeline); another sources of overhead
• Operation Start-up penalty (from CRAY-1)
• Vector load/store 12
• Vector multply 7
• Vector add 6
Assume convoys don't overlap; vector length =
Convoyn: Start 1st result last result
1. LV 0 12 11+n (12+n-1)
2. MULV, LV 12+n 12+n+12 23+2n Load start-up
3. ADDV 24+2n 24+2n+6 29+3n Wait convoy 2
5381/7381
4. SV 30+3n 30+3n+12 41+4n Wait convoy 3
Why startup time for each
vector instruction?
• Why not overlap startup time of back-to-back
vector instructions?
• Cray machines built from many ECL chips
operating at high clock rates; hard to do?
• Berkeley vector design (“T0”) didn’t know it
wasn’t supposed to do overlap, so no startup
times for functional units (except load)
5381/7381
Vector Load/Store Units & Memories
• Start-up overheads usually longer fo LSUs
• Memory system must sustain (# lanes x word) /clock cycle
• Many Vector Procs. use banks (vs. simple interleaving):
1) support multiple loads/stores per cycle
=> multiple banks & address banks independently
2) support non-sequential accesses (see soon)
• Note: No. memory banks > memory latency to avoid stalls
– m banks => m words per memory lantecy l clocks
– if m < l, then gap in memory pipeline:
clock: 0… ll+1 l+2…l+m- 1 l+m …2 l
word: --… 01 2… m-1 --… m
– may have 1024 banks in SRAM
5381/7381
Vector Length
5381/7381
Strip Mining
• Suppose Vector Length > Max. Vector Length (MVL)?
• Strip mining: generation of code such that each vector
operation is done for a size Š to the MVL
• 1st loop do short piece (n mod MVL), rest VL = MVL
low = 1
VL = (n mod MVL) /*find the odd size piece*/
do 1 j = 0,(n / MVL) /*outer loop*/
do 10 i = low,low+VL-1 /*runs for length VL*/
Y(i) = a*X(i) + Y(i) /*main operation*/
10 continue
low = low+VL /*start of next vector*/
VL = MVL /*reset the length to max*/
1 continue 5381/7381
Common Vector Metrics
• R: MFLOPS rate on an infinite-length vector
– vector “speed of light”
– Real problems do not have unlimited vector lengths, and the start-up penalties encountered in
real problems will be larger
– (Rn is the MFLOPS rate for a vector of length n)
• NV: The vector length needed to make vector mode faster than scalar mode
– measures both start-up and speed of scalars relative to vectors, quality of connection of scalar
unit to vector unit
5381/7381
Vector Stride
• Suppose adjacent elements not sequential in memory
do 10 i = 1,100
do 10 j = 1,100
A(i,j) = 0.0
do 10 k = 1,100
10 A(i,j) = A(i,j)+B(i,k)*C(k,j)
• Either B or C accesses not adjacent (800 bytes between)
• stride: distance separating elements that are to be
merged into a single vector (caches do unit stride)
=> LVWS (load vector with stride) instruction
• Strides => can cause bank conflicts
(e.g., stride = 32 and 16 banks) 5381/7381
5381/7381
Example Execution of Vector Code
Vector Vector Vector
Scalar Memory Pipeline Multiply Pipeline Adder Pipeline
5381/7381
Applications
Limited to scientific computing?
• Multimedia Processing (compress., graphics, audio synth, image proc.)
• Standard benchmark kernels (Matrix Multiply, FFT, Convolution, Sort)
• Lossy Compression (JPEG, MPEG video and audio)
• Lossless Compression (Zero removal, RLE, Differencing, LZW)
• Cryptography (RSA, DES/IDEA, SHA/MD5)
• Speech and handwriting recognition
• Operating systems/Networking (memcpy, memset, parity, checksum)
• Databases (hash/join, data mining, image/video serving)
• Language run-time support (stdlib, garbage collection)
• even SPECint95
5381/7381
Vector for Multimedia?
+
• Claim: overall speedup 1.5 to 2X for 2D/3D graphics, audio, video, speech,
comm., ...
– use in drivers or added to library routines; no compiler
5381/7381
MMX Instructions
5381/7381
Mediaprocesing:
Vectorizable? Vector Lengths?
Kernel Vector length
• Matrix transpose/multiply # vertices at once
• DCT (video, communication) image width
• FFT (audio) 256-1024
• Motion estimation (video) image width, iw/16
• Gamma correction (video) image width
• Haar transform (media mining) image width
• Median filter (image processing) image width
• Separable convolution (img. proc.) image width
5381/7381
Vector Pitfalls
• Pitfall: Concentrating on peak performance and ignoring start-up overhead: NV (length
faster than scalar) > 100!
• Pitfall: Increasing vector performance, without comparable increases in scalar
performance
(Amdahl's Law)
– failure of Cray competitor from his former company
• Pitfall: Good processor vector performance without providing good memory
bandwidth
– MMX?
5381/7381
Vector Advantages
• Easy to get high performance; N operations:
– are independent
– use same functional unit
– access disjoint registers
– access registers in same order as previous instructions
– access contiguous memory words or known pattern
– can exploit large memory bandwidth
– hide memory latency (and any other latency)
• Scalable (get higher performance as more HW resources available)
• Compact: Describe N operations with 1 short instruction (v. VLIW)
• Predictable (real-time) performance vs. statistical performance (cache)
• Multimedia ready: choose N * 64b, 2N * 32b, 4N * 16b, 8N * 8b
• Mature, developed compiler technology
• Vector Disadvantage: Out of Fashion
5381/7381
Vector Summary
• Alternate model accomodates long memory latency,
doesn’t rely on caches as does Out-Of-Order,
superscalar/VLIW designs
• Much easier for hardware: more powerful instructions,
more predictable memory accesses, fewer harzards,
fewer branches, fewer mispredicted branches, ...
• What % of computation is vectorizable?
• Is vector a good match to new apps such as
multidemia, DSP?
5381/7381