arne
/
cryptopp
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Switch to vec_mergeh and vec_mergel for PowerPC in binary curves

pull/795/head
Jeffrey Walton 2019-01-19 04:22:55 -05:00
parent ac49c970bc
commit e3bc76832c
No known key found for this signature in database
GPG Key ID: B36AB348921B1838
1 changed files with 5 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
gf2n_simd.cpp
View File

@ -361,26 +361,20 @@ inline uint64x2_p VMULL_11LE(const uint64x2_p& a, const uint64x2_p& b)
inline void inline void
F2N_Multiply_128x128_POWER8(uint64x2_p& c1, uint64x2_p& c0, const uint64x2_p& a, const uint64x2_p& b) F2N_Multiply_128x128_POWER8(uint64x2_p& c1, uint64x2_p& c0, const uint64x2_p& a, const uint64x2_p& b)
{ {
const uint8_t mb1[] = {8,9,10,11, 12,13,14,15, 8,9,10,11, 12,13,14,15};
const uint8_t mb2[] = {8,9,10,11, 12,13,14,15, 16,17,18,19, 20,21,22,23};
const uint8x16_p m1 = (uint8x16_p)VecLoad(mb1);
const uint8x16_p m2 = (uint8x16_p)VecLoad(mb2);
uint64x2_p t1, t2, z0={0}; uint64x2_p t1, t2, z0={0};
c0 = VMULL_00LE(a, b); c0 = VMULL_00LE(a, b);
c1 = VMULL_11LE(a, b); c1 = VMULL_11LE(a, b);
t1 = VecPermute(a, a, m1); t1 = vec_mergel(a, a);
t1 = VecXor(a, t1); t1 = VecXor(a, t1);
t2 = VecPermute(b, b, m1); t2 = vec_mergel(b, b);
t2 = VecXor(b, t2); t2 = VecXor(b, t2);
t1 = VMULL_00LE(t1, t2); t1 = VMULL_00LE(t1, t2);
t1 = VecXor(c0, t1); t1 = VecXor(c0, t1);
t1 = VecXor(c1, t1); t1 = VecXor(c1, t1);
t2 = t1; t2 = t1;
t1 = VecPermute(z0, t1, m2); t1 = vec_mergeh(z0, t1);
t2 = VecPermute(t2, z0, m2); t2 = vec_mergel(t2, z0);
c0 = VecXor(c0, t1); c0 = VecXor(c0, t1);
c1 = VecXor(c1, t2); c1 = VecXor(c1, t2);
} }
@ -389,13 +383,10 @@ F2N_Multiply_128x128_POWER8(uint64x2_p& c1, uint64x2_p& c0, const uint64x2_p& a,
template <unsigned int N> template <unsigned int N>
inline uint64x2_p ShiftLeft128_POWER8(uint64x2_p x) inline uint64x2_p ShiftLeft128_POWER8(uint64x2_p x)
{ {
const uint8_t mb[] = {0,1,2,3, 4,5,6,7, 16,17,18,19, 20,21,22,23};
const uint8x16_p m = (uint8x16_p)VecLoad(mb);
uint64x2_p u=x, v, z={0}; uint64x2_p u=x, v, z={0};
x = VecShiftLeft<N>(x); x = VecShiftLeft<N>(x);
u = VecShiftRight<64-N>(u); u = VecShiftRight<64-N>(u);
v = VecPermute(z, u, m); v = vec_mergeh(z, u);
x = VecOr(x, v); x = VecOr(x, v);
return x; return x;
} }