TWO
st->h[0].d = TWO(52)*TWO0;
st->h[1].d = TWO(52)*TWO32;
st->h[2].d = TWO(52)*TWO64;
st->h[3].d = TWO(52)*TWO96;
st->r[0] = r0.d - TWO(52) * TWO0;
st->r[2] = r1.d - TWO(52) * TWO32;
st->r[4] = r2.d - TWO(52) * TWO64;
st->r[6] = r3.d - TWO(52) * TWO96;
st->r[1] = (st->r[0] + TWO(52) * TWO(16) * TWO0) - TWO(52) * TWO(16) * TWO0;
st->r[3] = (st->r[2] + TWO(52) * TWO(16) * TWO32) - TWO(52) * TWO(16) * TWO32;
st->r[5] = (st->r[4] + TWO(52) * TWO(16) * TWO64) - TWO(52) * TWO(16) * TWO64;
st->r[7] = (st->r[6] + TWO(52) * TWO(16) * TWO96) - TWO(52) * TWO(16) * TWO96;
st->s[1] = (st->s[0] + TWO(52) * TWO(16) * TWO0 / TWO96) - TWO(52) * TWO(16) * TWO0 / TWO96;
st->s[3] = (st->s[2] + TWO(52) * TWO(16) * TWO32 / TWO96) - TWO(52) * TWO(16) * TWO32 / TWO96;
st->s[5] = (st->s[4] + TWO(52) * TWO(16) * TWO64 / TWO96) - TWO(52) * TWO(16) * TWO64 / TWO96;
h0lo = st->h[0].d - TWO(52) * TWO0;
h1lo = st->h[1].d - TWO(52) * TWO32;
h2lo = st->h[2].d - TWO(52) * TWO64;
h3lo = st->h[3].d - TWO(52) * TWO96;
x0 = in0.d - TWO(52) * TWO0;
x1 = in1.d - TWO(52) * TWO32;
x2 = in2.d - TWO(52) * TWO64;
x3 = in3.d - TWO(52) * TWO96;
x0 = in0.d - TWO(52) * TWO0;
x1 = in1.d - TWO(52) * TWO32;
x2 = in2.d - TWO(52) * TWO64;
x3 = in3.d - TWO(52) * TWO96;
c0lo = (h0lo + TWO(52) * TWO32) - TWO(52) * TWO32;
c1lo = (h1lo + TWO(52) * TWO64) - TWO(52) * TWO64;
c2lo = (h2lo + TWO(52) * TWO96) - TWO(52) * TWO96;
c3lo = (h3lo + TWO(52) * TWO130) - TWO(52) * TWO130;
c0hi = (h0hi + TWO(52) * TWO32) - TWO(52) * TWO32;
c1hi = (h1hi + TWO(52) * TWO64) - TWO(52) * TWO64;
c2hi = (h2hi + TWO(52) * TWO96) - TWO(52) * TWO96;
c3hi = (h3hi + TWO(52) * TWO130) - TWO(52) * TWO130;
c0lo = (h0lo + TWO(52) * TWO32) - TWO(52) * TWO32;
c1lo = (h1lo + TWO(52) * TWO64) - TWO(52) * TWO64;
c2lo = (h2lo + TWO(52) * TWO96) - TWO(52) * TWO96;
c3lo = (h3lo + TWO(52) * TWO130) - TWO(52) * TWO130;
c0hi = (h0hi + TWO(52) * TWO32) - TWO(52) * TWO32;
c1hi = (h1hi + TWO(52) * TWO64) - TWO(52) * TWO64;
c2hi = (h2hi + TWO(52) * TWO96) - TWO(52) * TWO96;
c3hi = (h3hi + TWO(52) * TWO130) - TWO(52) * TWO130;
st->h[1].d = x1 + TWO(52) * TWO32;
st->h[2].d = x2 + TWO(52) * TWO64;
st->h[3].d = x3 + TWO(52) * TWO96;
st->h[0].d = x0 + TWO(52) * TWO0;
#define TWO0 TWO(0)
#define TWO32 TWO(32)
#define TWO64 (TWO32 * TWO(32))
#define TWO96 (TWO64 * TWO(32))
#define TWO130 (TWO96 * TWO(34))