(no commit message)
[libreriscv.git] / openpower / sv / bitmanip.mdwn
index 7a6bf065a27e0af37b947ba599e69628cec77d31..25cd54eaf75624a47d7be87d06b811627ca5b578 100644 (file)
@@ -1,5 +1,15 @@
 [[!tag standards]]
 
+# bitmanipulation
+
+**DRAFT STATUS**
+
+this extension amalgamates bitnanipulation primitives from many sources, including RISC-V bitmanip, Packed SIMD, AVX-512 and OpenPOWER VSX.  Vectorisation and SIMD are removed: these are straight scalar (element) operations.  Vectorisation Context is provided by [[openpower/sv]].
+
+ternaryv is experimental and is the only operation that may be considered a "Packed SIMD".  It is added as a variant of the already well-justified ternary operation (done in AVX512 as an immediate only) "because it looks fun". As it is based on the LUT4 concept it will allow accelerated emulation of FPGAs.  Other vendors of ISAs are buying FPGA companies to achieve a similar objective.
+
+general-purpose Galois Field operations are added so as to avoid huge opcode proliferation across many areas of Computer Science.  however for convenience and also to avoid setup costs, some of the more common operations (clmul, crc32) are also added.  The expectation is that these operations would all be covered by the same pipeline.
+
 # summary
 
 minor opcode allocation
@@ -12,16 +22,14 @@ minor opcode allocation
     |  011   |Rc| gf*       |
     |  101   |1 | ternaryv  |
     |  101   |0 | ternarycr |
-    |  110   || 1/2-op    |
-    |  111   |Rc| reserved  |
+    |  110   |Rc| 1/2-op    |
+    |  111   |Rc| bitmaski  |
 
 1-op and variants
 
 | dest | src1 | subop | op       |
 | ---- | ---- | ----- | -------- |
 | RT   | RA   | ..    | bmatflip | 
-| RT   | RA   | size  | crc32    | 
-| RT   | RA   | size  | crc32c   | 
 
 2-op and variants
 
@@ -32,6 +40,7 @@ minor opcode allocation
 | RT   | RA   | RB   | bdep  | dep/ext  | 
 | RT   | RA   | RB   | bext  | dep/ext  | 
 | RT   | RA   | RB   |       | grev  |
+| RT   | RA   | RB   |       | clmul*  |
 | RT   | RA   | RB   |       | gorc |  
 | RT   | RA   | RB   | shuf  | shuffle | 
 | RT   | RA   | RB   | unshuf| shuffle | 
@@ -51,10 +60,10 @@ minor opcode allocation
 | -- | -- | --- | --- | ----- | -------- |--| ------ |
 | NN | RT | RA  | RB  | RC    | mode 001 |Rc| ternary |
 | NN | RT | RA  | RB  | im0-4 | im5-7 00 |Rc| ternaryi |
-| NN | RS | RA  | RB  | deg   | 00  011  |Rc| gfmul |
+| NN | RS | RA  | RB  | RC    | 00  011  |Rc| gfmul |
 | NN | RS | RA  | RB  | deg   | 01  011  |Rc| gfadd |
 | NN | RT | RA  | RB  | deg   | 10  011  |Rc| gfinv |
-| NN | RS | RA  | RB  | deg   | 11  011  |Rc| gf rsvd |
+| NN | RS | RA  | RB  | deg   | 11  011  |Rc| gfmuli |
 
 | 0.5|6.10|11.15| 16.23 |24.27 | 28.30 |31| name |
 | -- | -- | --- | ----- | ---- | ----- |--| ------ |
@@ -74,6 +83,9 @@ ops
 | NN | RA | RB  |     |       | 1  | 0100 110 |Rc| rsvd |
 | NN | RA | RB  | sh  | itype | SH | 1000 110 |Rc| bmopsi |
 | NN | RA | RB  |     |       |    | 1100 110 |Rc| rsvd |
+| NN | RA | RB  |     |       |    | 1100 110 |Rc| rsvd |
+| NN | RA | RB  |     |       |    | 1100 110 |Rc| rsvd |
+| NN | RA | RB  |     |       |    | 1100 110 |Rc| rsvd |
 | NN | RA | RB  |     |       | 0  | 0001 110 |Rc| rsvd |
 | NN | RA | RB  |     |       | 0  | 0101 110 |Rc| rsvd |
 | NN | RA | RB  | RC  | 00    | 0  | 0010 110 |Rc| gorc |
