Update NEON transpose functions.

Unify coding style.

Change-Id: I5826f40c02c882df7353391e0c9dd6cef6bd4b97
This commit is contained in:
Linfeng Zhang
2016-08-31 14:58:40 -07:00
parent 3dfba04dec
commit bee7d837ab

View File

@@ -117,7 +117,7 @@ static INLINE void transpose_u8_16x8(
uint8x8_t *o2, uint8x8_t *o3, uint8x8_t *o4, uint8x8_t *o5, uint8x8_t *o6,
uint8x8_t *o7, uint8x8_t *o8, uint8x8_t *o9, uint8x8_t *o10, uint8x8_t *o11,
uint8x8_t *o12, uint8x8_t *o13, uint8x8_t *o14, uint8x8_t *o15) {
// Input:
// Swap 8 bit elements. Goes from:
// i0: 00 01 02 03 04 05 06 07 08 09 0A 0B 0C 0D 0E 0F
// i1: 10 11 12 13 14 15 16 17 18 19 1A 1B 1C 1D 1E 1F
// i2: 20 21 22 23 24 25 26 27 28 29 2A 2B 2C 2D 2E 2F
@@ -126,56 +126,55 @@ static INLINE void transpose_u8_16x8(
// i5: 50 51 52 53 54 55 56 57 58 59 5A 5B 5C 5D 5E 5F
// i6: 60 61 62 63 64 65 66 67 68 69 6A 6B 6C 6D 6E 6F
// i7: 70 71 72 73 74 75 76 77 78 79 7A 7B 7C 7D 7E 7F
uint8x16x2_t b0, b1, b2, b3;
uint16x8x2_t c0, c1, c2, c3;
uint32x4x2_t d0, d1, d2, d3;
// to:
// b0.val[0]: 00 10 02 12 04 14 06 16 08 18 0A 1A 0C 1C 0E 1E
// b0.val[1]: 01 11 03 13 05 15 07 17 09 19 0B 1B 0D 1D 0F 1F
// b1.val[0]: 20 30 22 32 24 34 26 36 28 38 2A 3A 2C 3C 2E 3E
// b1.val[1]: 21 31 23 33 25 35 27 37 29 39 2B 3B 2D 3D 2F 3F
// b2.val[0]: 40 50 42 52 44 54 46 56 48 58 4A 5A 4C 5C 4E 5E
// b2.val[1]: 41 51 43 53 45 55 47 57 49 59 4B 5B 4D 5D 4F 5F
// b3.val[0]: 60 70 62 72 64 74 66 76 68 78 6A 7A 6C 7C 6E 7E
// b3.val[1]: 61 71 63 73 65 75 67 77 69 79 6B 7B 6D 7D 6F 7F
const uint8x16x2_t b0 = vtrnq_u8(i0, i1);
const uint8x16x2_t b1 = vtrnq_u8(i2, i3);
const uint8x16x2_t b2 = vtrnq_u8(i4, i5);
const uint8x16x2_t b3 = vtrnq_u8(i6, i7);
// b0: 00 10 02 12 04 14 06 16 08 18 0A 1A 0C 1C 0E 1E
// 01 11 03 13 05 15 07 17 09 19 0B 1B 0D 1D 0F 1F
// b1: 20 30 22 32 24 34 26 36 28 38 2A 3A 2C 3C 2E 3E
// 21 31 23 33 25 35 27 37 29 39 2B 3B 2D 3D 2F 3F
// b2: 40 50 42 52 44 54 46 56 48 58 4A 5A 4C 5C 4E 5E
// 41 51 43 53 45 55 47 57 49 59 4B 5B 4D 5D 4F 5F
// b3: 60 70 62 72 64 74 66 76 68 78 6A 7A 6C 7C 6E 7E
// 61 71 63 73 65 75 67 77 69 79 6B 7B 6D 7D 6F 7F
b0 = vtrnq_u8(i0, i1);
b1 = vtrnq_u8(i2, i3);
b2 = vtrnq_u8(i4, i5);
b3 = vtrnq_u8(i6, i7);
// Swap 16 bit elements resulting in:
// c0.val[0]: 00 10 20 30 04 14 24 34 08 18 28 38 0C 1C 2C 3C
// c0.val[1]: 02 12 22 32 06 16 26 36 0A 1A 2A 3A 0E 1E 2E 3E
// c1.val[0]: 01 11 21 31 05 15 25 35 09 19 29 39 0D 1D 2D 3D
// c1.val[1]: 03 13 23 33 07 17 27 37 0B 1B 2B 3B 0F 1F 2F 3F
// c2.val[0]: 40 50 60 70 44 54 64 74 48 58 68 78 4C 5C 6C 7C
// c2.val[1]: 42 52 62 72 46 56 66 76 4A 5A 6A 7A 4E 5E 6E 7E
// c3.val[0]: 41 51 61 71 45 55 65 75 49 59 69 79 4D 5D 6D 7D
// c3.val[1]: 43 53 63 73 47 57 67 77 4B 5B 6B 7B 4F 5F 6F 7F
const uint16x8x2_t c0 = vtrnq_u16(vreinterpretq_u16_u8(b0.val[0]),
vreinterpretq_u16_u8(b1.val[0]));
const uint16x8x2_t c1 = vtrnq_u16(vreinterpretq_u16_u8(b0.val[1]),
vreinterpretq_u16_u8(b1.val[1]));
const uint16x8x2_t c2 = vtrnq_u16(vreinterpretq_u16_u8(b2.val[0]),
vreinterpretq_u16_u8(b3.val[0]));
const uint16x8x2_t c3 = vtrnq_u16(vreinterpretq_u16_u8(b2.val[1]),
vreinterpretq_u16_u8(b3.val[1]));
// c0: 00 10 20 30 04 14 24 34 08 18 28 38 0C 1C 2C 3C
// 02 12 22 32 06 16 26 36 0A 1A 2A 3A 0E 1E 2E 3E
// c1: 01 11 21 31 05 15 25 35 09 19 29 39 0D 1D 2D 3D
// 03 13 23 33 07 17 27 37 0B 1B 2B 3B 0F 1F 2F 3F
// c2: 40 50 60 70 44 54 64 74 48 58 68 78 4C 5C 6C 7C
// 42 52 62 72 46 56 66 76 4A 5A 6A 7A 4E 5E 6E 7E
// c3: 41 51 61 71 45 55 65 75 49 59 69 79 4D 5D 6D 7D
// 43 53 63 73 47 57 67 77 4B 5B 6B 7B 4F 5F 6F 7F
c0 = vtrnq_u16(vreinterpretq_u16_u8(b0.val[0]),
vreinterpretq_u16_u8(b1.val[0]));
c1 = vtrnq_u16(vreinterpretq_u16_u8(b0.val[1]),
vreinterpretq_u16_u8(b1.val[1]));
c2 = vtrnq_u16(vreinterpretq_u16_u8(b2.val[0]),
vreinterpretq_u16_u8(b3.val[0]));
c3 = vtrnq_u16(vreinterpretq_u16_u8(b2.val[1]),
vreinterpretq_u16_u8(b3.val[1]));
// d0: 00 10 20 30 40 50 60 70 08 18 28 38 48 58 68 78
// 04 14 24 34 44 54 64 74 0C 1C 2C 3C 4C 5C 6C 7C
// d1: 02 12 22 32 42 52 62 72 0A 1A 2A 3A 4A 5A 6A 7A
// 06 16 26 36 46 56 66 76 0E 1E 2E 3E 4E 5E 6E 7E
// d2: 01 11 21 31 41 51 61 71 09 19 29 39 49 59 69 79
// 05 15 25 35 45 55 65 75 0D 1D 2D 3D 4D 5D 6D 7D
// d3: 03 13 23 33 43 53 63 73 0B 1B 2B 3B 4B 5B 6B 7B
// 07 17 27 37 47 57 67 77 0F 1F 2F 3F 4F 5F 6F 7F
d0 = vtrnq_u32(vreinterpretq_u32_u16(c0.val[0]),
vreinterpretq_u32_u16(c2.val[0]));
d1 = vtrnq_u32(vreinterpretq_u32_u16(c0.val[1]),
vreinterpretq_u32_u16(c2.val[1]));
d2 = vtrnq_u32(vreinterpretq_u32_u16(c1.val[0]),
vreinterpretq_u32_u16(c3.val[0]));
d3 = vtrnq_u32(vreinterpretq_u32_u16(c1.val[1]),
vreinterpretq_u32_u16(c3.val[1]));
// Swap 32 bit elements resulting in:
// d0.val[0]: 00 10 20 30 40 50 60 70 08 18 28 38 48 58 68 78
// d0.val[1]: 04 14 24 34 44 54 64 74 0C 1C 2C 3C 4C 5C 6C 7C
// d1.val[0]: 02 12 22 32 42 52 62 72 0A 1A 2A 3A 4A 5A 6A 7A
// d1.val[1]: 06 16 26 36 46 56 66 76 0E 1E 2E 3E 4E 5E 6E 7E
// d2.val[0]: 01 11 21 31 41 51 61 71 09 19 29 39 49 59 69 79
// d2.val[1]: 05 15 25 35 45 55 65 75 0D 1D 2D 3D 4D 5D 6D 7D
// d3.val[0]: 03 13 23 33 43 53 63 73 0B 1B 2B 3B 4B 5B 6B 7B
// d3.val[1]: 07 17 27 37 47 57 67 77 0F 1F 2F 3F 4F 5F 6F 7F
const uint32x4x2_t d0 = vtrnq_u32(vreinterpretq_u32_u16(c0.val[0]),
vreinterpretq_u32_u16(c2.val[0]));
const uint32x4x2_t d1 = vtrnq_u32(vreinterpretq_u32_u16(c0.val[1]),
vreinterpretq_u32_u16(c2.val[1]));
const uint32x4x2_t d2 = vtrnq_u32(vreinterpretq_u32_u16(c1.val[0]),
vreinterpretq_u32_u16(c3.val[0]));
const uint32x4x2_t d3 = vtrnq_u32(vreinterpretq_u32_u16(c1.val[1]),
vreinterpretq_u32_u16(c3.val[1]));
// Output:
// o0 : 00 10 20 30 40 50 60 70
@@ -221,7 +220,7 @@ static INLINE void transpose_u8_8x16(
const uint8x8_t i15, uint8x16_t *o0, uint8x16_t *o1, uint8x16_t *o2,
uint8x16_t *o3, uint8x16_t *o4, uint8x16_t *o5, uint8x16_t *o6,
uint8x16_t *o7) {
// Input:
// Combine 8 bit elements. Goes from:
// i0 : 00 01 02 03 04 05 06 07
// i1 : 10 11 12 13 14 15 16 17
// i2 : 20 21 22 23 24 25 26 27
@@ -238,73 +237,73 @@ static INLINE void transpose_u8_8x16(
// i13: D0 D1 D2 D3 D4 D5 D6 D7
// i14: E0 E1 E2 E3 E4 E5 E6 E7
// i15: F0 F1 F2 F3 F4 F5 F6 F7
uint8x16x2_t b0, b1, b2, b3;
uint16x8x2_t c0, c1, c2, c3;
uint32x4x2_t d0, d1, d2, d3;
// to:
// a0: 00 01 02 03 04 05 06 07 80 81 82 83 84 85 86 87
// a1: 10 11 12 13 14 15 16 17 90 91 92 93 94 95 96 97
// a2: 20 21 22 23 24 25 26 27 A0 A1 A2 A3 A4 A5 A6 A7
// a3: 30 31 32 33 34 35 36 37 B0 B1 B2 B3 B4 B5 B6 B7
// a4: 40 41 42 43 44 45 46 47 C0 C1 C2 C3 C4 C5 C6 C7
// a5: 50 51 52 53 54 55 56 57 D0 D1 D2 D3 D4 D5 D6 D7
// a6: 60 61 62 63 64 65 66 67 E0 E1 E2 E3 E4 E5 E6 E7
// a7: 70 71 72 73 74 75 76 77 F0 F1 F2 F3 F4 F5 F6 F7
const uint8x16_t a0 = vcombine_u8(i0, i8);
const uint8x16_t a1 = vcombine_u8(i1, i9);
const uint8x16_t a2 = vcombine_u8(i2, i10);
const uint8x16_t a3 = vcombine_u8(i3, i11);
const uint8x16_t a4 = vcombine_u8(i4, i12);
const uint8x16_t a5 = vcombine_u8(i5, i13);
const uint8x16_t a6 = vcombine_u8(i6, i14);
const uint8x16_t a7 = vcombine_u8(i7, i15);
// b0: 00 01 02 03 04 05 06 07 80 81 82 83 84 85 86 87
// 10 11 12 13 14 15 16 17 90 91 92 93 94 95 96 97
// b1: 20 21 22 23 24 25 26 27 A0 A1 A2 A3 A4 A5 A6 A7
// 30 31 32 33 34 35 36 37 B0 B1 B2 B3 B4 B5 B6 B7
// b2: 40 41 42 43 44 45 46 47 C0 C1 C2 C3 C4 C5 C6 C7
// 50 51 52 53 54 55 56 57 D0 D1 D2 D3 D4 D5 D6 D7
// b3: 60 61 62 63 64 65 66 67 E0 E1 E2 E3 E4 E5 E6 E7
// 70 71 72 73 74 75 76 77 F0 F1 F2 F3 F4 F5 F6 F7
b0.val[0] = vcombine_u8(i0, i8);
b0.val[1] = vcombine_u8(i1, i9);
b1.val[0] = vcombine_u8(i2, i10);
b1.val[1] = vcombine_u8(i3, i11);
b2.val[0] = vcombine_u8(i4, i12);
b2.val[1] = vcombine_u8(i5, i13);
b3.val[0] = vcombine_u8(i6, i14);
b3.val[1] = vcombine_u8(i7, i15);
// Swap 8 bit elements resulting in:
// b0.val[0]: 00 10 02 12 04 14 06 16 80 90 82 92 84 94 86 96
// b0.val[1]: 01 11 03 13 05 15 07 17 81 91 83 93 85 95 87 97
// b1.val[0]: 20 30 22 32 24 34 26 36 A0 B0 A2 B2 A4 B4 A6 B6
// b1.val[1]: 21 31 23 33 25 35 27 37 A1 B1 A3 B3 A5 B5 A7 B7
// b2.val[0]: 40 50 42 52 44 54 46 56 C0 D0 C2 D2 C4 D4 C6 D6
// b2.val[1]: 41 51 43 53 45 55 47 57 C1 D1 C3 D3 C5 D5 C7 D7
// b3.val[0]: 60 70 62 72 64 74 66 76 E0 F0 E2 F2 E4 F4 E6 F6
// b3.val[1]: 61 71 63 73 65 75 67 77 E1 F1 E3 F3 E5 F5 E7 F7
const uint8x16x2_t b0 = vtrnq_u8(a0, a1);
const uint8x16x2_t b1 = vtrnq_u8(a2, a3);
const uint8x16x2_t b2 = vtrnq_u8(a4, a5);
const uint8x16x2_t b3 = vtrnq_u8(a6, a7);
// b0: 00 10 02 12 04 14 06 16 80 90 82 92 84 94 86 96
// 01 11 03 13 05 15 07 17 81 91 83 93 85 95 87 97
// b1: 20 30 22 32 24 34 26 36 A0 B0 A2 B2 A4 B4 A6 B6
// 21 31 23 33 25 35 27 37 A1 B1 A3 B3 A5 B5 A7 B7
// b2: 40 50 42 52 44 54 46 56 C0 D0 C2 D2 C4 D4 C6 D6
// 41 51 43 53 45 55 47 57 C1 D1 C3 D3 C5 D5 C7 D7
// b3: 60 70 62 72 64 74 66 76 E0 F0 E2 F2 E4 F4 E6 F6
// 61 71 63 73 65 75 67 77 E1 F1 E3 F3 E5 F5 E7 F7
b0 = vtrnq_u8(b0.val[0], b0.val[1]);
b1 = vtrnq_u8(b1.val[0], b1.val[1]);
b2 = vtrnq_u8(b2.val[0], b2.val[1]);
b3 = vtrnq_u8(b3.val[0], b3.val[1]);
// Swap 16 bit elements resulting in:
// c0.val[0]: 00 10 20 30 04 14 24 34 80 90 A0 B0 84 94 A4 B4
// c0.val[1]: 02 12 22 32 06 16 26 36 82 92 A2 B2 86 96 A6 B6
// c1.val[0]: 01 11 21 31 05 15 25 35 81 91 A1 B1 85 95 A5 B5
// c1.val[1]: 03 13 23 33 07 17 27 37 83 93 A3 B3 87 97 A7 B7
// c2.val[0]: 40 50 60 70 44 54 64 74 C0 D0 E0 F0 C4 D4 E4 F4
// c2.val[1]: 42 52 62 72 46 56 66 76 C2 D2 E2 F2 C6 D6 E6 F6
// c3.val[0]: 41 51 61 71 45 55 65 75 C1 D1 E1 F1 C5 D5 E5 F5
// c3.val[1]: 43 53 63 73 47 57 67 77 C3 D3 E3 F3 C7 D7 E7 F7
const uint16x8x2_t c0 = vtrnq_u16(vreinterpretq_u16_u8(b0.val[0]),
vreinterpretq_u16_u8(b1.val[0]));
const uint16x8x2_t c1 = vtrnq_u16(vreinterpretq_u16_u8(b0.val[1]),
vreinterpretq_u16_u8(b1.val[1]));
const uint16x8x2_t c2 = vtrnq_u16(vreinterpretq_u16_u8(b2.val[0]),
vreinterpretq_u16_u8(b3.val[0]));
const uint16x8x2_t c3 = vtrnq_u16(vreinterpretq_u16_u8(b2.val[1]),
vreinterpretq_u16_u8(b3.val[1]));
// c0: 00 10 20 30 04 14 24 34 80 90 A0 B0 84 94 A4 B4
// 02 12 22 32 06 16 26 36 82 92 A2 B2 86 96 A6 B6
// c1: 01 11 21 31 05 15 25 35 81 91 A1 B1 85 95 A5 B5
// 03 13 23 33 07 17 27 37 83 93 A3 B3 87 97 A7 B7
// c2: 40 50 60 70 44 54 64 74 C0 D0 E0 F0 C4 D4 E4 F4
// 42 52 62 72 46 56 66 76 C2 D2 E2 F2 C6 D6 E6 F6
// c3: 41 51 61 71 45 55 65 75 C1 D1 E1 F1 C5 D5 E5 F5
// 43 53 63 73 47 57 67 77 C3 D3 E3 F3 C7 D7 E7 F7
c0 = vtrnq_u16(vreinterpretq_u16_u8(b0.val[0]),
vreinterpretq_u16_u8(b1.val[0]));
c1 = vtrnq_u16(vreinterpretq_u16_u8(b0.val[1]),
vreinterpretq_u16_u8(b1.val[1]));
c2 = vtrnq_u16(vreinterpretq_u16_u8(b2.val[0]),
vreinterpretq_u16_u8(b3.val[0]));
c3 = vtrnq_u16(vreinterpretq_u16_u8(b2.val[1]),
vreinterpretq_u16_u8(b3.val[1]));
// d0: 00 10 20 30 40 50 60 70 80 90 A0 B0 C0 D0 E0 F0
// 04 14 24 34 44 54 64 74 84 94 A4 B4 C4 D4 E4 F4
// d1: 02 12 22 32 42 52 62 72 82 92 A2 B2 C2 D2 E2 F2
// 06 16 26 36 46 56 66 76 86 96 A6 B6 C6 D6 E6 F6
// d2: 01 11 21 31 41 51 61 71 81 91 A1 B1 C1 D1 E1 F1
// 05 15 25 35 45 55 65 75 85 95 A5 B5 C5 D5 E5 F5
// d3: 03 13 23 33 43 53 63 73 83 93 A3 B3 C3 D3 E3 F3
// 07 17 27 37 47 57 67 77 87 97 A7 B7 C7 D7 E7 F7
d0 = vtrnq_u32(vreinterpretq_u32_u16(c0.val[0]),
vreinterpretq_u32_u16(c2.val[0]));
d1 = vtrnq_u32(vreinterpretq_u32_u16(c0.val[1]),
vreinterpretq_u32_u16(c2.val[1]));
d2 = vtrnq_u32(vreinterpretq_u32_u16(c1.val[0]),
vreinterpretq_u32_u16(c3.val[0]));
d3 = vtrnq_u32(vreinterpretq_u32_u16(c1.val[1]),
vreinterpretq_u32_u16(c3.val[1]));
// Swap 32 bit elements resulting in:
// d0.val[0]: 00 10 20 30 40 50 60 70 80 90 A0 B0 C0 D0 E0 F0
// d0.val[1]: 04 14 24 34 44 54 64 74 84 94 A4 B4 C4 D4 E4 F4
// d1.val[0]: 02 12 22 32 42 52 62 72 82 92 A2 B2 C2 D2 E2 F2
// d1.val[1]: 06 16 26 36 46 56 66 76 86 96 A6 B6 C6 D6 E6 F6
// d2.val[0]: 01 11 21 31 41 51 61 71 81 91 A1 B1 C1 D1 E1 F1
// d2.val[1]: 05 15 25 35 45 55 65 75 85 95 A5 B5 C5 D5 E5 F5
// d3.val[0]: 03 13 23 33 43 53 63 73 83 93 A3 B3 C3 D3 E3 F3
// d3.val[1]: 07 17 27 37 47 57 67 77 87 97 A7 B7 C7 D7 E7 F7
const uint32x4x2_t d0 = vtrnq_u32(vreinterpretq_u32_u16(c0.val[0]),
vreinterpretq_u32_u16(c2.val[0]));
const uint32x4x2_t d1 = vtrnq_u32(vreinterpretq_u32_u16(c0.val[1]),
vreinterpretq_u32_u16(c2.val[1]));
const uint32x4x2_t d2 = vtrnq_u32(vreinterpretq_u32_u16(c1.val[0]),
vreinterpretq_u32_u16(c3.val[0]));
const uint32x4x2_t d3 = vtrnq_u32(vreinterpretq_u32_u16(c1.val[1]),
vreinterpretq_u32_u16(c3.val[1]));
// Output:
// o0: 00 10 20 30 40 50 60 70 80 90 A0 B0 C0 D0 E0 F0