add slide
[libreriscv.git] / simple_v_extension.mdwn
1 # Variable-width Variable-packed SIMD / Simple-V / Parallelism Extension Proposal
2
3 Key insight: Simple-V is intended as an abstraction layer to provide
4 a consistent "API" to parallelisation of existing *and future* operations.
5 *Actual* internal hardware-level parallelism is *not* required, such
6 that Simple-V may be viewed as providing a "compact" or "consolidated"
7 means of issuing multiple near-identical arithmetic instructions to an
8 instruction queue (FIFO), pending execution.
9
10 *Actual* parallelism, if added independently of Simple-V in the form
11 of Out-of-order restructuring (including parallel ALU lanes) or VLIW
12 implementations, or SIMD, or anything else, would then benefit *if*
13 Simple-V was added on top.
14
15 [[!toc ]]
16
17 # Introduction
18
19 This proposal exists so as to be able to satisfy several disparate
20 requirements: power-conscious, area-conscious, and performance-conscious
21 designs all pull an ISA and its implementation in different conflicting
22 directions, as do the specific intended uses for any given implementation.
23
24 The existing P (SIMD) proposal and the V (Vector) proposals,
25 whilst each extremely powerful in their own right and clearly desirable,
26 are also:
27
28 * Clearly independent in their origins (Cray and AndesStar v3 respectively)
29 so need work to adapt to the RISC-V ethos and paradigm
30 * Are sufficiently large so as to make adoption (and exploration for
31 analysis and review purposes) prohibitively expensive
32 * Both contain partial duplication of pre-existing RISC-V instructions
33 (an undesirable characteristic)
34 * Both have independent, incompatible and disparate methods for introducing
35 parallelism at the instruction level
36 * Both require that their respective parallelism paradigm be implemented
37 along-side and integral to their respective functionality *or not at all*.
38 * Both independently have methods for introducing parallelism that
39 could, if separated, benefit
40 *other areas of RISC-V not just DSP or Floating-point respectively*.
41
42 There are also key differences between Vectorisation and SIMD (full
43 details outlined in the Appendix), the key points being:
44
45 * SIMD has an extremely seductively compelling ease of implementation argument:
46 each operation is passed to the ALU, which is where the parallelism
47 lies. There is *negligeable* (if any) impact on the rest of the core
48 (with life instead being made hell for compiler writers and applications
49 writers due to extreme ISA proliferation).
50 * By contrast, Vectorisation has quite some complexity (for considerable
51 flexibility, reduction in opcode proliferation and much more).
52 * Vectorisation typically includes much more comprehensive memory load
53 and store schemes (unit stride, constant-stride and indexed), which
54 in turn have ramifications: virtual memory misses (TLB cache misses)
55 and even multiple page-faults... all caused by a *single instruction*,
56 yet with a clear benefit that the regularisation of LOAD/STOREs can
57 be optimised for minimal impact on caches and maximised throughput.
58 * By contrast, SIMD can use "standard" memory load/stores (32-bit aligned
59 to pages), and these load/stores have absolutely nothing to do with the
60 SIMD / ALU engine, no matter how wide the operand. Simplicity but with
61 more impact on instruction and data caches.
62
63 Overall it makes a huge amount of sense to have a means and method
64 of introducing instruction parallelism in a flexible way that provides
65 implementors with the option to choose exactly where they wish to offer
66 performance improvements and where they wish to optimise for power
67 and/or area (and if that can be offered even on a per-operation basis that
68 would provide even more flexibility).
69
70 Additionally it makes sense to *split out* the parallelism inherent within
71 each of P and V, and to see if each of P and V then, in *combination* with
72 a "best-of-both" parallelism extension, could be added on *on top* of
73 this proposal, to topologically provide the exact same functionality of
74 each of P and V. Each of P and V then can focus on providing the best
75 operations possible for their respective target areas, without being
76 hugely concerned about the actual parallelism.
77
78 Furthermore, an additional goal of this proposal is to reduce the number
79 of opcodes utilised by each of P and V as they currently stand, leveraging
80 existing RISC-V opcodes where possible, and also potentially allowing
81 P and V to make use of Compressed Instructions as a result.
82
83 # Analysis and discussion of Vector vs SIMD
84
85 There are six combined areas between the two proposals that help with
86 parallelism (increased performance, reduced power / area) without
87 over-burdening the ISA with a huge proliferation of
88 instructions:
89
90 * Fixed vs variable parallelism (fixed or variable "M" in SIMD)
91 * Implicit vs fixed instruction bit-width (integral to instruction or not)
92 * Implicit vs explicit type-conversion (compounded on bit-width)
93 * Implicit vs explicit inner loops.
94 * Single-instruction LOAD/STORE.
95 * Masks / tagging (selecting/preventing certain indexed elements from execution)
96
97 The pros and cons of each are discussed and analysed below.
98
99 ## Fixed vs variable parallelism length
100
101 In David Patterson and Andrew Waterman's analysis of SIMD and Vector
102 ISAs, the analysis comes out clearly in favour of (effectively) variable
103 length SIMD. As SIMD is a fixed width, typically 4, 8 or in extreme cases
104 16 or 32 simultaneous operations, the setup, teardown and corner-cases of SIMD
105 are extremely burdensome except for applications whose requirements
106 *specifically* match the *precise and exact* depth of the SIMD engine.
107
108 Thus, SIMD, no matter what width is chosen, is never going to be acceptable
109 for general-purpose computation, and in the context of developing a
110 general-purpose ISA, is never going to satisfy 100 percent of implementors.
111
112 To explain this further: for increased workloads over time, as the
113 performance requirements increase for new target markets, implementors
114 choose to extend the SIMD width (so as to again avoid mixing parallelism
115 into the instruction issue phases: the primary "simplicity" benefit of
116 SIMD in the first place), with the result that the entire opcode space
117 effectively doubles with each new SIMD width that's added to the ISA.
118
119 That basically leaves "variable-length vector" as the clear *general-purpose*
120 winner, at least in terms of greatly simplifying the instruction set,
121 reducing the number of instructions required for any given task, and thus
122 reducing power consumption for the same.
123
124 ## Implicit vs fixed instruction bit-width
125
126 SIMD again has a severe disadvantage here, over Vector: huge proliferation
127 of specialist instructions that target 8-bit, 16-bit, 32-bit, 64-bit, and
128 have to then have operations *for each and between each*. It gets very
129 messy, very quickly.
130
131 The V-Extension on the other hand proposes to set the bit-width of
132 future instructions on a per-register basis, such that subsequent instructions
133 involving that register are *implicitly* of that particular bit-width until
134 otherwise changed or reset.
135
136 This has some extremely useful properties, without being particularly
137 burdensome to implementations, given that instruction decode already has
138 to direct the operation to a correctly-sized width ALU engine, anyway.
139
140 Not least: in places where an ISA was previously constrained (due for
141 whatever reason, including limitations of the available operand space),
142 implicit bit-width allows the meaning of certain operations to be
143 type-overloaded *without* pollution or alteration of frozen and immutable
144 instructions, in a fully backwards-compatible fashion.
145
146 ## Implicit and explicit type-conversion
147
148 The Draft 2.3 V-extension proposal has (deprecated) polymorphism to help
149 deal with over-population of instructions, such that type-casting from
150 integer (and floating point) of various sizes is automatically inferred
151 due to "type tagging" that is set with a special instruction. A register
152 will be *specifically* marked as "16-bit Floating-Point" and, if added
153 to an operand that is specifically tagged as "32-bit Integer" an implicit
154 type-conversion will take place *without* requiring that type-conversion
155 to be explicitly done with its own separate instruction.
156
157 However, implicit type-conversion is not only quite burdensome to
158 implement (explosion of inferred type-to-type conversion) but also is
159 never really going to be complete. It gets even worse when bit-widths
160 also have to be taken into consideration. Each new type results in
161 an increased O(N^2) conversion space that, as anyone who has examined
162 python's source code (which has built-in polymorphic type-conversion),
163 knows that the task is more complex than it first seems.
164
165 Overall, type-conversion is generally best to leave to explicit
166 type-conversion instructions, or in definite specific use-cases left to
167 be part of an actual instruction (DSP or FP)
168
169 ## Zero-overhead loops vs explicit loops
170
171 The initial Draft P-SIMD Proposal by Chuanhua Chang of Andes Technology
172 contains an extremely interesting feature: zero-overhead loops. This
173 proposal would basically allow an inner loop of instructions to be
174 repeated indefinitely, a fixed number of times.
175
176 Its specific advantage over explicit loops is that the pipeline in a DSP
177 can potentially be kept completely full *even in an in-order single-issue
178 implementation*. Normally, it requires a superscalar architecture and
179 out-of-order execution capabilities to "pre-process" instructions in
180 order to keep ALU pipelines 100% occupied.
181
182 By bringing that capability in, this proposal could offer a way to increase
183 pipeline activity even in simpler implementations in the one key area
184 which really matters: the inner loop.
185
186 However when looking at much more comprehensive schemes
187 "A portable specification of zero-overhead loop control hardware
188 applied to embedded processors" (ZOLC), optimising only the single
189 inner loop seems inadequate, tending to suggest that ZOLC may be
190 better off being proposed as an entirely separate Extension.
191
192 ## Single-instruction LOAD/STORE
193
194 In traditional Vector Architectures there are instructions which
195 result in multiple register-memory transfer operations resulting
196 from a single instruction. They're complicated to implement in hardware,
197 yet the benefits are a huge consistent regularisation of memory accesses
198 that can be highly optimised with respect to both actual memory and any
199 L1, L2 or other caches. In Hwacha EECS-2015-263 it is explicitly made
200 clear the consequences of getting this architecturally wrong:
201 L2 cache-thrashing at the very least.
202
203 Complications arise when Virtual Memory is involved: TLB cache misses
204 need to be dealt with, as do page faults. Some of the tradeoffs are
205 discussed in <http://people.eecs.berkeley.edu/~krste/thesis.pdf>, Section
206 4.6, and an article by Jeff Bush when faced with some of these issues
207 is particularly enlightening
208 <https://jbush001.github.io/2015/11/03/lost-in-translation.html>
209
210 Interestingly, none of this complexity is faced in SIMD architectures...
211 but then they do not get the opportunity to optimise for highly-streamlined
212 memory accesses either.
213
214 With the "bang-per-buck" ratio being so high and the indirect improvement
215 in L1 Instruction Cache usage (reduced instruction count), as well as
216 the opportunity to optimise L1 and L2 cache usage, the case for including
217 Vector LOAD/STORE is compelling.
218
219 ## Mask and Tagging (Predication)
220
221 Tagging (aka Masks aka Predication) is a pseudo-method of implementing
222 simplistic branching in a parallel fashion, by allowing execution on
223 elements of a vector to be switched on or off depending on the results
224 of prior operations in the same array position.
225
226 The reason for considering this is simple: by *definition* it
227 is not possible to perform individual parallel branches in a SIMD
228 (Single-Instruction, **Multiple**-Data) context. Branches (modifying
229 of the Program Counter) will result in *all* parallel data having
230 a different instruction executed on it: that's just the definition of
231 SIMD, and it is simply unavoidable.
232
233 So these are the ways in which conditional execution may be implemented:
234
235 * explicit compare and branch: BNE x, y -> offs would jump offs
236 instructions if x was not equal to y
237 * explicit store of tag condition: CMP x, y -> tagbit
238 * implicit (condition-code) such as ADD results in a carry, carry bit
239 implicitly (or sometimes explicitly) goes into a "tag" (mask) register
240
241 The first of these is a "normal" branch method, which is flat-out impossible
242 to parallelise without look-ahead and effectively rewriting instructions.
243 This would defeat the purpose of RISC.
244
245 The latter two are where parallelism becomes easy to do without complexity:
246 every operation is modified to be "conditionally executed" (in an explicit
247 way directly in the instruction format *or* implicitly).
248
249 RVV (Vector-Extension) proposes to have *explicit* storing of the compare
250 in a tag/mask register, and to *explicitly* have every vector operation
251 *require* that its operation be "predicated" on the bits within an
252 explicitly-named tag/mask register.
253
254 SIMD (P-Extension) has not yet published precise documentation on what its
255 schema is to be: there is however verbal indication at the time of writing
256 that:
257
258 > The "compare" instructions in the DSP/SIMD ISA proposed by Andes will
259 > be executed using the same compare ALU logic for the base ISA with some
260 > minor modifications to handle smaller data types. The function will not
261 > be duplicated.
262
263 This is an *implicit* form of predication as the base RV ISA does not have
264 condition-codes or predication. By adding a CSR it becomes possible
265 to also tag certain registers as "predicated if referenced as a destination".
266 Example:
267
268 // in future operations from now on, if r0 is the destination use r5 as
269 // the PREDICATION register
270 SET_IMPLICIT_CSRPREDICATE r0, r5
271 // store the compares in r5 as the PREDICATION register
272 CMPEQ8 r5, r1, r2
273 // r0 is used here. ah ha! that means it's predicated using r5!
274 ADD8 r0, r1, r3
275
276 With enough registers (and in RISC-V there are enough registers) some fairly
277 complex predication can be set up and yet still execute without significant
278 stalling, even in a simple non-superscalar architecture.
279
280 (For details on how Branch Instructions would be retro-fitted to indirectly
281 predicated equivalents, see Appendix)
282
283 ## Conclusions
284
285 In the above sections the five different ways where parallel instruction
286 execution has closely and loosely inter-related implications for the ISA and
287 for implementors, were outlined. The pluses and minuses came out as
288 follows:
289
290 * Fixed vs variable parallelism: <b>variable</b>
291 * Implicit (indirect) vs fixed (integral) instruction bit-width: <b>indirect</b>
292 * Implicit vs explicit type-conversion: <b>explicit</b>
293 * Implicit vs explicit inner loops: <b>implicit but best done separately</b>
294 * Single-instruction Vector LOAD/STORE: <b>Complex but highly beneficial</b>
295 * Tag or no-tag: <b>Complex but highly beneficial</b>
296
297 In particular:
298
299 * variable-length vectors came out on top because of the high setup, teardown
300 and corner-cases associated with the fixed width of SIMD.
301 * Implicit bit-width helps to extend the ISA to escape from
302 former limitations and restrictions (in a backwards-compatible fashion),
303 whilst also leaving implementors free to simmplify implementations
304 by using actual explicit internal parallelism.
305 * Implicit (zero-overhead) loops provide a means to keep pipelines
306 potentially 100% occupied in a single-issue in-order implementation
307 i.e. *without* requiring a super-scalar or out-of-order architecture,
308 but doing a proper, full job (ZOLC) is an entirely different matter.
309
310 Constructing a SIMD/Simple-Vector proposal based around four of these six
311 requirements would therefore seem to be a logical thing to do.
312
313 # Note on implementation of parallelism
314
315 One extremely important aspect of this proposal is to respect and support
316 implementors desire to focus on power, area or performance. In that regard,
317 it is proposed that implementors be free to choose whether to implement
318 the Vector (or variable-width SIMD) parallelism as sequential operations
319 with a single ALU, fully parallel (if practical) with multiple ALUs, or
320 a hybrid combination of both.
321
322 In Broadcom's Videocore-IV, they chose hybrid, and called it "Virtual
323 Parallelism". They achieve a 16-way SIMD at an **instruction** level
324 by providing a combination of a 4-way parallel ALU *and* an externally
325 transparent loop that feeds 4 sequential sets of data into each of the
326 4 ALUs.
327
328 Also in the same core, it is worth noting that particularly uncommon
329 but essential operations (Reciprocal-Square-Root for example) are
330 *not* part of the 4-way parallel ALU but instead have a *single* ALU.
331 Under the proposed Vector (varible-width SIMD) implementors would
332 be free to do precisely that: i.e. free to choose *on a per operation
333 basis* whether and how much "Virtual Parallelism" to deploy.
334
335 It is absolutely critical to note that it is proposed that such choices MUST
336 be **entirely transparent** to the end-user and the compiler. Whilst
337 a Vector (varible-width SIMD) may not precisely match the width of the
338 parallelism within the implementation, the end-user **should not care**
339 and in this way the performance benefits are gained but the ISA remains
340 straightforward. All that happens at the end of an instruction run is: some
341 parallel units (if there are any) would remain offline, completely
342 transparently to the ISA, the program, and the compiler.
343
344 To make that clear: should an implementor choose a particularly wide
345 SIMD-style ALU, each parallel unit *must* have predication so that
346 the parallel SIMD ALU may emulate variable-length parallel operations.
347 Thus the "SIMD considered harmful" trap of having huge complexity and extra
348 instructions to deal with corner-cases is thus avoided, and implementors
349 get to choose precisely where to focus and target the benefits of their
350 implementation efforts, without "extra baggage".
351
352 In addition, implementors will be free to choose whether to provide an
353 absolute bare minimum level of compliance with the "API" (software-traps
354 when vectorisation is detected), all the way up to full supercomputing
355 level all-hardware parallelism. Options are covered in the Appendix.
356
357 # CSRs <a name="csrs"></a>
358
359 There are a number of CSRs needed, which are used at the instruction
360 decode phase to re-interpret RV opcodes (a practice that has
361 precedent in the setting of MISA to enable / disable extensions).
362
363 * Integer Register N is Vector of length M: r(N) -> r(N..N+M-1)
364 * Integer Register N is of implicit bitwidth M (M=default,8,16,32,64)
365 * Floating-point Register N is Vector of length M: r(N) -> r(N..N+M-1)
366 * Floating-point Register N is of implicit bitwidth M (M=default,8,16,32,64)
367 * Integer Register N is a Predication Register (note: a key-value store)
368 * Vector Length CSR (VSETVL, VGETVL)
369
370 Also (see Appendix, "Context Switch Example") it may turn out to be important
371 to have a separate (smaller) set of CSRs for M-Mode (and S-Mode) so that
372 Vectorised LOAD / STORE may be used to load and store multiple registers:
373 something that is missing from the Base RV ISA.
374
375 Notes:
376
377 * for the purposes of LOAD / STORE, Integer Registers which are
378 marked as a Vector will result in a Vector LOAD / STORE.
379 * Vector Lengths are *not* the same as vsetl but are an integral part
380 of vsetl.
381 * Actual vector length is *multipled* by how many blocks of length
382 "bitwidth" may fit into an XLEN-sized register file.
383 * Predication is a key-value store due to the implicit referencing,
384 as opposed to having the predicate register explicitly in the instruction.
385 * Whilst the predication CSR is a key-value store it *generates* easier-to-use
386 state information.
387 * TODO: assess whether the same technique could be applied to the other
388 Vector CSRs, particularly as pointed out in Section 17.8 (Draft RV 0.4,
389 V2.3-Draft ISA Reference) it becomes possible to greatly reduce state
390 needed for context-switches (empty slots need never be stored).
391
392 ## Predication CSR
393
394 The Predication CSR is a key-value store indicating whether, if a given
395 destination register (integer or floating-point) is referred to in an
396 instruction, it is to be predicated. The first entry is whether predication
397 is enabled. The second entry is whether the register index refers to a
398 floating-point or an integer register. The third entry is the index
399 of that register which is to be predicated (if referred to). The fourth entry
400 is the integer register that is treated as a bitfield, indexable by the
401 vector element index.
402
403 | RegNo | 6 | 5 | (4..0) | (4..0) |
404 | ----- | - | - | ------- | ------- |
405 | r0 | pren0 | i/f | regidx | predidx |
406 | r1 | pren1 | i/f | regidx | predidx |
407 | .. | pren.. | i/f | regidx | predidx |
408 | r15 | pren15 | i/f | regidx | predidx |
409
410 The Predication CSR Table is a key-value store, so implementation-wise
411 it will be faster to turn the table around (maintain topologically
412 equivalent state):
413
414 fp_pred_enabled[32];
415 int_pred_enabled[32];
416 for (i = 0; i < 16; i++)
417 if CSRpred[i].pren:
418 idx = CSRpred[i].regidx
419 predidx = CSRpred[i].predidx
420 if CSRpred[i].type == 0: # integer
421 int_pred_enabled[idx] = 1
422 int_pred_reg[idx] = predidx
423 else:
424 fp_pred_enabled[idx] = 1
425 fp_pred_reg[idx] = predidx
426
427 So when an operation is to be predicated, it is the internal state that
428 is used. In Section 6.4.2 of Hwacha's Manual (EECS-2015-262) the following
429 pseudo-code for operations is given, where p is the explicit (direct)
430 reference to the predication register to be used:
431
432 for (int i=0; i<vl; ++i)
433 if ([!]preg[p][i])
434 (d ? vreg[rd][i] : sreg[rd]) =
435 iop(s1 ? vreg[rs1][i] : sreg[rs1],
436 s2 ? vreg[rs2][i] : sreg[rs2]); // for insts with 2 inputs
437
438 This instead becomes an *indirect* reference using the *internal* state
439 table generated from the Predication CSR key-value store:
440
441 if type(iop) == INT:
442 pred_enabled = int_pred_enabled
443 preg = int_pred_reg[rd]
444 else:
445 pred_enabled = fp_pred_enabled
446 preg = fp_pred_reg[rd]
447
448 for (int i=0; i<vl; ++i)
449 if (preg_enabled[rd] && [!]preg[i])
450 (d ? vreg[rd][i] : sreg[rd]) =
451 iop(s1 ? vreg[rs1][i] : sreg[rs1],
452 s2 ? vreg[rs2][i] : sreg[rs2]); // for insts with 2 inputs
453
454 ## MAXVECTORDEPTH
455
456 MAXVECTORDEPTH is the same concept as MVL in RVV. However in Simple-V,
457 given that its primary (base, unextended) purpose is for 3D, Video and
458 other purposes (not requiring supercomputing capability), it makes sense
459 to limit MAXVECTORDEPTH to the regfile bitwidth (32 for RV32, 64 for RV64
460 and so on).
461
462 The reason for setting this limit is so that predication registers, when
463 marked as such, may fit into a single register as opposed to fanning out
464 over several registers. This keeps the implementation a little simpler.
465 Note also (as also described in the VSETVL section) that the *minimum*
466 for MAXVECTORDEPTH must be the total number of registers (15 for RV32E
467 and 31 for RV32 or RV64).
468
469 Note that RVV on top of Simple-V may choose to over-ride this decision.
470
471 ## Vector-length CSRs
472
473 Vector lengths are interpreted as meaning "any instruction referring to
474 r(N) generates implicit identical instructions referring to registers
475 r(N+M-1) where M is the Vector Length". Vector Lengths may be set to
476 use up to 16 registers in the register file.
477
478 One separate CSR table is needed for each of the integer and floating-point
479 register files:
480
481 | RegNo | (3..0) |
482 | ----- | ------ |
483 | r0 | vlen0 |
484 | r1 | vlen1 |
485 | .. | vlen.. |
486 | r31 | vlen31 |
487
488 An array of 32 4-bit CSRs is needed (4 bits per register) to indicate
489 whether a register was, if referred to in any standard instructions,
490 implicitly to be treated as a vector.
491
492 Note:
493
494 * A vector length of 1 indicates that it is to be treated as a scalar.
495 Bitwidths (on the same register) are interpreted and meaningful.
496 * A vector length of 0 indicates that the parallelism is to be switched
497 off for this register (treated as a scalar). When length is 0,
498 the bitwidth CSR for the register is *ignored*.
499
500 Internally, implementations may choose to use the non-zero vector length
501 to set a bit-field per register, to be used in the instruction decode phase.
502 In this way any standard (current or future) operation involving
503 register operands may detect if the operation is to be vector-vector,
504 vector-scalar or scalar-scalar (standard) simply through a single
505 bit test.
506
507 Note that when using the "vsetl rs1, rs2" instruction (caveat: when the
508 bitwidth is specifically not set) it becomes:
509
510 CSRvlength = MIN(MIN(CSRvectorlen[rs1], MAXVECTORDEPTH), rs2)
511
512 This is in contrast to RVV:
513
514 CSRvlength = MIN(MIN(rs1, MAXVECTORDEPTH), rs2)
515
516 ## Element (SIMD) bitwidth CSRs
517
518 Element bitwidths may be specified with a per-register CSR, and indicate
519 how a register (integer or floating-point) is to be subdivided.
520
521 | RegNo | (2..0) |
522 | ----- | ------ |
523 | r0 | vew0 |
524 | r1 | vew1 |
525 | .. | vew.. |
526 | r31 | vew31 |
527
528 vew may be one of the following (giving a table "bytestable", used below):
529
530 | vew | bitwidth |
531 | --- | -------- |
532 | 000 | default |
533 | 001 | 8 |
534 | 010 | 16 |
535 | 011 | 32 |
536 | 100 | 64 |
537 | 101 | 128 |
538 | 110 | rsvd |
539 | 111 | rsvd |
540
541 Extending this table (with extra bits) is covered in the section
542 "Implementing RVV on top of Simple-V".
543
544 Note that when using the "vsetl rs1, rs2" instruction, taking bitwidth
545 into account, it becomes:
546
547 vew = CSRbitwidth[rs1]
548 if (vew == 0)
549 bytesperreg = (XLEN/8) # or FLEN as appropriate
550 else:
551 bytesperreg = bytestable[vew] # 1 2 4 8 16
552 simdmult = (XLEN/8) / bytesperreg # or FLEN as appropriate
553 vlen = CSRvectorlen[rs1] * simdmult
554 CSRvlength = MIN(MIN(vlen, MAXVECTORDEPTH), rs2)
555
556 The reason for multiplying the vector length by the number of SIMD elements
557 (in each individual register) is so that each SIMD element may optionally be
558 predicated.
559
560 An example of how to subdivide the register file when bitwidth != default
561 is given in the section "Bitwidth Virtual Register Reordering".
562
563 # Instructions
564
565 By being a topological remap of RVV concepts, the following RVV instructions
566 remain exactly the same: VMPOP, VMFIRST, VEXTRACT, VINSERT, VMERGE, VSELECT,
567 VSLIDE, VCLASS and VPOPC. Two instructions, VCLIP and VCLIPI, do not
568 have RV Standard equivalents, so are left out of Simple-V.
569 All other instructions from RVV are topologically re-mapped and retain
570 their complete functionality, intact.
571
572 ## Instruction Format
573
574 The instruction format for Simple-V does not actually have *any* explicit
575 compare operations, *any* arithmetic, floating point or *any*
576 memory instructions.
577 Instead it *overloads* pre-existing branch operations into predicated
578 variants, and implicitly overloads arithmetic operations and LOAD/STORE
579 depending on CSR configurations for vector length, bitwidth and
580 predication. *This includes Compressed instructions* as well as any
581 future instructions and Custom Extensions.
582
583 * For analysis of RVV see [[v_comparative_analysis]] which begins to
584 outline topologically-equivalent mappings of instructions
585 * Also see Appendix "Retro-fitting Predication into branch-explicit ISA"
586 for format of Branch opcodes.
587
588 **TODO**: *analyse and decide whether the implicit nature of predication
589 as proposed is or is not a lot of hassle, and if explicit prefixes are
590 a better idea instead. Parallelism therefore effectively may end up
591 as always being 64-bit opcodes (32 for the prefix, 32 for the instruction)
592 with some opportunities for to use Compressed bringing it down to 48.
593 Also to consider is whether one or both of the last two remaining Compressed
594 instruction codes in Quadrant 1 could be used as a parallelism prefix,
595 bringing parallelised opcodes down to 32-bit (when combined with C)
596 and having the benefit of being explicit.*
597
598 ## VSETVL
599
600 VSETVL is slightly different from RVV in that the minimum vector length
601 is required to be at least the number of registers in the register file,
602 and no more than XLEN. This allows vector LOAD/STORE to be used to switch
603 the entire bank of registers using a single instruction (see Appendix,
604 "Context Switch Example"). The reason for limiting VSETVL to XLEN is
605 down to the fact that predication bits fit into a single register of length
606 XLEN bits.
607
608 The second minor change is that when VSETVL is requested to be stored
609 into x0, it is *ignored* silently.
610
611 Unlike RVV, implementors *must* provide pseudo-parallelism (using sequential
612 loops in hardware) if actual hardware-parallelism in the ALUs is not deployed.
613 A hybrid is also permitted (as used in Broadcom's VideoCore-IV) however this
614 must be *entirely* transparent to the ISA.
615
616 ## Branch Instruction:
617
618 Branch operations use standard RV opcodes that are reinterpreted to be
619 "predicate variants" in the instance where either of the two src registers
620 have their corresponding CSRvectorlen[src] entry as non-zero. When this
621 reinterpretation is enabled the predicate target register rs3 is to be
622 treated as a bitfield (up to a maximum of XLEN bits corresponding to a
623 maximum of XLEN elements).
624
625 If either of src1 or src2 are scalars (CSRvectorlen[src] == 0) the comparison
626 goes ahead as vector-scalar or scalar-vector. Implementors should note that
627 this could require considerable multi-porting of the register file in order
628 to parallelise properly, so may have to involve the use of register cacheing
629 and transparent copying (see Multiple-Banked Register File Architectures
630 paper).
631
632 In instances where no vectorisation is detected on either src registers
633 the operation is treated as an absolutely standard scalar branch operation.
634
635 This is the overloaded table for Integer-base Branch operations. Opcode
636 (bits 6..0) is set in all cases to 1100011.
637
638 [[!table data="""
639 31 .. 25 |24 ... 20 | 19 15 | 14 12 | 11 .. 8 | 7 | 6 ... 0 |
640 imm[12,10:5]| rs2 | rs1 | funct3 | imm[4:1] | imm[11] | opcode |
641 7 | 5 | 5 | 3 | 4 | 1 | 7 |
642 reserved | src2 | src1 | BPR | predicate rs3 || BRANCH |
643 reserved | src2 | src1 | 000 | predicate rs3 || BEQ |
644 reserved | src2 | src1 | 001 | predicate rs3 || BNE |
645 reserved | src2 | src1 | 010 | predicate rs3 || rsvd |
646 reserved | src2 | src1 | 011 | predicate rs3 || rsvd |
647 reserved | src2 | src1 | 100 | predicate rs3 || BLE |
648 reserved | src2 | src1 | 101 | predicate rs3 || BGE |
649 reserved | src2 | src1 | 110 | predicate rs3 || BLTU |
650 reserved | src2 | src1 | 111 | predicate rs3 || BGEU |
651 """]]
652
653 Note that just as with the standard (scalar, non-predicated) branch
654 operations, BLT, BGT, BLEU and BTGU may be synthesised by inverting
655 src1 and src2.
656
657 Below is the overloaded table for Floating-point Predication operations.
658 Interestingly no change is needed to the instruction format because
659 FP Compare already stores a 1 or a zero in its "rd" integer register
660 target, i.e. it's not actually a Branch at all: it's a compare.
661 The target needs to simply change to be a predication bitfield (done
662 implicitly).
663
664 As with
665 Standard RVF/D/Q, Opcode (bits 6..0) is set in all cases to 1010011.
666 Likewise Single-precision, fmt bits 26..25) is still set to 00.
667 Double-precision is still set to 01, whilst Quad-precision
668 appears not to have a definition in V2.3-Draft (but should be unaffected).
669
670 It is however noted that an entry "FNE" (the opposite of FEQ) is missing,
671 and whilst in ordinary branch code this is fine because the standard
672 RVF compare can always be followed up with an integer BEQ or a BNE (or
673 a compressed comparison to zero or non-zero), in predication terms that
674 becomes more of an impact as an explicit (scalar) instruction is needed
675 to invert the predicate bitmask. An additional encoding funct3=011 is
676 therefore proposed to cater for this.
677
678 [[!table data="""
679 31 .. 27| 26 .. 25 |24 ... 20 | 19 15 | 14 12 | 11 .. 7 | 6 ... 0 |
680 funct5 | fmt | rs2 | rs1 | funct3 | rd | opcode |
681 5 | 2 | 5 | 5 | 3 | 4 | 7 |
682 10100 | 00/01/11 | src2 | src1 | 010 | pred rs3 | FEQ |
683 10100 | 00/01/11 | src2 | src1 | **011**| pred rs3 | FNE |
684 10100 | 00/01/11 | src2 | src1 | 001 | pred rs3 | FLT |
685 10100 | 00/01/11 | src2 | src1 | 000 | pred rs3 | FLE |
686 """]]
687
688 Note (**TBD**): floating-point exceptions will need to be extended
689 to cater for multiple exceptions (and statuses of the same). The
690 usual approach is to have an array of status codes and bit-fields,
691 and one exception, rather than throw separate exceptions for each
692 Vector element.
693
694 In Hwacha EECS-2015-262 Section 6.7.2 the following pseudocode is given
695 for predicated compare operations of function "cmp":
696
697 for (int i=0; i<vl; ++i)
698 if ([!]preg[p][i])
699 preg[pd][i] = cmp(s1 ? vreg[rs1][i] : sreg[rs1],
700 s2 ? vreg[rs2][i] : sreg[rs2]);
701
702 With associated predication, vector-length adjustments and so on,
703 and temporarily ignoring bitwidth (which makes the comparisons more
704 complex), this becomes:
705
706 if I/F == INT: # integer type cmp
707 pred_enabled = int_pred_enabled # TODO: exception if not set!
708 preg = int_pred_reg[rd]
709 reg = int_regfile
710 else:
711 pred_enabled = fp_pred_enabled # TODO: exception if not set!
712 preg = fp_pred_reg[rd]
713 reg = fp_regfile
714
715 s1 = CSRvectorlen[src1] > 1;
716 s2 = CSRvectorlen[src2] > 1;
717 for (int i=0; i<vl; ++i)
718 preg[rs3][i] = cmp(s1 ? reg[src1+i] : reg[src1],
719 s2 ? reg[src2+i] : reg[src2]);
720
721 Notes:
722
723 * Predicated SIMD comparisons would break src1 and src2 further down
724 into bitwidth-sized chunks (see Appendix "Bitwidth Virtual Register
725 Reordering") setting Vector-Length times (number of SIMD elements) bits
726 in Predicate Register rs3 as opposed to just Vector-Length bits.
727 * Predicated Branches do not actually have an adjustment to the Program
728 Counter, so all of bits 25 through 30 in every case are not needed.
729 * There are plenty of reserved opcodes for which bits 25 through 30 could
730 be put to good use if there is a suitable use-case.
731 * FEQ and FNE (and BEQ and BNE) are included in order to save one
732 instruction having to invert the resultant predicate bitfield.
733 FLT and FLE may be inverted to FGT and FGE if needed by swapping
734 src1 and src2 (likewise the integer counterparts).
735
736 ## Compressed Branch Instruction:
737
738 [[!table data="""
739 15..13 | 12...10 | 9..7 | 6..5 | 4..2 | 1..0 | name |
740 funct3 | imm | rs10 | imm | | op | |
741 3 | 3 | 3 | 2 | 3 | 2 | |
742 C.BPR | pred rs3 | src1 | I/F B | src2 | C1 | |
743 110 | pred rs3 | src1 | I/F 0 | src2 | C1 | P.EQ |
744 111 | pred rs3 | src1 | I/F 0 | src2 | C1 | P.NE |
745 110 | pred rs3 | src1 | I/F 1 | src2 | C1 | P.LT |
746 111 | pred rs3 | src1 | I/F 1 | src2 | C1 | P.LE |
747 """]]
748
749 Notes:
750
751 * Bits 5 13 14 and 15 make up the comparator type
752 * Bit 6 indicates whether to use integer or floating-point comparisons
753 * In both floating-point and integer cases there are four predication
754 comparators: EQ/NEQ/LT/LE (with GT and GE being synthesised by inverting
755 src1 and src2).
756
757 ## LOAD / STORE Instructions <a name="load_store"></a>
758
759 For full analysis of topological adaptation of RVV LOAD/STORE
760 see [[v_comparative_analysis]]. All three types (LD, LD.S and LD.X)
761 may be implicitly overloaded into the one base RV LOAD instruction,
762 and likewise for STORE.
763
764 Revised LOAD:
765
766 [[!table data="""
767 31 | 30 | 29 25 | 24 20 | 19 15 | 14 12 | 11 7 | 6 0 |
768 imm[11:0] |||| rs1 | funct3 | rd | opcode |
769 1 | 1 | 5 | 5 | 5 | 3 | 5 | 7 |
770 ? | s | rs2 | imm[4:0] | base | width | dest | LOAD |
771 """]]
772
773 The exact same corresponding adaptation is also carried out on the single,
774 double and quad precision floating-point LOAD-FP and STORE-FP operations,
775 which fit the exact same instruction format. Thus all three types
776 (unit, stride and indexed) may be fitted into FLW, FLD and FLQ,
777 as well as FSW, FSD and FSQ.
778
779 Notes:
780
781 * LOAD remains functionally (topologically) identical to RVV LOAD
782 (for both integer and floating-point variants).
783 * Predication CSR-marking register is not explicitly shown in instruction, it's
784 implicit based on the CSR predicate state for the rd (destination) register
785 * rs2, the source, may *also be marked as a vector*, which implicitly
786 is taken to indicate "Indexed Load" (LD.X)
787 * Bit 30 indicates "element stride" or "constant-stride" (LD or LD.S)
788 * Bit 31 is reserved (ideas under consideration: auto-increment)
789 * **TODO**: include CSR SIMD bitwidth in the pseudo-code below.
790 * **TODO**: clarify where width maps to elsize
791
792 Pseudo-code (excludes CSR SIMD bitwidth for simplicity):
793
794 if (unit-strided) stride = elsize;
795 else stride = areg[as2]; // constant-strided
796
797 pred_enabled = int_pred_enabled
798 preg = int_pred_reg[rd]
799
800 for (int i=0; i<vl; ++i)
801 if (preg_enabled[rd] && [!]preg[i])
802 for (int j=0; j<seglen+1; j++)
803 {
804 if CSRvectorised[rs2])
805 offs = vreg[rs2][i]
806 else
807 offs = i*(seglen+1)*stride;
808 vreg[rd+j][i] = mem[sreg[base] + offs + j*stride];
809 }
810
811 Taking CSR (SIMD) bitwidth into account involves using the vector
812 length and register encoding according to the "Bitwidth Virtual Register
813 Reordering" scheme shown in the Appendix (see function "regoffs").
814
815 A similar instruction exists for STORE, with identical topological
816 translation of all features. **TODO**
817
818 ## Compressed LOAD / STORE Instructions
819
820 Compressed LOAD and STORE are of the same format, where bits 2-4 are
821 a src register instead of dest:
822
823 [[!table data="""
824 15 13 | 12 10 | 9 7 | 6 5 | 4 2 | 1 0 |
825 funct3 | imm | rs10 | imm | rd0 | op |
826 3 | 3 | 3 | 2 | 3 | 2 |
827 C.LW | offset[5:3] | base | offset[2|6] | dest | C0 |
828 """]]
829
830 Unfortunately it is not possible to fit the full functionality
831 of vectorised LOAD / STORE into C.LD / C.ST: the "X" variants (Indexed)
832 require another operand (rs2) in addition to the operand width
833 (which is also missing), offset, base, and src/dest.
834
835 However a close approximation may be achieved by taking the top bit
836 of the offset in each of the five types of LD (and ST), reducing the
837 offset to 4 bits and utilising the 5th bit to indicate whether "stride"
838 is to be enabled. In this way it is at least possible to introduce
839 that functionality.
840
841 (**TODO**: *assess whether the loss of one bit from offset is worth having
842 "stride" capability.*)
843
844 We also assume (including for the "stride" variant) that the "width"
845 parameter, which is missing, is derived and implicit, just as it is
846 with the standard Compressed LOAD/STORE instructions. For C.LW, C.LD
847 and C.LQ, the width is implicitly 4, 8 and 16 respectively, whilst for
848 C.FLW and C.FLD the width is implicitly 4 and 8 respectively.
849
850 Interestingly we note that the Vectorised Simple-V variant of
851 LOAD/STORE (Compressed and otherwise), due to it effectively using the
852 standard register file(s), is the direct functional equivalent of
853 standard load-multiple and store-multiple instructions found in other
854 processors.
855
856 In Section 12.3 riscv-isa manual V2.3-draft it is noted the comments on
857 page 76, "For virtual memory systems some data accesses could be resident
858 in physical memory and some not". The interesting question then arises:
859 how does RVV deal with the exact same scenario?
860 Expired U.S. Patent 5895501 (Filing Date Sep 3 1996) describes a method
861 of detecting early page / segmentation faults and adjusting the TLB
862 in advance, accordingly: other strategies are explored in the Appendix
863 Section "Virtual Memory Page Faults".
864
865 # Exceptions
866
867 > What does an ADD of two different-sized vectors do in simple-V?
868
869 * if the two source operands are not the same, throw an exception.
870 * if the destination operand is also a vector, and the source is longer
871 than the destination, throw an exception.
872
873 > And what about instructions like JALR? 
874 > What does jumping to a vector do?
875
876 * Throw an exception. Whether that actually results in spawning threads
877 as part of the trap-handling remains to be seen.
878
879 # Impementing V on top of Simple-V
880
881 With Simple-V converting the original RVV draft concept-for-concept
882 from explicit opcodes to implicit overloading of existing RV Standard
883 Extensions, certain features were (deliberately) excluded that need
884 to be added back in for RVV to reach its full potential. This is
885 made slightly complicated by the fact that RVV itself has two
886 levels: Base and reserved future functionality.
887
888 * Representation Encoding is entirely left out of Simple-V in favour of
889 implicitly taking the exact (explicit) meaning from RV Standard Extensions.
890 * VCLIP and VCLIPI do not have corresponding RV Standard Extension
891 opcodes (and are the only such operations).
892 * Extended Element bitwidths (1 through to 24576 bits) were left out
893 of Simple-V as, again, there is no corresponding RV Standard Extension
894 that covers anything even below 32-bit operands.
895 * Polymorphism was entirely left out of Simple-V due to the inherent
896 complexity of automatic type-conversion.
897 * Vector Register files were specifically left out of Simple-V in favour
898 of fitting on top of the integer and floating-point files. An
899 "RVV re-retro-fit" needs to be able to mark (implicitly marked)
900 registers as being actually in a separate *vector* register file.
901 * Fortunately in RVV (Draft 0.4, V2.3-Draft), the "base" vector
902 register file size is 5 bits (32 registers), whilst the "Extended"
903 variant of RVV specifies 8 bits (256 registers) and has yet to
904 be published.
905 * One big difference: Sections 17.12 and 17.17, there are only two possible
906 predication registers in RVV "Base". Through the "indirect" method,
907 Simple-V provides a key-value CSR table that allows (arbitrarily)
908 up to 16 (TBD) of either the floating-point or integer registers to
909 be marked as "predicated" (key), and if so, which integer register to
910 use as the predication mask (value).
911
912 **TODO**
913
914 # Implementing P (renamed to DSP) on top of Simple-V
915
916 * Implementors indicate chosen bitwidth support in Vector-bitwidth CSR
917 (caveat: anything not specified drops through to software-emulation / traps)
918 * TODO
919
920 # Appendix
921
922 ## V-Extension to Simple-V Comparative Analysis
923
924 This section has been moved to its own page [[v_comparative_analysis]]
925
926 ## P-Ext ISA
927
928 This section has been moved to its own page [[p_comparative_analysis]]
929
930 ## Comparison of "Traditional" SIMD, Alt-RVP, Simple-V and RVV Proposals <a name="parallelism_comparisons"></a>
931
932 This section compares the various parallelism proposals as they stand,
933 including traditional SIMD, in terms of features, ease of implementation,
934 complexity, flexibility, and die area.
935
936 ### [[harmonised_rvv_rvp]]
937
938 This is an interesting proposal under development to retro-fit the AndesStar
939 P-Ext into V-Ext.
940
941 ### [[alt_rvp]]
942
943 Primary benefit of Alt-RVP is the simplicity with which parallelism
944 may be introduced (effective multiplication of regfiles and associated ALUs).
945
946 * plus: the simplicity of the lanes (combined with the regularity of
947 allocating identical opcodes multiple independent registers) meaning
948 that SRAM or 2R1W can be used for entire regfile (potentially).
949 * minus: a more complex instruction set where the parallelism is much
950 more explicitly directly specified in the instruction and
951 * minus: if you *don't* have an explicit instruction (opcode) and you
952 need one, the only place it can be added is... in the vector unit and
953 * minus: opcode functions (and associated ALUs) duplicated in Alt-RVP are
954 not useable or accessible in other Extensions.
955 * plus-and-minus: Lanes may be utilised for high-speed context-switching
956 but with the down-side that they're an all-or-nothing part of the Extension.
957 No Alt-RVP: no fast register-bank switching.
958 * plus: Lane-switching would mean that complex operations not suited to
959 parallelisation can be carried out, followed by further parallel Lane-based
960 work, without moving register contents down to memory (and back)
961 * minus: Access to registers across multiple lanes is challenging. "Solution"
962 is to drop data into memory and immediately back in again (like MMX).
963
964 ### Simple-V
965
966 Primary benefit of Simple-V is the OO abstraction of parallel principles
967 from actual (internal) parallel hardware. It's an API in effect that's
968 designed to be slotted in to an existing implementation (just after
969 instruction decode) with minimum disruption and effort.
970
971 * minus: the complexity (if full parallelism is to be exploited)
972 of having to use register renames, OoO, VLIW, register file cacheing,
973 all of which has been done before but is a pain
974 * plus: transparent re-use of existing opcodes as-is just indirectly
975 saying "this register's now a vector" which
976 * plus: means that future instructions also get to be inherently
977 parallelised because there's no "separate vector opcodes"
978 * plus: Compressed instructions may also be (indirectly) parallelised
979 * minus: the indirect nature of Simple-V means that setup (setting
980 a CSR register to indicate vector length, a separate one to indicate
981 that it is a predicate register and so on) means a little more setup
982 time than Alt-RVP or RVV's "direct and within the (longer) instruction"
983 approach.
984 * plus: shared register file meaning that, like Alt-RVP, complex
985 operations not suited to parallelisation may be carried out interleaved
986 between parallelised instructions *without* requiring data to be dropped
987 down to memory and back (into a separate vectorised register engine).
988 * plus-and-maybe-minus: re-use of integer and floating-point 32-wide register
989 files means that huge parallel workloads would use up considerable
990 chunks of the register file. However in the case of RV64 and 32-bit
991 operations, that effectively means 64 slots are available for parallel
992 operations.
993 * plus: inherent parallelism (actual parallel ALUs) doesn't actually need to
994 be added, yet the instruction opcodes remain unchanged (and still appear
995 to be parallel). consistent "API" regardless of actual internal parallelism:
996 even an in-order single-issue implementation with a single ALU would still
997 appear to have parallel vectoristion.
998 * hard-to-judge: if actual inherent underlying ALU parallelism is added it's
999 hard to say if there would be pluses or minuses (on die area). At worse it
1000 would be "no worse" than existing register renaming, OoO, VLIW and register
1001 file cacheing schemes.
1002
1003 ### RVV (as it stands, Draft 0.4 Section 17, RISC-V ISA V2.3-Draft)
1004
1005 RVV is extremely well-designed and has some amazing features, including
1006 2D reorganisation of memory through LOAD/STORE "strides".
1007
1008 * plus: regular predictable workload means that implementations may
1009 streamline effects on L1/L2 Cache.
1010 * plus: regular and clear parallel workload also means that lanes
1011 (similar to Alt-RVP) may be used as an implementation detail,
1012 using either SRAM or 2R1W registers.
1013 * plus: separate engine with no impact on the rest of an implementation
1014 * minus: separate *complex* engine with no RTL (ALUs, Pipeline stages) reuse
1015 really feasible.
1016 * minus: no ISA abstraction or re-use either: additions to other Extensions
1017 do not gain parallelism, resulting in prolific duplication of functionality
1018 inside RVV *and out*.
1019 * minus: when operations require a different approach (scalar operations
1020 using the standard integer or FP regfile) an entire vector must be
1021 transferred out to memory, into standard regfiles, then back to memory,
1022 then back to the vector unit, this to occur potentially multiple times.
1023 * minus: will never fit into Compressed instruction space (as-is. May
1024 be able to do so if "indirect" features of Simple-V are partially adopted).
1025 * plus-and-slight-minus: extended variants may address up to 256
1026 vectorised registers (requires 48/64-bit opcodes to do it).
1027 * minus-and-partial-plus: separate engine plus complexity increases
1028 implementation time and die area, meaning that adoption is likely only
1029 to be in high-performance specialist supercomputing (where it will
1030 be absolutely superb).
1031
1032 ### Traditional SIMD
1033
1034 The only really good things about SIMD are how easy it is to implement and
1035 get good performance. Unfortunately that makes it quite seductive...
1036
1037 * plus: really straightforward, ALU basically does several packed operations
1038 at once. Parallelism is inherent at the ALU, making the addition of
1039 SIMD-style parallelism an easy decision that has zero significant impact
1040 on the rest of any given architectural design and layout.
1041 * plus (continuation): SIMD in simple in-order single-issue designs can
1042 therefore result in superb throughput, easily achieved even with a very
1043 simple execution model.
1044 * minus: ridiculously complex setup and corner-cases that disproportionately
1045 increase instruction count on what would otherwise be a "simple loop",
1046 should the number of elements in an array not happen to exactly match
1047 the SIMD group width.
1048 * minus: getting data usefully out of registers (if separate regfiles
1049 are used) means outputting to memory and back.
1050 * minus: quite a lot of supplementary instructions for bit-level manipulation
1051 are needed in order to efficiently extract (or prepare) SIMD operands.
1052 * minus: MASSIVE proliferation of ISA both in terms of opcodes in one
1053 dimension and parallelism (width): an at least O(N^2) and quite probably
1054 O(N^3) ISA proliferation that often results in several thousand
1055 separate instructions. all requiring separate and distinct corner-case
1056 algorithms!
1057 * minus: EVEN BIGGER proliferation of SIMD ISA if the functionality of
1058 8, 16, 32 or 64-bit reordering is built-in to the SIMD instruction.
1059 For example: add (high|low) 16-bits of r1 to (low|high) of r2 requires
1060 four separate and distinct instructions: one for (r1:low r2:high),
1061 one for (r1:high r2:low), one for (r1:high r2:high) and one for
1062 (r1:low r2:low) *per function*.
1063 * minus: EVEN BIGGER proliferation of SIMD ISA if there is a mismatch
1064 between operand and result bit-widths. In combination with high/low
1065 proliferation the situation is made even worse.
1066 * minor-saving-grace: some implementations *may* have predication masks
1067 that allow control over individual elements within the SIMD block.
1068
1069 ## Comparison *to* Traditional SIMD: Alt-RVP, Simple-V and RVV Proposals <a name="simd_comparison"></a>
1070
1071 This section compares the various parallelism proposals as they stand,
1072 *against* traditional SIMD as opposed to *alongside* SIMD. In other words,
1073 the question is asked "How can each of the proposals effectively implement
1074 (or replace) SIMD, and how effective would they be"?
1075
1076 ### [[alt_rvp]]
1077
1078 * Alt-RVP would not actually replace SIMD but would augment it: just as with
1079 a SIMD architecture where the ALU becomes responsible for the parallelism,
1080 Alt-RVP ALUs would likewise be so responsible... with *additional*
1081 (lane-based) parallelism on top.
1082 * Thus at least some of the downsides of SIMD ISA O(N^5) proliferation by
1083 at least one dimension are avoided (architectural upgrades introducing
1084 128-bit then 256-bit then 512-bit variants of the exact same 64-bit
1085 SIMD block)
1086 * Thus, unfortunately, Alt-RVP would suffer the same inherent proliferation
1087 of instructions as SIMD, albeit not quite as badly (due to Lanes).
1088 * In the same discussion for Alt-RVP, an additional proposal was made to
1089 be able to subdivide the bits of each register lane (columns) down into
1090 arbitrary bit-lengths (RGB 565 for example).
1091 * A recommendation was given instead to make the subdivisions down to 32-bit,
1092 16-bit or even 8-bit, effectively dividing the registerfile into
1093 Lane0(H), Lane0(L), Lane1(H) ... LaneN(L) or further. If inter-lane
1094 "swapping" instructions were then introduced, some of the disadvantages
1095 of SIMD could be mitigated.
1096
1097 ### RVV
1098
1099 * RVV is designed to replace SIMD with a better paradigm: arbitrary-length
1100 parallelism.
1101 * However whilst SIMD is usually designed for single-issue in-order simple
1102 DSPs with a focus on Multimedia (Audio, Video and Image processing),
1103 RVV's primary focus appears to be on Supercomputing: optimisation of
1104 mathematical operations that fit into the OpenCL space.
1105 * Adding functions (operations) that would normally fit (in parallel)
1106 into a SIMD instruction requires an equivalent to be added to the
1107 RVV Extension, if one does not exist. Given the specialist nature of
1108 some SIMD instructions (8-bit or 16-bit saturated or halving add),
1109 this possibility seems extremely unlikely to occur, even if the
1110 implementation overhead of RVV were acceptable (compared to
1111 normal SIMD/DSP-style single-issue in-order simplicity).
1112
1113 ### Simple-V
1114
1115 * Simple-V borrows hugely from RVV as it is intended to be easy to
1116 topologically transplant every single instruction from RVV (as
1117 designed) into Simple-V equivalents, with *zero loss of functionality
1118 or capability*.
1119 * With the "parallelism" abstracted out, a hypothetical SIMD-less "DSP"
1120 Extension which contained the basic primitives (non-parallelised
1121 8, 16 or 32-bit SIMD operations) inherently *become* parallel,
1122 automatically.
1123 * Additionally, standard operations (ADD, MUL) that would normally have
1124 to have special SIMD-parallel opcodes added need no longer have *any*
1125 of the length-dependent variants (2of 32-bit ADDs in a 64-bit register,
1126 4of 32-bit ADDs in a 128-bit register) because Simple-V takes the
1127 *standard* RV opcodes (present and future) and automatically parallelises
1128 them.
1129 * By inheriting the RVV feature of arbitrary vector-length, then just as
1130 with RVV the corner-cases and ISA proliferation of SIMD is avoided.
1131 * Whilst not entirely finalised, registers are expected to be
1132 capable of being subdivided down to an implementor-chosen bitwidth
1133 in the underlying hardware (r1 becomes r1[31..24] r1[23..16] r1[15..8]
1134 and r1[7..0], or just r1[31..16] r1[15..0]) where implementors can
1135 choose to have separate independent 8-bit ALUs or dual-SIMD 16-bit
1136 ALUs that perform twin 8-bit operations as they see fit, or anything
1137 else including no subdivisions at all.
1138 * Even though implementors have that choice even to have full 64-bit
1139 (with RV64) SIMD, they *must* provide predication that transparently
1140 switches off appropriate units on the last loop, thus neatly fitting
1141 underlying SIMD ALU implementations *into* the arbitrary vector-length
1142 RVV paradigm, keeping the uniform consistent API that is a key strategic
1143 feature of Simple-V.
1144 * With Simple-V fitting into the standard register files, certain classes
1145 of SIMD operations such as High/Low arithmetic (r1[31..16] + r2[15..0])
1146 can be done by applying *Parallelised* Bit-manipulation operations
1147 followed by parallelised *straight* versions of element-to-element
1148 arithmetic operations, even if the bit-manipulation operations require
1149 changing the bitwidth of the "vectors" to do so. Predication can
1150 be utilised to skip high words (or low words) in source or destination.
1151 * In essence, the key downside of SIMD - massive duplication of
1152 identical functions over time as an architecture evolves from 32-bit
1153 wide SIMD all the way up to 512-bit, is avoided with Simple-V, through
1154 vector-style parallelism being dropped on top of 8-bit or 16-bit
1155 operations, all the while keeping a consistent ISA-level "API" irrespective
1156 of implementor design choices (or indeed actual implementations).
1157
1158 ### Example Instruction translation: <a name="example_translation"></a>
1159
1160 Instructions "ADD r2 r4 r4" would result in three instructions being
1161 generated and placed into the FIFO:
1162
1163 * ADD r2 r4 r4
1164 * ADD r2 r5 r5
1165 * ADD r2 r6 r6
1166
1167 ## Example of vector / vector, vector / scalar, scalar / scalar => vector add
1168
1169 register CSRvectorlen[XLEN][4]; # not quite decided yet about this one...
1170 register CSRpredicate[XLEN][4]; # 2^4 is max vector length
1171 register CSRreg_is_vectorised[XLEN]; # just for fun support scalars as well
1172 register x[32][XLEN];
1173
1174 function op_add(rd, rs1, rs2, predr)
1175 {
1176    /* note that this is ADD, not PADD */
1177    int i, id, irs1, irs2;
1178    # checks CSRvectorlen[rd] == CSRvectorlen[rs] etc. ignored
1179    # also destination makes no sense as a scalar but what the hell...
1180    for (i = 0, id=0, irs1=0, irs2=0; i<CSRvectorlen[rd]; i++)
1181       if (CSRpredicate[predr][i]) # i *think* this is right...
1182          x[rd+id] <= x[rs1+irs1] + x[rs2+irs2];
1183       # now increment the idxs
1184       if (CSRreg_is_vectorised[rd]) # bitfield check rd, scalar/vector?
1185          id += 1;
1186       if (CSRreg_is_vectorised[rs1]) # bitfield check rs1, scalar/vector?
1187          irs1 += 1;
1188       if (CSRreg_is_vectorised[rs2]) # bitfield check rs2, scalar/vector?
1189          irs2 += 1;
1190 }
1191
1192 ## Retro-fitting Predication into branch-explicit ISA <a name="predication_retrofit"></a>
1193
1194 One of the goals of this parallelism proposal is to avoid instruction
1195 duplication. However, with the base ISA having been designed explictly
1196 to *avoid* condition-codes entirely, shoe-horning predication into it
1197 bcomes quite challenging.
1198
1199 However what if all branch instructions, if referencing a vectorised
1200 register, were instead given *completely new analogous meanings* that
1201 resulted in a parallel bit-wise predication register being set? This
1202 would have to be done for both C.BEQZ and C.BNEZ, as well as BEQ, BNE,
1203 BLT and BGE.
1204
1205 We might imagine that FEQ, FLT and FLT would also need to be converted,
1206 however these are effectively *already* in the precise form needed and
1207 do not need to be converted *at all*! The difference is that FEQ, FLT
1208 and FLE *specifically* write a 1 to an integer register if the condition
1209 holds, and 0 if not. All that needs to be done here is to say, "if
1210 the integer register is tagged with a bit that says it is a predication
1211 register, the **bit** in the integer register is set based on the
1212 current vector index" instead.
1213
1214 There is, in the standard Conditional Branch instruction, more than
1215 adequate space to interpret it in a similar fashion:
1216
1217 [[!table data="""
1218 31 |30 ..... 25 |24..20|19..15| 14...12| 11.....8 | 7 | 6....0 |
1219 imm[12] | imm[10:5] |rs2 | rs1 | funct3 | imm[4:1] | imm[11] | opcode |
1220 1 | 6 | 5 | 5 | 3 | 4 | 1 | 7 |
1221 offset[12,10:5] || src2 | src1 | BEQ | offset[11,4:1] || BRANCH |
1222 """]]
1223
1224 This would become:
1225
1226 [[!table data="""
1227 31 | 30 .. 25 |24 ... 20 | 19 15 | 14 12 | 11 .. 8 | 7 | 6 ... 0 |
1228 imm[12] | imm[10:5]| rs2 | rs1 | funct3 | imm[4:1] | imm[11] | opcode |
1229 1 | 6 | 5 | 5 | 3 | 4 | 1 | 7 |
1230 reserved || src2 | src1 | BEQ | predicate rs3 || BRANCH |
1231 """]]
1232
1233 Similarly the C.BEQZ and C.BNEZ instruction format may be retro-fitted,
1234 with the interesting side-effect that there is space within what is presently
1235 the "immediate offset" field to reinterpret that to add in not only a bit
1236 field to distinguish between floating-point compare and integer compare,
1237 not only to add in a second source register, but also use some of the bits as
1238 a predication target as well.
1239
1240 [[!table data="""
1241 15..13 | 12 ....... 10 | 9...7 | 6 ......... 2 | 1 .. 0 |
1242 funct3 | imm | rs10 | imm | op |
1243 3 | 3 | 3 | 5 | 2 |
1244 C.BEQZ | offset[8,4:3] | src | offset[7:6,2:1,5] | C1 |
1245 """]]
1246
1247 Now uses the CS format:
1248
1249 [[!table data="""
1250 15..13 | 12 . 10 | 9 .. 7 | 6 .. 5 | 4..2 | 1 .. 0 |
1251 funct3 | imm | rs10 | imm | | op |
1252 3 | 3 | 3 | 2 | 3 | 2 |
1253 C.BEQZ | pred rs3 | src1 | I/F B | src2 | C1 |
1254 """]]
1255
1256 Bit 6 would be decoded as "operation refers to Integer or Float" including
1257 interpreting src1 and src2 accordingly as outlined in Table 12.2 of the
1258 "C" Standard, version 2.0,
1259 whilst Bit 5 would allow the operation to be extended, in combination with
1260 funct3 = 110 or 111: a combination of four distinct (predicated) comparison
1261 operators. In both floating-point and integer cases those could be
1262 EQ/NEQ/LT/LE (with GT and GE being synthesised by inverting src1 and src2).
1263
1264 ## Register reordering <a name="register_reordering"></a>
1265
1266 ### Register File
1267
1268 | Reg Num | Bits |
1269 | ------- | ---- |
1270 | r0 | (32..0) |
1271 | r1 | (32..0) |
1272 | r2 | (32..0) |
1273 | r3 | (32..0) |
1274 | r4 | (32..0) |
1275 | r5 | (32..0) |
1276 | r6 | (32..0) |
1277 | r7 | (32..0) |
1278 | .. | (32..0) |
1279 | r31| (32..0) |
1280
1281 ### Vectorised CSR
1282
1283 May not be an actual CSR: may be generated from Vector Length CSR:
1284 single-bit is less burdensome on instruction decode phase.
1285
1286 | 7 | 6 | 5 | 4 | 3 | 2 | 1 | 0 |
1287 | - | - | - | - | - | - | - | - |
1288 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 1 |
1289
1290 ### Vector Length CSR
1291
1292 | Reg Num | (3..0) |
1293 | ------- | ---- |
1294 | r0 | 2 |
1295 | r1 | 0 |
1296 | r2 | 1 |
1297 | r3 | 1 |
1298 | r4 | 3 |
1299 | r5 | 0 |
1300 | r6 | 0 |
1301 | r7 | 1 |
1302
1303 ### Virtual Register Reordering
1304
1305 This example assumes the above Vector Length CSR table
1306
1307 | Reg Num | Bits (0) | Bits (1) | Bits (2) |
1308 | ------- | -------- | -------- | -------- |
1309 | r0 | (32..0) | (32..0) |
1310 | r2 | (32..0) |
1311 | r3 | (32..0) |
1312 | r4 | (32..0) | (32..0) | (32..0) |
1313 | r7 | (32..0) |
1314
1315 ### Bitwidth Virtual Register Reordering
1316
1317 This example goes a little further and illustrates the effect that a
1318 bitwidth CSR has been set on a register. Preconditions:
1319
1320 * RV32 assumed
1321 * CSRintbitwidth[2] = 010 # integer r2 is 16-bit
1322 * CSRintvlength[2] = 3 # integer r2 is a vector of length 3
1323 * vsetl rs1, 5 # set the vector length to 5
1324
1325 This is interpreted as follows:
1326
1327 * Given that the context is RV32, ELEN=32.
1328 * With ELEN=32 and bitwidth=16, the number of SIMD elements is 2
1329 * Therefore the actual vector length is up to *six* elements
1330 * However vsetl sets a length 5 therefore the last "element" is skipped
1331
1332 So when using an operation that uses r2 as a source (or destination)
1333 the operation is carried out as follows:
1334
1335 * 16-bit operation on r2(15..0) - vector element index 0
1336 * 16-bit operation on r2(31..16) - vector element index 1
1337 * 16-bit operation on r3(15..0) - vector element index 2
1338 * 16-bit operation on r3(31..16) - vector element index 3
1339 * 16-bit operation on r4(15..0) - vector element index 4
1340 * 16-bit operation on r4(31..16) **NOT** carried out due to length being 5
1341
1342 Predication has been left out of the above example for simplicity, however
1343 predication is ANDed with the latter stages (vsetl not equal to maximum
1344 capacity).
1345
1346 Note also that it is entirely an implementor's choice as to whether to have
1347 actual separate ALUs down to the minimum bitwidth, or whether to have something
1348 more akin to traditional SIMD (at any level of subdivision: 8-bit SIMD
1349 operations carried out 32-bits at a time is perfectly acceptable, as is
1350 8-bit SIMD operations carried out 16-bits at a time requiring two ALUs).
1351 Regardless of the internal parallelism choice, *predication must
1352 still be respected*, making Simple-V in effect the "consistent public API".
1353
1354 vew may be one of the following (giving a table "bytestable", used below):
1355
1356 | vew | bitwidth | bytestable |
1357 | --- | -------- | ---------- |
1358 | 000 | default | XLEN/8 |
1359 | 001 | 8 | 1 |
1360 | 010 | 16 | 2 |
1361 | 011 | 32 | 4 |
1362 | 100 | 64 | 8 |
1363 | 101 | 128 | 16 |
1364 | 110 | rsvd | rsvd |
1365 | 111 | rsvd | rsvd |
1366
1367 Pseudocode for vector length taking CSR SIMD-bitwidth into account:
1368
1369 vew = CSRbitwidth[rs1]
1370 if (vew == 0)
1371 bytesperreg = (XLEN/8) # or FLEN as appropriate
1372 else:
1373 bytesperreg = bytestable[vew] # 1 2 4 8 16
1374 simdmult = (XLEN/8) / bytesperreg # or FLEN as appropriate
1375 vlen = CSRvectorlen[rs1] * simdmult
1376
1377 To index an element in a register rnum where the vector element index is i:
1378
1379 function regoffs(rnum, i):
1380 regidx = floor(i / simdmult) # integer-div rounded down
1381 byteidx = i % simdmult # integer-remainder
1382 return rnum + regidx, # actual real register
1383 byteidx * 8, # low
1384 byteidx * 8 + (vew-1), # high
1385
1386 ### Insights
1387
1388 SIMD register file splitting still to consider. For RV64, benefits of doubling
1389 (quadrupling in the case of Half-Precision IEEE754 FP) the apparent
1390 size of the floating point register file to 64 (128 in the case of HP)
1391 seem pretty clear and worth the complexity.
1392
1393 64 virtual 32-bit F.P. registers and given that 32-bit FP operations are
1394 done on 64-bit registers it's not so conceptually difficult.  May even
1395 be achieved by *actually* splitting the regfile into 64 virtual 32-bit
1396 registers such that a 64-bit FP scalar operation is dropped into (r0.H
1397 r0.L) tuples.  Implementation therefore hidden through register renaming.
1398
1399 Implementations intending to introduce VLIW, OoO and parallelism
1400 (even without Simple-V) would then find that the instructions are
1401 generated quicker (or in a more compact fashion that is less heavy
1402 on caches). Interestingly we observe then that Simple-V is about
1403 "consolidation of instruction generation", where actual parallelism
1404 of underlying hardware is an implementor-choice that could just as
1405 equally be applied *without* Simple-V even being implemented.
1406
1407 ## Analysis of CSR decoding on latency <a name="csr_decoding_analysis"></a>
1408
1409 It could indeed have been logically deduced (or expected), that there
1410 would be additional decode latency in this proposal, because if
1411 overloading the opcodes to have different meanings, there is guaranteed
1412 to be some state, some-where, directly related to registers.
1413
1414 There are several cases:
1415
1416 * All operands vector-length=1 (scalars), all operands
1417 packed-bitwidth="default": instructions are passed through direct as if
1418 Simple-V did not exist.  Simple-V is, in effect, completely disabled.
1419 * At least one operand vector-length > 1, all operands
1420 packed-bitwidth="default": any parallel vector ALUs placed on "alert",
1421 virtual parallelism looping may be activated.
1422 * All operands vector-length=1 (scalars), at least one
1423 operand packed-bitwidth != default: degenerate case of SIMD,
1424 implementation-specific complexity here (packed decode before ALUs or
1425 *IN* ALUs)
1426 * At least one operand vector-length > 1, at least one operand
1427 packed-bitwidth != default: parallel vector ALUs (if any)
1428 placed on "alert", virtual parallelsim looping may be activated,
1429 implementation-specific SIMD complexity kicks in (packed decode before
1430 ALUs or *IN* ALUs).
1431
1432 Bear in mind that the proposal includes that the decision whether
1433 to parallelise in hardware or whether to virtual-parallelise (to
1434 dramatically simplify compilers and also not to run into the SIMD
1435 instruction proliferation nightmare) *or* a transprent combination
1436 of both, be done on a *per-operand basis*, so that implementors can
1437 specifically choose to create an application-optimised implementation
1438 that they believe (or know) will sell extremely well, without having
1439 "Extra Standards-Mandated Baggage" that would otherwise blow their area
1440 or power budget completely out the window.
1441
1442 Additionally, two possible CSR schemes have been proposed, in order to
1443 greatly reduce CSR space:
1444
1445 * per-register CSRs (vector-length and packed-bitwidth)
1446 * a smaller number of CSRs with the same information but with an *INDEX*
1447 specifying WHICH register in one of three regfiles (vector, fp, int)
1448 the length and bitwidth applies to.
1449
1450 (See "CSR vector-length and CSR SIMD packed-bitwidth" section for details)
1451
1452 In addition, LOAD/STORE has its own associated proposed CSRs that
1453 mirror the STRIDE (but not yet STRIDE-SEGMENT?) functionality of
1454 V (and Hwacha).
1455
1456 Also bear in mind that, for reasons of simplicity for implementors,
1457 I was coming round to the idea of permitting implementors to choose
1458 exactly which bitwidths they would like to support in hardware and which
1459 to allow to fall through to software-trap emulation.
1460
1461 So the question boils down to:
1462
1463 * whether either (or both) of those two CSR schemes have significant
1464 latency that could even potentially require an extra pipeline decode stage
1465 * whether there are implementations that can be thought of which do *not*
1466 introduce significant latency
1467 * whether it is possible to explicitly (through quite simply
1468 disabling Simple-V-Ext) or implicitly (detect the case all-vlens=1,
1469 all-simd-bitwidths=default) switch OFF any decoding, perhaps even to
1470 the extreme of skipping an entire pipeline stage (if one is needed)
1471 * whether packed bitwidth and associated regfile splitting is so complex
1472 that it should definitely, definitely be made mandatory that implementors
1473 move regfile splitting into the ALU, and what are the implications of that
1474 * whether even if that *is* made mandatory, is software-trapped
1475 "unsupported bitwidths" still desirable, on the basis that SIMD is such
1476 a complete nightmare that *even* having a software implementation is
1477 better, making Simple-V have more in common with a software API than
1478 anything else.
1479
1480 Whilst the above may seem to be severe minuses, there are some strong
1481 pluses:
1482
1483 * Significant reduction of V's opcode space: over 95%.
1484 * Smaller reduction of P's opcode space: around 10%.
1485 * The potential to use Compressed instructions in both Vector and SIMD
1486 due to the overloading of register meaning (implicit vectorisation,
1487 implicit packing)
1488 * Not only present but also future extensions automatically gain parallelism.
1489 * Already mentioned but worth emphasising: the simplification to compiler
1490 writers and assembly-level writers of having the same consistent ISA
1491 regardless of whether the internal level of parallelism (number of
1492 parallel ALUs) is only equal to one ("virtual" parallelism), or is
1493 greater than one, should not be underestimated.
1494
1495 ## Reducing Register Bank porting
1496
1497 This looks quite reasonable.
1498 <https://www.princeton.edu/~rblee/ELE572Papers/MultiBankRegFile_ISCA2000.pdf>
1499
1500 The main details are outlined on page 4.  They propose a 2-level register
1501 cache hierarchy, note that registers are typically only read once, that
1502 you never write back from upper to lower cache level but always go in a
1503 cycle lower -> upper -> ALU -> lower, and at the top of page 5 propose
1504 a scheme where you look ahead by only 2 instructions to determine which
1505 registers to bring into the cache.
1506
1507 The nice thing about a vector architecture is that you *know* that
1508 *even more* registers are going to be pulled in: Hwacha uses this fact
1509 to optimise L1/L2 cache-line usage (avoid thrashing), strangely enough
1510 by *introducing* deliberate latency into the execution phase.
1511
1512 ## Overflow registers in combination with predication
1513
1514 **TODO**: propose overflow registers be actually one of the integer regs
1515 (flowing to multiple regs).
1516
1517 **TODO**: propose "mask" (predication) registers likewise. combination with
1518 standard RV instructions and overflow registers extremely powerful, see
1519 Aspex ASP.
1520
1521 When integer overflow is stored in an easily-accessible bit (or another
1522 register), parallelisation turns this into a group of bits which can
1523 potentially be interacted with in predication, in interesting and powerful
1524 ways. For example, by taking the integer-overflow result as a predication
1525 field and shifting it by one, a predicated vectorised "add one" can emulate
1526 "carry" on arbitrary (unlimited) length addition.
1527
1528 However despite RVV having made room for floating-point exceptions, neither
1529 RVV nor base RV have taken integer-overflow (carry) into account, which
1530 makes proposing it quite challenging given that the relevant (Base) RV
1531 sections are frozen. Consequently it makes sense to forgo this feature.
1532
1533 ## Context Switch Example <a name="context_switch"></a>
1534
1535 An unusual side-effect of Simple-V mapping onto the standard register files
1536 is that LOAD-multiple and STORE-multiple are accidentally available, as long
1537 as it is acceptable that the register(s) to be loaded/stored are contiguous
1538 (per instruction). An additional accidental benefit is that Compressed LD/ST
1539 may also be used.
1540
1541 To illustrate how this works, here is some example code from FreeRTOS
1542 (GPLv2 licensed, portasm.S):
1543
1544 /* Macro for saving task context */
1545 .macro portSAVE_CONTEXT
1546 .global pxCurrentTCB
1547 /* make room in stack */
1548 addi sp, sp, -REGBYTES * 32
1549
1550 /* Save Context */
1551 STORE x1, 0x0(sp)
1552 STORE x2, 1 * REGBYTES(sp)
1553 STORE x3, 2 * REGBYTES(sp)
1554 ...
1555 ...
1556 STORE x30, 29 * REGBYTES(sp)
1557 STORE x31, 30 * REGBYTES(sp)
1558
1559 /* Store current stackpointer in task control block (TCB) */
1560 LOAD t0, pxCurrentTCB //pointer
1561 STORE sp, 0x0(t0)
1562 .endm
1563
1564 /* Saves current error program counter (EPC) as task program counter */
1565 .macro portSAVE_EPC
1566 csrr t0, mepc
1567 STORE t0, 31 * REGBYTES(sp)
1568 .endm
1569
1570 /* Saves current return adress (RA) as task program counter */
1571 .macro portSAVE_RA
1572 STORE ra, 31 * REGBYTES(sp)
1573 .endm
1574
1575 /* Macro for restoring task context */
1576 .macro portRESTORE_CONTEXT
1577
1578 .global pxCurrentTCB
1579 /* Load stack pointer from the current TCB */
1580 LOAD sp, pxCurrentTCB
1581 LOAD sp, 0x0(sp)
1582
1583 /* Load task program counter */
1584 LOAD t0, 31 * REGBYTES(sp)
1585 csrw mepc, t0
1586
1587 /* Run in machine mode */
1588 li t0, MSTATUS_PRV1
1589 csrs mstatus, t0
1590
1591 /* Restore registers,
1592 Skip global pointer because that does not change */
1593 LOAD x1, 0x0(sp)
1594 LOAD x4, 3 * REGBYTES(sp)
1595 LOAD x5, 4 * REGBYTES(sp)
1596 ...
1597 ...
1598 LOAD x30, 29 * REGBYTES(sp)
1599 LOAD x31, 30 * REGBYTES(sp)
1600
1601 addi sp, sp, REGBYTES * 32
1602 mret
1603 .endm
1604
1605 The important bits are the Load / Save context, which may be replaced
1606 with firstly setting up the Vectors and secondly using a *single* STORE
1607 (or LOAD) including using C.ST or C.LD, to indicate that the entire
1608 bank of registers is to be loaded/saved:
1609
1610 /* a few things are assumed here: (a) that when switching to
1611 M-Mode an entirely different set of CSRs is used from that
1612 which is used in U-Mode and (b) that the M-Mode x1 and x4
1613 vectors are also not used anywhere else in M-Mode, consequently
1614 only need to be set up just the once.
1615 */
1616 .macroVectorSetup
1617 MVECTORCSRx1 = 31, defaultlen
1618 MVECTORCSRx4 = 28, defaultlen
1619
1620 /* Save Context */
1621 SETVL x0, x0, 31 /* x0 ignored silently */
1622 STORE x1, 0x0(sp) // x1 marked as 31-long vector of default bitwidth
1623
1624 /* Restore registers,
1625 Skip global pointer because that does not change */
1626 LOAD x1, 0x0(sp)
1627 SETVL x0, x0, 28 /* x0 ignored silently */
1628 LOAD x4, 3 * REGBYTES(sp) // x4 marked as 28-long default bitwidth
1629
1630 Note that although it may just be a bug in portasm.S, x2 and x3 appear not
1631 to be being restored. If however this is a bug and they *do* need to be
1632 restored, then the SETVL call may be moved to *outside* the Save / Restore
1633 Context assembly code, into the macroVectorSetup, as long as vectors are
1634 never used anywhere else (i.e. VL is never altered by M-Mode).
1635
1636 In effect the entire bank of repeated LOAD / STORE instructions is replaced
1637 by one single (compressed if it is available) instruction.
1638
1639 ## Virtual Memory page-faults on LOAD/STORE
1640
1641
1642 ### Notes from conversations
1643
1644 > I was going through the C.LOAD / C.STORE section 12.3 of V2.3-Draft
1645 > riscv-isa-manual in order to work out how to re-map RVV onto the standard
1646 > ISA, and came across an interesting comments at the bottom of pages 75
1647 > and 76:
1648
1649 > " A common mechanism used in other ISAs to further reduce save/restore
1650 > code size is load- multiple and store-multiple instructions. "
1651
1652 > Fascinatingly, due to Simple-V proposing to use the *standard* register
1653 > file, both C.LOAD / C.STORE *and* LOAD / STORE would in effect be exactly
1654 > that: load-multiple and store-multiple instructions. Which brings us
1655 > on to this comment:
1656
1657 > "For virtual memory systems, some data accesses could be resident in
1658 > physical memory and
1659 > some could not, which requires a new restart mechanism for partially
1660 > executed instructions."
1661
1662 > Which then of course brings us to the interesting question: how does RVV
1663 > cope with the scenario when, particularly with LD.X (Indexed / indirect
1664 > loads), part-way through the loading a page fault occurs?
1665
1666 > Has this been noted or discussed before?
1667
1668 For applications-class platforms, the RVV exception model is
1669 element-precise (that is, if an exception occurs on element j of a
1670 vector instruction, elements 0..j-1 have completed execution and elements
1671 j+1..vl-1 have not executed).
1672
1673 Certain classes of embedded platforms where exceptions are always fatal
1674 might choose to offer resumable/swappable interrupts but not precise
1675 exceptions.
1676
1677
1678 > Is RVV designed in any way to be re-entrant?
1679
1680 Yes.
1681
1682
1683 > What would the implications be for instructions that were in a FIFO at
1684 > the time, in out-of-order and VLIW implementations, where partial decode
1685 > had taken place?
1686
1687 The usual bag of tricks for maintaining precise exceptions applies to
1688 vector machines as well. Register renaming makes the job easier, and
1689 it's relatively cheaper for vectors, since the control cost is amortized
1690 over longer registers.
1691
1692
1693 > Would it be reasonable at least to say *bypass* (and freeze) the
1694 > instruction FIFO (drop down to a single-issue execution model temporarily)
1695 > for the purposes of executing the instructions in the interrupt (whilst
1696 > setting up the VM page), then re-continue the instruction with all
1697 > state intact?
1698
1699 This approach has been done successfully, but it's desirable to be
1700 able to swap out the vector unit state to support context switches on
1701 exceptions that result in long-latency I/O.
1702
1703
1704 > Or would it be better to switch to an entirely separate secondary
1705 > hyperthread context?
1706
1707 > Does anyone have any ideas or know if there is any academic literature
1708 > on solutions to this problem?
1709
1710 The Vector VAX offered imprecise but restartable and swappable exceptions:
1711 http://mprc.pku.edu.cn/~liuxianhua/chn/corpus/Notes/articles/isca/1990/VAX%20vector%20architecture.pdf
1712
1713 Sec. 4.6 of Krste's dissertation assesses some of
1714 the tradeoffs and references a bunch of related work:
1715 http://people.eecs.berkeley.edu/~krste/thesis.pdf
1716
1717
1718 ----
1719
1720 Started reading section 4.6 of Krste's thesis, noted the "IEE85 F.P
1721 exceptions" and thought, "hmmm that could go into a CSR, must re-read
1722 the section on FP state CSRs in RVV 0.4-Draft again" then i suddenly
1723 thought, "ah ha! what if the memory exceptions were, instead of having
1724 an immediate exception thrown, were simply stored in a type of predication
1725 bit-field with a flag "error this element failed"?
1726
1727 Then, *after* the vector load (or store, or even operation) was
1728 performed, you could *then* raise an exception, at which point it
1729 would be possible (yes in software... I know....) to go "hmmm, these
1730 indexed operations didn't work, let's get them into memory by triggering
1731 page-loads", then *re-run the entire instruction* but this time with a
1732 "memory-predication CSR" that stops the already-performed operations
1733 (whether they be loads, stores or an arithmetic / FP operation) from
1734 being carried out a second time.
1735
1736 This theoretically could end up being done multiple times in an SMP
1737 environment, and also for LD.X there would be the remote outside annoying
1738 possibility that the indexed memory address could end up being modified.
1739
1740 The advantage would be that the order of execution need not be
1741 sequential, which potentially could have some big advantages.
1742 Am still thinking through the implications as any dependent operations
1743 (particularly ones already decoded and moved into the execution FIFO)
1744 would still be there (and stalled). hmmm.
1745
1746 ----
1747
1748 > > # assume internal parallelism of 8 and MAXVECTORLEN of 8
1749 > > VSETL r0, 8
1750 > > FADD x1, x2, x3
1751 >
1752 > > x3[0]: ok
1753 > > x3[1]: exception
1754 > > x3[2]: ok
1755 > > ...
1756 > > ...
1757 > > x3[7]: ok
1758 >
1759 > > what happens to result elements 2-7?  those may be *big* results
1760 > > (RV128)
1761 > > or in the RVV-Extended may be arbitrary bit-widths far greater.
1762 >
1763 >  (you replied:)
1764 >
1765 > Thrown away.
1766
1767 discussion then led to the question of OoO architectures
1768
1769 > The costs of the imprecise-exception model are greater than the benefit.
1770 > Software doesn't want to cope with it.  It's hard to debug.  You can't
1771 > migrate state between different microarchitectures--unless you force all
1772 > implementations to support the same imprecise-exception model, which would
1773 > greatly limit implementation flexibility.  (Less important, but still
1774 > relevant, is that the imprecise model increases the size of the context
1775 > structure, as the microarchitectural guts have to be spilled to memory.)
1776
1777 -----
1778
1779 >> >  it just occurred to me that there's another reason why the data
1780 >> > should be left instead of zeroed.  if the standard register file is
1781 >> > used, such that vectorised operations are translated to mean "please
1782 >> > insert multiple register-contiguous operations into the instruction
1783 >> > FIFO" and predication is used to *skip* some of those, then if the
1784 >> > next "vector" operation uses the (standard) registers that were masked
1785 >> > *out* of the previous operation it may proceed without blocking.
1786 >> >
1787 >> >  if however zeroing is made mandatory then that optimisation becomes
1788 >> > flat-out impossible to deploy.
1789 >> >
1790 >> >  whilst i haven't fully thought through the full implications, i
1791 >> > suspect RVV might also be able to benefit by being able to fit more
1792 >> > overlapping operations into the available SRAM by doing something
1793 >> > similar.
1794 >
1795 >
1796 > Luke, this is called density time masking. It doesn’t apply to only your
1797 > model with the “standard register file” is used. it applies to any
1798 > architecture that attempts to speed up by skipping computation and writeback
1799 > of masked elements.
1800 >
1801 > That said, the writing of zeros need not be explicit. It is possible to add
1802 > a “zero bit” per element that, when set, forces a zero to be read from the
1803 > vector (although the underlying storage may have old data). In this case,
1804 > there may be a way to implement DTM as well.
1805
1806
1807
1808 ## Implementation Paradigms <a name="implementation_paradigms"></a>
1809
1810 TODO: assess various implementation paradigms. These are listed roughly
1811 in order of simplicity (minimum compliance, for ultra-light-weight
1812 embedded systems or to reduce design complexity and the burden of
1813 design implementation and compliance, in non-critical areas), right the
1814 way to high-performance systems.
1815
1816 * Full (or partial) software-emulated (via traps): full support for CSRs
1817 required, however when a register is used that is detected (in hardware)
1818 to be vectorised, an exception is thrown.
1819 * Single-issue In-order, reduced pipeline depth (traditional SIMD / DSP)
1820 * In-order 5+ stage pipelines with instruction FIFOs and mild register-renaming
1821 * Out-of-order with instruction FIFOs and aggressive register-renaming
1822 * VLIW
1823
1824 Also to be taken into consideration:
1825
1826 * "Virtual" vectorisation: single-issue loop, no internal ALU parallelism
1827 * Comphrensive vectorisation: FIFOs and internal parallelism
1828 * Hybrid Parallelism
1829
1830 ### Full or partial software-emulation
1831
1832 The absolute, absolute minimal implementation is to provide the full
1833 set of CSRs and detection logic for when any of the source or destination
1834 registers are vectorised. On detection, a trap is thrown, whether it's
1835 a branch, LOAD, STORE, or an arithmetic operation.
1836
1837 Implementors are entirely free to choose whether to allow absolutely every
1838 single operation to be software-emulated, or whether to provide some emulation
1839 and some hardware support. In particular, for an RV32E implementation
1840 where fast context-switching is a requirement (see "Context Switch Example"),
1841 it makes no sense to allow Vectorised-LOAD/STORE to be implemented as an
1842 exception, as every context-switch will result in double-traps.
1843
1844 # TODO Research
1845
1846 > For great floating point DSPs check TI’s C3x, C4X, and C6xx DSPs
1847
1848 Idea: basic simple butterfly swap on a few element indices, primarily targetted
1849 at SIMD / DSP. High-byte low-byte swapping, high-word low-word swapping,
1850 perhaps allow reindexing of permutations up to 4 elements? 8? Reason:
1851 such operations are less costly than a full indexed-shuffle, which requires
1852 a separate instruction cycle.
1853
1854 Predication "all zeros" needs to be "leave alone". Detection of
1855 ADD r1, rs1, rs0 cases result in nop on predication index 0, whereas
1856 ADD r0, rs1, rs2 is actually a desirable copy from r2 into r0.
1857 Destruction of destination indices requires a copy of the entire vector
1858 in advance to avoid.
1859
1860 TBD: floating-point compare and other exception handling
1861
1862 # References
1863
1864 * SIMD considered harmful <https://www.sigarch.org/simd-instructions-considered-harmful/>
1865 * Link to first proposal <https://groups.google.com/a/groups.riscv.org/forum/#!topic/isa-dev/GuukrSjgBH8>
1866 * Recommendation by Jacob Bachmeyer to make zero-overhead loop an
1867 "implicit program-counter" <https://groups.google.com/a/groups.riscv.org/d/msg/isa-dev/vYVi95gF2Mo/SHz6a4_lAgAJ>
1868 * Re-continuing P-Extension proposal <https://groups.google.com/a/groups.riscv.org/forum/#!msg/isa-dev/IkLkQn3HvXQ/SEMyC9IlAgAJ>
1869 * First Draft P-SIMD (DSP) proposal <https://groups.google.com/a/groups.riscv.org/forum/#!topic/isa-dev/vYVi95gF2Mo>
1870 * B-Extension discussion <https://groups.google.com/a/groups.riscv.org/forum/#!topic/isa-dev/zi_7B15kj6s>
1871 * Broadcom VideoCore-IV <https://docs.broadcom.com/docs/12358545>
1872 Figure 2 P17 and Section 3 on P16.
1873 * Hwacha <https://www2.eecs.berkeley.edu/Pubs/TechRpts/2015/EECS-2015-262.html>
1874 * Hwacha <https://www2.eecs.berkeley.edu/Pubs/TechRpts/2015/EECS-2015-263.html>
1875 * Vector Workshop <http://riscv.org/wp-content/uploads/2015/06/riscv-vector-workshop-june2015.pdf>
1876 * Predication <https://groups.google.com/a/groups.riscv.org/forum/#!topic/isa-dev/XoP4BfYSLXA>
1877 * Branch Divergence <https://jbush001.github.io/2014/12/07/branch-divergence-in-parallel-kernels.html>
1878 * Life of Triangles (3D) <https://jbush001.github.io/2016/02/27/life-of-triangle.html>
1879 * Videocore-IV <https://github.com/hermanhermitage/videocoreiv/wiki/VideoCore-IV-3d-Graphics-Pipeline>
1880 * Discussion proposing CSRs that change ISA definition
1881 <https://groups.google.com/a/groups.riscv.org/forum/#!topic/isa-dev/InzQ1wr_3Ak>
1882 * Zero-overhead loops <https://pdfs.semanticscholar.org/dbaa/66985cc730d4b44d79f519e96ec9c43ab5b7.pdf>
1883 * Multi-ported VLIW Register File Implementation <https://ce-publications.et.tudelft.nl/publications/1517_multiple_contexts_in_a_multiported_vliw_register_file_impl.pdf>
1884 * Fast context save/restore proposal <https://groups.google.com/a/groups.riscv.org/d/msgid/isa-dev/57F823FA.6030701%40gmail.com>
1885 * Register File Bank Cacheing <https://www.princeton.edu/~rblee/ELE572Papers/MultiBankRegFile_ISCA2000.pdf>
1886 * Expired Patent on Vector Virtual Memory solutions
1887 <https://patentimages.storage.googleapis.com/fc/f6/e2/2cbee92fcd8743/US5895501.pdf>
1888 * Discussion on RVV "re-entrant" capabilities allowing operations to be
1889 restarted if an exception occurs (VM page-table miss)
1890 <https://groups.google.com/a/groups.riscv.org/d/msg/isa-dev/IuNFitTw9fM/CCKBUlzsAAAJ>
1891 * Dot Product Vector <https://people.eecs.berkeley.edu/~biancolin/papers/arith17.pdf>
1892 * RVV slides 2017 <https://content.riscv.org/wp-content/uploads/2017/12/Wed-1330-RISCVRogerEspasaVEXT-v4.pdf>
1893 * Wavefront skipping using BRAMS <http://www.ece.ubc.ca/~lemieux/publications/severance-fpga2015.pdf>
1894 * Streaming Pipelines <http://www.ece.ubc.ca/~lemieux/publications/severance-fpga2014.pdf>
1895 * Barcelona SIMD Presentation <https://content.riscv.org/wp-content/uploads/2018/05/09.05.2018-9.15-9.30am-RISCV201805-Andes-proposed-P-extension.pdf>