@@ -82,6 +94,7 @@ ops
 | NN | RA | RB  | sh  | 00    | 0  | 1110 110 |Rc| gorcwi |
 | NN | RA | RB  | RC  | 00    | 1  | 1110 110 |Rc| bmator  |
 | NN | RA | RB  | RC  | 01    | 0  | 0010 110 |Rc| grev |
+| NN | RA | RB  | RC  | 01    | 1  | 0010 110 |Rc| clmul |
 | NN | RA | RB  | sh  | 01    | SH | 1010 110 |Rc| grevi |
 | NN | RA | RB  | RC  | 01    | 0  | 0110 110 |Rc| grevw |
 | NN | RA | RB  | sh  | 01    | 0  | 1110 110 |Rc| grevwi |
@@ -91,7 +104,8 @@ ops
 | NN | RA | RB  | RC  | 10    | 0  | 0110 110 |Rc| shflw |
 | NN | RA | RB  | RC  | 10    | 0  | 1110 110 |Rc| bdep   |
 | NN | RA | RB  | RC  | 10    | 1  | 1110 110 |Rc| bext  |
-| NN | RA | RB  |     | 11    |    | 1110 110 |Rc| rsvd  |
+| NN | RA | RB  | RC  | 11    | 0  | 1110 110 |Rc| clmulr  |
+| NN | RA | RB  | RC  | 11    | 1  | 1110 110 |Rc| clmulh  |
 | NN | RA | RB  |     |       |    | NN11 110 |Rc| rsvd  |
 
 # bit to byte permute
@@ -154,6 +168,22 @@ signed and unsigned min/max for integer.  this is sort-of partly synthesiseable
 
 signed/unsigned min/max gives more flexibility.
 
+```
+uint_xlen_t min(uint_xlen_t rs1, uint_xlen_t rs2)
+{ return (int_xlen_t)rs1 < (int_xlen_t)rs2 ? rs1 : rs2;
+}
+uint_xlen_t max(uint_xlen_t rs1, uint_xlen_t rs2)
+{ return (int_xlen_t)rs1 > (int_xlen_t)rs2 ? rs1 : rs2;
+}
+uint_xlen_t minu(uint_xlen_t rs1, uint_xlen_t rs2)
+{ return rs1 < rs2 ? rs1 : rs2;
+}
+uint_xlen_t maxu(uint_xlen_t rs1, uint_xlen_t rs2)
+{ return rs1 > rs2 ? rs1 : rs2;
+}
+```
+
+
 # ternary bitops
 
 Similar to FPGA LUTs: for every bit perform a lookup into a table using an 8bit immediate, or in another register
@@ -209,12 +239,17 @@ another mode selection would be CRs not Ints.
 # bitmask set
 
 based on RV bitmanip singlebit set, instruction format similar to shift
-[[isa/fixedshift]].  bmext is actually covered already (shift-with-mask).
-however bitmask-invert is not, and set/clr are not covered, although they can ise the same Shift ALU.
+[[isa/fixedshift]].  bmext is actually covered already (shift-with-mask rldicl but only immediate version).
+however bitmask-invert is not, and set/clr are not covered, although they can use the same Shift ALU.
+
+bmext (RB) version is not the same as rldicl because bmext is a right shift by RC, where rldicl is a left rotate.  for the immediate version this does not matter, so a bmexti is not required.
+bmrev however there is no direct equivalent and consequently a bmrevi is required.
+
+| 0.5|6.10|11.15|16.20|21.25| 26..30  |31| name  |
+| -- | -- | --- | --- | --- | ------- |--| ----- |
+| NN | RT | RA  | RB  | RC  | mode 010 |Rc| bm*   |
+| NN | RT | RA  | RB  | RC  | 0 1  111 |Rc| bmrev |
 
-| 0.5|6.10|11.15|16.20|21.25| 26..30  |31|
-| -- | -- | --- | --- | --- | ------- |--|
-| NN | RT | RA  | RB  | RC  | mode 010 |Rc|
 
 ```
 uint_xlen_t bmset(RA, RB, sh)
@@ -246,6 +281,29 @@ uint_xlen_t bmext(RA, RB, sh)
 }
 ```
 
