whitespace
[libreriscv.git] / openpower / sv / rfc / ls010.mdwn
1 # RFC ls010 SVP64 Zero-Overhead Loop Prefix Subsystem
2
3 Credits and acknowledgements:
4
5 * Luke Leighton
6 * Jacob Lifshay
7 * Hendrik Boom
8 * Richard Wilbur
9 * Alexandre Oliva
10 * Cesar Strauss
11 * NLnet Foundation, for funding
12 * OpenPOWER Foundation
13 * Paul Mackerras
14 * Toshaan Bharvani
15 * IBM for the Power ISA itself
16
17 Links:
18
19 * <https://bugs.libre-soc.org/show_bug.cgi?id=1045>
20
21 # Introduction
22
23 Simple-V is a type of Vectorisation best described as a "Prefix Loop
24 Subsystem" similar to the 5 decades-old Zilog Z80 `LDIR` instruction and
25 to the 8086 `REP` Prefix instruction. More advanced features are similar
26 to the Z80 `CPIR` instruction. If viewed one-dimensionally as an actual
27 Vector ISA it introduces over 1.5 million 64-bit Vector instructions.
28 SVP64, the instruction format used by Simple-V, is therefore best viewed
29 as an orthogonal RISC-paradigm "Prefixing" subsystem instead.
30
31 Except where explicitly stated all bit numbers remain as in the rest of
32 the Power ISA: in MSB0 form (the bits are numbered from 0 at the MSB on
33 the left and counting up as you move rightwards to the LSB end). All bit
34 ranges are inclusive (so `4:6` means bits 4, 5, and 6, in MSB0 order).
35 **All register numbering and element numbering however is LSB0 ordering**
36 which is a different convention from that used elsewhere in the Power ISA.
37
38 The SVP64 prefix always comes before the suffix in PC order and must be
39 considered an independent "Defined word" that augments the behaviour of
40 the following instruction, but does **not** change the actual Decoding
41 of that following instruction. **All prefixed instructions retain their
42 non-prefixed encoding and definition**.
43
44 Two apparent exceptions to the above hard rule exist: SV Branch-Conditional
45 operations and LD/ST-update "Post-Increment" Mode. Post-Increment
46 was considered sufficiently high priority (significantly reducing hot-loop
47 instruction count) that one bit in the Prefix is reserved for it.
48 Vectorised Branch-Conditional operations "embed" the original Scalar
49 Branch-Conditional behaviour into a much more advanced variant that
50 is highly suited to High-Performance Computation (HPC), Supercomputing,
51 and parallel GPU Workloads.
52
53 *Architectural Resource Allocation note: it is prohibited to accept RFCs
54 which fundamentally violate this hard requirement. Under no circumstances
55 must the Suffix space have an alternate instruction encoding allocated
56 within SVP64 that is entirely different from the non-prefixed Defined
57 Word. Hardware Implementors critically rely on this inviolate guarantee
58 to implement High-Performance Multi-Issue micro-architectures that can
59 sustain 100% throughput*
60
61 Subset implementations in hardware are permitted, as long as certain
62 rules are followed, allowing for full soft-emulation including future
63 revisions. Compliancy Subsets exist to ensure minimum levels of binary
64 interoperability expectations within certain environments.
65
66 ## SVP64 encoding features
67
68 A number of features need to be compacted into a very small space of
69 only 24 bits:
70
71 * Independent per-register Scalar/Vector tagging and range extension on
72 every register
73 * Element width overrides on both source and destination
74 * Predication on both source and destination
75 * Two different sources of predication: INT and CR Fields
76 * SV Modes including saturation (for Audio, Video and DSP), mapreduce,
77 fail-first and predicate-result mode.
78
79 Different classes of operations require different formats. The earlier
80 sections cover the common formats and the four separate modes follow:
81 CR operations (crops), Arithmetic/Logical (termed "normal"), Load/Store
82 and Branch-Conditional.
83
84 ## Definition of Reserved in this spec.
85
86 For the new fields added in SVP64, instructions that have any of their
87 fields set to a reserved value must cause an illegal instruction trap,
88 to allow emulation of future instruction sets, or for subsets of SVP64 to
89 be implemented in hardware and the rest emulated. This includes SVP64
90 SPRs: reading or writing values which are not supported in hardware
91 must also raise illegal instruction traps in order to allow emulation.
92 Unless otherwise stated, reserved values are always all zeros.
93
94 This is unlike OpenPower ISA v3.1, which in many instances does not
95 require a trap if reserved fields are nonzero. Where the standard Power
96 ISA definition is intended the red keyword `RESERVED` is used.
97
98 ## Definition of "UnVectoriseable"
99
100 Any operation that inherently makes no sense if repeated is termed
101 "UnVectoriseable" or "UnVectorised". Examples include `sc` or `sync`
102 which have no registers. `mtmsr` is also classed as UnVectoriseable
103 because there is only one `MSR`.
104
105 ## Register files, elements, and Element-width Overrides
106
107 In the Upper Compliancy Levels of SVP64 the size of the GPR and FPR
108 Register files are expanded from 32 to 128 entries, and the number of
109 CR Fields expanded from CR0-CR7 to CR0-CR127. (Note: A future version
110 of SVP64 is anticipated to extend the VSR register file).
111
112 Memory access remains exactly the same: the effects of `MSR.LE` remain
113 exactly the same, affecting as they already do and remain **only**
114 on the Load and Store memory-register operation byte-order, and having
115 nothing to do with the ordering of the contents of register files or
116 register-register operations.
117
118 To be absolutely clear:
119
120 ```
121 There are no conceptual arithmetic ordering or other changes over the
122 Scalar Power ISA definitions to registers or register files or to
123 arithmetic or Logical Operations beyond element-width subdivision
124 ```
125
126 Element offset
127 numbering is naturally **LSB0-sequentially-incrementing from zero, not
128 MSB0-incrementing** including when element-width overrides are used,
129 at which point the elements progress through each register
130 sequentially from the LSB end
131 (confusingly numbered the highest in MSB0 ordering) and progress
132 incrementally to the MSB end (confusingly numbered the lowest in
133 MSB0 ordering).
134
135 When exclusively using MSB0-numbering, SVP64
136 becomes unnecessarily complex to both express and subsequently understand:
137 the required conditional subtractions from 63,
138 31, 15 and 7 needed to express the fact that elements are LSB0-sequential
139 unfortunately become a hostile minefield, obscuring both
140 intent and meaning. Therefore for the
141 purposes of this section the more natural **LSB0 numbering is assumed**
142 and it is left to the reader to translate to MSB0 numbering.
143
144 The Canonical specification for how element-sequential numbering and
145 element-width overrides is defined is expressed in the following c
146 structure, assuming a Little-Endian system, and naturally using LSB0
147 numbering everywhere because the ANSI c specification is inherently LSB0:
148
149 ```
150 #pragma pack
151 typedef union {
152 uint8_t b[]; // elwidth 8
153 uint16_t s[]; // elwidth 16
154 uint32_t i[]; // elwidth 32
155 uint64_t l[]; // elwidth 64
156 uint8_t actual_bytes[8];
157 } el_reg_t;
158
159 elreg_t int_regfile[128];
160
161 void get_register_element(el_reg_t* el, int gpr, int element, int width) {
162 switch (width) {
163 case 64: el->l = int_regfile[gpr].l[element];
164 case 32: el->i = int_regfile[gpr].i[element];
165 case 16: el->s = int_regfile[gpr].s[element];
166 case 8 : el->b = int_regfile[gpr].b[element];
167 }
168 }
169 void set_register_element(el_reg_t* el, int gpr, int element, int width) {
170 switch (width) {
171 case 64: int_regfile[gpr].l[element] = el->l;
172 case 32: int_regfile[gpr].i[element] = el->i;
173 case 16: int_regfile[gpr].s[element] = el->s;
174 case 8 : int_regfile[gpr].b[element] = el->b;
175 }
176 }
177 ```
178
179 Example Vector-looped add operation implementation when elwidths are 64-bit:
180
181 ```
182 # add RT, RA,RB using the "uint64_t" union member, "l"
183 for i in range(VL):
184 int_regfile[RT].l[i] = int_regfile[RA].l[i] + int_regfile[RB].l[i]
185 ```
186
187 However if elwidth overrides are set to 16 for both source and destination:
188
189 ```
190 # add RT, RA, RB using the "uint64_t" union member "s"
191 for i in range(VL):
192 int_regfile[RT].s[i] = int_regfile[RA].s[i] + int_regfile[RB].s[i]
193 ```
194
195 Hardware Architectural note: to avoid a Read-Modify-Write at the register
196 file it is strongly recommended to implement byte-level write-enable lines
197 exactly as has been implemented in DRAM ICs for many decades. Additionally
198 the predicate mask bit is advised to be associated with the element
199 operation and alongside the result ultimately passed to the register file.
200 When element-width is set to 64-bit the relevant predicate mask bit
201 may be repeated eight times and pull all eight write-port byte-level
202 lines HIGH. Clearly when element-width is set to 8-bit the relevant
203 predicate mask bit corresponds directly with one single byte-level
204 write-enable line. It is up to the Hardware Architect to then amortise
205 (merge) elements together into both PredicatedSIMD Pipelines as well
206 as simultaneous non-overlapping Register File writes, to achieve High
207 Performance designs.
208
209 ## Scalar Identity Behaviour
210
211 SVP64 is designed so that when the prefix is all zeros, and
212 VL=1, no effect or
213 influence occurs (no augmentation) such that all standard Power ISA
214 v3.0/v3 1 instructions covered by the prefix are "unaltered". This
215 is termed `scalar identity behaviour` (based on the mathematical
216 definition for "identity", as in, "identity matrix" or better "identity
217 transformation").
218
219 Note that this is completely different from when VL=0. VL=0 turns all
220 operations under its influence into `nops` (regardless of the prefix)
221 whereas when VL=1 and the SV prefix is all zeros, the operation simply
222 acts as if SV had not been applied at all to the instruction (an
223 "identity transformation").
224
225 The fact that `VL` is dynamic and can be set to any value at runtime based
226 on program conditions and behaviour means very specifically that
227 `scalar identity behaviour` is **not** a redundant encoding. If the
228 only means by which VL could be set was by way of static-compiled
229 immediates then this assertion would be false. VL should not
230 be confused with MAXVL when understanding this key aspect of SimpleV.
231
232 ## Register Naming and size
233
234 As indicated above SV Registers are simply the GPR, FPR and CR
235 register files extended linearly to larger sizes; SV Vectorisation
236 iterates sequentially through these registers (LSB0 sequential ordering
237 from 0 to VL-1).
238
239 Where the integer regfile in standard scalar Power ISA v3.0B/v3.1B is
240 r0 to r31, SV extends this as r0 to r127. Likewise FP registers are
241 extended to 128 (fp0 to fp127), and CR Fields are extended to 128 entries,
242 CR0 thru CR127.
243
244 The names of the registers therefore reflects a simple linear extension
245 of the Power ISA v3.0B / v3.1B register naming, and in hardware this
246 would be reflected by a linear increase in the size of the underlying
247 SRAM used for the regfiles.
248
249 Note: when an EXTRA field (defined below) is zero, SV is deliberately
250 designed so that the register fields are identical to as if SV was not in
251 effect i.e. under these circumstances (EXTRA=0) the register field names
252 RA, RB etc. are interpreted and treated as v3.0B / v3.1B scalar registers.
253 This is part of `scalar identity behaviour` described above.
254
255 **Condition Register(s)**
256
257 The Scalar Power ISA Condition Register is a 64 bit register where the top
258 32 MSBs (numbered 0:31 in MSB0 numbering) are not used. This convention is
259 *preserved*
260 in SVP64 and an additional 15 Condition Registers provided in
261 order to store the new CR Fields, CR8-CR15, CR16-CR23 etc. sequentially.
262 The top 32 MSBs in each new SVP64 Condition Register are *also* not used:
263 only the bottom 32 bits (numbered 32:63 in MSB0 numbering).
264
265 *Programmer's note: using `sv.mfcr` without element-width overrides
266 to take into account the fact that the top 32 MSBs are zero and thus
267 effectively doubling the number of GPR registers required to hold all 128
268 CR Fields would seem the only option because normally elwidth overrides
269 would halve the capacity of the instruction. However in this case it
270 is possible to use destination element-width overrides (for `sv.mfcr`.
271 source overrides would be used on the GPR of `sv.mtocrf`), whereupon
272 truncation of the 64-bit Condition Register(s) occurs, throwing away
273 the zeros and storing the remaining (valid, desired) 32-bit values
274 sequentially into (LSB0-convention) lower-numbered and upper-numbered
275 halves of GPRs respectively. The programmer is expected to be aware
276 however that the full width of the entire 64-bit Condition Register
277 is considered to be "an element". This is **not** like any other
278 Condition-Register instructions because all other CR instructions,
279 on closer investigation, will be observed to all be CR-bit or CR-Field
280 related. Thus a `VL` of 16 must be used*
281
282 ## Future expansion.
283
284 With the way that EXTRA fields are defined and applied to register fields,
285 future versions of SV may involve 256 or greater registers. Backwards
286 binary compatibility may be achieved with a PCR bit (Program Compatibility
287 Register) or an MSR bit analogous to SF.
288 Further discussion is out of scope for this version of SVP64.
289
290 Additionally, a future variant of SVP64 will be applied to the Scalar
291 (Quad-precision and 128-bit) VSX instructions. Element-width overrides
292 are an opportunity to expand a future version of the Power ISA
293 to 256-bit, 512-bit and
294 1024-bit operations, as well as doubling or quadrupling the number
295 of VSX registers to 128 or 256. Again further discussion is out of
296 scope for this version of SVP64.
297
298 --------
299
300 \newpage{}
301
302 # New 64-bit Instruction Encoding spaces
303
304 The following seven new areas are defined within Primary Opcode 9 (EXT009)
305 as a new 64-bit encoding space, alongside EXT1xx.
306
307 | 0-5 | 6 | 7 | 8-31 | 32| Description |
308 |-----|---|---|-------|---|------------------------------------|
309 | PO | 0 | x | xxxx | 0 | `RESERVED2` (56-bit) |
310 | PO | 0 | 0 | !zero | 1 | SVP64Single:EXT232-263, or `RESERVED3` |
311 | PO | 0 | 0 | 0000 | 1 | Scalar EXT232-263 |
312 | PO | 0 | 1 | nnnn | 1 | SVP64:EXT232-263 |
313 | PO | 1 | 0 | 0000 | x | `RESERVED1` (32-bit) |
314 | PO | 1 | 0 | !zero | n | SVP64Single:EXT000-063 or `RESERVED4` |
315 | PO | 1 | 1 | nnnn | n | SVP64:EXT000-063 |
316
317 Note that for the future SVP64Single Encoding (currently RESERVED3 and 4)
318 it is prohibited to have bits 8-31 be zero, unlike for SVP64 Vector space,
319 for which bits 8-31 can be zero (termed `scalar identity behaviour`). This
320 prohibition allows SVP64Single to share its Encoding space with Scalar
321 Ext232-263 and Scalar EXT300-363.
322
323 Also that RESERVED1 and 2 are candidates for future Major opcode
324 areas EXT200-231 and EXT300-363 respectively, however as RESERVED areas
325 they may equally be allocated entirely differently.
326
327 *Architectural Resource Allocation Note: **under no circumstances** must
328 different Defined Words be allocated within any `EXT{z}` prefixed
329 or unprefixed space for a given value of `z`. Even if UnVectoriseable
330 an instruction Defined Word space must have the exact same Instruction
331 and exact same Instruction Encoding in all spaces (including
332 being RESERVED if UnVectoriseable) or not be allocated at all.
333 This is required as an inviolate hard rule governing Primary Opcode 9
334 that may not be revoked under any circumstances. A useful way to think
335 of this is that the Prefix Encoding is, like the 8086 REP instruction,
336 an independent 32-bit Defined Word. The only semi-exceptions are
337 the Post-Increment Mode of LD/ST-Update and Vectorised Branch-Conditional.*
338
339 Ecoding spaces and their potential are illustrated:
340
341 | Encoding | Available bits | Scalar | Vectoriseable | SVP64Single |
342 |----------|----------------|--------|---------------|--------------|
343 |EXT000-063| 32 | yes | yes |yes |
344 |EXT100-163| 64 (?) | yes | no |no |
345 |R3SERVED2 | 56 | N/A |not applicable |not applicable|
346 |EXT232-263| 32 | yes | yes |yes |
347 |RESERVED1 | 32 | N/A | no |no |
348
349 Prefixed-Prefixed (96-bit) instructions are prohibited. RESERVED2 presently
350 remains unallocated as of yet and therefore its potential is not yet defined
351 (Not Applicable). RESERVED1 is also unallocated at present, but it is
352 known in advance that the area is UnVectoriseable and also cannot be
353 Prefixed with SVP64Single.
354
355 # Remapped Encoding (`RM[0:23]`)
356
357 In the SVP64 Vector Prefix spaces, the 24 bits 8-31 are termed `RM`. Bits
358 32-37 are the Primary Opcode of the Suffix "Defined Word". 38-63 are the
359 remainder of the Defined Word. Note that the new EXT232-263 SVP64 area
360 it is obviously mandatory that bit 32 is required to be set to 1.
361
362 | 0-5 | 6 | 7 | 8-31 | 32-37 | 38-64 |Description |
363 |-----|---|---|----------|--------|----------|-----------------------|
364 | PO | 0 | 1 | RM[0:23] | 1nnnnn | xxxxxxxx | SVP64:EXT232-263 |
365 | PO | 1 | 1 | RM[0:23] | nnnnnn | xxxxxxxx | SVP64:EXT000-063 |
366
367 It is important to note that unlike v3.1 64-bit prefixed instructions
368 there is insufficient space in `RM` to provide identification of
369 any SVP64 Fields without first partially decoding the 32-bit suffix.
370 Similar to the "Forms" (X-Form, D-Form) the `RM` format is individually
371 associated with every instruction. However this still does not adversely
372 affect Multi-Issue Decoding because the identification of the *length*
373 of anything in the 64-bit space has been kept brutally simple (EXT009),
374 and further decoding of any number of 64-bit Encodings in parallel at
375 that point is fully independent.
376
377 Extreme caution and care must be taken when extending SVP64
378 in future, to not create unnecessary relationships between prefix and
379 suffix that could complicate decoding, adding latency.
380
381 ## Common RM fields
382
383 The following fields are common to all Remapped Encodings:
384
385 | Field Name | Field bits | Description |
386 |------------|------------|----------------------------------------|
387 | MASKMODE | `0` | Execution (predication) Mask Kind |
388 | MASK | `1:3` | Execution Mask |
389 | SUBVL | `8:9` | Sub-vector length |
390
391 The following fields are optional or encoded differently depending
392 on context after decoding of the Scalar suffix:
393
394 | Field Name | Field bits | Description |
395 |------------|------------|----------------------------------------|
396 | ELWIDTH | `4:5` | Element Width |
397 | ELWIDTH_SRC | `6:7` | Element Width for Source |
398 | EXTRA | `10:18` | Register Extra encoding |
399 | MODE | `19:23` | changes Vector behaviour |
400
401 * MODE changes the behaviour of the SV operation (result saturation,
402 mapreduce)
403 * SUBVL groups elements together into vec2, vec3, vec4 for use in 3D
404 and Audio/Video DSP work
405 * ELWIDTH and ELWIDTH_SRC overrides the instruction's destination and
406 source operand width
407 * MASK (and MASK_SRC) and MASKMODE provide predication (two types of
408 sources: scalar INT and Vector CR).
409 * Bits 10 to 18 (EXTRA) are further decoded depending on the RM category
410 for the instruction, which is determined only by decoding the Scalar 32
411 bit suffix.
412
413 Similar to Power ISA `X-Form` etc. EXTRA bits are given designations,
414 such as `RM-1P-3S1D` which indicates for this example that the operation
415 is to be single-predicated and that there are 3 source operand EXTRA
416 tags and one destination operand tag.
417
418 Note that if ELWIDTH != ELWIDTH_SRC this may result in reduced performance
419 or increased latency in some implementations due to lane-crossing.
420
421 ## Mode
422
423 Mode is an augmentation of SV behaviour. Different types of instructions
424 have different needs, similar to Power ISA v3.1 64 bit prefix 8LS and MTRR
425 formats apply to different instruction types. Modes include Reduction,
426 Iteration, arithmetic saturation, and Fail-First. More specific details
427 in each section and in the SVP64 appendix
428
429 * For condition register operations see [[sv/cr_ops]]
430 * For LD/ST Modes, see [[sv/ldst]].
431 * For Branch modes, see [[sv/branches]]
432 * For arithmetic and logical, see [[sv/normal]]
433
434 ## ELWIDTH Encoding
435
436 Default behaviour is set to 0b00 so that zeros follow the convention
437 of `scalar identity behaviour`. In this case it means that elwidth
438 overrides are not applicable. Thus if a 32 bit instruction operates
439 on 32 bit, `elwidth=0b00` specifies that this behaviour is unmodified.
440 Likewise when a processor is switched from 64 bit to 32 bit mode,
441 `elwidth=0b00` states that, again, the behaviour is not to be modified.
442
443 Only when elwidth is nonzero is the element width overridden to the
444 explicitly required value.
445
446 ### Elwidth for Integers:
447
448 | Value | Mnemonic | Description |
449 |-------|----------------|------------------------------------|
450 | 00 | DEFAULT | default behaviour for operation |
451 | 01 | `ELWIDTH=w` | Word: 32-bit integer |
452 | 10 | `ELWIDTH=h` | Halfword: 16-bit integer |
453 | 11 | `ELWIDTH=b` | Byte: 8-bit integer |
454
455 This encoding is chosen such that the byte width may be computed as
456 `8<<(3-ew)`
457
458 ### Elwidth for FP Registers:
459
460 | Value | Mnemonic | Description |
461 |-------|----------------|------------------------------------|
462 | 00 | DEFAULT | default behaviour for FP operation |
463 | 01 | `ELWIDTH=f32` | 32-bit IEEE 754 Single floating-point |
464 | 10 | `ELWIDTH=f16` | 16-bit IEEE 754 Half floating-point |
465 | 11 | `ELWIDTH=bf16` | Reserved for `bf16` |
466
467 Note:
468 [`bf16`](https://en.wikipedia.org/wiki/Bfloat16_floating-point_format)
469 is reserved for a future implementation of SV
470
471 Note that any IEEE754 FP operation in Power ISA ending in "s" (`fadds`)
472 shall perform its operation at **half** the ELWIDTH then padded back out
473 to ELWIDTH. `sv.fadds/ew=f32` shall perform an IEEE754 FP16 operation
474 that is then "padded" to fill out to an IEEE754 FP32. When ELWIDTH=DEFAULT
475 clearly the behaviour of `sv.fadds` is performed at 32-bit accuracy
476 then padded back out to fit in IEEE754 FP64, exactly as for Scalar
477 v3.0B "single" FP. Any FP operation ending in "s" where ELWIDTH=f16 or
478 ELWIDTH=bf16 is reserved and must raise an illegal instruction (IEEE754
479 FP8 or BF8 are not defined).
480
481 ### Elwidth for CRs (no meaning)
482
483 Element-width overrides for CR Fields has no meaning. The bits
484 are therefore used for other purposes, or when Rc=1, the Elwidth
485 applies to the result being tested (a GPR or FPR), but not to the
486 Vector of CR Fields.
487
488 ## SUBVL Encoding
489
490 The default for SUBVL is 1 and its encoding is 0b00 to indicate that
491 SUBVL is effectively disabled (a SUBVL for-loop of only one element). this
492 lines up in combination with all other "default is all zeros" behaviour.
493
494 | Value | Mnemonic | Subvec | Description |
495 |-------|-----------|---------|------------------------|
496 | 00 | `SUBVL=1` | single | Sub-vector length of 1 |
497 | 01 | `SUBVL=2` | vec2 | Sub-vector length of 2 |
498 | 10 | `SUBVL=3` | vec3 | Sub-vector length of 3 |
499 | 11 | `SUBVL=4` | vec4 | Sub-vector length of 4 |
500
501 The SUBVL encoding value may be thought of as an inclusive range of a
502 sub-vector. SUBVL=2 represents a vec2, its encoding is 0b01, therefore
503 this may be considered to be elements 0b00 to 0b01 inclusive.
504
505 ## MASK/MASK_SRC & MASKMODE Encoding
506
507 One bit (`MASKMODE`) indicates the mode: CR or Int predication. The two
508 types may not be mixed.
509
510 Special note: to disable predication this field must be set to zero in
511 combination with Integer Predication also being set to 0b000. this has the
512 effect of enabling "all 1s" in the predicate mask, which is equivalent to
513 "not having any predication at all" and consequently, in combination with
514 all other default zeros, fully disables SV (`scalar identity behaviour`).
515
516 `MASKMODE` may be set to one of 2 values:
517
518 | Value | Description |
519 |-----------|------------------------------------------------------|
520 | 0 | MASK/MASK_SRC are encoded using Integer Predication |
521 | 1 | MASK/MASK_SRC are encoded using CR-based Predication |
522
523 Integer Twin predication has a second set of 3 bits that uses the same
524 encoding thus allowing either the same register (r3, r10 or r31) to be
525 used for both src and dest, or different regs (one for src, one for dest).
526
527 Likewise CR based twin predication has a second set of 3 bits, allowing
528 a different test to be applied.
529
530 Note that it is assumed that Predicate Masks (whether INT or CR) are
531 read *before* the operations proceed. In practice (for CR Fields)
532 this creates an unnecessary block on parallelism. Therefore, it is up
533 to the programmer to ensure that the CR fields used as Predicate Masks
534 are not being written to by any parallel Vector Loop. Doing so results
535 in **UNDEFINED** behaviour, according to the definition outlined in the
536 Power ISA v3.0B Specification.
537
538 Hardware Implementations are therefore free and clear to delay reading
539 of individual CR fields until the actual predicated element operation
540 needs to take place, safe in the knowledge that no programmer will have
541 issued a Vector Instruction where previous elements could have overwritten
542 (destroyed) not-yet-executed CR-Predicated element operations.
543
544 ### Integer Predication (MASKMODE=0)
545
546 When the predicate mode bit is zero the 3 bits are interpreted as below.
547 Twin predication has an identical 3 bit field similarly encoded.
548
549 `MASK` and `MASK_SRC` may be set to one of 8 values, to provide the
550 following meaning:
551
552 | Value | Mnemonic | Element `i` enabled if: |
553 |-------|----------|------------------------------|
554 | 000 | ALWAYS | predicate effectively all 1s |
555 | 001 | 1 << R3 | `i == R3` |
556 | 010 | R3 | `R3 & (1 << i)` is non-zero |
557 | 011 | ~R3 | `R3 & (1 << i)` is zero |
558 | 100 | R10 | `R10 & (1 << i)` is non-zero |
559 | 101 | ~R10 | `R10 & (1 << i)` is zero |
560 | 110 | R30 | `R30 & (1 << i)` is non-zero |
561 | 111 | ~R30 | `R30 & (1 << i)` is zero |
562
563 r10 and r30 are at the high end of temporary and unused registers,
564 so as not to interfere with register allocation from ABIs.
565
566 ### CR-based Predication (MASKMODE=1)
567
568 When the predicate mode bit is one the 3 bits are interpreted as below.
569 Twin predication has an identical 3 bit field similarly encoded.
570
571 `MASK` and `MASK_SRC` may be set to one of 8 values, to provide the
572 following meaning:
573
574 | Value | Mnemonic | Element `i` is enabled if |
575 |-------|----------|--------------------------|
576 | 000 | lt | `CR[offs+i].LT` is set |
577 | 001 | nl/ge | `CR[offs+i].LT` is clear |
578 | 010 | gt | `CR[offs+i].GT` is set |
579 | 011 | ng/le | `CR[offs+i].GT` is clear |
580 | 100 | eq | `CR[offs+i].EQ` is set |
581 | 101 | ne | `CR[offs+i].EQ` is clear |
582 | 110 | so/un | `CR[offs+i].FU` is set |
583 | 111 | ns/nu | `CR[offs+i].FU` is clear |
584
585 `offs` is defined as CR32 (4x8) so as to mesh cleanly with Vectorised
586 Rc=1 operations (see below). Rc=1 operations start from CR8 (TBD).
587
588 The CR Predicates chosen must start on a boundary that Vectorised CR
589 operations can access cleanly, in full. With EXTRA2 restricting starting
590 points to multiples of 8 (CR0, CR8, CR16...) both Vectorised Rc=1 and
591 CR Predicate Masks have to be adapted to fit on these boundaries as well.
592
593 ## Extra Remapped Encoding <a name="extra_remap"> </a>
594
595 Shows all instruction-specific fields in the Remapped Encoding
596 `RM[10:18]` for all instruction variants. Note that due to the very
597 tight space, the encoding mode is *not* included in the prefix itself.
598 The mode is "applied", similar to Power ISA "Forms" (X-Form, D-Form)
599 on a per-instruction basis, and, like "Forms" are given a designation
600 (below) of the form `RM-nP-nSnD`. The full list of which instructions
601 use which remaps is here [[opcode_regs_deduped]].
602
603 **Please note the following**:
604
605 ```
606 Machine-readable CSV files have been provided which will make the task
607 of creating SV-aware ISA decoders, documentation, assembler tools
608 compiler tools Simulators documentation all aspects of SVP64 easier
609 and less prone to mistakes. Please avoid manual re-creation of
610 information from the written specification wording, and use the
611 CSV files or use the Canonical tool which creates the CSV files,
612 named sv_analysis.py. The information contained within sv_analysis.py
613 is considered to be part of this Specification, even encoded as it
614 is in python3.
615 ```
616
617 The mappings are part of the SVP64 Specification in exactly the same
618 way as X-Form, D-Form. New Scalar instructions added to the Power ISA
619 will need a corresponding SVP64 Mapping, which can be derived by-rote
620 from examining the Register "Profile" of the instruction.
621
622 There are two categories: Single and Twin Predication. Due to space
623 considerations further subdivision of Single Predication is based on
624 whether the number of src operands is 2 or 3. With only 9 bits available
625 some compromises have to be made.
626
627 * `RM-1P-3S1D` Single Predication dest/src1/2/3, applies to 4-operand
628 instructions (fmadd, isel, madd).
629 * `RM-1P-2S1D` Single Predication dest/src1/2 applies to 3-operand
630 instructions (src1 src2 dest)
631 * `RM-2P-1S1D` Twin Predication (src=1, dest=1)
632 * `RM-2P-2S1D` Twin Predication (src=2, dest=1) primarily for LDST (Indexed)
633 * `RM-2P-1S2D` Twin Predication (src=1, dest=2) primarily for LDST Update
634
635 ### RM-1P-3S1D
636
637 | Field Name | Field bits | Description |
638 |------------|------------|----------------------------------------|
639 | Rdest\_EXTRA2 | `10:11` | extends Rdest (R\*\_EXTRA2 Encoding) |
640 | Rsrc1\_EXTRA2 | `12:13` | extends Rsrc1 (R\*\_EXTRA2 Encoding) |
641 | Rsrc2\_EXTRA2 | `14:15` | extends Rsrc2 (R\*\_EXTRA2 Encoding) |
642 | Rsrc3\_EXTRA2 | `16:17` | extends Rsrc3 (R\*\_EXTRA2 Encoding) |
643 | EXTRA2_MODE | `18` | used by `divmod2du` and `maddedu` for RS |
644
645 These are for 3 operand in and either 1 or 2 out instructions.
646 3-in 1-out includes `madd RT,RA,RB,RC`. (DRAFT) instructions
647 such as `maddedu` have an implicit second destination, RS, the
648 selection of which is determined by bit 18.
649
650 ### RM-1P-2S1D
651
652 | Field Name | Field bits | Description |
653 |------------|------------|-------------------------------------------|
654 | Rdest\_EXTRA3 | `10:12` | extends Rdest |
655 | Rsrc1\_EXTRA3 | `13:15` | extends Rsrc1 |
656 | Rsrc2\_EXTRA3 | `16:18` | extends Rsrc3 |
657
658 These are for 2 operand 1 dest instructions, such as `add RT, RA,
659 RB`. However also included are unusual instructions with an implicit
660 dest that is identical to its src reg, such as `rlwinmi`.
661
662 Normally, with instructions such as `rlwinmi`, the scalar v3.0B ISA would
663 not have sufficient bit fields to allow an alternative destination.
664 With SV however this becomes possible. Therefore, the fact that the
665 dest is implicitly also a src should not mislead: due to the *prefix*
666 they are different SV regs.
667
668 * `rlwimi RA, RS, ...`
669 * Rsrc1_EXTRA3 applies to RS as the first src
670 * Rsrc2_EXTRA3 applies to RA as the secomd src
671 * Rdest_EXTRA3 applies to RA to create an **independent** dest.
672
673 With the addition of the EXTRA bits, the three registers
674 each may be *independently* made vector or scalar, and be independently
675 augmented to 7 bits in length.
676
677 ### RM-2P-1S1D/2S
678
679 | Field Name | Field bits | Description |
680 |------------|------------|----------------------------|
681 | Rdest_EXTRA3 | `10:12` | extends Rdest |
682 | Rsrc1_EXTRA3 | `13:15` | extends Rsrc1 |
683 | MASK_SRC | `16:18` | Execution Mask for Source |
684
685 `RM-2P-2S` is for `stw` etc. and is Rsrc1 Rsrc2.
686
687 ### RM-1P-2S1D
688
689 single-predicate, three registers (2 read, 1 write)
690
691 | Field Name | Field bits | Description |
692 |------------|------------|----------------------------|
693 | Rdest_EXTRA3 | `10:12` | extends Rdest |
694 | Rsrc1_EXTRA3 | `13:15` | extends Rsrc1 |
695 | Rsrc2_EXTRA3 | `16:18` | extends Rsrc2 |
696
697 ### RM-2P-2S1D/1S2D/3S
698
699 The primary purpose for this encoding is for Twin Predication on LOAD
700 and STORE operations. see [[sv/ldst]] for detailed anslysis.
701
702 **RM-2P-2S1D:**
703
704 | Field Name | Field bits | Description |
705 |------------|------------|----------------------------|
706 | Rdest_EXTRA2 | `10:11` | extends Rdest (R\*\_EXTRA2 Encoding) |
707 | Rsrc1_EXTRA2 | `12:13` | extends Rsrc1 (R\*\_EXTRA2 Encoding) |
708 | Rsrc2_EXTRA2 | `14:15` | extends Rsrc2 (R\*\_EXTRA2 Encoding) |
709 | MASK_SRC | `16:18` | Execution Mask for Source |
710
711 **RM-2P-1S2D:**
712
713 For RM-2P-1S2D the EXTRA2 dest and src names are switched (Rsrc_EXTRA2
714 is in bits 10:11, Rdest1_EXTRA2 in 12:13)
715
716 | Field Name | Field bits | Description |
717 |------------|------------|----------------------------|
718 | Rsrc2_EXTRA2 | `10:11` | extends Rsrc2 (R\*\_EXTRA2 Encoding) |
719 | Rsrc1_EXTRA2 | `12:13` | extends Rsrc1 (R\*\_EXTRA2 Encoding) |
720 | Rdest_EXTRA2 | `14:15` | extends Rdest (R\*\_EXTRA2 Encoding) |
721 | MASK_SRC | `16:18` | Execution Mask for Source |
722
723 **RM-2P-3S:**
724
725 Also that for RM-2P-3S (to cover `stdx` etc.) the names are switched to 3 src:
726 Rsrc1_EXTRA2, Rsrc2_EXTRA2, Rsrc3_EXTRA2.
727
728 | Field Name | Field bits | Description |
729 |------------|------------|----------------------------|
730 | Rsrc1_EXTRA2 | `10:11` | extends Rsrc1 (R\*\_EXTRA2 Encoding) |
731 | Rsrc2_EXTRA2 | `12:13` | extends Rsrc2 (R\*\_EXTRA2 Encoding) |
732 | Rsrc3_EXTRA2 | `14:15` | extends Rsrc3 (R\*\_EXTRA2 Encoding) |
733 | MASK_SRC | `16:18` | Execution Mask for Source |
734
735 Note also that LD with update indexed, which takes 2 src and
736 creates 2 dest registers (e.g. `lhaux RT,RA,RB`), does not have room
737 for 4 registers and also Twin Predication. Therefore these are treated as
738 RM-2P-2S1D and the src spec for RA is also used for the same RA as a dest.
739
740 Note that if ELWIDTH != ELWIDTH_SRC this may result in reduced performance
741 or increased latency in some implementations due to lane-crossing.
742
743 ## R\*\_EXTRA2/3
744
745 EXTRA is the means by which two things are achieved:
746
747 1. Registers are marked as either Vector *or Scalar*
748 2. Register field numbers (limited typically to 5 bit)
749 are extended in range, both for Scalar and Vector.
750
751 The register files are therefore extended:
752
753 * INT (GPR) is extended from r0-31 to r0-127
754 * FP (FPR) is extended from fp0-32 to fp0-fp127
755 * CR Fields are extended from CR0-7 to CR0-127
756
757 However due to pressure in `RM.EXTRA` not all these registers
758 are accessible by all instructions, particularly those with
759 a large number of operands (`madd`, `isel`).
760
761 In the following tables register numbers are constructed from the
762 standard v3.0B / v3.1B 32 bit register field (RA, FRA) and the EXTRA2 or
763 EXTRA3 field from the SV Prefix, determined by the specific RM-xx-yyyy
764 designation for a given instruction. The prefixing is arranged so that
765 interoperability between prefixing and nonprefixing of scalar registers
766 is direct and convenient (when the EXTRA field is all zeros).
767
768 A pseudocode algorithm explains the relationship, for INT/FP (see
769 SVP64 appendix for CRs)
770
771 ```
772 if extra3_mode:
773 spec = EXTRA3
774 else:
775 spec = EXTRA2 << 1 # same as EXTRA3, shifted
776 if spec[0]: # vector
777 return (RA << 2) | spec[1:2]
778 else: # scalar
779 return (spec[1:2] << 5) | RA
780 ```
781
782 Future versions may extend to 256 by shifting Vector numbering up.
783 Scalar will not be altered.
784
785 Note that in some cases the range of starting points for Vectors
786 is limited.
787
788 ### INT/FP EXTRA3
789
790 If EXTRA3 is zero, maps to "scalar identity" (scalar Power ISA field
791 naming).
792
793 Fields are as follows:
794
795 * Value: R_EXTRA3
796 * Mode: register is tagged as scalar or vector
797 * Range/Inc: the range of registers accessible from this EXTRA
798 encoding, and the "increment" (accessibility). "/4" means
799 that this EXTRA encoding may only give access (starting point)
800 every 4th register.
801 * MSB..LSB: the bit field showing how the register opcode field
802 combines with EXTRA to give (extend) the register number (GPR)
803
804 | Value | Mode | Range/Inc | 6..0 |
805 |-----------|-------|---------------|---------------------|
806 | 000 | Scalar | `r0-r31`/1 | `0b00 RA` |
807 | 001 | Scalar | `r32-r63`/1 | `0b01 RA` |
808 | 010 | Scalar | `r64-r95`/1 | `0b10 RA` |
809 | 011 | Scalar | `r96-r127`/1 | `0b11 RA` |
810 | 100 | Vector | `r0-r124`/4 | `RA 0b00` |
811 | 101 | Vector | `r1-r125`/4 | `RA 0b01` |
812 | 110 | Vector | `r2-r126`/4 | `RA 0b10` |
813 | 111 | Vector | `r3-r127`/4 | `RA 0b11` |
814
815 ### INT/FP EXTRA2
816
817 If EXTRA2 is zero will map to
818 "scalar identity behaviour" i.e Scalar Power ISA register naming:
819
820 | Value | Mode | Range/inc | 6..0 |
821 |----------|-------|---------------|-----------|
822 | 00 | Scalar | `r0-r31`/1 | `0b00 RA` |
823 | 01 | Scalar | `r32-r63`/1 | `0b01 RA` |
824 | 10 | Vector | `r0-r124`/4 | `RA 0b00` |
825 | 11 | Vector | `r2-r126`/4 | `RA 0b10` |
826
827 **Note that unlike in EXTRA3, in EXTRA2**:
828
829 * the GPR Vectors may only start from
830 `r0, r2, r4, r6, r8` and likewise FPR Vectors.
831 * the GPR Scalars may only go from `r0, r1, r2.. r63` and likewise FPR Scalars.
832
833 as there is insufficient bits to cover the full range.
834
835 ### CR Field EXTRA3
836
837 CR Field encoding is essentially the same but made more complex due to CRs
838 being bit-based, because the application of SVP64 element-numbering applies
839 to the CR *Field* numbering not the CR register *bit* numbering.
840 Note that Vectors may only start from `CR0, CR4, CR8, CR12, CR16, CR20`...
841 and Scalars may only go from `CR0, CR1, ... CR31`
842
843 Encoding shown MSB down to LSB
844
845 For a 5-bit operand (BA, BB, BT):
846
847 | Value | Mode | Range/Inc | 8..5 | 4..2 | 1..0 |
848 |-------|------|---------------|-----------| --------|---------|
849 | 000 | Scalar | `CR0-CR7`/1 | 0b0000 | BA[0:2] | BA[3:4] |
850 | 001 | Scalar | `CR8-CR15`/1 | 0b0001 | BA[0:2] | BA[3:4] |
851 | 010 | Scalar | `CR16-CR23`/1 | 0b0010 | BA[0:2] | BA[3:4] |
852 | 011 | Scalar | `CR24-CR31`/1 | 0b0011 | BA[0:2] | BA[3:4] |
853 | 100 | Vector | `CR0-CR112`/16 | BA[0:2] 0 | 0b000 | BA[3:4] |
854 | 101 | Vector | `CR4-CR116`/16 | BA[0:2] 0 | 0b100 | BA[3:4] |
855 | 110 | Vector | `CR8-CR120`/16 | BA[0:2] 1 | 0b000 | BA[3:4] |
856 | 111 | Vector | `CR12-CR124`/16 | BA[0:2] 1 | 0b100 | BA[3:4] |
857
858 For a 3-bit operand (e.g. BFA):
859
860 | Value | Mode | Range/Inc | 6..3 | 2..0 |
861 |-------|------|---------------|-----------| --------|
862 | 000 | Scalar | `CR0-CR7`/1 | 0b0000 | BFA |
863 | 001 | Scalar | `CR8-CR15`/1 | 0b0001 | BFA |
864 | 010 | Scalar | `CR16-CR23`/1 | 0b0010 | BFA |
865 | 011 | Scalar | `CR24-CR31`/1 | 0b0011 | BFA |
866 | 100 | Vector | `CR0-CR112`/16 | BFA 0 | 0b000 |
867 | 101 | Vector | `CR4-CR116`/16 | BFA 0 | 0b100 |
868 | 110 | Vector | `CR8-CR120`/16 | BFA 1 | 0b000 |
869 | 111 | Vector | `CR12-CR124`/16 | BFA 1 | 0b100 |
870
871 ### CR EXTRA2
872
873 CR encoding is essentially the same but made more complex due to CRs
874 being bit-based, because the application of SVP64 element-numbering applies
875 to the CR *Field* numbering not the CR register *bit* numbering.
876 See separate section for explanation and pseudocode.
877 Note that Vectors may only start from CR0, CR8, CR16, CR24, CR32...
878
879 Encoding shown MSB down to LSB
880
881 For a 5-bit operand (BA, BB, BC):
882
883 | Value | Mode | Range/Inc | 8..5 | 4..2 | 1..0 |
884 |-------|--------|----------------|---------|---------|---------|
885 | 00 | Scalar | `CR0-CR7`/1 | 0b0000 | BA[0:2] | BA[3:4] |
886 | 01 | Scalar | `CR8-CR15`/1 | 0b0001 | BA[0:2] | BA[3:4] |
887 | 10 | Vector | `CR0-CR112`/16 | BA[0:2] 0 | 0b000 | BA[3:4] |
888 | 11 | Vector | `CR8-CR120`/16 | BA[0:2] 1 | 0b000 | BA[3:4] |
889
890 For a 3-bit operand (e.g. BFA):
891
892 | Value | Mode | Range/Inc | 6..3 | 2..0 |
893 |-------|------|---------------|-----------| --------|
894 | 00 | Scalar | `CR0-CR7`/1 | 0b0000 | BFA |
895 | 01 | Scalar | `CR8-CR15`/1 | 0b0001 | BFA |
896 | 10 | Vector | `CR0-CR112`/16 | BFA 0 | 0b000 |
897 | 11 | Vector | `CR8-CR120`/16 | BFA 1 | 0b000 |
898
899 --------
900
901 \newpage{}
902
903
904 # Normal SVP64 Modes, for Arithmetic and Logical Operations
905
906 Normal SVP64 Mode covers Arithmetic and Logical operations
907 to provide suitable additional behaviour. The Mode
908 field is bits 19-23 of the [[svp64]] RM Field.
909
910 ## Mode
911
912 Mode is an augmentation of SV behaviour, providing additional
913 functionality. Some of these alterations are element-based (saturation),
914 others involve post-analysis (predicate result) and others are
915 Vector-based (mapreduce, fail-on-first).
916
917 [[sv/ldst]], [[sv/cr_ops]] and [[sv/branches]] are covered separately:
918 the following Modes apply to Arithmetic and Logical SVP64 operations:
919
920 * **simple** mode is straight vectorisation. no augmentations: the
921 vector comprises an array of independently created results.
922 * **ffirst** or data-dependent fail-on-first: see separate section.
923 the vector may be truncated depending on certain criteria.
924 *VL is altered as a result*.
925 * **sat mode** or saturation: clamps each element result to a min/max
926 rather than overflows / wraps. allows signed and unsigned clamping
927 for both INT and FP.
928 * **reduce mode**. if used correctly, a mapreduce (or a prefix sum)
929 is performed. see [[svp64/appendix]].
930 note that there are comprehensive caveats when using this mode.
931 * **pred-result** will test the result (CR testing selects a bit of CR
932 and inverts it, just like branch conditional testing) and if the
933 test fails it is as if the *destination* predicate bit was zero even
934 before starting the operation. When Rc=1 the CR element however is
935 still stored in the CR regfile, even if the test failed. See appendix
936 for details.
937
938 Note that ffirst and reduce modes are not anticipated to be
939 high-performance in some implementations. ffirst due to interactions
940 with VL, and reduce due to it requiring additional operations to produce
941 a result. simple, saturate and pred-result are however inter-element
942 independent and may easily be parallelised to give high performance,
943 regardless of the value of VL.
944
945 The Mode table for Arithmetic and Logical operations is laid out as
946 follows:
947
948 | 0-1 | 2 | 3 4 | description |
949 | --- | --- |---------|-------------------------- |
950 | 00 | 0 | dz sz | simple mode |
951 | 00 | 1 | 0 RG | scalar reduce mode (mapreduce) |
952 | 00 | 1 | 1 / | reserved |
953 | 01 | inv | CR-bit | Rc=1: ffirst CR sel |
954 | 01 | inv | VLi RC1 | Rc=0: ffirst z/nonz |
955 | 10 | N | dz sz | sat mode: N=0/1 u/s |
956 | 11 | inv | CR-bit | Rc=1: pred-result CR sel |
957 | 11 | inv | zz RC1 | Rc=0: pred-result z/nonz |
958
959 Fields:
960
961 * **sz / dz** if predication is enabled will put zeros into the dest
962 (or as src in the case of twin pred) when the predicate bit is zero.
963 otherwise the element is ignored or skipped, depending on context.
964 * **zz**: both sz and dz are set equal to this flag
965 * **inv CR bit** just as in branches (BO) these bits allow testing of
966 a CR bit and whether it is set (inv=0) or unset (inv=1)
967 * **RG** inverts the Vector Loop order (VL-1 downto 0) rather
968 than the normal 0..VL-1
969 * **N** sets signed/unsigned saturation.
970 * **RC1** as if Rc=1, enables access to `VLi`.
971 * **VLi** VL inclusive: in fail-first mode, the truncation of
972 VL *includes* the current element at the failure point rather
973 than excludes it from the count.
974
975 For LD/ST Modes, see [[sv/ldst]]. For Condition Registers see
976 [[sv/cr_ops]]. For Branch modes, see [[sv/branches]].
977
978 ## Rounding, clamp and saturate
979
980 To help ensure for example that audio quality is not compromised by
981 overflow, "saturation" is provided, as well as a way to detect when
982 saturation occurred if desired (Rc=1). When Rc=1 there will be a *vector*
983 of CRs, one CR per element in the result (Note: this is different from
984 VSX which has a single CR per block).
985
986 When N=0 the result is saturated to within the maximum range of an
987 unsigned value. For integer ops this will be 0 to 2^elwidth-1. Similar
988 logic applies to FP operations, with the result being saturated to
989 maximum rather than returning INF, and the minimum to +0.0
990
991 When N=1 the same occurs except that the result is saturated to the min
992 or max of a signed result, and for FP to the min and max value rather
993 than returning +/- INF.
994
995 When Rc=1, the CR "overflow" bit is set on the CR associated with
996 the element, to indicate whether saturation occurred. Note that
997 due to the hugely detrimental effect it has on parallel processing,
998 XER.SO is **ignored** completely and is **not** brought into play here.
999 The CR overflow bit is therefore simply set to zero if saturation did
1000 not occur, and to one if it did. This behaviour (ignoring XER.SO) is
1001 actually optional in the SFFS Compliancy Subset: for SVP64 it is made
1002 mandatory *but only on Vectorised instructions*.
1003
1004 Note also that saturate on operations that set OE=1 must raise an Illegal
1005 Instruction due to the conflicting use of the CR.so bit for storing
1006 if saturation occurred. Vectorised Integer Operations that produce a
1007 Carry-Out (CA, CA32): these two bits will be `UNDEFINED` if saturation
1008 is also requested.
1009
1010 Note that the operation takes place at the maximum bitwidth (max of
1011 src and dest elwidth) and that truncation occurs to the range of the
1012 dest elwidth.
1013
1014 *Programmer's Note: Post-analysis of the Vector of CRs to find out if any
1015 given element hit saturation may be done using a mapreduced CR op (cror),
1016 or by using the new crrweird instruction with Rc=1, which will transfer
1017 the required CR bits to a scalar integer and update CR0, which will allow
1018 testing the scalar integer for nonzero. see [[sv/cr_int_predication]].
1019 Alternatively, a Data-Dependent Fail-First may be used to truncate the
1020 Vector Length to non-saturated elements, greatly increasing the productivity
1021 of parallelised inner hot-loops.*
1022
1023 ## Reduce mode
1024
1025 Reduction in SVP64 is similar in essence to other Vector Processing ISAs,
1026 but leverages the underlying scalar Base v3.0B operations. Thus it is
1027 more a convention that the programmer may utilise to give the appearance
1028 and effect of a Horizontal Vector Reduction. Due to the unusual decoupling
1029 it is also possible to perform prefix-sum (Fibonacci Series) in certain
1030 circumstances. Details are in the SVP64 appendix
1031
1032 Reduce Mode should not be confused with Parallel Reduction [[sv/remap]].
1033 As explained in the [[sv/appendix]] Reduce Mode switches off the check
1034 which would normally stop looping if the result register is scalar.
1035 Thus, the result scalar register, if also used as a source scalar,
1036 may be used to perform sequential accumulation. This *deliberately*
1037 sets up a chain of Register Hazard Dependencies, whereas Parallel Reduce
1038 [[sv/remap]] deliberately issues a Tree-Schedule of operations that may
1039 be parallelised.
1040
1041 ## Data-dependent Fail-on-first
1042
1043 Data-dependent fail-on-first is very different from LD/ST Fail-First
1044 (also known as Fault-First) and is actually CR-field-driven.
1045 Vector elements are required to appear
1046 to be executed in sequential Program Order. When REMAP is not active,
1047 element 0 would be the first.
1048
1049 Data-driven (CR-driven) fail-on-first activates when Rc=1 or other
1050 CR-creating operation produces a result (including cmp). Similar to
1051 branch, an analysis of the CR is performed and if the test fails, the
1052 vector operation terminates and discards all element operations **at and
1053 above the current one**, and VL is truncated to either the *previous*
1054 element or the current one, depending on whether VLi (VL "inclusive")
1055 is clear or set, respectively.
1056
1057 Thus the new VL comprises a contiguous vector of results, all of which
1058 pass the testing criteria (equal to zero, less than zero etc as defined
1059 by the CR-bit test).
1060
1061 *Note: when VLi is clear, the behaviour at first seems counter-intuitive.
1062 A result is calculated but if the test fails it is prohibited from being
1063 actually written. This becomes intuitive again when it is remembered
1064 that the length that VL is set to is the number of *written* elements, and
1065 only when VLI is set will the current element be included in that count.*
1066
1067 The CR-based data-driven fail-on-first is "new" and not found in ARM SVE
1068 or RVV. At the same time it is "old" because it is almost identical to
1069 a generalised form of Z80's `CPIR` instruction. It is extremely useful
1070 for reducing instruction count, however requires speculative execution
1071 involving modifications of VL to get high performance implementations.
1072 An additional mode (RC1=1) effectively turns what would otherwise be an
1073 arithmetic operation into a type of `cmp`. The CR is stored (and the
1074 CR.eq bit tested against the `inv` field). If the CR.eq bit is equal to
1075 `inv` then the Vector is truncated and the loop ends.
1076
1077 VLi is only available as an option when `Rc=0` (or for instructions
1078 which do not have Rc). When set, the current element is always also
1079 included in the count (the new length that VL will be set to). This may
1080 be useful in combination with "inv" to truncate the Vector to *exclude*
1081 elements that fail a test, or, in the case of implementations of strncpy,
1082 to include the terminating zero.
1083
1084 In CR-based data-driven fail-on-first there is only the option to select
1085 and test one bit of each CR (just as with branch BO). For more complex
1086 tests this may be insufficient. If that is the case, a vectorised crop
1087 such as crand, cror or [[sv/cr_int_predication]] crweirder may be used,
1088 and ffirst applied to the crop instead of to the arithmetic vector. Note
1089 that crops are covered by the [[sv/cr_ops]] Mode format.
1090
1091 Use of Fail-on-first with Vertical-First Mode is not prohibited but is
1092 not really recommended. The effect of truncating VL
1093 may have unintended and unexpected consequences on subsequent instructions.
1094 VLi set will be fine: it is when VLi is clear that problems may be faced.
1095
1096 *Programmer's note: `VLi` is only accessible in normal operations which in
1097 turn limits the CR field bit-testing to only `EQ/NE`. [[sv/cr_ops]] are
1098 not so limited. Thus it is possible to use for example `sv.cror/ff=gt/vli
1099 *0,*0,*0`, which is not a `nop` because it allows Fail-First Mode to
1100 perform a test and truncate VL.*
1101
1102 *Hardware implementor's note: effective Sequential Program Order must
1103 be preserved. Speculative Execution is perfectly permitted as long as
1104 the speculative elements are held back from writing to register files
1105 (kept in Resevation Stations), until such time as the relevant CR Field
1106 bit(s) has been analysed. All Speculative elements sequentially beyond
1107 the test-failure point **MUST** be cancelled. This is no different from
1108 standard Out-of-Order Execution and the modification effort to efficiently
1109 support Data-Dependent Fail-First within a pre-existing Multi-Issue
1110 Out-of-Order Engine is anticipated to be minimal. In-Order systems on
1111 the other hand are expected, unavoidably, to be low-performance*.
1112
1113 Two extremely important aspects of ffirst are:
1114
1115 * LDST ffirst may never set VL equal to zero. This because on the first
1116 element an exception must be raised "as normal".
1117 * CR-based data-dependent ffirst on the other hand **can** set VL equal
1118 to zero. This is the only means in the entirety of SV that VL may be set
1119 to zero (with the exception of via the SV.STATE SPR). When VL is set
1120 zero due to the first element failing the CR bit-test, all subsequent
1121 vectorised operations are effectively `nops` which is
1122 *precisely the desired and intended behaviour*.
1123
1124 The second crucial aspect, compared to LDST Ffirst:
1125
1126 * LD/ST Failfirst may (beyond the initial first element
1127 conditions) truncate VL for any architecturally suitable reason. Beyond
1128 the first element LD/ST Failfirst is arbitrarily speculative and 100%
1129 non-deterministic.
1130 * CR-based data-dependent first on the other hand MUST NOT truncate VL
1131 arbitrarily to a length decided by the hardware: VL MUST only be
1132 truncated based explicitly on whether a test fails. This because it is
1133 a precise Deterministic test on which algorithms can and will will rely.
1134
1135 **Floating-point Exceptions**
1136
1137 When Floating-point exceptions are enabled VL must be truncated at
1138 the point where the Exception appears not to have occurred. If `VLi`
1139 is set then VL must include the faulting element, and thus the faulting
1140 element will always raise its exception. If however `VLi` is clear then
1141 VL **excludes** the faulting element and thus the exception will **never**
1142 be raised.
1143
1144 Although very strongly discouraged the Exception Mode that permits
1145 Floating Point Exception notification to arrive too late to unwind
1146 is permitted (under protest, due it violating the otherwise 100%
1147 Deterministic nature of Data-dependent Fail-first) and is `UNDEFINED`
1148 behaviour.
1149
1150 **Use of lax FP Exception Notification Mode could result in parallel
1151 computations proceeding with invalid results that have to be explicitly
1152 detected, whereas with the strict FP Execption Mode enabled, FFirst
1153 truncates VL, allows subsequent parallel computation to avoid the
1154 exceptions entirely**
1155
1156 ## Data-dependent fail-first on CR operations (crand etc)
1157
1158 Operations that actually produce or alter CR Field as a result have
1159 their own SVP64 Mode, described in [[sv/cr_ops]].
1160
1161 ## pred-result mode
1162
1163 This mode merges common CR testing with predication, saving on instruction
1164 count. Below is the pseudocode excluding predicate zeroing and elwidth
1165 overrides. Note that the pseudocode for SVP64 CR-ops is slightly different.
1166
1167 ```
1168 for i in range(VL):
1169 # predication test, skip all masked out elements.
1170 if predicate_masked_out(i):
1171 continue
1172 result = op(iregs[RA+i], iregs[RB+i])
1173 CRnew = analyse(result) # calculates eq/lt/gt
1174 # Rc=1 always stores the CR field
1175 if Rc=1 or RC1:
1176 CR.field[offs+i] = CRnew
1177 # now test CR, similar to branch
1178 if RC1 or CR.field[BO[0:1]] != BO[2]:
1179 continue # test failed: cancel store
1180 # result optionally stored but CR always is
1181 iregs[RT+i] = result
1182 ```
1183
1184 The reason for allowing the CR element to be stored is so that
1185 post-analysis of the CR Vector may be carried out. For example:
1186 Saturation may have occurred (and been prevented from updating, by the
1187 test) but it is desirable to know *which* elements fail saturation.
1188
1189 Note that RC1 Mode basically turns all operations into `cmp`. The
1190 calculation is performed but it is only the CR that is written. The
1191 element result is *always* discarded, never written (just like `cmp`).
1192
1193 Note that predication is still respected: predicate zeroing is slightly
1194 different: elements that fail the CR test *or* are masked out are zero'd.
1195
1196 --------
1197
1198 \newpage{}
1199
1200 # SV Load and Store
1201
1202 **Rationale**
1203
1204 All Vector ISAs dating back fifty years have extensive and comprehensive
1205 Load and Store operations that go far beyond the capabilities of Scalar
1206 RISC and most CISC processors, yet at their heart on an individual element
1207 basis may be found to be no different from RISC Scalar equivalents.
1208
1209 The resource savings from Vector LD/ST are significant and stem
1210 from the fact that one single instruction can trigger a dozen (or in
1211 some microarchitectures such as Cray or NEC SX Aurora) hundreds of
1212 element-level Memory accesses.
1213
1214 Additionally, and simply: if the Arithmetic side of an ISA supports
1215 Vector Operations, then in order to keep the ALUs 100% occupied the
1216 Memory infrastructure (and the ISA itself) correspondingly needs Vector
1217 Memory Operations as well.
1218
1219 Vectorised Load and Store also presents an extra dimension (literally)
1220 which creates scenarios unique to Vector applications, that a Scalar (and
1221 even a SIMD) ISA simply never encounters. SVP64 endeavours to add the
1222 modes typically found in *all* Scalable Vector ISAs, without changing the
1223 behaviour of the underlying Base (Scalar) v3.0B operations in any way.
1224 (The sole apparent exception is Post-Increment Mode on LD/ST-update
1225 instructions)
1226
1227 ## Modes overview
1228
1229 Vectorisation of Load and Store requires creation, from scalar operations,
1230 a number of different modes:
1231
1232 * **fixed aka "unit" stride** - contiguous sequence with no gaps
1233 * **element strided** - sequential but regularly offset, with gaps
1234 * **vector indexed** - vector of base addresses and vector of offsets
1235 * **Speculative fail-first** - where it makes sense to do so
1236 * **Structure Packing** - covered in SV by [[sv/remap]] and Pack/Unpack Mode.
1237
1238 *Despite being constructed from Scalar LD/ST none of these Modes exist
1239 or make sense in any Scalar ISA. They **only** exist in Vector ISAs*
1240
1241 Also included in SVP64 LD/ST is both signed and unsigned Saturation,
1242 as well as Element-width overrides and Twin-Predication.
1243
1244 Note also that Indexed [[sv/remap]] mode may be applied to both v3.0
1245 LD/ST Immediate instructions *and* v3.0 LD/ST Indexed instructions.
1246 LD/ST-Indexed should not be conflated with Indexed REMAP mode:
1247 clarification is provided below.
1248
1249 **Determining the LD/ST Modes**
1250
1251 A minor complication (caused by the retro-fitting of modern Vector
1252 features to a Scalar ISA) is that certain features do not exactly make
1253 sense or are considered a security risk. Fail-first on Vector Indexed
1254 would allow attackers to probe large numbers of pages from userspace,
1255 where strided fail-first (by creating contiguous sequential LDs) does not.
1256
1257 In addition, reduce mode makes no sense. Realistically we need an
1258 alternative table definition for [[sv/svp64]] `RM.MODE`. The following
1259 modes make sense:
1260
1261 * saturation
1262 * predicate-result (mostly for cache-inhibited LD/ST)
1263 * simple (no augmentation)
1264 * fail-first (where Vector Indexed is banned)
1265 * Signed Effective Address computation (Vector Indexed only)
1266
1267 More than that however it is necessary to fit the usual Vector ISA
1268 capabilities onto both Power ISA LD/ST with immediate and to LD/ST
1269 Indexed. They present subtly different Mode tables, which, due to lack
1270 of space, have the following quirks:
1271
1272 * LD/ST Immediate has no individual control over src/dest zeroing,
1273 whereas LD/ST Indexed does.
1274 * LD/ST Indexed has limited zeroing on pred-result, LD/ST Immediate has
1275 *no* option to select zeroing on pred-result.
1276
1277 ## Format and fields
1278
1279 Fields used in tables below:
1280
1281 * **sz / dz** if predication is enabled will put zeros into the dest
1282 (or as src in the case of twin pred) when the predicate bit is zero.
1283 otherwise the element is ignored or skipped, depending on context.
1284 * **zz**: both sz and dz are set equal to this flag.
1285 * **inv CR bit** just as in branches (BO) these bits allow testing of
1286 a CR bit and whether it is set (inv=0) or unset (inv=1)
1287 * **N** sets signed/unsigned saturation.
1288 * **RC1** as if Rc=1, stores CRs *but not the result*
1289 * **SEA** - Signed Effective Address, if enabled performs sign-extension on
1290 registers that have been reduced due to elwidth overrides
1291 * **PI** - post-increment mode (applies to LD/ST with update only).
1292 the Effective Address utilised is always just RA, i.e. the computation of
1293 EA is stored in RA **after** it is actually used.
1294 * **LF** - Load/Store Fail or Fault First: for any reason Load or Store Vectors
1295 may be truncated to (at least) one element, and VL altered to indicate such.
1296
1297 **LD/ST immediate**
1298
1299 The table for [[sv/svp64]] for `immed(RA)` which is `RM.MODE`
1300 (bits 19:23 of `RM`) is:
1301
1302 | 0-1 | 2 | 3 4 | description |
1303 | --- | --- |---------|--------------------------- |
1304 | 00 | 0 | zz els | simple mode |
1305 | 00 | 1 | PI LF | post-increment and Fault-First |
1306 | 01 | inv | CR-bit | Rc=1: ffirst CR sel |
1307 | 01 | inv | els RC1 | Rc=0: ffirst z/nonz |
1308 | 10 | N | zz els | sat mode: N=0/1 u/s |
1309 | 11 | inv | CR-bit | Rc=1: pred-result CR sel |
1310 | 11 | inv | els RC1 | Rc=0: pred-result z/nonz |
1311
1312 The `els` bit is only relevant when `RA.isvec` is clear: this indicates
1313 whether stride is unit or element:
1314
1315 ```
1316 if RA.isvec:
1317 svctx.ldstmode = indexed
1318 elif els == 0:
1319 svctx.ldstmode = unitstride
1320 elif immediate != 0:
1321 svctx.ldstmode = elementstride
1322 ```
1323
1324 An immediate of zero is a safety-valve to allow `LD-VSPLAT`: in effect
1325 the multiplication of the immediate-offset by zero results in reading from
1326 the exact same memory location, *even with a Vector register*. (Normally
1327 this type of behaviour is reserved for the mapreduce modes)
1328
1329 For `LD-VSPLAT`, on non-cache-inhibited Loads, the read can occur just
1330 the once and be copied, rather than hitting the Data Cache multiple
1331 times with the same memory read at the same location. The benefit of
1332 Cache-inhibited LD-splats is that it allows for memory-mapped peripherals
1333 to have multiple data values read in quick succession and stored in
1334 sequentially numbered registers (but, see Note below).
1335
1336 For non-cache-inhibited ST from a vector source onto a scalar destination:
1337 with the Vector loop effectively creating multiple memory writes to
1338 the same location, we can deduce that the last of these will be the
1339 "successful" one. Thus, implementations are free and clear to optimise
1340 out the overwriting STs, leaving just the last one as the "winner".
1341 Bear in mind that predicate masks will skip some elements (in source
1342 non-zeroing mode). Cache-inhibited ST operations on the other hand
1343 **MUST** write out a Vector source multiple successive times to the exact
1344 same Scalar destination. Just like Cache-inhibited LDs, multiple values
1345 may be written out in quick succession to a memory-mapped peripheral
1346 from sequentially-numbered registers.
1347
1348 Note that any memory location may be Cache-inhibited
1349 (Power ISA v3.1, Book III, 1.6.1, p1033)
1350
1351 *Programmer's Note: an immediate also with a Scalar source as a "VSPLAT"
1352 mode is simply not possible: there are not enough Mode bits. One single
1353 Scalar Load operation may be used instead, followed by any arithmetic
1354 operation (including a simple mv) in "Splat" mode.*
1355
1356 **LD/ST Indexed**
1357
1358 The modes for `RA+RB` indexed version are slightly different
1359 but are the same `RM.MODE` bits (19:23 of `RM`):
1360
1361 | 0-1 | 2 | 3 4 | description |
1362 | --- | --- |---------|-------------------------- |
1363 | 00 | SEA | dz sz | simple mode |
1364 | 01 | SEA | dz sz | Strided (scalar only source) |
1365 | 10 | N | dz sz | sat mode: N=0/1 u/s |
1366 | 11 | inv | CR-bit | Rc=1: pred-result CR sel |
1367 | 11 | inv | zz RC1 | Rc=0: pred-result z/nonz |
1368
1369 Vector Indexed Strided Mode is qualified as follows:
1370
1371 if mode = 0b01 and !RA.isvec and !RB.isvec:
1372 svctx.ldstmode = elementstride
1373
1374 A summary of the effect of Vectorisation of src or dest:
1375
1376 ```
1377 imm(RA) RT.v RA.v no stride allowed
1378 imm(RA) RT.s RA.v no stride allowed
1379 imm(RA) RT.v RA.s stride-select allowed
1380 imm(RA) RT.s RA.s not vectorised
1381 RA,RB RT.v {RA|RB}.v Standard Indexed
1382 RA,RB RT.s {RA|RB}.v Indexed but single LD (no VSPLAT)
1383 RA,RB RT.v {RA&RB}.s VSPLAT possible. stride selectable
1384 RA,RB RT.s {RA&RB}.s not vectorised (scalar identity)
1385 ```
1386
1387 Signed Effective Address computation is only relevant for Vector Indexed
1388 Mode, when elwidth overrides are applied. The source override applies to
1389 RB, and before adding to RA in order to calculate the Effective Address,
1390 if SEA is set RB is sign-extended from elwidth bits to the full 64 bits.
1391 For other Modes (ffirst, saturate), all EA computation with elwidth
1392 overrides is unsigned.
1393
1394 Note that cache-inhibited LD/ST when VSPLAT is activated will perform
1395 **multiple** LD/ST operations, sequentially. Even with scalar src
1396 a Cache-inhibited LD will read the same memory location *multiple
1397 times*, storing the result in successive Vector destination registers.
1398 This because the cache-inhibit instructions are typically used to read
1399 and write memory-mapped peripherals. If a genuine cache-inhibited
1400 LD-VSPLAT is required then a single *scalar* cache-inhibited LD should
1401 be performed, followed by a VSPLAT-augmented mv, copying the one *scalar*
1402 value into multiple register destinations.
1403
1404 Note also that cache-inhibited VSPLAT with Predicate-result is possible.
1405 This allows for example to issue a massive batch of memory-mapped
1406 peripheral reads, stopping at the first NULL-terminated character and
1407 truncating VL to that point. No branch is needed to issue that large
1408 burst of LDs, which may be valuable in Embedded scenarios.
1409
1410 ## Vectorisation of Scalar Power ISA v3.0B
1411
1412 Scalar Power ISA Load/Store operations may be seen from their
1413 pseudocode to be of the form:
1414
1415 ```
1416 lbux RT, RA, RB
1417 EA <- (RA) + (RB)
1418 RT <- MEM(EA)
1419 ```
1420
1421 and for immediate variants:
1422
1423 ```
1424 lb RT,D(RA)
1425 EA <- RA + EXTS(D)
1426 RT <- MEM(EA)
1427 ```
1428
1429 Thus in the first example, the source registers may each be independently
1430 marked as scalar or vector, and likewise the destination; in the second
1431 example only the one source and one dest may be marked as scalar or
1432 vector.
1433
1434 Thus we can see that Vector Indexed may be covered, and, as demonstrated
1435 with the pseudocode below, the immediate can be used to give unit
1436 stride or element stride. With there being no way to tell which from
1437 the Power v3.0B Scalar opcode alone, the choice is provided instead by
1438 the SV Context.
1439
1440 ```
1441 # LD not VLD! format - ldop RT, immed(RA)
1442 # op_width: lb=1, lh=2, lw=4, ld=8
1443 op_load(RT, RA, op_width, immed, svctx, RAupdate):
1444  ps = get_pred_val(FALSE, RA); # predication on src
1445  pd = get_pred_val(FALSE, RT); # ... AND on dest
1446  for (i=0, j=0, u=0; i < VL && j < VL;):
1447 # skip nonpredicates elements
1448 if (RA.isvec) while (!(ps & 1<<i)) i++;
1449 if (RAupdate.isvec) while (!(ps & 1<<u)) u++;
1450 if (RT.isvec) while (!(pd & 1<<j)) j++;
1451 if postinc:
1452 offs = 0; # added afterwards
1453 if RA.isvec: srcbase = ireg[RA+i]
1454 else srcbase = ireg[RA]
1455 elif svctx.ldstmode == elementstride:
1456 # element stride mode
1457 srcbase = ireg[RA]
1458 offs = i * immed # j*immed for a ST
1459 elif svctx.ldstmode == unitstride:
1460 # unit stride mode
1461 srcbase = ireg[RA]
1462 offs = immed + (i * op_width) # j*op_width for ST
1463 elif RA.isvec:
1464 # quirky Vector indexed mode but with an immediate
1465 srcbase = ireg[RA+i]
1466 offs = immed;
1467 else
1468 # standard scalar mode (but predicated)
1469 # no stride multiplier means VSPLAT mode
1470 srcbase = ireg[RA]
1471 offs = immed
1472
1473 # compute EA
1474 EA = srcbase + offs
1475 # load from memory
1476 ireg[RT+j] <= MEM[EA];
1477 # check post-increment of EA
1478 if postinc: EA = srcbase + immed;
1479 # update RA?
1480 if RAupdate: ireg[RAupdate+u] = EA;
1481 if (!RT.isvec)
1482 break # destination scalar, end now
1483 if (RA.isvec) i++;
1484 if (RAupdate.isvec) u++;
1485 if (RT.isvec) j++;
1486 ```
1487
1488 Indexed LD is:
1489
1490 ```
1491 # format: ldop RT, RA, RB
1492 function op_ldx(RT, RA, RB, RAupdate=False) # LD not VLD!
1493  ps = get_pred_val(FALSE, RA); # predication on src
1494  pd = get_pred_val(FALSE, RT); # ... AND on dest
1495  for (i=0, j=0, k=0, u=0; i < VL && j < VL && k < VL):
1496 # skip nonpredicated RA, RB and RT
1497 if (RA.isvec) while (!(ps & 1<<i)) i++;
1498 if (RAupdate.isvec) while (!(ps & 1<<u)) u++;
1499 if (RB.isvec) while (!(ps & 1<<k)) k++;
1500 if (RT.isvec) while (!(pd & 1<<j)) j++;
1501 if svctx.ldstmode == elementstride:
1502 EA = ireg[RA] + ireg[RB]*j # register-strided
1503 else
1504 EA = ireg[RA+i] + ireg[RB+k] # indexed address
1505 if RAupdate: ireg[RAupdate+u] = EA
1506 ireg[RT+j] <= MEM[EA];
1507 if (!RT.isvec)
1508 break # destination scalar, end immediately
1509 if (RA.isvec) i++;
1510 if (RAupdate.isvec) u++;
1511 if (RB.isvec) k++;
1512 if (RT.isvec) j++;
1513 ```
1514
1515 Note that Element-Strided uses the Destination Step because with both
1516 sources being Scalar as a prerequisite condition of activation of
1517 Element-Stride Mode, the source step (being Scalar) would never advance.
1518
1519 Note in both cases that [[sv/svp64]] allows RA-as-a-dest in "update"
1520 mode (`ldux`) to be effectively a *completely different* register from
1521 RA-as-a-source. This because there is room in svp64 to extend RA-as-src
1522 as well as RA-as-dest, both independently as scalar or vector *and*
1523 independently extending their range.
1524
1525 *Programmer's note: being able to set RA-as-a-source as separate from
1526 RA-as-a-destination as Scalar is **extremely valuable** once it is
1527 remembered that Simple-V element operations must be in Program Order,
1528 especially in loops, for saving on multiple address computations. Care
1529 does have to be taken however that RA-as-src is not overwritten by
1530 RA-as-dest unless intentionally desired, especially in element-strided
1531 Mode.*
1532
1533 ## LD/ST Indexed vs Indexed REMAP
1534
1535 Unfortunately the word "Indexed" is used twice in completely different
1536 contexts, potentially causing confusion.
1537
1538 * There has existed instructions in the Power ISA `ld RT,RA,RB` since
1539 its creation: these are called "LD/ST Indexed" instructions and their
1540 name and meaning is well-established.
1541 * There now exists, in Simple-V, a REMAP mode called "Indexed"
1542 Mode that can be applied to *any* instruction **including those
1543 named LD/ST Indexed**.
1544
1545 Whilst it may be costly in terms of register reads to allow REMAP Indexed
1546 Mode to be applied to any Vectorised LD/ST Indexed operation such as
1547 `sv.ld *RT,RA,*RB`, or even misleadingly labelled as redundant, firstly
1548 the strict application of the RISC Paradigm that Simple-V follows makes
1549 it awkward to consider *preventing* the application of Indexed REMAP to
1550 such operations, and secondly they are not actually the same at all.
1551
1552 Indexed REMAP, as applied to RB in the instruction `sv.ld *RT,RA,*RB`
1553 effectively performs an *in-place* re-ordering of the offsets, RB.
1554 To achieve the same effect without Indexed REMAP would require taking
1555 a *copy* of the Vector of offsets starting at RB, manually explicitly
1556 reordering them, and finally using the copy of re-ordered offsets in a
1557 non-REMAP'ed `sv.ld`. Using non-strided LD as an example, pseudocode
1558 showing what actually occurs, where the pseudocode for `indexed_remap`
1559 may be found in [[sv/remap]]:
1560
1561 ```
1562 # sv.ld *RT,RA,*RB with Index REMAP applied to RB
1563 for i in 0..VL-1:
1564 if remap.indexed:
1565 rb_idx = indexed_remap(i) # remap
1566 else:
1567 rb_idx = i # use the index as-is
1568 EA = GPR(RA) + GPR(RB+rb_idx)
1569 GPR(RT+i) = MEM(EA, 8)
1570 ```
1571
1572 Thus it can be seen that the use of Indexed REMAP saves copying
1573 and manual reordering of the Vector of RB offsets.
1574
1575 ## LD/ST ffirst
1576
1577 LD/ST ffirst treats the first LD/ST in a vector (element 0 if REMAP
1578 is not active) as an ordinary one, with all behaviour with respect to
1579 Interrupts Exceptions Page Faults Memory Management being identical
1580 in every regard to Scalar v3.0 Power ISA LD/ST. However for elements
1581 1 and above, if an exception would occur, then VL is **truncated**
1582 to the previous element: the exception is **not** then raised because
1583 the LD/ST that would otherwise have caused an exception is *required*
1584 to be cancelled. Additionally an implementor may choose to truncate VL
1585 for any arbitrary reason *except for the very first*.
1586
1587 ffirst LD/ST to multiple pages via a Vectorised Index base is
1588 considered a security risk due to the abuse of probing multiple
1589 pages in rapid succession and getting speculative feedback on which
1590 pages would fail. Therefore Vector Indexed LD/ST is prohibited
1591 entirely, and the Mode bit instead used for element-strided LD/ST.
1592
1593 ```
1594 for(i = 0; i < VL; i++)
1595 reg[rt + i] = mem[reg[ra] + i * reg[rb]];
1596 ```
1597
1598 High security implementations where any kind of speculative probing of
1599 memory pages is considered a risk should take advantage of the fact
1600 that implementations may truncate VL at any point, without requiring
1601 software to be rewritten and made non-portable. Such implementations may
1602 choose to *always* set VL=1 which will have the effect of terminating
1603 any speculative probing (and also adversely affect performance), but
1604 will at least not require applications to be rewritten.
1605
1606 Low-performance simpler hardware implementations may also choose (always)
1607 to also set VL=1 as the bare minimum compliant implementation of LD/ST
1608 Fail-First. It is however critically important to remember that the first
1609 element LD/ST **MUST** be treated as an ordinary LD/ST, i.e. **MUST**
1610 raise exceptions exactly like an ordinary LD/ST.
1611
1612 For ffirst LD/STs, VL may be truncated arbitrarily to a nonzero value
1613 for any implementation-specific reason. For example: it is perfectly
1614 reasonable for implementations to alter VL when ffirst LD or ST operations
1615 are initiated on a nonaligned boundary, such that within a loop the
1616 subsequent iteration of that loop begins the following ffirst LD/ST
1617 operations on an aligned boundary such as the beginning of a cache line,
1618 or beginning of a Virtual Memory page. Likewise, to reduce workloads or
1619 balance resources.
1620
1621 Vertical-First Mode is slightly strange in that only one element at a time
1622 is ever executed anyway. Given that programmers may legitimately choose
1623 to alter srcstep and dststep in non-sequential order as part of explicit
1624 loops, it is neither possible nor safe to make speculative assumptions
1625 about future LD/STs. Therefore, Fail-First LD/ST in Vertical-First is
1626 `UNDEFINED`. This is very different from Arithmetic (Data-dependent)
1627 FFirst where Vertical-First Mode is fully deterministic, not speculative.
1628
1629 ## LOAD/STORE Elwidths <a name="elwidth"></a>
1630
1631 Loads and Stores are almost unique in that the Power Scalar ISA
1632 provides a width for the operation (lb, lh, lw, ld). Only `extsb` and
1633 others like it provide an explicit operation width. There are therefore
1634 *three* widths involved:
1635
1636 * operation width (lb=8, lh=16, lw=32, ld=64)
1637 * src element width override (8/16/32/default)
1638 * destination element width override (8/16/32/default)
1639
1640 Some care is therefore needed to express and make clear the transformations,
1641 which are expressly in this order:
1642
1643 * Calculate the Effective Address from RA at full width
1644 but (on Indexed Load) allow srcwidth overrides on RB
1645 * Load at the operation width (lb/lh/lw/ld) as usual
1646 * byte-reversal as usual
1647 * Non-saturated mode:
1648 - zero-extension or truncation from operation width to dest elwidth
1649 - place result in destination at dest elwidth
1650 * Saturated mode:
1651 - Sign-extension or truncation from operation width to dest width
1652 - signed/unsigned saturation down to dest elwidth
1653
1654 In order to respect Power v3.0B Scalar behaviour the memory side
1655 is treated effectively as completely separate and distinct from SV
1656 augmentation. This is primarily down to quirks surrounding LE/BE and
1657 byte-reversal.
1658
1659 It is rather unfortunately possible to request an elwidth override on
1660 the memory side which does not mesh with the overridden operation width:
1661 these result in `UNDEFINED` behaviour. The reason is that the effect
1662 of attempting a 64-bit `sv.ld` operation with a source elwidth override
1663 of 8/16/32 would result in overlapping memory requests, particularly
1664 on unit and element strided operations. Thus it is `UNDEFINED` when
1665 the elwidth is smaller than the memory operation width. Examples include
1666 `sv.lw/sw=16/els` which requests (overlapping) 4-byte memory reads offset
1667 from each other at 2-byte intervals. Store likewise is also `UNDEFINED`
1668 where the dest elwidth override is less than the operation width.
1669
1670 Note the following regarding the pseudocode to follow:
1671
1672 * `scalar identity behaviour` SV Context parameter conditions turn this
1673 into a straight absolute fully-compliant Scalar v3.0B LD operation
1674 * `brev` selects whether the operation is the byte-reversed variant (`ldbrx`
1675 rather than `ld`)
1676 * `op_width` specifies the operation width (`lb`, `lh`, `lw`, `ld`) as
1677 a "normal" part of Scalar v3.0B LD
1678 * `imm_offs` specifies the immediate offset `ld r3, imm_offs(r5)`, again
1679 as a "normal" part of Scalar v3.0B LD
1680 * `svctx` specifies the SV Context and includes VL as well as
1681 source and destination elwidth overrides.
1682
1683 Below is the pseudocode for Unit-Strided LD (which includes Vector
1684 capability). Observe in particular that RA, as the base address in both
1685 Immediate and Indexed LD/ST, does not have element-width overriding
1686 applied to it.
1687
1688 Note that predication, predication-zeroing, and other modes except
1689 saturation have all been removed, for clarity and simplicity:
1690
1691 ```
1692 # LD not VLD!
1693 # this covers unit stride mode and a type of vector offset
1694 function op_ld(RT, RA, op_width, imm_offs, svctx)
1695 for (int i = 0, int j = 0; i < svctx.VL && j < svctx.VL):
1696 if not svctx.unit/el-strided:
1697 # strange vector mode, compute 64 bit address which is
1698 # not polymorphic! elwidth hardcoded to 64 here
1699 srcbase = get_polymorphed_reg(RA, 64, i)
1700 else:
1701 # unit / element stride mode, compute 64 bit address
1702 srcbase = get_polymorphed_reg(RA, 64, 0)
1703 # adjust for unit/el-stride
1704 srcbase += ....
1705
1706 # read the underlying memory
1707 memread <= MEM(srcbase + imm_offs, op_width)
1708
1709 # check saturation.
1710 if svpctx.saturation_mode:
1711 # ... saturation adjustment...
1712 memread = clamp(memread, op_width, svctx.dest_elwidth)
1713 else:
1714 # truncate/extend to over-ridden dest width.
1715 memread = adjust_wid(memread, op_width, svctx.dest_elwidth)
1716
1717 # takes care of inserting memory-read (now correctly byteswapped)
1718 # into regfile underlying LE-defined order, into the right place
1719 # within the NEON-like register, respecting destination element
1720 # bitwidth, and the element index (j)
1721 set_polymorphed_reg(RT, svctx.dest_elwidth, j, memread)
1722
1723 # increments both src and dest element indices (no predication here)
1724 i++;
1725 j++;
1726 ```
1727
1728 Note above that the source elwidth is *not used at all* in LD-immediate.
1729
1730 For LD/Indexed, the key is that in the calculation of the Effective Address,
1731 RA has no elwidth override but RB does. Pseudocode below is simplified
1732 for clarity: predication and all modes except saturation are removed:
1733
1734 ```
1735 # LD not VLD! ld*rx if brev else ld*
1736 function op_ld(RT, RA, RB, op_width, svctx, brev)
1737 for (int i = 0, int j = 0; i < svctx.VL && j < svctx.VL):
1738 if not svctx.el-strided:
1739 # RA not polymorphic! elwidth hardcoded to 64 here
1740 srcbase = get_polymorphed_reg(RA, 64, i)
1741 else:
1742 # element stride mode, again RA not polymorphic
1743 srcbase = get_polymorphed_reg(RA, 64, 0)
1744 # RB *is* polymorphic
1745 offs = get_polymorphed_reg(RB, svctx.src_elwidth, i)
1746 # sign-extend
1747 if svctx.SEA: offs = sext(offs, svctx.src_elwidth, 64)
1748
1749 # takes care of (merges) processor LE/BE and ld/ldbrx
1750 bytereverse = brev XNOR MSR.LE
1751
1752 # read the underlying memory
1753 memread <= MEM(srcbase + offs, op_width)
1754
1755 # optionally performs byteswap at op width
1756 if (bytereverse):
1757 memread = byteswap(memread, op_width)
1758
1759 if svpctx.saturation_mode:
1760 # ... saturation adjustment...
1761 memread = clamp(memread, op_width, svctx.dest_elwidth)
1762 else:
1763 # truncate/extend to over-ridden dest width.
1764 memread = adjust_wid(memread, op_width, svctx.dest_elwidth)
1765
1766 # takes care of inserting memory-read (now correctly byteswapped)
1767 # into regfile underlying LE-defined order, into the right place
1768 # within the NEON-like register, respecting destination element
1769 # bitwidth, and the element index (j)
1770 set_polymorphed_reg(RT, svctx.dest_elwidth, j, memread)
1771
1772 # increments both src and dest element indices (no predication here)
1773 i++;
1774 j++;
1775 ```
1776
1777 ## Remapped LD/ST
1778
1779 In the [[sv/remap]] page the concept of "Remapping" is described. Whilst
1780 it is expensive to set up (2 64-bit opcodes minimum) it provides a way to
1781 arbitrarily perform 1D, 2D and 3D "remapping" of up to 64 elements worth
1782 of LDs or STs. The usual interest in such re-mapping is for example in
1783 separating out 24-bit RGB channel data into separate contiguous registers.
1784
1785 REMAP easily covers this capability, and with dest elwidth overrides
1786 and saturation may do so with built-in conversion that would normally
1787 require additional width-extension, sign-extension and min/max Vectorised
1788 instructions as post-processing stages.
1789
1790 Thus we do not need to provide specialist LD/ST "Structure Packed" opcodes
1791 because the generic abstracted concept of "Remapping", when applied to
1792 LD/ST, will give that same capability, with far more flexibility.
1793
1794 It is worth noting that Pack/Unpack Modes of SVSTATE, which may be
1795 established through `svstep`, are also an easy way to perform regular
1796 Structure Packing, at the vec2/vec3/vec4 granularity level. Beyond that,
1797 REMAP will need to be used.
1798
1799 --------
1800
1801 \newpage{}
1802
1803 # Condition Register SVP64 Operations
1804
1805 Condition Register Fields are only 4 bits wide: this presents some
1806 interesting conceptual challenges for SVP64, which was designed
1807 primarily for vectors of arithmetic and logical operations. However
1808 if predicates may be bits of CR Fields it makes sense to extend
1809 Simple-V to cover CR Operations, especially given that Vectorised Rc=1
1810 may be processed by Vectorised CR Operations tbat usefully in turn
1811 may become Predicate Masks to yet more Vector operations, like so:
1812
1813 ```
1814 sv.cmpi/ew=8 *B,*ra,0 # compare bytes against zero
1815 sv.cmpi/ew=8 *B2,*ra,13. # and against newline
1816 sv.cror PM.EQ,B.EQ,B2.EQ # OR compares to create mask
1817 sv.stb/sm=EQ ... # store only nonzero/newline
1818 ```
1819
1820 Element width however is clearly meaningless for a 4-bit collation of
1821 Conditions, EQ LT GE SO. Likewise, arithmetic saturation (an important
1822 part of Arithmetic SVP64) has no meaning. An alternative Mode Format is
1823 required, and given that elwidths are meaningless for CR Fields the bits
1824 in SVP64 `RM` may be used for other purposes.
1825
1826 This alternative mapping **only** applies to instructions that **only**
1827 reference a CR Field or CR bit as the sole exclusive result. This section
1828 **does not** apply to instructions which primarily produce arithmetic
1829 results that also, as an aside, produce a corresponding CR Field (such as
1830 when Rc=1). Instructions that involve Rc=1 are definitively arithmetic
1831 in nature, where the corresponding Condition Register Field can be
1832 considered to be a "co-result". Such CR Field "co-result" arithmeric
1833 operations are firmly out of scope for this section, being covered fully
1834 by [[sv/normal]].
1835
1836 * Examples of v3.0B instructions to which this section does
1837 apply is
1838 - `mfcr` and `cmpi` (3 bit operands) and
1839 - `crnor` and `crand` (5 bit operands).
1840 * Examples to which this section does **not** apply include
1841 `fadds.` and `subf.` which both produce arithmetic results
1842 (and a CR Field co-result).
1843
1844 The CR Mode Format still applies to `sv.cmpi` because despite
1845 taking a GPR as input, the output from the Base Scalar v3.0B `cmpi`
1846 instruction is purely to a Condition Register Field.
1847
1848 Other modes are still applicable and include:
1849
1850 * **Data-dependent fail-first**.
1851 useful to truncate VL based on analysis of a Condition Register result bit.
1852 * **Reduction**.
1853 Reduction is useful for analysing a Vector of Condition Register Fields
1854 and reducing it to one single Condition Register Field.
1855
1856 Predicate-result does not make any sense because when Rc=1 a co-result
1857 is created (a CR Field). Testing the co-result allows the decision to
1858 be made to store or not store the main result, and for CR Ops the CR
1859 Field result *is* the main result.
1860
1861 ## Format
1862
1863 SVP64 RM `MODE` (includes `ELWIDTH_SRC` bits) for CR-based operations:
1864
1865 |6 | 7 |19-20| 21 | 22 23 | description |
1866 |--|---|-----| --- |---------|----------------- |
1867 |/ | / |0 RG | 0 | dz sz | simple mode |
1868 |/ | / |0 RG | 1 | dz sz | scalar reduce mode (mapreduce) |
1869 |zz|SNZ|1 VLI| inv | CR-bit | Ffirst 3-bit mode |
1870 |/ |SNZ|1 VLI| inv | dz sz | Ffirst 5-bit mode (implies CR-bit from result) |
1871
1872 Fields:
1873
1874 * **sz / dz** if predication is enabled will put zeros into the dest
1875 (or as src in the case of twin pred) when the predicate bit is zero.
1876 otherwise the element is ignored or skipped, depending on context.
1877 * **zz** set both sz and dz equal to this flag
1878 * **SNZ** In fail-first mode, on the bit being tested, when sz=1 and
1879 SNZ=1 a value "1" is put in place of "0".
1880 * **inv CR-bit** just as in branches (BO) these bits allow testing of
1881 a CR bit and whether it is set (inv=0) or unset (inv=1)
1882 * **RG** inverts the Vector Loop order (VL-1 downto 0) rather
1883 than the normal 0..VL-1
1884 * **SVM** sets "subvector" reduce mode
1885 * **VLi** VL inclusive: in fail-first mode, the truncation of
1886 VL *includes* the current element at the failure point rather
1887 than excludes it from the count.
1888
1889 ## Data-dependent fail-first on CR operations
1890
1891 The principle of data-dependent fail-first is that if, during the course
1892 of sequentially evaluating an element's Condition Test, one such test
1893 is encountered which fails, then VL (Vector Length) is truncated (set)
1894 at that point. In the case of Arithmetic SVP64 Operations the Condition
1895 Register Field generated from Rc=1 is used as the basis for the truncation
1896 decision. However with CR-based operations that CR Field result to be
1897 tested is provided *by the operation itself*.
1898
1899 Data-dependent SVP64 Vectorised Operations involving the creation
1900 or modification of a CR can require an extra two bits, which are not
1901 available in the compact space of the SVP64 RM `MODE` Field. With the
1902 concept of element width overrides being meaningless for CR Fields it
1903 is possible to use the `ELWIDTH` field for alternative purposes.
1904
1905 Condition Register based operations such as `sv.mfcr` and `sv.crand`
1906 can thus be made more flexible. However the rules that apply in this
1907 section also apply to future CR-based instructions.
1908
1909 There are two primary different types of CR operations:
1910
1911 * Those which have a 3-bit operand field (referring to a CR Field)
1912 * Those which have a 5-bit operand (referring to a bit within the
1913 whole 32-bit CR)
1914
1915 Examining these two types it is observed that the difference may
1916 be considered to be that the 5-bit variant *already* provides the
1917 prerequisite information about which CR Field bit (EQ, GE, LT, SO) is
1918 to be operated on by the instruction. Thus, logically, we may set the
1919 following rule:
1920
1921 * When a 5-bit CR Result field is used in an instruction, the
1922 5-bit variant of Data-Dependent Fail-First
1923 must be used. i.e. the bit of the CR field to be tested is
1924 the one that has just been modified (created) by the operation.
1925 * When a 3-bit CR Result field is used the 3-bit variant
1926 must be used, providing as it does the missing `CRbit` field
1927 in order to select which CR Field bit of the result shall
1928 be tested (EQ, LE, GE, SO)
1929
1930 The reason why the 3-bit CR variant needs the additional CR-bit field
1931 should be obvious from the fact that the 3-bit CR Field from the base
1932 Power ISA v3.0B operation clearly does not contain and is missing the
1933 two CR Field Selector bits. Thus, these two bits (to select EQ, LE,
1934 GE or SO) must be provided in another way.
1935
1936 Examples of the former type:
1937
1938 * crand, cror, crnor. These all are 5-bit (BA, BB, BT). The bit
1939 to be tested against `inv` is the one selected by `BT`
1940 * mcrf. This has only 3-bit (BF, BFA). In order to select the
1941 bit to be tested, the alternative encoding must be used.
1942 With `CRbit` coming from the SVP64 RM bits 22-23 the bit
1943 of BF to be tested is identified.
1944
1945 Just as with SVP64 [[sv/branches]] there is the option to truncate
1946 VL to include the element being tested (`VLi=1`) and to exclude it
1947 (`VLi=0`).
1948
1949 Also exactly as with [[sv/normal]] fail-first, VL cannot, unlike
1950 [[sv/ldst]], be set to an arbitrary value. Deterministic behaviour
1951 is *required*.
1952
1953 ## Reduction and Iteration
1954
1955 Bearing in mind as described in the svp64 Appendix, SVP64 Horizontal
1956 Reduction is a deterministic schedule on top of base Scalar v3.0
1957 operations, the same rules apply to CR Operations, i.e. that programmers
1958 must follow certain conventions in order for an *end result* of a
1959 reduction to be achieved. Unlike other Vector ISAs *there are no explicit
1960 reduction opcodes* in SVP64: Schedules however achieve the same effect.
1961
1962 Due to these conventions only reduction on operations such as `crand`
1963 and `cror` are meaningful because these have Condition Register Fields
1964 as both input and output. Meaningless operations are not prohibited
1965 because the cost in hardware of doing so is prohibitive, but neither
1966 are they `UNDEFINED`. Implementations are still required to execute them
1967 but are at liberty to optimise out any operations that would ultimately
1968 be overwritten, as long as Strict Program Order is still obvservable by
1969 the programmer.
1970
1971 Also bear in mind that 'Reverse Gear' may be enabled, which can be
1972 used in combination with overlapping CR operations to iteratively
1973 accumulate results. Issuing a `sv.crand` operation for example with
1974 `BA` differing from `BB` by one Condition Register Field would result
1975 in a cascade effect, where the first-encountered CR Field would set the
1976 result to zero, and also all subsequent CR Field elements thereafter:
1977
1978 ```
1979 # sv.crand/mr/rg CR4.ge.v, CR5.ge.v, CR4.ge.v
1980 for i in VL-1 downto 0 # reverse gear
1981 CR.field[4+i].ge &= CR.field[5+i].ge
1982 ```
1983
1984 `sv.crxor` with reduction would be particularly useful for parity
1985 calculation for example, although there are many ways in which the same
1986 calculation could be carried out after transferring a vector of CR Fields
1987 to a GPR using crweird operations.
1988
1989 Implementations are free and clear to optimise these reductions in any way
1990 they see fit, as long as the end-result is compatible with Strict Program
1991 Order being observed, and Interrupt latency is not adversely impacted.
1992
1993 ## Unusual and quirky CR operations
1994
1995 **cmp and other compare ops**
1996
1997 `cmp` and `cmpi` etc take GPRs as sources and create a CR Field as a result.
1998
1999 cmpli BF,L,RA,UI
2000 cmpeqb BF,RA,RB
2001
2002 With `ELWIDTH` applying to the source GPR operands this is perfectly fine.
2003
2004 **crweird operations**
2005
2006 There are 4 weird CR-GPR operations and one reasonable one in
2007 the [[cr_int_predication]] set:
2008
2009 * crrweird
2010 * mtcrweird
2011 * crweirder
2012 * crweird
2013 * mcrfm - reasonably normal and referring to CR Fields for src and dest.
2014
2015 The "weird" operations have a non-standard behaviour, being able to
2016 treat *individual bits* of a GPR effectively as elements. They are
2017 expected to be Micro-coded by most Hardware implementations.
2018
2019
2020 --------
2021
2022 \newpage{}
2023
2024 # SVP64 Branch Conditional behaviour
2025
2026 Please note: although similar, SVP64 Branch instructions should be
2027 considered completely separate and distinct from standard scalar
2028 OpenPOWER-approved v3.0B branches. **v3.0B branches are in no way
2029 impacted, altered, changed or modified in any way, shape or form by the
2030 SVP64 Vectorised Variants**.
2031
2032 It is also extremely important to note that Branches are the sole
2033 pseudo-exception in SVP64 to `Scalar Identity Behaviour`. SVP64 Branches
2034 contain additional modes that are useful for scalar operations (i.e. even
2035 when VL=1 or when using single-bit predication).
2036
2037 **Rationale**
2038
2039 Scalar 3.0B Branch Conditional operations, `bc`, `bctar` etc. test
2040 a Condition Register. However for parallel processing it is simply
2041 impossible to perform multiple independent branches: the Program
2042 Counter simply cannot branch to multiple destinations based on multiple
2043 conditions. The best that can be done is to test multiple Conditions
2044 and make a decision of a *single* branch, based on analysis of a *Vector*
2045 of CR Fields which have just been calculated from a *Vector* of results.
2046
2047 In 3D Shader binaries, which are inherently parallelised and predicated,
2048 testing all or some results and branching based on multiple tests is
2049 extremely common, and a fundamental part of Shader Compilers. Example:
2050 without such multi-condition test-and-branch, if a predicate mask is
2051 all zeros a large batch of instructions may be masked out to `nop`,
2052 and it would waste CPU cycles to run them. 3D GPU ISAs can test for
2053 this scenario and, with the appropriate predicate-analysis instruction,
2054 jump over fully-masked-out operations, by spotting that *all* Conditions
2055 are false.
2056
2057 Unless Branches are aware and capable of such analysis, additional
2058 instructions would be required which perform Horizontal Cumulative
2059 analysis of Vectorised Condition Register Fields, in order to reduce
2060 the Vector of CR Fields down to one single yes or no decision that a
2061 Scalar-only v3.0B Branch-Conditional could cope with. Such instructions
2062 would be unavoidable, required, and costly by comparison to a single
2063 Vector-aware Branch. Therefore, in order to be commercially competitive,
2064 `sv.bc` and other Vector-aware Branch Conditional instructions are a
2065 high priority for 3D GPU (and OpenCL-style) workloads.
2066
2067 Given that Power ISA v3.0B is already quite powerful, particularly
2068 the Condition Registers and their interaction with Branches, there are
2069 opportunities to create extremely flexible and compact Vectorised Branch
2070 behaviour. In addition, the side-effects (updating of CTR, truncation
2071 of VL, described below) make it a useful instruction even if the branch
2072 points to the next instruction (no actual branch).
2073
2074 ## Overview
2075
2076 When considering an "array" of branch-tests, there are four
2077 primarily-useful modes: AND, OR, NAND and NOR of all Conditions.
2078 NAND and NOR may be synthesised from AND and OR by inverting `BO[1]`
2079 which just leaves two modes:
2080
2081 * Branch takes place on the **first** CR Field test to succeed
2082 (a Great Big OR of all condition tests). Exit occurs
2083 on the first **successful** test.
2084 * Branch takes place only if **all** CR field tests succeed:
2085 a Great Big AND of all condition tests. Exit occurs
2086 on the first **failed** test.
2087
2088 Early-exit is enacted such that the Vectorised Branch does not
2089 perform needless extra tests, which will help reduce reads on
2090 the Condition Register file.
2091
2092 *Note: Early-exit is **MANDATORY** (required) behaviour. Branches
2093 **MUST** exit at the first sequentially-encountered failure point,
2094 for exactly the same reasons for which it is mandatory in programming
2095 languages doing early-exit: to avoid damaging side-effects and to provide
2096 deterministic behaviour. Speculative testing of Condition Register
2097 Fields is permitted, as is speculative calculation of CTR, as long as,
2098 as usual in any Out-of-Order microarchitecture, that speculative testing
2099 is cancelled should an early-exit occur. i.e. the speculation must be
2100 "precise": Program Order must be preserved*
2101
2102 Also note that when early-exit occurs in Horizontal-first Mode, srcstep,
2103 dststep etc. are all reset, ready to begin looping from the beginning
2104 for the next instruction. However for Vertical-first Mode srcstep
2105 etc. are incremented "as usual" i.e. an early-exit has no special impact,
2106 regardless of whether the branch occurred or not. This can leave srcstep
2107 etc. in what may be considered an unusual state on exit from a loop and
2108 it is up to the programmer to reset srcstep, dststep etc. to known-good
2109 values *(easily achieved with `setvl`)*.
2110
2111 Additional useful behaviour involves two primary Modes (both of which
2112 may be enabled and combined):
2113
2114 * **VLSET Mode**: identical to Data-Dependent Fail-First Mode
2115 for Arithmetic SVP64 operations, with more
2116 flexibility and a close interaction and integration into the
2117 underlying base Scalar v3.0B Branch instruction.
2118 Truncation of VL takes place around the early-exit point.
2119 * **CTR-test Mode**: gives much more flexibility over when and why
2120 CTR is decremented, including options to decrement if a Condition
2121 test succeeds *or if it fails*.
2122
2123 With these side-effects, basic Boolean Logic Analysis advises that it
2124 is important to provide a means to enact them each based on whether
2125 testing succeeds *or fails*. This results in a not-insignificant number
2126 of additional Mode Augmentation bits, accompanying VLSET and CTR-test
2127 Modes respectively.
2128
2129 Predicate skipping or zeroing may, as usual with SVP64, be controlled by
2130 `sz`. Where the predicate is masked out and zeroing is enabled, then in
2131 such circumstances the same Boolean Logic Analysis dictates that rather
2132 than testing only against zero, the option to test against one is also
2133 prudent. This introduces a new immediate field, `SNZ`, which works in
2134 conjunction with `sz`.
2135
2136 Vectorised Branches can be used in either SVP64 Horizontal-First or
2137 Vertical-First Mode. Essentially, at an element level, the behaviour
2138 is identical in both Modes, although the `ALL` bit is meaningless in
2139 Vertical-First Mode.
2140
2141 It is also important to bear in mind that, fundamentally, Vectorised
2142 Branch-Conditional is still extremely close to the Scalar v3.0B
2143 Branch-Conditional instructions, and that the same v3.0B Scalar
2144 Branch-Conditional instructions are still *completely separate and
2145 independent*, being unaltered and unaffected by their SVP64 variants in
2146 every conceivable way.
2147
2148 *Programming note: One important point is that SVP64 instructions are
2149 64 bit. (8 bytes not 4). This needs to be taken into consideration
2150 when computing branch offsets: the offset is relative to the start of
2151 the instruction, which **includes** the SVP64 Prefix*
2152
2153 ## Format and fields
2154
2155 With element-width overrides being meaningless for Condition Register
2156 Fields, bits 4 thru 7 of SVP64 RM may be used for additional Mode bits.
2157
2158 SVP64 RM `MODE` (includes repurposing `ELWIDTH` bits 4:5, and
2159 `ELWIDTH_SRC` bits 6-7 for *alternate* uses) for Branch Conditional:
2160
2161 | 4 | 5 | 6 | 7 | 17 | 18 | 19 | 20 | 21 | 22 23 | description |
2162 | - | - | - | - | -- | -- | -- | -- | --- |--------|----------------- |
2163 |ALL|SNZ| / | / | SL |SLu | 0 | 0 | / | LRu sz | simple mode |
2164 |ALL|SNZ| / |VSb| SL |SLu | 0 | 1 | VLI | LRu sz | VLSET mode |
2165 |ALL|SNZ|CTi| / | SL |SLu | 1 | 0 | / | LRu sz | CTR-test mode |
2166 |ALL|SNZ|CTi|VSb| SL |SLu | 1 | 1 | VLI | LRu sz | CTR-test+VLSET mode |
2167
2168 Brief description of fields:
2169
2170 * **sz=1** if predication is enabled and `sz=1` and a predicate
2171 element bit is zero, `SNZ` will
2172 be substituted in place of the CR bit selected by `BI`,
2173 as the Condition tested.
2174 Contrast this with
2175 normal SVP64 `sz=1` behaviour, where *only* a zero is put in
2176 place of masked-out predicate bits.
2177 * **sz=0** When `sz=0` skipping occurs as usual on
2178 masked-out elements, but unlike all
2179 other SVP64 behaviour which entirely skips an element with
2180 no related side-effects at all, there are certain
2181 special circumstances where CTR
2182 may be decremented. See CTR-test Mode, below.
2183 * **ALL** when set, all branch conditional tests must pass in order for
2184 the branch to succeed. When clear, it is the first sequentially
2185 encountered successful test that causes the branch to succeed.
2186 This is identical behaviour to how programming languages perform
2187 early-exit on Boolean Logic chains.
2188 * **VLI** VLSET is identical to Data-dependent Fail-First mode.
2189 In VLSET mode, VL *may* (depending on `VSb`) be truncated.
2190 If VLI (Vector Length Inclusive) is clear,
2191 VL is truncated to *exclude* the current element, otherwise it is
2192 included. SVSTATE.MVL is not altered: only VL.
2193 * **SL** identical to `LR` except applicable to SVSTATE. If `SL`
2194 is set, SVSTATE is transferred to SVLR (conditionally on
2195 whether `SLu` is set).
2196 * **SLu**: SVSTATE Link Update, like `LRu` except applies to SVSTATE.
2197 * **LRu**: Link Register Update, used in conjunction with LK=1
2198 to make LR update conditional
2199 * **VSb** In VLSET Mode, after testing,
2200 if VSb is set, VL is truncated if the test succeeds. If VSb is clear,
2201 VL is truncated if a test *fails*. Masked-out (skipped)
2202 bits are not considered
2203 part of testing when `sz=0`
2204 * **CTi** CTR inversion. CTR-test Mode normally decrements per element
2205 tested. CTR inversion decrements if a test *fails*. Only relevant
2206 in CTR-test Mode.
2207
2208 LRu and CTR-test modes are where SVP64 Branches subtly differ from
2209 Scalar v3.0B Branches. `sv.bcl` for example will always update LR,
2210 whereas `sv.bcl/lru` will only update LR if the branch succeeds.
2211
2212 Of special interest is that when using ALL Mode (Great Big AND of all
2213 Condition Tests), if `VL=0`, which is rare but can occur in Data-Dependent
2214 Modes, the Branch will always take place because there will be no failing
2215 Condition Tests to prevent it. Likewise when not using ALL Mode (Great
2216 Big OR of all Condition Tests) and `VL=0` the Branch is guaranteed not
2217 to occur because there will be no *successful* Condition Tests to make
2218 it happen.
2219
2220 ## Vectorised CR Field numbering, and Scalar behaviour
2221
2222 It is important to keep in mind that just like all SVP64 instructions,
2223 the `BI` field of the base v3.0B Branch Conditional instruction may be
2224 extended by SVP64 EXTRA augmentation, as well as be marked as either
2225 Scalar or Vector. It is also crucially important to keep in mind that for
2226 CRs, SVP64 sequentially increments the CR *Field* numbers. CR *Fields*
2227 are treated as elements, not bit-numbers of the CR *register*.
2228
2229 The `BI` operand of Branch Conditional operations is five bits, in scalar
2230 v3.0B this would select one bit of the 32 bit CR, comprising eight CR
2231 Fields of 4 bits each. In SVP64 there are 16 32 bit CRs, containing
2232 128 4-bit CR Fields. Therefore, the 2 LSBs of `BI` select the bit from
2233 the CR Field (EQ LT GT SO), and the top 3 bits are extended to either
2234 scalar or vector and to select CR Fields 0..127 as specified in SVP64
2235 [[sv/svp64/appendix]].
2236
2237 When the CR Fields selected by SVP64-Augmented `BI` is marked as scalar,
2238 then as the usual SVP64 rules apply: the Vector loop ends at the first
2239 element tested (the first CR *Field*), after taking predication into
2240 consideration. Thus, also as usual, when a predicate mask is given, and
2241 `BI` marked as scalar, and `sz` is zero, srcstep skips forward to the
2242 first non-zero predicated element, and only that one element is tested.
2243
2244 In other words, the fact that this is a Branch Operation (instead of an
2245 arithmetic one) does not result, ultimately, in significant changes as
2246 to how SVP64 is fundamentally applied, except with respect to:
2247
2248 * the unique properties associated with conditionally
2249 changing the Program Counter (aka "a Branch"), resulting in early-out
2250 opportunities
2251 * CTR-testing
2252
2253 Both are outlined below, in later sections.
2254
2255 ## Horizontal-First and Vertical-First Modes
2256
2257 In SVP64 Horizontal-First Mode, the first failure in ALL mode (Great Big
2258 AND) results in early exit: no more updates to CTR occur (if requested);
2259 no branch occurs, and LR is not updated (if requested). Likewise for
2260 non-ALL mode (Great Big Or) on first success early exit also occurs,
2261 however this time with the Branch proceeding. In both cases the testing
2262 of the Vector of CRs should be done in linear sequential order (or in
2263 REMAP re-sequenced order): such that tests that are sequentially beyond
2264 the exit point are *not* carried out. (*Note: it is standard practice
2265 in Programming languages to exit early from conditional tests, however a
2266 little unusual to consider in an ISA that is designed for Parallel Vector
2267 Processing. The reason is to have strictly-defined guaranteed behaviour*)
2268
2269 In Vertical-First Mode, setting the `ALL` bit results in `UNDEFINED`
2270 behaviour. Given that only one element is being tested at a time in
2271 Vertical-First Mode, a test designed to be done on multiple bits is
2272 meaningless.
2273
2274 ## Description and Modes
2275
2276 Predication in both INT and CR modes may be applied to `sv.bc` and other
2277 SVP64 Branch Conditional operations, exactly as they may be applied to
2278 other SVP64 operations. When `sz` is zero, any masked-out Branch-element
2279 operations are not included in condition testing, exactly like all other
2280 SVP64 operations, *including* side-effects such as potentially updating
2281 LR or CTR, which will also be skipped. There is *one* exception here,
2282 which is when `BO[2]=0, sz=0, CTR-test=0, CTi=1` and the relevant element
2283 predicate mask bit is also zero: under these special circumstances CTR
2284 will also decrement.
2285
2286 When `sz` is non-zero, this normally requests insertion of a zero in
2287 place of the input data, when the relevant predicate mask bit is zero.
2288 This would mean that a zero is inserted in place of `CR[BI+32]` for
2289 testing against `BO`, which may not be desirable in all circumstances.
2290 Therefore, an extra field is provided `SNZ`, which, if set, will insert
2291 a **one** in place of a masked-out element, instead of a zero.
2292
2293 (*Note: Both options are provided because it is useful to deliberately
2294 cause the Branch-Conditional Vector testing to fail at a specific point,
2295 controlled by the Predicate mask. This is particularly useful in `VLSET`
2296 mode, which will truncate SVSTATE.VL at the point of the first failed
2297 test.*)
2298
2299 Normally, CTR mode will decrement once per Condition Test, resulting under
2300 normal circumstances that CTR reduces by up to VL in Horizontal-First
2301 Mode. Just as when v3.0B Branch-Conditional saves at least one instruction
2302 on tight inner loops through auto-decrementation of CTR, likewise it
2303 is also possible to save instruction count for SVP64 loops in both
2304 Vertical-First and Horizontal-First Mode, particularly in circumstances
2305 where there is conditional interaction between the element computation
2306 and testing, and the continuation (or otherwise) of a given loop. The
2307 potential combinations of interactions is why CTR testing options have
2308 been added.
2309
2310 Also, the unconditional bit `BO[0]` is still relevant when Predication
2311 is applied to the Branch because in `ALL` mode all nonmasked bits have
2312 to be tested, and when `sz=0` skipping occurs. Even when VLSET mode is
2313 not used, CTR may still be decremented by the total number of nonmasked
2314 elements, acting in effect as either a popcount or cntlz depending
2315 on which mode bits are set. In short, Vectorised Branch becomes an
2316 extremely powerful tool.
2317
2318 **Micro-Architectural Implementation Note**: *when implemented on top
2319 of a Multi-Issue Out-of-Order Engine it is possible to pass a copy of
2320 the predicate and the prerequisite CR Fields to all Branch Units, as
2321 well as the current value of CTR at the time of multi-issue, and for
2322 each Branch Unit to compute how many times CTR would be subtracted,
2323 in a fully-deterministic and parallel fashion. A SIMD-based Branch
2324 Unit, receiving and processing multiple CR Fields covered by multiple
2325 predicate bits, would do the exact same thing. Obviously, however, if
2326 CTR is modified within any given loop (mtctr) the behaviour of CTR is
2327 no longer deterministic.*
2328
2329 ### Link Register Update
2330
2331 For a Scalar Branch, unconditional updating of the Link Register LR
2332 is useful and practical. However, if a loop of CR Fields is tested,
2333 unconditional updating of LR becomes problematic.
2334
2335 For example when using `bclr` with `LRu=1,LK=0` in Horizontal-First Mode,
2336 LR's value will be unconditionally overwritten after the first element,
2337 such that for execution (testing) of the second element, LR has the value
2338 `CIA+8`. This is covered in the `bclrl` example, in a later section.
2339
2340 The addition of a LRu bit modifies behaviour in conjunction with LK,
2341 as follows:
2342
2343 * `sv.bc` When LRu=0,LK=0, Link Register is not updated
2344 * `sv.bcl` When LRu=0,LK=1, Link Register is updated unconditionally
2345 * `sv.bcl/lru` When LRu=1,LK=1, Link Register will
2346 only be updated if the Branch Condition fails.
2347 * `sv.bc/lru` When LRu=1,LK=0, Link Register will only be updated if
2348 the Branch Condition succeeds.
2349
2350 This avoids destruction of LR during loops (particularly Vertical-First
2351 ones).
2352
2353 **SVLR and SVSTATE**
2354
2355 For precisely the reasons why `LK=1` was added originally to the Power
2356 ISA, with SVSTATE being a peer of the Program Counter it becomes necessary
2357 to also add an SVLR (SVSTATE Link Register) and corresponding control bits
2358 `SL` and `SLu`.
2359
2360 ### CTR-test
2361
2362 Where a standard Scalar v3.0B branch unconditionally decrements CTR when
2363 `BO[2]` is clear, CTR-test Mode introduces more flexibility which allows
2364 CTR to be used for many more types of Vector loops constructs.
2365
2366 CTR-test mode and CTi interaction is as follows: note that `BO[2]`
2367 is still required to be clear for CTR decrements to be considered,
2368 exactly as is the case in Scalar Power ISA v3.0B
2369
2370 * **CTR-test=0, CTi=0**: CTR decrements on a per-element basis
2371 if `BO[2]` is zero. Masked-out elements when `sz=0` are
2372 skipped (i.e. CTR is *not* decremented when the predicate
2373 bit is zero and `sz=0`).
2374 * **CTR-test=0, CTi=1**: CTR decrements on a per-element basis
2375 if `BO[2]` is zero and a masked-out element is skipped
2376 (`sz=0` and predicate bit is zero). This one special case is the
2377 **opposite** of other combinations, as well as being
2378 completely different from normal SVP64 `sz=0` behaviour)
2379 * **CTR-test=1, CTi=0**: CTR decrements on a per-element basis
2380 if `BO[2]` is zero and the Condition Test succeeds.
2381 Masked-out elements when `sz=0` are skipped (including
2382 not decrementing CTR)
2383 * **CTR-test=1, CTi=1**: CTR decrements on a per-element basis
2384 if `BO[2]` is zero and the Condition Test *fails*.
2385 Masked-out elements when `sz=0` are skipped (including
2386 not decrementing CTR)
2387
2388 `CTR-test=0, CTi=1, sz=0` requires special emphasis because it is the
2389 only time in the entirety of SVP64 that has side-effects when
2390 a predicate mask bit is clear. **All** other SVP64 operations
2391 entirely skip an element when sz=0 and a predicate mask bit is zero.
2392 It is also critical to emphasise that in this unusual mode,
2393 no other side-effects occur: **only** CTR is decremented, i.e. the
2394 rest of the Branch operation is skipped.
2395
2396 ### VLSET Mode
2397
2398 VLSET Mode truncates the Vector Length so that subsequent instructions
2399 operate on a reduced Vector Length. This is similar to Data-dependent
2400 Fail-First and LD/ST Fail-First, where for VLSET the truncation occurs
2401 at the Branch decision-point.
2402
2403 Interestingly, due to the side-effects of `VLSET` mode it is actually
2404 useful to use Branch Conditional even to perform no actual branch
2405 operation, i.e to point to the instruction after the branch. Truncation of
2406 VL would thus conditionally occur yet control flow alteration would not.
2407
2408 `VLSET` mode with Vertical-First is particularly unusual. Vertical-First
2409 is designed to be used for explicit looping, where an explicit call to
2410 `svstep` is required to move both srcstep and dststep on to the next
2411 element, until VL (or other condition) is reached. Vertical-First Looping
2412 is expected (required) to terminate if the end of the Vector, VL, is
2413 reached. If however that loop is terminated early because VL is truncated,
2414 VLSET with Vertical-First becomes meaningless. Resolving this would
2415 require two branches: one Conditional, the other branching unconditionally
2416 to create the loop, where the Conditional one jumps over it.
2417
2418 Therefore, with `VSb`, the option to decide whether truncation should
2419 occur if the branch succeeds *or* if the branch condition fails allows
2420 for the flexibility required. This allows a Vertical-First Branch to
2421 *either* be used as a branch-back (loop) *or* as part of a conditional
2422 exit or function call from *inside* a loop, and for VLSET to be integrated
2423 into both types of decision-making.
2424
2425 In the case of a Vertical-First branch-back (loop), with `VSb=0` the
2426 branch takes place if success conditions are met, but on exit from that
2427 loop (branch condition fails), VL will be truncated. This is extremely
2428 useful.
2429
2430 `VLSET` mode with Horizontal-First when `VSb=0` is still useful, because
2431 it can be used to truncate VL to the first predicated (non-masked-out)
2432 element.
2433
2434 The truncation point for VL, when VLi is clear, must not include skipped
2435 elements that preceded the current element being tested. Example:
2436 `sz=0, VLi=0, predicate mask = 0b110010` and the Condition Register
2437 failure point is at CR Field element 4.
2438
2439 * Testing at element 0 is skipped because its predicate bit is zero
2440 * Testing at element 1 passed
2441 * Testing elements 2 and 3 are skipped because their
2442 respective predicate mask bits are zero
2443 * Testing element 4 fails therefore VL is truncated to **2**
2444 not 4 due to elements 2 and 3 being skipped.
2445
2446 If `sz=1` in the above example *then* VL would have been set to 4 because
2447 in non-zeroing mode the zero'd elements are still effectively part of the
2448 Vector (with their respective elements set to `SNZ`)
2449
2450 If `VLI=1` then VL would be set to 5 regardless of sz, due to being inclusive
2451 of the element actually being tested.
2452
2453 ### VLSET and CTR-test combined
2454
2455 If both CTR-test and VLSET Modes are requested, it is important to
2456 observe the correct order. What occurs depends on whether VLi is enabled,
2457 because VLi affects the length, VL.
2458
2459 If VLi (VL truncate inclusive) is set:
2460
2461 1. compute the test including whether CTR triggers
2462 2. (optionally) decrement CTR
2463 3. (optionally) truncate VL (VSb inverts the decision)
2464 4. decide (based on step 1) whether to terminate looping
2465 (including not executing step 5)
2466 5. decide whether to branch.
2467
2468 If VLi is clear, then when a test fails that element
2469 and any following it
2470 should **not** be considered part of the Vector. Consequently:
2471
2472 1. compute the branch test including whether CTR triggers
2473 2. if the test fails against VSb, truncate VL to the *previous*
2474 element, and terminate looping. No further steps executed.
2475 3. (optionally) decrement CTR
2476 4. decide whether to branch.
2477
2478 ## Boolean Logic combinations
2479
2480 In a Scalar ISA, Branch-Conditional testing even of vector results may be
2481 performed through inversion of tests. NOR of all tests may be performed
2482 by inversion of the scalar condition and branching *out* from the scalar
2483 loop around elements, using scalar operations.
2484
2485 In a parallel (Vector) ISA it is the ISA itself which must perform
2486 the prerequisite logic manipulation. Thus for SVP64 there are an
2487 extraordinary number of nesessary combinations which provide completely
2488 different and useful behaviour. Available options to combine:
2489
2490 * `BO[0]` to make an unconditional branch would seem irrelevant if
2491 it were not for predication and for side-effects (CTR Mode
2492 for example)
2493 * Enabling CTR-test Mode and setting `BO[2]` can still result in the
2494 Branch
2495 taking place, not because the Condition Test itself failed, but
2496 because CTR reached zero **because**, as required by CTR-test mode,
2497 CTR was decremented as a **result** of Condition Tests failing.
2498 * `BO[1]` to select whether the CR bit being tested is zero or nonzero
2499 * `R30` and `~R30` and other predicate mask options including CR and
2500 inverted CR bit testing
2501 * `sz` and `SNZ` to insert either zeros or ones in place of masked-out
2502 predicate bits
2503 * `ALL` or `ANY` behaviour corresponding to `AND` of all tests and
2504 `OR` of all tests, respectively.
2505 * Predicate Mask bits, which combine in effect with the CR being
2506 tested.
2507 * Inversion of Predicate Masks (`~r3` instead of `r3`, or using
2508 `NE` rather than `EQ`) which results in an additional
2509 level of possible ANDing, ORing etc. that would otherwise
2510 need explicit instructions.
2511
2512 The most obviously useful combinations here are to set `BO[1]` to zero
2513 in order to turn `ALL` into Great-Big-NAND and `ANY` into Great-Big-NOR.
2514 Other Mode bits which perform behavioural inversion then have to work
2515 round the fact that the Condition Testing is NOR or NAND. The alternative
2516 to not having additional behavioural inversion (`SNZ`, `VSb`, `CTi`)
2517 would be to have a second (unconditional) branch directly after the first,
2518 which the first branch jumps over. This contrivance is avoided by the
2519 behavioural inversion bits.
2520
2521 ## Pseudocode and examples
2522
2523 Please see the SVP64 appendix regarding CR bit ordering and for
2524 the definition of `CR{n}`
2525
2526 For comparative purposes this is a copy of the v3.0B `bc` pseudocode
2527
2528 ```
2529 if (mode_is_64bit) then M <- 0
2530 else M <- 32
2531 if ¬BO[2] then CTR <- CTR - 1
2532 ctr_ok <- BO[2] | ((CTR[M:63] != 0) ^ BO[3])
2533 cond_ok <- BO[0] | ¬(CR[BI+32] ^ BO[1])
2534 if ctr_ok & cond_ok then
2535 if AA then NIA <-iea EXTS(BD || 0b00)
2536 else NIA <-iea CIA + EXTS(BD || 0b00)
2537 if LK then LR <-iea CIA + 4
2538 ```
2539
2540 Simplified pseudocode including LRu and CTR skipping, which illustrates
2541 clearly that SVP64 Scalar Branches (VL=1) are **not** identical to
2542 v3.0B Scalar Branches. The key areas where differences occur are the
2543 inclusion of predication (which can still be used when VL=1), in when and
2544 why CTR is decremented (CTRtest Mode) and whether LR is updated (which
2545 is unconditional in v3.0B when LK=1, and conditional in SVP64 when LRu=1).
2546
2547 Inline comments highlight the fact that the Scalar Branch behaviour and
2548 pseudocode is still clearly visible and embedded within the Vectorised
2549 variant:
2550
2551 ```
2552 if (mode_is_64bit) then M <- 0
2553 else M <- 32
2554 # the bit of CR to test, if the predicate bit is zero,
2555 # is overridden
2556 testbit = CR[BI+32]
2557 if ¬predicate_bit then testbit = SVRMmode.SNZ
2558 # otherwise apart from the override ctr_ok and cond_ok
2559 # are exactly the same
2560 ctr_ok <- BO[2] | ((CTR[M:63] != 0) ^ BO[3])
2561 cond_ok <- BO[0] | ¬(testbit ^ BO[1])
2562 if ¬predicate_bit & ¬SVRMmode.sz then
2563 # this is entirely new: CTR-test mode still decrements CTR
2564 # even when predicate-bits are zero
2565 if ¬BO[2] & CTRtest & ¬CTi then
2566 CTR = CTR - 1
2567 # instruction finishes here
2568 else
2569 # usual BO[2] CTR-mode now under CTR-test mode as well
2570 if ¬BO[2] & ¬(CTRtest & (cond_ok ^ CTi)) then CTR <- CTR - 1
2571 # new VLset mode, conditional test truncates VL
2572 if VLSET and VSb = (cond_ok & ctr_ok) then
2573 if SVRMmode.VLI then SVSTATE.VL = srcstep+1
2574 else SVSTATE.VL = srcstep
2575 # usual LR is now conditional, but also joined by SVLR
2576 lr_ok <- LK
2577 svlr_ok <- SVRMmode.SL
2578 if ctr_ok & cond_ok then
2579 if AA then NIA <-iea EXTS(BD || 0b00)
2580 else NIA <-iea CIA + EXTS(BD || 0b00)
2581 if SVRMmode.LRu then lr_ok <- ¬lr_ok
2582 if SVRMmode.SLu then svlr_ok <- ¬svlr_ok
2583 if lr_ok then LR <-iea CIA + 4
2584 if svlr_ok then SVLR <- SVSTATE
2585 ```
2586
2587 Below is the pseudocode for SVP64 Branches, which is a little less
2588 obvious but identical to the above. The lack of obviousness is down to
2589 the early-exit opportunities.
2590
2591 Effective pseudocode for Horizontal-First Mode:
2592
2593 ```
2594 if (mode_is_64bit) then M <- 0
2595 else M <- 32
2596 cond_ok = not SVRMmode.ALL
2597 for srcstep in range(VL):
2598 # select predicate bit or zero/one
2599 if predicate[srcstep]:
2600 # get SVP64 extended CR field 0..127
2601 SVCRf = SVP64EXTRA(BI>>2)
2602 CRbits = CR{SVCRf}
2603 testbit = CRbits[BI & 0b11]
2604 # testbit = CR[BI+32+srcstep*4]
2605 else if not SVRMmode.sz:
2606 # inverted CTR test skip mode
2607 if ¬BO[2] & CTRtest & ¬CTI then
2608 CTR = CTR - 1
2609 continue # skip to next element
2610 else
2611 testbit = SVRMmode.SNZ
2612 # actual element test here
2613 ctr_ok <- BO[2] | ((CTR[M:63] != 0) ^ BO[3])
2614 el_cond_ok <- BO[0] | ¬(testbit ^ BO[1])
2615 # check if CTR dec should occur
2616 ctrdec = ¬BO[2]
2617 if CTRtest & (el_cond_ok ^ CTi) then
2618 ctrdec = 0b0
2619 if ctrdec then CTR <- CTR - 1
2620 # merge in the test
2621 if SVRMmode.ALL:
2622 cond_ok &= (el_cond_ok & ctr_ok)
2623 else
2624 cond_ok |= (el_cond_ok & ctr_ok)
2625 # test for VL to be set (and exit)
2626 if VLSET and VSb = (el_cond_ok & ctr_ok) then
2627 if SVRMmode.VLI then SVSTATE.VL = srcstep+1
2628 else SVSTATE.VL = srcstep
2629 break
2630 # early exit?
2631 if SVRMmode.ALL != (el_cond_ok & ctr_ok):
2632 break
2633 # SVP64 rules about Scalar registers still apply!
2634 if SVCRf.scalar:
2635 break
2636 # loop finally done, now test if branch (and update LR)
2637 lr_ok <- LK
2638 svlr_ok <- SVRMmode.SL
2639 if cond_ok then
2640 if AA then NIA <-iea EXTS(BD || 0b00)
2641 else NIA <-iea CIA + EXTS(BD || 0b00)
2642 if SVRMmode.LRu then lr_ok <- ¬lr_ok
2643 if SVRMmode.SLu then svlr_ok <- ¬svlr_ok
2644 if lr_ok then LR <-iea CIA + 4
2645 if svlr_ok then SVLR <- SVSTATE
2646 ```
2647
2648 Pseudocode for Vertical-First Mode:
2649
2650 ```
2651 # get SVP64 extended CR field 0..127
2652 SVCRf = SVP64EXTRA(BI>>2)
2653 CRbits = CR{SVCRf}
2654 # select predicate bit or zero/one
2655 if predicate[srcstep]:
2656 if BRc = 1 then # CR0 vectorised
2657 CR{SVCRf+srcstep} = CRbits
2658 testbit = CRbits[BI & 0b11]
2659 else if not SVRMmode.sz:
2660 # inverted CTR test skip mode
2661 if ¬BO[2] & CTRtest & ¬CTI then
2662 CTR = CTR - 1
2663 SVSTATE.srcstep = new_srcstep
2664 exit # no branch testing
2665 else
2666 testbit = SVRMmode.SNZ
2667 # actual element test here
2668 cond_ok <- BO[0] | ¬(testbit ^ BO[1])
2669 # test for VL to be set (and exit)
2670 if VLSET and cond_ok = VSb then
2671 if SVRMmode.VLI
2672 SVSTATE.VL = new_srcstep+1
2673 else
2674 SVSTATE.VL = new_srcstep
2675 ```
2676
2677 ### Example Shader code
2678
2679 ```
2680 // assume f() g() or h() modify a and/or b
2681 while(a > 2) {
2682 if(b < 5)
2683 f();
2684 else
2685 g();
2686 h();
2687 }
2688 ```
2689
2690 which compiles to something like:
2691
2692 ```
2693 vec<i32> a, b;
2694 // ...
2695 pred loop_pred = a > 2;
2696 // loop continues while any of a elements greater than 2
2697 while(loop_pred.any()) {
2698 // vector of predicate bits
2699 pred if_pred = loop_pred & (b < 5);
2700 // only call f() if at least 1 bit set
2701 if(if_pred.any()) {
2702 f(if_pred);
2703 }
2704 label1:
2705 // loop mask ANDs with inverted if-test
2706 pred else_pred = loop_pred & ~if_pred;
2707 // only call g() if at least 1 bit set
2708 if(else_pred.any()) {
2709 g(else_pred);
2710 }
2711 h(loop_pred);
2712 }
2713 ```
2714
2715 which will end up as:
2716
2717 ```
2718 # start from while loop test point
2719 b looptest
2720 while_loop:
2721 sv.cmpi CR80.v, b.v, 5 # vector compare b into CR64 Vector
2722 sv.bc/m=r30/~ALL/sz CR80.v.LT skip_f # skip when none
2723 # only calculate loop_pred & pred_b because needed in f()
2724 sv.crand CR80.v.SO, CR60.v.GT, CR80.V.LT # if = loop & pred_b
2725 f(CR80.v.SO)
2726 skip_f:
2727 # illustrate inversion of pred_b. invert r30, test ALL
2728 # rather than SOME, but masked-out zero test would FAIL,
2729 # therefore masked-out instead is tested against 1 not 0
2730 sv.bc/m=~r30/ALL/SNZ CR80.v.LT skip_g
2731 # else = loop & ~pred_b, need this because used in g()
2732 sv.crternari(A&~B) CR80.v.SO, CR60.v.GT, CR80.V.LT
2733 g(CR80.v.SO)
2734 skip_g:
2735 # conditionally call h(r30) if any loop pred set
2736 sv.bclr/m=r30/~ALL/sz BO[1]=1 h()
2737 looptest:
2738 sv.cmpi CR60.v a.v, 2 # vector compare a into CR60 vector
2739 sv.crweird r30, CR60.GT # transfer GT vector to r30
2740 sv.bc/m=r30/~ALL/sz BO[1]=1 while_loop
2741 end:
2742 ```
2743
2744 ### LRu example
2745
2746 show why LRu would be useful in a loop. Imagine the following
2747 c code:
2748
2749 ```
2750 for (int i = 0; i < 8; i++) {
2751 if (x < y) break;
2752 }
2753 ```
2754
2755 Under these circumstances exiting from the loop is not only based on
2756 CTR it has become conditional on a CR result. Thus it is desirable that
2757 NIA *and* LR only be modified if the conditions are met
2758
2759 v3.0 pseudocode for `bclrl`:
2760
2761 ```
2762 if (mode_is_64bit) then M <- 0
2763 else M <- 32
2764 if ¬BO[2] then CTR <- CTR - 1
2765 ctr_ok <- BO[2] | ((CTR[M:63] != 0) ^ BO[3])
2766 cond_ok <- BO[0] | ¬(CR[BI+32] ^ BO[1])
2767 if ctr_ok & cond_ok then NIA <-iea LR[0:61] || 0b00
2768 if LK then LR <-iea CIA + 4
2769 ```
2770
2771 the latter part for SVP64 `bclrl` becomes:
2772
2773 ```
2774 for i in 0 to VL-1:
2775 ...
2776 ...
2777 cond_ok <- BO[0] | ¬(CR[BI+32] ^ BO[1])
2778 lr_ok <- LK
2779 if ctr_ok & cond_ok then
2780 NIA <-iea LR[0:61] || 0b00
2781 if SVRMmode.LRu then lr_ok <- ¬lr_ok
2782 if lr_ok then LR <-iea CIA + 4
2783 # if NIA modified exit loop
2784 ```
2785
2786 The reason why should be clear from this being a Vector loop:
2787 unconditional destruction of LR when LK=1 makes `sv.bclrl` ineffective,
2788 because the intention going into the loop is that the branch should be to
2789 the copy of LR set at the *start* of the loop, not half way through it.
2790 However if the change to LR only occurs if the branch is taken then it
2791 becomes a useful instruction.
2792
2793 The following pseudocode should **not** be implemented because it
2794 violates the fundamental principle of SVP64 which is that SVP64 looping
2795 is a thin wrapper around Scalar Instructions. The pseducode below is
2796 more an actual Vector ISA Branch and as such is not at all appropriate:
2797
2798 ```
2799 for i in 0 to VL-1:
2800 ...
2801 ...
2802 cond_ok <- BO[0] | ¬(CR[BI+32] ^ BO[1])
2803 if ctr_ok & cond_ok then NIA <-iea LR[0:61] || 0b00
2804 # only at the end of looping is LK checked.
2805 # this completely violates the design principle of SVP64
2806 # and would actually need to be a separate (scalar)
2807 # instruction "set LR to CIA+4 but retrospectively"
2808 # which is clearly impossible
2809 if LK then LR <-iea CIA + 4
2810 ```
2811
2812 [[!tag opf_rfc]]