+bitmask extract with reverse.  can be done by bitinverting all of RA and getting bits of RA from the opposite end.
+
+```
+msb = rb[5:0];
+rev[0:msb] = ra[msb:0];
+rt = ZE(rev[msb:0]);
+
+uint_xlen_t bmextrev(RA, RB, sh)
+{
+    int shamt = (RB & (XLEN - 1));
+    shamt = (XLEN-1)-shamt;  # shift other end
+    bra = bitreverse(RA)     # swap LSB-MSB
+    mask = (2<<sh)-1;
+    return mask & (bra >> shamt);
+}
+```
+
+| 0.5|6.10|11.15|16.20|21.26| 27..30  |31| name   |
+| -- | -- | --- | --- | --- | ------- |--| ------ |
+| NN | RT | RA  | RB  | sh  | 0   111 |Rc| bmrevi |
+
+
+
 # grev
 
 based on RV bitmanip
@@ -445,6 +503,8 @@ uint_xlen_t clmulr(uint_xlen_t RA, uint_xlen_t RB)
 ```
 # Galois Field
 
+see <https://courses.csail.mit.edu/6.857/2016/files/ffield.py>
+
 ## Multiply
 
 this requires 3 parameters and a "degree"
@@ -454,10 +514,12 @@ this requires 3 parameters and a "degree"
 realistically with the degree also needing to be an immediate it should be brought down to an overwrite version:
 
     RS = GFMUL(RS, RA, gfdegree, modulo=RB)
+    RS = GFMUL(RS, RA, gfdegree=RC, modulo=RB)
 
 | 0.5|6.10|11.15|16.20|21.25| 26..30  |31|
 | -- | -- | --- | --- | --- | ------- |--|
 | NN | RS | RA  | RB  | deg | 00  011 |Rc|
+| NN | RS | RA  | RB  | RC  | 11  011 |Rc|
 
 where the SimpleV variant may override RS-as-src differently from RS-as-dest
 
@@ -514,10 +576,14 @@ if __name__ == "__main__":
 ## GF add
 
     RS = GFADD(RS, RA|0, gfdegree, modulo=RB)
+    RS = GFADDI(RS, RA|0, gfdegree=RC, modulo=RB)
 
-| 0.5|6.10|11.15|16.20|21.25| 26..30  |31|
-| -- | -- | --- | --- | --- | ------- |--|
-| NN | RS | RA  | RB  | deg | 01  011 |Rc|
+| 0.5|6.10|11.15|16.20|21.25| 26..30  |31| name  |
+| -- | -- | --- | --- | --- | ------- |--| ----- |
+| NN | RS | RA  | RB  | RC  | 0 1  011 |Rc| gfadd |
+| NN | RS | RA  | RB  | RC  | 1 1  111 |Rc| gfaddi |
+
+GFMOD is a pseudo-op where RA=0
 
 ## gf invert
 
@@ -553,36 +619,6 @@ def gf_invert(a, mod=0x1B) :
   return g1
 ```
 
-# crc
-
-* <https://stackoverflow.com/questions/21171733/calculating-constants-for-crc32-using-pclmulqdq>
-* <https://en.wikipedia.org/wiki/Cyclic_redundancy_check>
-
-```
-uint_xlen_t crc32(uint_xlen_t x, int nbits)
-{
-    for (int i = 0; i < nbits; i++)
-        x = (x >> 1) ^ (0xEDB88320 & ~((x&1)-1));
-    return x;
-}
-uint_xlen_t crc32c(uint_xlen_t x, int nbits)
-{
-    for (int i = 0; i < nbits; i++)
-        x = (x >> 1) ^ (0x82F63B78 & ~((x&1)-1));
-    return x;
-}
-uint_xlen_t crc32_b(uint_xlen_t RA) { return crc32(RA, 8); }
-uint_xlen_t crc32_h(uint_xlen_t RA) { return crc32(RA, 16); }
-uint_xlen_t crc32_w(uint_xlen_t RA) { return crc32(RA, 32); }
-uint_xlen_t crc32c_b(uint_xlen_t RA) { return crc32c(RA, 8); }
-uint_xlen_t crc32c_h(uint_xlen_t RA) { return crc32c(RA, 16); }
-uint_xlen_t crc32c_w(uint_xlen_t RA) { return crc32c(RA, 32); }
-#if XLEN > 32
-uint_xlen_t crc32_d (uint_xlen_t RA) { return crc32 (RA, 64); }
-uint_xlen_t crc32c_d(uint_xlen_t RA) { return crc32c(RA, 64); }
-#endif
-```
-
 # bitmatrix
 
 ```