#include "av1/common/daala_tx.h" #include "av1/common/odintrin.h" /* clang-format off */ # define OD_DCT_RSHIFT(_a, _b) OD_UNBIASED_RSHIFT32(_a, _b) /* TODO: Daala DCT overflow checks need to be ported as a later test */ # if defined(OD_DCT_CHECK_OVERFLOW) # else # define OD_DCT_OVERFLOW_CHECK(val, scale, offset, idx) # endif #define OD_FDCT_2(p0, p1) \ /* Embedded 2-point orthonormal Type-II fDCT. */ \ do { \ /* 13573/32768 ~= Tan[pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(p1, 13573, 16384, 100); \ p0 -= (p1*13573 + 16384) >> 15; \ /* 5793/8192 ~= Sin[pi/4] ~= 0.707106781186547 */ \ OD_DCT_OVERFLOW_CHECK(p0, 5793, 4096, 101); \ p1 += (p0*5793 + 4096) >> 13; \ /* 3393/8192 ~= Tan[pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(p1, 3393, 4096, 102); \ p0 -= (p1*3393 + 4096) >> 13; \ } \ while (0) #define OD_IDCT_2(p0, p1) \ /* Embedded 2-point orthonormal Type-II iDCT. */ \ do { \ /* 3393/8192 ~= Tan[pi/8] ~= 0.414213562373095 */ \ p0 += (p1*3393 + 4096) >> 13; \ /* 5793/8192 ~= Sin[pi/4] ~= 0.707106781186547 */ \ p1 -= (p0*5793 + 4096) >> 13; \ /* 13573/32768 ~= Tan[pi/8] ~= 0.414213562373095 */ \ p0 += (p1*13573 + 16384) >> 15; \ } \ while (0) #define OD_FDCT_2_ASYM(p0, p1, p1h) \ /* Embedded 2-point asymmetric Type-II fDCT. */ \ do { \ p0 += p1h; \ p1 = p0 - p1; \ } \ while (0) #define OD_IDCT_2_ASYM(p0, p1, p1h) \ /* Embedded 2-point asymmetric Type-II iDCT. */ \ do { \ p1 = p0 - p1; \ p1h = OD_DCT_RSHIFT(p1, 1); \ p0 -= p1h; \ } \ while (0) #define OD_FDST_2(p0, p1) \ /* Embedded 2-point orthonormal Type-IV fDST. */ \ do { \ /* 10947/16384 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ OD_DCT_OVERFLOW_CHECK(p1, 10947, 8192, 103); \ p0 -= (p1*10947 + 8192) >> 14; \ /* 473/512 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ OD_DCT_OVERFLOW_CHECK(p0, 473, 256, 104); \ p1 += (p0*473 + 256) >> 9; \ /* 10947/16384 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ OD_DCT_OVERFLOW_CHECK(p1, 10947, 8192, 105); \ p0 -= (p1*10947 + 8192) >> 14; \ } \ while (0) #define OD_IDST_2(p0, p1) \ /* Embedded 2-point orthonormal Type-IV iDST. */ \ do { \ /* 10947/16384 ~= Tan[3*Pi/16]) ~= 0.668178637919299 */ \ p0 += (p1*10947 + 8192) >> 14; \ /* 473/512 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ p1 -= (p0*473 + 256) >> 9; \ /* 10947/16384 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ p0 += (p1*10947 + 8192) >> 14; \ } \ while (0) #define OD_FDST_2_ASYM(p0, p1) \ /* Embedded 2-point asymmetric Type-IV fDST. */ \ do { \ /* 11507/16384 ~= 4*Sin[Pi/8] - 2*Tan[Pi/8] ~= 0.702306604714169 */ \ OD_DCT_OVERFLOW_CHECK(p1, 11507, 8192, 187); \ p0 -= (p1*11507 + 8192) >> 14; \ /* 669/1024 ~= Cos[Pi/8]/Sqrt[2] ~= 0.653281482438188 */ \ OD_DCT_OVERFLOW_CHECK(p0, 669, 512, 188); \ p1 += (p0*669 + 512) >> 10; \ /* 4573/4096 ~= 4*Sin[Pi/8] - Tan[Pi/8] ~= 1.11652016708726 */ \ OD_DCT_OVERFLOW_CHECK(p1, 4573, 2048, 189); \ p0 -= (p1*4573 + 2048) >> 12; \ } \ while (0) #define OD_IDST_2_ASYM(p0, p1) \ /* Embedded 2-point asymmetric Type-IV iDST. */ \ do { \ /* 4573/4096 ~= 4*Sin[Pi/8] - Tan[Pi/8] ~= 1.11652016708726 */ \ p0 += (p1*4573 + 2048) >> 12; \ /* 669/1024 ~= Cos[Pi/8]/Sqrt[2] ~= 0.653281482438188 */ \ p1 -= (p0*669 + 512) >> 10; \ /* 11507/16384 ~= 4*Sin[Pi/8] - 2*Tan[Pi/8] ~= 0.702306604714169 */ \ p0 += (p1*11507 + 8192) >> 14; \ } \ while (0) #define OD_FDCT_4(q0, q2, q1, q3) \ /* Embedded 4-point orthonormal Type-II fDCT. */ \ do { \ int q2h; \ int q3h; \ q3 = q0 - q3; \ q3h = OD_DCT_RSHIFT(q3, 1); \ q0 -= q3h; \ q2 += q1; \ q2h = OD_DCT_RSHIFT(q2, 1); \ q1 = q2h - q1; \ OD_FDCT_2_ASYM(q0, q2, q2h); \ OD_FDST_2_ASYM(q3, q1); \ } \ while (0) #define OD_IDCT_4(q0, q2, q1, q3) \ /* Embedded 4-point orthonormal Type-II iDCT. */ \ do { \ int q1h; \ int q3h; \ OD_IDST_2_ASYM(q3, q2); \ OD_IDCT_2_ASYM(q0, q1, q1h); \ q3h = OD_DCT_RSHIFT(q3, 1); \ q0 += q3h; \ q3 = q0 - q3; \ q2 = q1h - q2; \ q1 -= q2; \ } \ while (0) #define OD_FDCT_4_ASYM(q0, q2, q2h, q1, q3, q3h) \ /* Embedded 4-point asymmetric Type-II fDCT. */ \ do { \ q0 += q3h; \ q3 = q0 - q3; \ q1 = q2h - q1; \ q2 = q1 - q2; \ OD_FDCT_2(q0, q2); \ OD_FDST_2(q3, q1); \ } \ while (0) #define OD_IDCT_4_ASYM(q0, q2, q1, q1h, q3, q3h) \ /* Embedded 4-point asymmetric Type-II iDCT. */ \ do { \ OD_IDST_2(q3, q2); \ OD_IDCT_2(q0, q1); \ q1 = q2 - q1; \ q1h = OD_DCT_RSHIFT(q1, 1); \ q2 = q1h - q2; \ q3 = q0 - q3; \ q3h = OD_DCT_RSHIFT(q3, 1); \ q0 -= q3h; \ } \ while (0) #define OD_FDST_4(q0, q2, q1, q3) \ /* Embedded 4-point orthonormal Type-IV fDST. */ \ do { \ int q0h; \ int q1h; \ /* 13573/32768 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(q1, 13573, 16384, 190); \ q2 += (q1*13573 + 16384) >> 15; \ /* 5793/8192 ~= Sin[Pi/4] ~= 0.707106781186547 */ \ OD_DCT_OVERFLOW_CHECK(q2, 5793, 4096, 191); \ q1 -= (q2*5793 + 4096) >> 13; \ /* 3393/8192 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(q1, 3393, 4096, 192); \ q2 += (q1*3393 + 4096) >> 13; \ q0 += q2; \ q0h = OD_DCT_RSHIFT(q0, 1); \ q2 = q0h - q2; \ q1 += q3; \ q1h = OD_DCT_RSHIFT(q1, 1); \ q3 -= q1h; \ /* 537/1024 ~= (1/Sqrt[2] - Cos[3*Pi/16]/2)/Sin[3*Pi/16] ~= 0.524455699240090 */ \ OD_DCT_OVERFLOW_CHECK(q1, 537, 512, 193); \ q2 -= (q1*537 + 512) >> 10; \ /* 1609/2048 ~= Sqrt[2]*Sin[3*Pi/16] ~= 0.785694958387102 */ \ OD_DCT_OVERFLOW_CHECK(q2, 1609, 1024, 194); \ q1 += (q2*1609 + 1024) >> 11; \ /* 7335/32768 ~= (1/Sqrt[2] - Cos[3*Pi/16])/Sin[3*Pi/16] ~= 0.223847182092655 */ \ OD_DCT_OVERFLOW_CHECK(q1, 7335, 16384, 195); \ q2 += (q1*7335 + 16384) >> 15; \ /* 5091/8192 ~= (1/Sqrt[2] - Cos[7*Pi/16]/2)/Sin[7*Pi/16] ~= 0.6215036383171189 */ \ OD_DCT_OVERFLOW_CHECK(q0, 5091, 4096, 196); \ q3 += (q0*5091 + 4096) >> 13; \ /* 5681/4096 ~= Sqrt[2]*Sin[7*Pi/16] ~= 1.38703984532215 */ \ OD_DCT_OVERFLOW_CHECK(q3, 5681, 2048, 197); \ q0 -= (q3*5681 + 2048) >> 12; \ /* 4277/8192 ~= (1/Sqrt[2] - Cos[7*Pi/16])/Sin[7*Pi/16] ~= 0.52204745462729 */ \ OD_DCT_OVERFLOW_CHECK(q0, 4277, 4096, 198); \ q3 += (q0*4277 + 4096) >> 13; \ } \ while (0) #define OD_IDST_4(q0, q2, q1, q3) \ /* Embedded 4-point orthonormal Type-IV iDST. */ \ do { \ int q0h; \ int q2h; \ /* 4277/8192 ~= (1/Sqrt[2] - Cos[7*Pi/16])/Sin[7*Pi/16] ~= 0.52204745462729 */ \ q3 -= (q0*4277 + 4096) >> 13; \ /* 5681/4096 ~= Sqrt[2]*Sin[7*Pi/16] ~= 1.38703984532215 */ \ q0 += (q3*5681 + 2048) >> 12; \ /* 5091/8192 ~= (1/Sqrt[2] - Cos[7*Pi/16]/2)/Sin[7*Pi/16] ~= 0.6215036383171189 */ \ q3 -= (q0*5091 + 4096) >> 13; \ /* 7335/32768 ~= (1/Sqrt[2] - Cos[3*Pi/16])/Sin[3*Pi/16] ~= 0.223847182092655 */ \ q1 -= (q2*7335 + 16384) >> 15; \ /* 1609/2048 ~= Sqrt[2]*Sin[3*Pi/16] ~= 0.785694958387102 */ \ q2 -= (q1*1609 + 1024) >> 11; \ /* 537/1024 ~= (1/Sqrt[2] - Cos[3*Pi/16]/2)/Sin[3*Pi/16] ~= 0.524455699240090 */ \ q1 += (q2*537 + 512) >> 10; \ q2h = OD_DCT_RSHIFT(q2, 1); \ q3 += q2h; \ q2 -= q3; \ q0h = OD_DCT_RSHIFT(q0, 1); \ q1 = q0h - q1; \ q0 -= q1; \ /* 3393/8192 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ q1 -= (q2*3393 + 4096) >> 13; \ /* 5793/8192 ~= Sin[Pi/4] ~= 0.707106781186547 */ \ q2 += (q1*5793 + 4096) >> 13; \ /* 13573/32768 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ q1 -= (q2*13573 + 16384) >> 15; \ } \ while (0) #define OD_FDST_4_ASYM(t0, t0h, t2, t1, t3) \ /* Embedded 4-point asymmetric Type-IV fDST. */ \ do { \ /* 7489/8192 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ OD_DCT_OVERFLOW_CHECK(t1, 7489, 4096, 106); \ t2 -= (t1*7489 + 4096) >> 13; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ OD_DCT_OVERFLOW_CHECK(t1, 11585, 8192, 107); \ t1 += (t2*11585 + 8192) >> 14; \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ OD_DCT_OVERFLOW_CHECK(t1, 19195, 16384, 108); \ t2 += (t1*19195 + 16384) >> 15; \ t3 += OD_DCT_RSHIFT(t2, 1); \ t2 -= t3; \ t1 = t0h - t1; \ t0 -= t1; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.820678790828660 */ \ OD_DCT_OVERFLOW_CHECK(t0, 6723, 4096, 109); \ t3 += (t0*6723 + 4096) >> 13; \ /* 8035/8192 ~= Sin[7*Pi/16] ~= 0.980785280403230 */ \ OD_DCT_OVERFLOW_CHECK(t3, 8035, 4096, 110); \ t0 -= (t3*8035 + 4096) >> 13; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.820678790828660 */ \ OD_DCT_OVERFLOW_CHECK(t0, 6723, 4096, 111); \ t3 += (t0*6723 + 4096) >> 13; \ /* 8757/16384 ~= Tan[5*Pi/32] ~= 0.534511135950792 */ \ OD_DCT_OVERFLOW_CHECK(t1, 8757, 8192, 112); \ t2 += (t1*8757 + 8192) >> 14; \ /* 6811/8192 ~= Sin[5*Pi/16] ~= 0.831469612302545 */ \ OD_DCT_OVERFLOW_CHECK(t2, 6811, 4096, 113); \ t1 -= (t2*6811 + 4096) >> 13; \ /* 8757/16384 ~= Tan[5*Pi/32] ~= 0.534511135950792 */ \ OD_DCT_OVERFLOW_CHECK(t1, 8757, 8192, 114); \ t2 += (t1*8757 + 8192) >> 14; \ } \ while (0) #define OD_IDST_4_ASYM(t0, t0h, t2, t1, t3) \ /* Embedded 4-point asymmetric Type-IV iDST. */ \ do { \ /* 8757/16384 ~= Tan[5*Pi/32] ~= 0.534511135950792 */ \ t1 -= (t2*8757 + 8192) >> 14; \ /* 6811/8192 ~= Sin[5*Pi/16] ~= 0.831469612302545 */ \ t2 += (t1*6811 + 4096) >> 13; \ /* 8757/16384 ~= Tan[5*Pi/32] ~= 0.534511135950792 */ \ t1 -= (t2*8757 + 8192) >> 14; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.820678790828660 */ \ t3 -= (t0*6723 + 4096) >> 13; \ /* 8035/8192 ~= Sin[7*Pi/16] ~= 0.980785280403230 */ \ t0 += (t3*8035 + 4096) >> 13; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.820678790828660 */ \ t3 -= (t0*6723 + 4096) >> 13; \ t0 += t2; \ t0h = OD_DCT_RSHIFT(t0, 1); \ t2 = t0h - t2; \ t1 += t3; \ t3 -= OD_DCT_RSHIFT(t1, 1); \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ t1 -= (t2*19195 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ t2 -= (t1*11585 + 8192) >> 14; \ /* 7489/8192 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ t1 += (t2*7489 + 4096) >> 13; \ } \ while (0) #define OD_FDCT_8(r0, r4, r2, r6, r1, r5, r3, r7) \ /* Embedded 8-point orthonormal Type-II fDCT. */ \ do { \ int r4h; \ int r5h; \ int r6h; \ int r7h; \ r7 = r0 - r7; \ r7h = OD_DCT_RSHIFT(r7, 1); \ r0 -= r7h; \ r6 += r1; \ r6h = OD_DCT_RSHIFT(r6, 1); \ r1 = r6h - r1; \ r5 = r2 - r5; \ r5h = OD_DCT_RSHIFT(r5, 1); \ r2 -= r5h; \ r4 += r3; \ r4h = OD_DCT_RSHIFT(r4, 1); \ r3 = r4h - r3; \ OD_FDCT_4_ASYM(r0, r4, r4h, r2, r6, r6h); \ OD_FDST_4_ASYM(r7, r7h, r3, r5, r1); \ } \ while (0) #define OD_IDCT_8(r0, r4, r2, r6, r1, r5, r3, r7) \ /* Embedded 8-point orthonormal Type-II iDCT. */ \ do { \ int r1h; \ int r3h; \ int r5h; \ int r7h; \ OD_IDST_4_ASYM(r7, r7h, r5, r6, r4); \ OD_IDCT_4_ASYM(r0, r2, r1, r1h, r3, r3h); \ r0 += r7h; \ r7 = r0 - r7; \ r6 = r1h - r6; \ r1 -= r6; \ r5h = OD_DCT_RSHIFT(r5, 1); \ r2 += r5h; \ r5 = r2 - r5; \ r4 = r3h - r4; \ r3 -= r4; \ } \ while (0) #define OD_FDCT_8_ASYM(r0, r4, r4h, r2, r6, r6h, r1, r5, r5h, r3, r7, r7h) \ /* Embedded 8-point asymmetric Type-II fDCT. */ \ do { \ r0 += r7h; \ r7 = r0 - r7; \ r1 = r6h - r1; \ r6 -= r1; \ r2 += r5h; \ r5 = r2 - r5; \ r3 = r4h - r3; \ r4 -= r3; \ OD_FDCT_4(r0, r4, r2, r6); \ OD_FDST_4(r7, r3, r5, r1); \ } \ while (0) #define OD_IDCT_8_ASYM(r0, r4, r2, r6, r1, r1h, r5, r5h, r3, r3h, r7, r7h) \ /* Embedded 8-point asymmetric Type-II iDCT. */ \ do { \ OD_IDST_4(r7, r5, r6, r4); \ OD_IDCT_4(r0, r2, r1, r3); \ r7 = r0 - r7; \ r7h = OD_DCT_RSHIFT(r7, 1); \ r0 -= r7h; \ r1 += r6; \ r1h = OD_DCT_RSHIFT(r1, 1); \ r6 = r1h - r6; \ r5 = r2 - r5; \ r5h = OD_DCT_RSHIFT(r5, 1); \ r2 -= r5h; \ r3 += r4; \ r3h = OD_DCT_RSHIFT(r3, 1); \ r4 = r3h - r4; \ } \ while (0) #define OD_FDST_8(t0, t4, t2, t6, t1, t5, t3, t7) \ /* Embedded 8-point orthonormal Type-IV fDST. */ \ do { \ int t0h; \ int t2h; \ int t5h; \ int t7h; \ /* 13573/32768 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(t1, 13573, 16384, 115); \ t6 -= (t1*13573 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186547 */ \ OD_DCT_OVERFLOW_CHECK(t6, 11585, 8192, 116); \ t1 += (t6*11585 + 8192) >> 14; \ /* 13573/32768 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(t1, 13573, 16384, 117); \ t6 -= (t1*13573 + 16384) >> 15; \ /* 21895/32768 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ OD_DCT_OVERFLOW_CHECK(t2, 21895, 16384, 118); \ t5 -= (t2*21895 + 16384) >> 15; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ OD_DCT_OVERFLOW_CHECK(t5, 15137, 8192, 119); \ t2 += (t5*15137 + 8192) >> 14; \ /* 10947/16384 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ OD_DCT_OVERFLOW_CHECK(t2, 10947, 8192, 120); \ t5 -= (t2*10947 + 8192) >> 14; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ OD_DCT_OVERFLOW_CHECK(t3, 3259, 8192, 121); \ t4 -= (t3*3259 + 8192) >> 14; \ /* 3135/8192 ~= Sin[Pi/8] ~= 0.382683432365090 */ \ OD_DCT_OVERFLOW_CHECK(t4, 3135, 4096, 122); \ t3 += (t4*3135 + 4096) >> 13; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ OD_DCT_OVERFLOW_CHECK(t3, 3259, 8192, 123); \ t4 -= (t3*3259 + 8192) >> 14; \ t7 += t1; \ t7h = OD_DCT_RSHIFT(t7, 1); \ t1 -= t7h; \ t2 = t3 - t2; \ t2h = OD_DCT_RSHIFT(t2, 1); \ t3 -= t2h; \ t0 -= t6; \ t0h = OD_DCT_RSHIFT(t0, 1); \ t6 += t0h; \ t5 = t4 - t5; \ t5h = OD_DCT_RSHIFT(t5, 1); \ t4 -= t5h; \ t1 += t5h; \ t5 = t1 - t5; \ t4 += t0h; \ t0 -= t4; \ t6 -= t2h; \ t2 += t6; \ t3 -= t7h; \ t7 += t3; \ /* TODO: Can we move this into another operation */ \ t7 = -t7; \ /* 7425/8192 ~= Tan[15*Pi/64] ~= 0.906347169019147 */ \ OD_DCT_OVERFLOW_CHECK(t7, 7425, 4096, 124); \ t0 -= (t7*7425 + 4096) >> 13; \ /* 8153/8192 ~= Sin[15*Pi/32] ~= 0.995184726672197 */ \ OD_DCT_OVERFLOW_CHECK(t0, 8153, 4096, 125); \ t7 += (t0*8153 + 4096) >> 13; \ /* 7425/8192 ~= Tan[15*Pi/64] ~= 0.906347169019147 */ \ OD_DCT_OVERFLOW_CHECK(t7, 7425, 4096, 126); \ t0 -= (t7*7425 + 4096) >> 13; \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.148335987538347 */ \ OD_DCT_OVERFLOW_CHECK(t1, 4861, 16384, 127); \ t6 -= (t1*4861 + 16384) >> 15; \ /* 1189/4096 ~= Sin[3*Pi/32] ~= 0.290284677254462 */ \ OD_DCT_OVERFLOW_CHECK(t6, 1189, 2048, 128); \ t1 += (t6*1189 + 2048) >> 12; \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.148335987538347 */ \ OD_DCT_OVERFLOW_CHECK(t1, 4861, 16384, 129); \ t6 -= (t1*4861 + 16384) >> 15; \ /* 2455/4096 ~= Tan[11*Pi/64] ~= 0.599376933681924 */ \ OD_DCT_OVERFLOW_CHECK(t5, 2455, 2048, 130); \ t2 -= (t5*2455 + 2048) >> 12; \ /* 7225/8192 ~= Sin[11*Pi/32] ~= 0.881921264348355 */ \ OD_DCT_OVERFLOW_CHECK(t2, 7225, 4096, 131); \ t5 += (t2*7225 + 4096) >> 13; \ /* 2455/4096 ~= Tan[11*Pi/64] ~= 0.599376933681924 */ \ OD_DCT_OVERFLOW_CHECK(t5, 2455, 2048, 132); \ t2 -= (t5*2455 + 2048) >> 12; \ /* 11725/32768 ~= Tan[7*Pi/64] ~= 0.357805721314524 */ \ OD_DCT_OVERFLOW_CHECK(t3, 11725, 16384, 133); \ t4 -= (t3*11725 + 16384) >> 15; \ /* 5197/8192 ~= Sin[7*Pi/32] ~= 0.634393284163645 */ \ OD_DCT_OVERFLOW_CHECK(t4, 5197, 4096, 134); \ t3 += (t4*5197 + 4096) >> 13; \ /* 11725/32768 ~= Tan[7*Pi/64] ~= 0.357805721314524 */ \ OD_DCT_OVERFLOW_CHECK(t3, 11725, 16384, 135); \ t4 -= (t3*11725 + 16384) >> 15; \ } \ while (0) #define OD_IDST_8(t0, t4, t2, t6, t1, t5, t3, t7) \ /* Embedded 8-point orthonormal Type-IV iDST. */ \ do { \ int t0h; \ int t2h; \ int t5h_; \ int t7h_; \ /* 11725/32768 ~= Tan[7*Pi/64] ~= 0.357805721314524 */ \ t1 += (t6*11725 + 16384) >> 15; \ /* 5197/8192 ~= Sin[7*Pi/32] ~= 0.634393284163645 */ \ t6 -= (t1*5197 + 4096) >> 13; \ /* 11725/32768 ~= Tan[7*Pi/64] ~= 0.357805721314524 */ \ t1 += (t6*11725 + 16384) >> 15; \ /* 2455/4096 ~= Tan[11*Pi/64] ~= 0.599376933681924 */ \ t2 += (t5*2455 + 2048) >> 12; \ /* 7225/8192 ~= Sin[11*Pi/32] ~= 0.881921264348355 */ \ t5 -= (t2*7225 + 4096) >> 13; \ /* 2455/4096 ~= Tan[11*Pi/64] ~= 0.599376933681924 */ \ t2 += (t5*2455 + 2048) >> 12; \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.148335987538347 */ \ t3 += (t4*4861 + 16384) >> 15; \ /* 1189/4096 ~= Sin[3*Pi/32] ~= 0.290284677254462 */ \ t4 -= (t3*1189 + 2048) >> 12; \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.148335987538347 */ \ t3 += (t4*4861 + 16384) >> 15; \ /* 7425/8192 ~= Tan[15*Pi/64] ~= 0.906347169019147 */ \ t0 += (t7*7425 + 4096) >> 13; \ /* 8153/8192 ~= Sin[15*Pi/32] ~= 0.995184726672197 */ \ t7 -= (t0*8153 + 4096) >> 13; \ /* 7425/8192 ~= Tan[15*Pi/64] ~= 0.906347169019147 */ \ t0 += (t7*7425 + 4096) >> 13; \ /* TODO: Can we move this into another operation */ \ t7 = -t7; \ t7 -= t6; \ t7h_ = OD_DCT_RSHIFT(t7, 1); \ t6 += t7h_; \ t2 -= t3; \ t2h = OD_DCT_RSHIFT(t2, 1); \ t3 += t2h; \ t0 += t1; \ t0h = OD_DCT_RSHIFT(t0, 1); \ t1 -= t0h; \ t5 = t4 - t5; \ t5h_ = OD_DCT_RSHIFT(t5, 1); \ t4 -= t5h_; \ t1 += t5h_; \ t5 = t1 - t5; \ t3 -= t0h; \ t0 += t3; \ t6 += t2h; \ t2 = t6 - t2; \ t4 += t7h_; \ t7 -= t4; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ t1 += (t6*3259 + 8192) >> 14; \ /* 3135/8192 ~= Sin[Pi/8] ~= 0.382683432365090 */ \ t6 -= (t1*3135 + 4096) >> 13; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ t1 += (t6*3259 + 8192) >> 14; \ /* 10947/16384 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ t5 += (t2*10947 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ t2 -= (t5*15137 + 8192) >> 14; \ /* 21895/32768 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ t5 += (t2*21895 + 16384) >> 15; \ /* 13573/32768 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ t3 += (t4*13573 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186547 */ \ t4 -= (t3*11585 + 8192) >> 14; \ /* 13573/32768 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ t3 += (t4*13573 + 16384) >> 15; \ } \ while (0) /* Rewrite this so that t0h can be passed in. */ #define OD_FDST_8_ASYM(t0, t4, t2, t6, t1, t5, t3, t7) \ /* Embedded 8-point asymmetric Type-IV fDST. */ \ do { \ int t0h; \ int t2h; \ int t5h; \ int t7h; \ /* 1035/2048 ~= (Sqrt[2] - Cos[7*Pi/32])/(2*Sin[7*Pi/32]) */ \ OD_DCT_OVERFLOW_CHECK(t1, 1035, 1024, 199); \ t6 += (t1*1035 + 1024) >> 11; \ /* 3675/4096 ~= Sqrt[2]*Sin[7*Pi/32] */ \ OD_DCT_OVERFLOW_CHECK(t6, 3675, 2048, 200); \ t1 -= (t6*3675 + 2048) >> 12; \ /* 851/8192 ~= (Cos[7*Pi/32] - 1/Sqrt[2])/Sin[7*Pi/32] */ \ OD_DCT_OVERFLOW_CHECK(t1, 851, 4096, 201); \ t6 -= (t1*851 + 4096) >> 13; \ /* 4379/8192 ~= (Sqrt[2] - Sin[5*Pi/32])/(2*Cos[5*Pi/32]) */ \ OD_DCT_OVERFLOW_CHECK(t2, 4379, 4096, 202); \ t5 += (t2*4379 + 4096) >> 13; \ /* 10217/8192 ~= Sqrt[2]*Cos[5*Pi/32] */ \ OD_DCT_OVERFLOW_CHECK(t5, 10217, 4096, 203); \ t2 -= (t5*10217 + 4096) >> 13; \ /* 4379/16384 ~= (1/Sqrt[2] - Sin[5*Pi/32])/Cos[5*Pi/32] */ \ OD_DCT_OVERFLOW_CHECK(t2, 4379, 8192, 204); \ t5 += (t2*4379 + 8192) >> 14; \ /* 12905/16384 ~= (Sqrt[2] - Cos[3*Pi/32])/(2*Sin[3*Pi/32]) */ \ OD_DCT_OVERFLOW_CHECK(t3, 12905, 8192, 205); \ t4 += (t3*12905 + 8192) >> 14; \ /* 3363/8192 ~= Sqrt[2]*Sin[3*Pi/32] */ \ OD_DCT_OVERFLOW_CHECK(t4, 3363, 4096, 206); \ t3 -= (t4*3363 + 4096) >> 13; \ /* 3525/4096 ~= (Cos[3*Pi/32] - 1/Sqrt[2])/Sin[3*Pi/32] */ \ OD_DCT_OVERFLOW_CHECK(t3, 3525, 2048, 207); \ t4 -= (t3*3525 + 2048) >> 12; \ /* 5417/8192 ~= (Sqrt[2] - Sin[Pi/32])/(2*Cos[Pi/32]) */ \ OD_DCT_OVERFLOW_CHECK(t0, 5417, 4096, 208); \ t7 += (t0*5417 + 4096) >> 13; \ /* 5765/4096 ~= Sqrt[2]*Cos[Pi/32] */ \ OD_DCT_OVERFLOW_CHECK(t7, 5765, 2048, 209); \ t0 -= (t7*5765 + 2048) >> 12; \ /* 2507/4096 ~= (1/Sqrt[2] - Sin[Pi/32])/Cos[Pi/32] */ \ OD_DCT_OVERFLOW_CHECK(t0, 2507, 2048, 210); \ t7 += (t0*2507 + 2048) >> 12; \ t0 += t1; \ t0h = OD_DCT_RSHIFT(t0, 1); \ t1 -= t0h; \ t2 -= t3; \ t2h = OD_DCT_RSHIFT(t2, 1); \ t3 += t2h; \ t5 -= t4; \ t5h = OD_DCT_RSHIFT(t5, 1); \ t4 += t5h; \ t7 += t6; \ t7h = OD_DCT_RSHIFT(t7, 1); \ t6 = t7h - t6; \ t4 = t7h - t4; \ t7 -= t4; \ t1 += t5h; \ t5 = t1 - t5; \ t6 += t2h; \ t2 = t6 - t2; \ t3 -= t0h; \ t0 += t3; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ OD_DCT_OVERFLOW_CHECK(t6, 3259, 8192, 211); \ t1 += (t6*3259 + 8192) >> 14; \ /* 3135/8192 ~= Sin[Pi/8] ~= 0.382683432365090 */ \ OD_DCT_OVERFLOW_CHECK(t1, 3135, 4096, 212); \ t6 -= (t1*3135 + 4096) >> 13; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ OD_DCT_OVERFLOW_CHECK(t6, 3259, 8192, 213); \ t1 += (t6*3259 + 8192) >> 14; \ /* 2737/4096 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ OD_DCT_OVERFLOW_CHECK(t2, 2737, 2048, 214); \ t5 += (t2*2737 + 2048) >> 12; \ /* 473/512 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ OD_DCT_OVERFLOW_CHECK(t5, 473, 256, 215); \ t2 -= (t5*473 + 256) >> 9; \ /* 2737/4096 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ OD_DCT_OVERFLOW_CHECK(t2, 2737, 2048, 216); \ t5 += (t2*2737 + 2048) >> 12; \ /* 3393/8192 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(t4, 3393, 4096, 217); \ t3 += (t4*3393 + 4096) >> 13; \ /* 5793/8192 ~= Sin[Pi/4] ~= 0.707106781186547 */ \ OD_DCT_OVERFLOW_CHECK(t3, 5793, 4096, 218); \ t4 -= (t3*5793 + 4096) >> 13; \ /* 3393/8192 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(t4, 3393, 4096, 219); \ t3 += (t4*3393 + 4096) >> 13; \ } \ while (0) #define OD_IDST_8_ASYM(t0, t4, t2, t6, t1, t5, t3, t7) \ /* Embedded 8-point asymmetric Type-IV iDST. */ \ do { \ int t0h; \ int t2h; \ int t5h__; \ int t7h__; \ /* 3393/8192 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ t6 -= (t1*3393 + 4096) >> 13; \ /* 5793/8192 ~= Sin[Pi/4] ~= 0.707106781186547 */ \ t1 += (t6*5793 + 4096) >> 13; \ /* 3393/8192 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ t6 -= (t1*3393 + 4096) >> 13; \ /* 2737/4096 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ t5 -= (t2*2737 + 2048) >> 12; \ /* 473/512 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ t2 += (t5*473 + 256) >> 9; \ /* 2737/4096 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ t5 -= (t2*2737 + 2048) >> 12; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ t4 -= (t3*3259 + 8192) >> 14; \ /* 3135/8192 ~= Sin[Pi/8] ~= 0.382683432365090 */ \ t3 += (t4*3135 + 4096) >> 13; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ t4 -= (t3*3259 + 8192) >> 14; \ t0 -= t6; \ t0h = OD_DCT_RSHIFT(t0, 1); \ t6 += t0h; \ t2 = t3 - t2; \ t2h = OD_DCT_RSHIFT(t2, 1); \ t3 -= t2h; \ t5 = t4 - t5; \ t5h__ = OD_DCT_RSHIFT(t5, 1); \ t4 -= t5h__; \ t7 += t1; \ t7h__ = OD_DCT_RSHIFT(t7, 1); \ t1 = t7h__ - t1; \ t3 = t7h__ - t3; \ t7 -= t3; \ t1 -= t5h__; \ t5 += t1; \ t6 -= t2h; \ t2 += t6; \ t4 += t0h; \ t0 -= t4; \ /* 2507/4096 ~= (1/Sqrt[2] - Sin[Pi/32])/Cos[Pi/32] */ \ t7 -= (t0*2507 + 2048) >> 12; \ /* 5765/4096 ~= Sqrt[2]*Cos[Pi/32] */ \ t0 += (t7*5765 + 2048) >> 12; \ /* 5417/8192 ~= (Sqrt[2] - Sin[Pi/32])/(2*Cos[Pi/32]) */ \ t7 -= (t0*5417 + 4096) >> 13; \ /* 3525/4096 ~= (Cos[3*Pi/32] - 1/Sqrt[2])/Sin[3*Pi/32] */ \ t1 += (t6*3525 + 2048) >> 12; \ /* 3363/8192 ~= Sqrt[2]*Sin[3*Pi/32] */ \ t6 += (t1*3363 + 4096) >> 13; \ /* 12905/16384 ~= (1/Sqrt[2] - Cos[3*Pi/32]/1)/Sin[3*Pi/32] */ \ t1 -= (t6*12905 + 8192) >> 14; \ /* 4379/16384 ~= (1/Sqrt[2] - Sin[5*Pi/32])/Cos[5*Pi/32] */ \ t5 -= (t2*4379 + 8192) >> 14; \ /* 10217/8192 ~= Sqrt[2]*Cos[5*Pi/32] */ \ t2 += (t5*10217 + 4096) >> 13; \ /* 4379/8192 ~= (Sqrt[2] - Sin[5*Pi/32])/(2*Cos[5*Pi/32]) */ \ t5 -= (t2*4379 + 4096) >> 13; \ /* 851/8192 ~= (Cos[7*Pi/32] - 1/Sqrt[2])/Sin[7*Pi/32] */ \ t3 += (t4*851 + 4096) >> 13; \ /* 3675/4096 ~= Sqrt[2]*Sin[7*Pi/32] */ \ t4 += (t3*3675 + 2048) >> 12; \ /* 1035/2048 ~= (Sqrt[2] - Cos[7*Pi/32])/(2*Sin[7*Pi/32]) */ \ t3 -= (t4*1035 + 1024) >> 11; \ } \ while (0) #define OD_FDCT_16(s0, s8, s4, sc, s2, sa, s6, se, \ s1, s9, s5, sd, s3, sb, s7, sf) \ /* Embedded 16-point orthonormal Type-II fDCT. */ \ do { \ int s8h; \ int sah; \ int sch; \ int seh; \ int sfh; \ sf = s0 - sf; \ sfh = OD_DCT_RSHIFT(sf, 1); \ s0 -= sfh; \ se += s1; \ seh = OD_DCT_RSHIFT(se, 1); \ s1 = seh - s1; \ sd = s2 - sd; \ s2 -= OD_DCT_RSHIFT(sd, 1); \ sc += s3; \ sch = OD_DCT_RSHIFT(sc, 1); \ s3 = sch - s3; \ sb = s4 - sb; \ s4 -= OD_DCT_RSHIFT(sb, 1); \ sa += s5; \ sah = OD_DCT_RSHIFT(sa, 1); \ s5 = sah - s5; \ s9 = s6 - s9; \ s6 -= OD_DCT_RSHIFT(s9, 1); \ s8 += s7; \ s8h = OD_DCT_RSHIFT(s8, 1); \ s7 = s8h - s7; \ OD_FDCT_8_ASYM(s0, s8, s8h, s4, sc, sch, s2, sa, sah, s6, se, seh); \ OD_FDST_8_ASYM(sf, s7, sb, s3, sd, s5, s9, s1); \ } \ while (0) #define OD_IDCT_16(s0, s8, s4, sc, s2, sa, s6, se, \ s1, s9, s5, sd, s3, sb, s7, sf) \ /* Embedded 16-point orthonormal Type-II iDCT. */ \ do { \ int s1h; \ int s3h; \ int s5h; \ int s7h; \ int sfh; \ OD_IDST_8_ASYM(sf, sb, sd, s9, se, sa, sc, s8); \ OD_IDCT_8_ASYM(s0, s4, s2, s6, s1, s1h, s5, s5h, s3, s3h, s7, s7h); \ sfh = OD_DCT_RSHIFT(sf, 1); \ s0 += sfh; \ sf = s0 - sf; \ se = s1h - se; \ s1 -= se; \ s2 += OD_DCT_RSHIFT(sd, 1); \ sd = s2 - sd; \ sc = s3h - sc; \ s3 -= sc; \ s4 += OD_DCT_RSHIFT(sb, 1); \ sb = s4 - sb; \ sa = s5h - sa; \ s5 -= sa; \ s6 += OD_DCT_RSHIFT(s9, 1); \ s9 = s6 - s9; \ s8 = s7h - s8; \ s7 -= s8; \ } \ while (0) #define OD_FDCT_16_ASYM(t0, t8, t8h, t4, tc, tch, t2, ta, tah, t6, te, teh, \ t1, t9, t9h, t5, td, tdh, t3, tb, tbh, t7, tf, tfh) \ /* Embedded 16-point asymmetric Type-II fDCT. */ \ do { \ t0 += tfh; \ tf = t0 - tf; \ t1 -= teh; \ te += t1; \ t2 += tdh; \ td = t2 - td; \ t3 -= tch; \ tc += t3; \ t4 += tbh; \ tb = t4 - tb; \ t5 -= tah; \ ta += t5; \ t6 += t9h; \ t9 = t6 - t9; \ t7 -= t8h; \ t8 += t7; \ OD_FDCT_8(t0, t8, t4, tc, t2, ta, t6, te); \ OD_FDST_8(tf, t7, tb, t3, td, t5, t9, t1); \ } \ while (0) #define OD_IDCT_16_ASYM(t0, t8, t4, tc, t2, ta, t6, te, \ t1, t1h, t9, t9h, t5, t5h, td, tdh, t3, t3h, tb, tbh, t7, t7h, tf, tfh) \ /* Embedded 16-point asymmetric Type-II iDCT. */ \ do { \ OD_IDST_8(tf, tb, td, t9, te, ta, tc, t8); \ OD_IDCT_8(t0, t4, t2, t6, t1, t5, t3, t7); \ t1 -= te; \ t1h = OD_DCT_RSHIFT(t1, 1); \ te += t1h; \ t9 = t6 - t9; \ t9h = OD_DCT_RSHIFT(t9, 1); \ t6 -= t9h; \ t5 -= ta; \ t5h = OD_DCT_RSHIFT(t5, 1); \ ta += t5h; \ td = t2 - td; \ tdh = OD_DCT_RSHIFT(td, 1); \ t2 -= tdh; \ t3 -= tc; \ t3h = OD_DCT_RSHIFT(t3, 1); \ tc += t3h; \ tb = t4 - tb; \ tbh = OD_DCT_RSHIFT(tb, 1); \ t4 -= tbh; \ t7 -= t8; \ t7h = OD_DCT_RSHIFT(t7, 1); \ t8 += t7h; \ tf = t0 - tf; \ tfh = OD_DCT_RSHIFT(tf, 1); \ t0 -= tfh; \ } \ while (0) #define OD_FDST_16(s0, s8, s4, sc, s2, sa, s6, se, \ s1, s9, s5, sd, s3, sb, s7, sf) \ /* Embedded 16-point orthonormal Type-IV fDST. */ \ do { \ int s0h; \ int s2h; \ int sdh; \ int sfh; \ /* 13573/32768 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(s3, 13573, 16384, 220); \ s1 += (se*13573 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186547 */ \ OD_DCT_OVERFLOW_CHECK(s1, 11585, 8192, 221); \ se -= (s1*11585 + 8192) >> 14; \ /* 13573/32768 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(s3, 13573, 16384, 222); \ s1 += (se*13573 + 16384) >> 15; \ /* 21895/32768 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ OD_DCT_OVERFLOW_CHECK(s2, 21895, 16384, 223); \ sd += (s2*21895 + 16384) >> 15; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ OD_DCT_OVERFLOW_CHECK(sd, 15137, 16384, 224); \ s2 -= (sd*15137 + 8192) >> 14; \ /* 21895/32768 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ OD_DCT_OVERFLOW_CHECK(s2, 21895, 16384, 225); \ sd += (s2*21895 + 16384) >> 15; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ OD_DCT_OVERFLOW_CHECK(s3, 3259, 8192, 226); \ sc += (s3*3259 + 8192) >> 14; \ /* 3135/8192 ~= Sin[Pi/8] ~= 0.382683432365090 */ \ OD_DCT_OVERFLOW_CHECK(sc, 3135, 4096, 227); \ s3 -= (sc*3135 + 4096) >> 13; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ OD_DCT_OVERFLOW_CHECK(s3, 3259, 8192, 228); \ sc += (s3*3259 + 8192) >> 14; \ /* 13573/32768 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(s5, 13573, 16384, 229); \ sa += (s5*13573 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186547 */ \ OD_DCT_OVERFLOW_CHECK(sa, 11585, 8192, 230); \ s5 -= (sa*11585 + 8192) >> 14; \ /* 13573/32768 ~= Tan[Pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(s5, 13573, 16384, 231); \ sa += (s5*13573 + 16384) >> 15; \ /* 13573/32768 ~= Tan[pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(s9, 13573, 16384, 232); \ s6 += (s9*13573 + 16384) >> 15; \ /* 11585/16384 ~= Sin[pi/4] ~= 0.707106781186547 */ \ OD_DCT_OVERFLOW_CHECK(s6, 11585, 8192, 233); \ s9 -= (s6*11585 + 8192) >> 14; \ /* 13573/32768 ~= Tan[pi/8] ~= 0.414213562373095 */ \ OD_DCT_OVERFLOW_CHECK(s9, 13573, 16384, 234); \ s6 += (s9*13573 + 16384) >> 15; \ sf += se; \ sfh = OD_DCT_RSHIFT(sf, 1); \ se = sfh - se; \ s0 += s1; \ s0h = OD_DCT_RSHIFT(s0, 1); \ s1 = s0h - s1; \ s2 = s3 - s2; \ s2h = OD_DCT_RSHIFT(s2, 1); \ s3 -= s2h; \ sd -= sc; \ sdh = OD_DCT_RSHIFT(sd, 1); \ sc += sdh; \ sa = s4 - sa; \ s4 -= OD_DCT_RSHIFT(sa, 1); \ s5 += sb; \ sb = OD_DCT_RSHIFT(s5, 1) - sb; \ s8 += s6; \ s6 -= OD_DCT_RSHIFT(s8, 1); \ s7 = s9 - s7; \ s9 -= OD_DCT_RSHIFT(s7, 1); \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.820678790828660 */ \ OD_DCT_OVERFLOW_CHECK(sb, 6723, 4096, 235); \ s4 += (sb*6723 + 4096) >> 13; \ /* 16069/16384 ~= Sin[7*Pi/16] ~= 0.980785280403230 */ \ OD_DCT_OVERFLOW_CHECK(s4, 16069, 8192, 236); \ sb -= (s4*16069 + 8192) >> 14; \ /* 6723/8192 ~= Tan[7*Pi/32]) ~= 0.820678790828660 */ \ OD_DCT_OVERFLOW_CHECK(sb, 6723, 4096, 237); \ s4 += (sb*6723 + 4096) >> 13; \ /* 8757/16384 ~= Tan[5*Pi/32]) ~= 0.534511135950792 */ \ OD_DCT_OVERFLOW_CHECK(s5, 8757, 8192, 238); \ sa += (s5*8757 + 8192) >> 14; \ /* 6811/8192 ~= Sin[5*Pi/16] ~= 0.831469612302545 */ \ OD_DCT_OVERFLOW_CHECK(sa, 6811, 4096, 239); \ s5 -= (sa*6811 + 4096) >> 13; \ /* 8757/16384 ~= Tan[5*Pi/32] ~= 0.534511135950792 */ \ OD_DCT_OVERFLOW_CHECK(s5, 8757, 8192, 240); \ sa += (s5*8757 + 8192) >> 14; \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ OD_DCT_OVERFLOW_CHECK(s9, 2485, 4096, 241); \ s6 += (s9*2485 + 4096) >> 13; \ /* 4551/8192 ~= Sin[3*Pi/16] ~= 0.555570233019602 */ \ OD_DCT_OVERFLOW_CHECK(s6, 4551, 4096, 242); \ s9 -= (s6*4551 + 4096) >> 13; \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ OD_DCT_OVERFLOW_CHECK(s9, 2485, 4096, 243); \ s6 += (s9*2485 + 4096) >> 13; \ /* 3227/32768 ~= Tan[Pi/32] ~= 0.09849140335716425 */ \ OD_DCT_OVERFLOW_CHECK(s8, 3227, 16384, 244); \ s7 += (s8*3227 + 16384) >> 15; \ /* 6393/32768 ~= Sin[Pi/16] ~= 0.19509032201612825 */ \ OD_DCT_OVERFLOW_CHECK(s7, 6393, 16384, 245); \ s8 -= (s7*6393 + 16384) >> 15; \ /* 3227/32768 ~= Tan[Pi/32] ~= 0.09849140335716425 */ \ OD_DCT_OVERFLOW_CHECK(s8, 3227, 16384, 246); \ s7 += (s8*3227 + 16384) >> 15; \ s1 -= s2h; \ s2 += s1; \ se += sdh; \ sd = se - sd; \ s3 += sfh; \ sf -= s3; \ sc = s0h - sc; \ s0 -= sc; \ sb += OD_DCT_RSHIFT(s8, 1); \ s8 = sb - s8; \ s4 += OD_DCT_RSHIFT(s7, 1); \ s7 -= s4; \ s6 += OD_DCT_RSHIFT(s5, 1); \ s5 = s6 - s5; \ s9 -= OD_DCT_RSHIFT(sa, 1); \ sa += s9; \ s8 += s0; \ s0 -= OD_DCT_RSHIFT(s8, 1); \ sf += s7; \ s7 = OD_DCT_RSHIFT(sf, 1) - s7; \ s1 -= s6; \ s6 += OD_DCT_RSHIFT(s1, 1); \ s9 += se; \ se = OD_DCT_RSHIFT(s9, 1) - se; \ s2 += sa; \ sa = OD_DCT_RSHIFT(s2, 1) - sa; \ s5 += sd; \ sd -= OD_DCT_RSHIFT(s5, 1); \ s4 = sc - s4; \ sc -= OD_DCT_RSHIFT(s4, 1); \ s3 -= sb; \ sb += OD_DCT_RSHIFT(s3, 1); \ /* 2799/4096 ~= (1/Sqrt[2] - Cos[31*Pi/64]/2)/Sin[31*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(sf, 2799, 2048, 247); \ s0 -= (sf*2799 + 2048) >> 12; \ /* 2893/2048 ~= Sqrt[2]*Sin[31*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s0, 2893, 1024, 248); \ sf += (s0*2893 + 1024) >> 11; \ /* 5397/8192 ~= (Cos[Pi/4] - Cos[31*Pi/64])/Sin[31*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(sf, 5397, 4096, 249); \ s0 -= (sf*5397 + 4096) >> 13; \ /* 41/64 ~= (1/Sqrt[2] - Cos[29*Pi/64]/2)/Sin[29*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s1, 41, 32, 250); \ se += (s1*41 + 32) >> 6; \ /* 2865/2048 ~= Sqrt[2]*Sin[29*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(se, 2865, 1024, 251); \ s1 -= (se*2865 + 1024) >> 11; \ /* 4641/8192 ~= (1/Sqrt[2] - Cos[29*Pi/64])/Sin[29*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s1, 4641, 4096, 252); \ se += (s1*4641 + 4096) >> 13; \ /* 2473/4096 ~= (1/Sqrt[2] - Cos[27*Pi/64]/2)/Sin[27*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s2, 2473, 2048, 253); \ sd += (s2*2473 + 2048) >> 12; \ /* 5619/4096 ~= Sqrt[2]*Sin[27*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(sd, 5619, 2048, 254); \ s2 -= (sd*5619 + 2048) >> 12; \ /* 7839/16384 ~= (1/Sqrt[2] - Cos[27*Pi/64])/Sin[27*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s2, 7839, 8192, 255); \ sd += (s2*7839 + 8192) >> 14; \ /* 5747/8192 ~= (1/Sqrt[2] - Cos[7*Pi/64]/2)/Sin[7*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s3, 5747, 4096, 256); \ sc -= (s3*5747 + 4096) >> 13; \ /* 3903/8192 ~= Sqrt[2]*Sin[7*Pi/64] ~= */ \ OD_DCT_OVERFLOW_CHECK(sc, 3903, 4096, 257); \ s3 += (sc*3903 + 4096) >> 13; \ /* 5701/8192 ~= (1/Sqrt[2] - Cos[7*Pi/64])/Sin[7*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s3, 5701, 4096, 258); \ sc += (s3*5701 + 4096) >> 13; \ /* 4471/8192 ~= (1/Sqrt[2] - Cos[23*Pi/64]/2)/Sin[23*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s4, 4471, 4096, 259); \ sb += (s4*4471 + 4096) >> 13; \ /* 1309/1024 ~= Sqrt[2]*Sin[23*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(sb, 1309, 512, 260); \ s4 -= (sb*1309 + 512) >> 10; \ /* 5067/16384 ~= (1/Sqrt[2] - Cos[23*Pi/64])/Sin[23*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s4, 5067, 8192, 261); \ sb += (s4*5067 + 8192) >> 14; \ /* 2217/4096 ~= (1/Sqrt[2] - Cos[11*Pi/64]/2)/Sin[11*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s5, 2217, 2048, 262); \ sa -= (s5*2217 + 2048) >> 12; \ /* 1489/2048 ~= Sqrt[2]*Sin[11*Pi/64] ~= 0.72705107329128 */ \ OD_DCT_OVERFLOW_CHECK(sa, 1489, 1024, 263); \ s5 += (sa*1489 + 1024) >> 11; \ /* 75/256 ~= (1/Sqrt[2] - Cos[11*Pi/64])/Sin[11*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s5, 75, 128, 264); \ sa += (s5*75 + 128) >> 8; \ /* 2087/4096 ~= (1/Sqrt[2] - Cos[19*Pi/64]/2)/Sin[19*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s9, 2087, 2048, 265); \ s6 -= (s9*2087 + 2048) >> 12; \ /* 4653/4096 ~= Sqrt[2]*Sin[19*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s6, 4653, 2048, 266); \ s9 += (s6*4653 + 2048) >> 12; \ /* 4545/32768 ~= (1/Sqrt[2] - Cos[19*Pi/64])/Sin[19*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s9, 4545, 16384, 267); \ s6 -= (s9*4545 + 16384) >> 15; \ /* 2053/4096 ~= (1/Sqrt[2] - Cos[15*Pi/64]/2)/Sin[15*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s8, 2053, 2048, 268); \ s7 += (s8*2053 + 2048) >> 12; \ /* 1945/2048 ~= Sqrt[2]*Sin[15*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s7, 1945, 1024, 269); \ s8 -= (s7*1945 + 1024) >> 11; \ /* 1651/32768 ~= (1/Sqrt[2] - Cos[15*Pi/64])/Sin[15*Pi/64] */ \ OD_DCT_OVERFLOW_CHECK(s8, 1651, 16384, 270); \ s7 -= (s8*1651 + 16384) >> 15; \ } \ while (0) #define OD_IDST_16(s0, s8, s4, sc, s2, sa, s6, se, \ s1, s9, s5, sd, s3, sb, s7, sf) \ /* Embedded 16-point orthonormal Type-IV iDST. */ \ do { \ int s0h; \ int s4h; \ int sbh; \ int sfh; \ /* 1651/32768 ~= (1/Sqrt[2] - Cos[15*Pi/64])/Sin[15*Pi/64] */ \ se += (s1*1651 + 16384) >> 15; \ /* 1945/2048 ~= Sqrt[2]*Sin[15*Pi/64] */ \ s1 += (se*1945 + 1024) >> 11; \ /* 2053/4096 ~= (1/Sqrt[2] - Cos[15*Pi/64]/2)/Sin[15*Pi/64] */ \ se -= (s1*2053 + 2048) >> 12; \ /* 4545/32768 ~= (1/Sqrt[2] - Cos[19*Pi/64])/Sin[19*Pi/64] */ \ s6 += (s9*4545 + 16384) >> 15; \ /* 4653/32768 ~= Sqrt[2]*Sin[19*Pi/64] */ \ s9 -= (s6*4653 + 2048) >> 12; \ /* 2087/4096 ~= (1/Sqrt[2] - Cos[19*Pi/64]/2)/Sin[19*Pi/64] */ \ s6 += (s9*2087 + 2048) >> 12; \ /* 75/256 ~= (1/Sqrt[2] - Cos[11*Pi/64])/Sin[11*Pi/64] */ \ s5 -= (sa*75 + 128) >> 8; \ /* 1489/2048 ~= Sqrt[2]*Sin[11*Pi/64] */ \ sa -= (s5*1489 + 1024) >> 11; \ /* 2217/4096 ~= (1/Sqrt[2] - Cos[11*Pi/64]/2)/Sin[11*Pi/64] */ \ s5 += (sa*2217 + 2048) >> 12; \ /* 5067/16384 ~= (1/Sqrt[2] - Cos[23*Pi/64])/Sin[23*Pi/64] */ \ sd -= (s2*5067 + 8192) >> 14; \ /* 1309/1024 ~= Sqrt[2]*Sin[23*Pi/64] */ \ s2 += (sd*1309 + 512) >> 10; \ /* 4471/8192 ~= (1/Sqrt[2] - Cos[23*Pi/64]/2)/Sin[23*Pi/64] */ \ sd -= (s2*4471 + 4096) >> 13; \ /* 5701/8192 ~= (1/Sqrt[2] - Cos[7*Pi/64])/Sin[7*Pi/64] */ \ s3 -= (sc*5701 + 4096) >> 13; \ /* 3903/8192 ~= Sqrt[2]*Sin[7*Pi/64] */ \ sc -= (s3*3903 + 4096) >> 13; \ /* 5747/8192 ~= (1/Sqrt[2] - Cos[7*Pi/64]/2)/Sin[7*Pi/64] */ \ s3 += (sc*5747 + 4096) >> 13; \ /* 7839/16384 ~= (1/Sqrt[2] - Cos[27*Pi/64])/Sin[27*Pi/64] */ \ sb -= (s4*7839 + 8192) >> 14; \ /* 5619/4096 ~= Sqrt[2]*Sin[27*Pi/64] */ \ s4 += (sb*5619 + 2048) >> 12; \ /* 2473/4096 ~= (1/Sqrt[2] - Cos[27*Pi/64]/2)/Sin[27*Pi/64] */ \ sb -= (s4*2473 + 2048) >> 12; \ /* 4641/8192 ~= (1/Sqrt[2] - Cos[29*Pi/64])/Sin[29*Pi/64] */ \ s7 -= (s8*4641 + 4096) >> 13; \ /* 2865/2048 ~= Sqrt[2]*Sin[29*Pi/64] */ \ s8 += (s7*2865 + 1024) >> 11; \ /* 41/64 ~= (1/Sqrt[2] - Cos[29*Pi/64]/2)/Sin[29*Pi/64] */ \ s7 -= (s8*41 + 32) >> 6; \ /* 5397/8192 ~= (Cos[Pi/4] - Cos[31*Pi/64])/Sin[31*Pi/64] */ \ s0 += (sf*5397 + 4096) >> 13; \ /* 2893/2048 ~= Sqrt[2]*Sin[31*Pi/64] */ \ sf -= (s0*2893 + 1024) >> 11; \ /* 2799/4096 ~= (1/Sqrt[2] - Cos[31*Pi/64]/2)/Sin[31*Pi/64] */ \ s0 += (sf*2799 + 2048) >> 12; \ sd -= OD_DCT_RSHIFT(sc, 1); \ sc += sd; \ s3 += OD_DCT_RSHIFT(s2, 1); \ s2 = s3 - s2; \ sb += OD_DCT_RSHIFT(sa, 1); \ sa -= sb; \ s5 = OD_DCT_RSHIFT(s4, 1) - s5; \ s4 -= s5; \ s7 = OD_DCT_RSHIFT(s9, 1) - s7; \ s9 -= s7; \ s6 -= OD_DCT_RSHIFT(s8, 1); \ s8 += s6; \ se = OD_DCT_RSHIFT(sf, 1) - se; \ sf -= se; \ s0 += OD_DCT_RSHIFT(s1, 1); \ s1 -= s0; \ s5 -= s9; \ s9 += OD_DCT_RSHIFT(s5, 1); \ sa = s6 - sa; \ s6 -= OD_DCT_RSHIFT(sa, 1); \ se += s2; \ s2 -= OD_DCT_RSHIFT(se, 1); \ s1 = sd - s1; \ sd -= OD_DCT_RSHIFT(s1, 1); \ s0 += s3; \ s0h = OD_DCT_RSHIFT(s0, 1); \ s3 = s0h - s3; \ sf += sc; \ sfh = OD_DCT_RSHIFT(sf, 1); \ sc -= sfh; \ sb = s7 - sb; \ sbh = OD_DCT_RSHIFT(sb, 1); \ s7 -= sbh; \ s4 -= s8; \ s4h = OD_DCT_RSHIFT(s4, 1); \ s8 += s4h; \ /* 3227/32768 ~= Tan[Pi/32] ~= 0.09849140335716425 */ \ se -= (s1*3227 + 16384) >> 15; \ /* 6393/32768 ~= Sin[Pi/16] ~= 0.19509032201612825 */ \ s1 += (se*6393 + 16384) >> 15; \ /* 3227/32768 ~= Tan[Pi/32] ~= 0.09849140335716425 */ \ se -= (s1*3227 + 16384) >> 15; \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ s6 -= (s9*2485 + 4096) >> 13; \ /* 4551/8192 ~= Sin[3*Pi/16] ~= 0.555570233019602 */ \ s9 += (s6*4551 + 4096) >> 13; \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ s6 -= (s9*2485 + 4096) >> 13; \ /* 8757/16384 ~= Tan[5*Pi/32] ~= 0.534511135950792 */ \ s5 -= (sa*8757 + 8192) >> 14; \ /* 6811/8192 ~= Sin[5*Pi/16] ~= 0.831469612302545 */ \ sa += (s5*6811 + 4096) >> 13; \ /* 8757/16384 ~= Tan[5*Pi/32]) ~= 0.534511135950792 */ \ s5 -= (sa*8757 + 8192) >> 14; \ /* 6723/8192 ~= Tan[7*Pi/32]) ~= 0.820678790828660 */ \ s2 -= (sd*6723 + 4096) >> 13; \ /* 16069/16384 ~= Sin[7*Pi/16] ~= 0.980785280403230 */ \ sd += (s2*16069 + 8192) >> 14; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.820678790828660 */ \ s2 -= (sd*6723 + 4096) >> 13; \ s9 += OD_DCT_RSHIFT(se, 1); \ se = s9 - se; \ s6 += OD_DCT_RSHIFT(s1, 1); \ s1 -= s6; \ sd = OD_DCT_RSHIFT(sa, 1) - sd; \ sa -= sd; \ s2 += OD_DCT_RSHIFT(s5, 1); \ s5 = s2 - s5; \ s3 -= sbh; \ sb += s3; \ sc += s4h; \ s4 = sc - s4; \ s8 = s0h - s8; \ s0 -= s8; \ s7 = sfh - s7; \ sf -= s7; \ /* 13573/32768 ~= Tan[pi/8] ~= 0.414213562373095 */ \ s6 -= (s9*13573 + 16384) >> 15; \ /* 11585/16384 ~= Sin[pi/4] ~= 0.707106781186547 */ \ s9 += (s6*11585 + 8192) >> 14; \ /* 13573/32768 ~= Tan[pi/8] ~= 0.414213562373095 */ \ s6 -= (s9*13573 + 16384) >> 15; \ /* 13573/32768 ~= Tan[pi/8] ~= 0.414213562373095 */ \ s5 -= (sa*13573 + 16384) >> 15; \ /* 11585/16384 ~= Sin[pi/4] ~= 0.707106781186547 */ \ sa += (s5*11585 + 8192) >> 14; \ /* 13573/32768 ~= Tan[pi/8] ~= 0.414213562373095 */ \ s5 -= (sa*13573 + 16384) >> 15; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ s3 -= (sc*3259 + 8192) >> 14; \ /* 3135/8192 ~= Sin[Pi/8] ~= 0.382683432365090 */ \ sc += (s3*3135 + 4096) >> 13; \ /* 3259/16384 ~= Tan[Pi/16] ~= 0.198912367379658 */ \ s3 -= (sc*3259 + 8192) >> 14; \ /* 21895/32768 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ sb -= (s4*21895 + 16384) >> 15; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ s4 += (sb*15137 + 8192) >> 14; \ /* 21895/32768 ~= Tan[3*Pi/16] ~= 0.668178637919299 */ \ sb -= (s4*21895 + 16384) >> 15; \ /* 13573/32768 ~= Tan[pi/8] ~= 0.414213562373095 */ \ s8 -= (s7*13573 + 16384) >> 15; \ /* 11585/16384 ~= Sin[pi/4] ~= 0.707106781186547 */ \ s7 += (s8*11585 + 8192) >> 14; \ /* 13573/32768 ~= Tan[pi/8] ~= 0.414213562373095 */ \ s8 -= (s7*13573 + 16384) >> 15; \ } \ while (0) /* TODO: rewrite this to match OD_FDST_16. */ #define OD_FDST_16_ASYM(t0, t0h, t8, t4, t4h, tc, t2, ta, t6, te, \ t1, t9, t5, td, t3, tb, t7, t7h, tf) \ /* Embedded 16-point asymmetric Type-IV fDST. */ \ do { \ int t2h; \ int t3h; \ int t6h; \ int t8h; \ int t9h; \ int tch; \ int tdh; \ /* TODO: Can we move these into another operation */ \ t8 = -t8; \ t9 = -t9; \ ta = -ta; \ tb = -tb; \ td = -td; \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ OD_DCT_OVERFLOW_CHECK(te, 13573, 8192, 136); \ t1 -= (te*13573 + 8192) >> 14; \ /* 11585/32768 ~= Sin[Pi/4]/2 ~= 0.353553390593274 */ \ OD_DCT_OVERFLOW_CHECK(t1, 11585, 16384, 137); \ te += (t1*11585 + 16384) >> 15; \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ OD_DCT_OVERFLOW_CHECK(te, 13573, 8192, 138); \ t1 -= (te*13573 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ OD_DCT_OVERFLOW_CHECK(td, 4161, 8192, 139); \ t2 += (td*4161 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ OD_DCT_OVERFLOW_CHECK(t2, 15137, 8192, 140); \ td -= (t2*15137 + 8192) >> 14; \ /* 14341/16384 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ OD_DCT_OVERFLOW_CHECK(td, 14341, 8192, 141); \ t2 += (td*14341 + 8192) >> 14; \ /* 14341/16384 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ OD_DCT_OVERFLOW_CHECK(t3, 14341, 8192, 142); \ tc -= (t3*14341 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ OD_DCT_OVERFLOW_CHECK(tc, 15137, 8192, 143); \ t3 += (tc*15137 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ OD_DCT_OVERFLOW_CHECK(t3, 4161, 8192, 144); \ tc -= (t3*4161 + 8192) >> 14; \ te = t0h - te; \ t0 -= te; \ tf = OD_DCT_RSHIFT(t1, 1) - tf; \ t1 -= tf; \ /* TODO: Can we move this into another operation */ \ tc = -tc; \ t2 = OD_DCT_RSHIFT(tc, 1) - t2; \ tc -= t2; \ t3 = OD_DCT_RSHIFT(td, 1) - t3; \ td = t3 - td; \ /* 7489/8192 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ OD_DCT_OVERFLOW_CHECK(t6, 7489, 4096, 145); \ t9 -= (t6*7489 + 4096) >> 13; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ OD_DCT_OVERFLOW_CHECK(t9, 11585, 8192, 146); \ t6 += (t9*11585 + 8192) >> 14; \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ OD_DCT_OVERFLOW_CHECK(t6, 19195, 16384, 147); \ t9 += (t6*19195 + 16384) >> 15; \ t8 += OD_DCT_RSHIFT(t9, 1); \ t9 -= t8; \ t6 = t7h - t6; \ t7 -= t6; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.820678790828660 */ \ OD_DCT_OVERFLOW_CHECK(t7, 6723, 4096, 148); \ t8 += (t7*6723 + 4096) >> 13; \ /* 16069/16384 ~= Sin[7*Pi/16] ~= 0.980785280403230 */ \ OD_DCT_OVERFLOW_CHECK(t8, 16069, 8192, 149); \ t7 -= (t8*16069 + 8192) >> 14; \ /* 6723/8192 ~= Tan[7*Pi/32]) ~= 0.820678790828660 */ \ OD_DCT_OVERFLOW_CHECK(t7, 6723, 4096, 150); \ t8 += (t7*6723 + 4096) >> 13; \ /* 17515/32768 ~= Tan[5*Pi/32]) ~= 0.534511135950792 */ \ OD_DCT_OVERFLOW_CHECK(t6, 17515, 16384, 151); \ t9 += (t6*17515 + 16384) >> 15; \ /* 13623/16384 ~= Sin[5*Pi/16] ~= 0.831469612302545 */ \ OD_DCT_OVERFLOW_CHECK(t9, 13623, 8192, 152); \ t6 -= (t9*13623 + 8192) >> 14; \ /* 17515/32768 ~= Tan[5*Pi/32] ~= 0.534511135950792 */ \ OD_DCT_OVERFLOW_CHECK(t6, 17515, 16384, 153); \ t9 += (t6*17515 + 16384) >> 15; \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ OD_DCT_OVERFLOW_CHECK(ta, 13573, 8192, 154); \ t5 += (ta*13573 + 8192) >> 14; \ /* 11585/32768 ~= Sin[Pi/4]/2 ~= 0.353553390593274 */ \ OD_DCT_OVERFLOW_CHECK(t5, 11585, 16384, 155); \ ta -= (t5*11585 + 16384) >> 15; \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ OD_DCT_OVERFLOW_CHECK(ta, 13573, 8192, 156); \ t5 += (ta*13573 + 8192) >> 14; \ tb += OD_DCT_RSHIFT(t5, 1); \ t5 = tb - t5; \ ta += t4h; \ t4 -= ta; \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ OD_DCT_OVERFLOW_CHECK(t5, 2485, 4096, 157); \ ta += (t5*2485 + 4096) >> 13; \ /* 18205/32768 ~= Sin[3*Pi/16] ~= 0.555570233019602 */ \ OD_DCT_OVERFLOW_CHECK(ta, 18205, 16384, 158); \ t5 -= (ta*18205 + 16384) >> 15; \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ OD_DCT_OVERFLOW_CHECK(t5, 2485, 4096, 159); \ ta += (t5*2485 + 4096) >> 13; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.820678790828660 */ \ OD_DCT_OVERFLOW_CHECK(t4, 6723, 4096, 160); \ tb -= (t4*6723 + 4096) >> 13; \ /* 16069/16384 ~= Sin[7*Pi/16] ~= 0.980785280403230 */ \ OD_DCT_OVERFLOW_CHECK(tb, 16069, 8192, 161); \ t4 += (tb*16069 + 8192) >> 14; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.820678790828660 */ \ OD_DCT_OVERFLOW_CHECK(t4, 6723, 4096, 162); \ tb -= (t4*6723 + 4096) >> 13; \ /* TODO: Can we move this into another operation */ \ t5 = -t5; \ tc -= tf; \ tch = OD_DCT_RSHIFT(tc, 1); \ tf += tch; \ t3 += t0; \ t3h = OD_DCT_RSHIFT(t3, 1); \ t0 -= t3h; \ td -= t1; \ tdh = OD_DCT_RSHIFT(td, 1); \ t1 += tdh; \ t2 += te; \ t2h = OD_DCT_RSHIFT(t2, 1); \ te -= t2h; \ t8 += t4; \ t8h = OD_DCT_RSHIFT(t8, 1); \ t4 = t8h - t4; \ t7 = tb - t7; \ t7h = OD_DCT_RSHIFT(t7, 1); \ tb = t7h - tb; \ t6 -= ta; \ t6h = OD_DCT_RSHIFT(t6, 1); \ ta += t6h; \ t9 = t5 - t9; \ t9h = OD_DCT_RSHIFT(t9, 1); \ t5 -= t9h; \ t0 -= t7h; \ t7 += t0; \ tf += t8h; \ t8 -= tf; \ te -= t6h; \ t6 += te; \ t1 += t9h; \ t9 -= t1; \ tb -= tch; \ tc += tb; \ t4 += t3h; \ t3 -= t4; \ ta -= tdh; \ td += ta; \ t5 = t2h - t5; \ t2 -= t5; \ /* TODO: Can we move these into another operation */ \ t8 = -t8; \ t9 = -t9; \ ta = -ta; \ tb = -tb; \ tc = -tc; \ td = -td; \ tf = -tf; \ /* 7799/8192 ~= Tan[31*Pi/128] ~= 0.952079146700925 */ \ OD_DCT_OVERFLOW_CHECK(tf, 7799, 4096, 163); \ t0 -= (tf*7799 + 4096) >> 13; \ /* 4091/4096 ~= Sin[31*Pi/64] ~= 0.998795456205172 */ \ OD_DCT_OVERFLOW_CHECK(t0, 4091, 2048, 164); \ tf += (t0*4091 + 2048) >> 12; \ /* 7799/8192 ~= Tan[31*Pi/128] ~= 0.952079146700925 */ \ OD_DCT_OVERFLOW_CHECK(tf, 7799, 4096, 165); \ t0 -= (tf*7799 + 4096) >> 13; \ /* 2417/32768 ~= Tan[3*Pi/128] ~= 0.0737644315224493 */ \ OD_DCT_OVERFLOW_CHECK(te, 2417, 16384, 166); \ t1 += (te*2417 + 16384) >> 15; \ /* 601/4096 ~= Sin[3*Pi/64] ~= 0.146730474455362 */ \ OD_DCT_OVERFLOW_CHECK(t1, 601, 2048, 167); \ te -= (t1*601 + 2048) >> 12; \ /* 2417/32768 ~= Tan[3*Pi/128] ~= 0.0737644315224493 */ \ OD_DCT_OVERFLOW_CHECK(te, 2417, 16384, 168); \ t1 += (te*2417 + 16384) >> 15; \ /* 14525/32768 ~= Tan[17*Pi/128] ~= 0.443269513890864 */ \ OD_DCT_OVERFLOW_CHECK(t8, 14525, 16384, 169); \ t7 -= (t8*14525 + 16384) >> 15; \ /* 3035/4096 ~= Sin[17*Pi/64] ~= 0.740951125354959 */ \ OD_DCT_OVERFLOW_CHECK(t7, 3035, 2048, 170); \ t8 += (t7*3035 + 2048) >> 12; \ /* 7263/16384 ~= Tan[17*Pi/128] ~= 0.443269513890864 */ \ OD_DCT_OVERFLOW_CHECK(t8, 7263, 8192, 171); \ t7 -= (t8*7263 + 8192) >> 14; \ /* 6393/8192 ~= Tan[27*Pi/128] ~= 0.780407659653944 */ \ OD_DCT_OVERFLOW_CHECK(td, 6393, 4096, 172); \ t2 -= (td*6393 + 4096) >> 13; \ /* 3973/4096 ~= Sin[27*Pi/64] ~= 0.970031253194544 */ \ OD_DCT_OVERFLOW_CHECK(t2, 3973, 2048, 173); \ td += (t2*3973 + 2048) >> 12; \ /* 6393/8192 ~= Tan[27*Pi/128] ~= 0.780407659653944 */ \ OD_DCT_OVERFLOW_CHECK(td, 6393, 4096, 174); \ t2 -= (td*6393 + 4096) >> 13; \ /* 9281/16384 ~= Tan[21*Pi/128] ~= 0.566493002730344 */ \ OD_DCT_OVERFLOW_CHECK(ta, 9281, 8192, 175); \ t5 -= (ta*9281 + 8192) >> 14; \ /* 7027/8192 ~= Sin[21*Pi/64] ~= 0.857728610000272 */ \ OD_DCT_OVERFLOW_CHECK(t5, 7027, 4096, 176); \ ta += (t5*7027 + 4096) >> 13; \ /* 9281/16384 ~= Tan[21*Pi/128] ~= 0.566493002730344 */ \ OD_DCT_OVERFLOW_CHECK(ta, 9281, 8192, 177); \ t5 -= (ta*9281 + 8192) >> 14; \ /* 11539/16384 ~= Tan[25*Pi/128] ~= 0.704279460865044 */ \ OD_DCT_OVERFLOW_CHECK(tc, 11539, 8192, 178); \ t3 -= (tc*11539 + 8192) >> 14; \ /* 7713/8192 ~= Sin[25*Pi/64] ~= 0.941544065183021 */ \ OD_DCT_OVERFLOW_CHECK(t3, 7713, 4096, 179); \ tc += (t3*7713 + 4096) >> 13; \ /* 11539/16384 ~= Tan[25*Pi/128] ~= 0.704279460865044 */ \ OD_DCT_OVERFLOW_CHECK(tc, 11539, 8192, 180); \ t3 -= (tc*11539 + 8192) >> 14; \ /* 10375/16384 ~= Tan[23*Pi/128] ~= 0.633243016177569 */ \ OD_DCT_OVERFLOW_CHECK(tb, 10375, 8192, 181); \ t4 -= (tb*10375 + 8192) >> 14; \ /* 7405/8192 ~= Sin[23*Pi/64] ~= 0.903989293123443 */ \ OD_DCT_OVERFLOW_CHECK(t4, 7405, 4096, 182); \ tb += (t4*7405 + 4096) >> 13; \ /* 10375/16384 ~= Tan[23*Pi/128] ~= 0.633243016177569 */ \ OD_DCT_OVERFLOW_CHECK(tb, 10375, 8192, 183); \ t4 -= (tb*10375 + 8192) >> 14; \ /* 8247/16384 ~= Tan[19*Pi/128] ~= 0.503357699799294 */ \ OD_DCT_OVERFLOW_CHECK(t9, 8247, 8192, 184); \ t6 -= (t9*8247 + 8192) >> 14; \ /* 1645/2048 ~= Sin[19*Pi/64] ~= 0.803207531480645 */ \ OD_DCT_OVERFLOW_CHECK(t6, 1645, 1024, 185); \ t9 += (t6*1645 + 1024) >> 11; \ /* 8247/16384 ~= Tan[19*Pi/128] ~= 0.503357699799294 */ \ OD_DCT_OVERFLOW_CHECK(t9, 8247, 8192, 186); \ t6 -= (t9*8247 + 8192) >> 14; \ } \ while (0) #define OD_IDST_16_ASYM(t0, t0h, t8, t4, tc, t2, t2h, ta, t6, te, teh, \ t1, t9, t5, td, t3, tb, t7, tf) \ /* Embedded 16-point asymmetric Type-IV iDST. */ \ do { \ int t1h_; \ int t3h_; \ int t4h; \ int t6h; \ int t9h_; \ int tbh_; \ int tch; \ /* 8247/16384 ~= Tan[19*Pi/128] ~= 0.503357699799294 */ \ t6 += (t9*8247 + 8192) >> 14; \ /* 1645/2048 ~= Sin[19*Pi/64] ~= 0.803207531480645 */ \ t9 -= (t6*1645 + 1024) >> 11; \ /* 8247/16384 ~= Tan[19*Pi/128] ~= 0.503357699799294 */ \ t6 += (t9*8247 + 8192) >> 14; \ /* 10375/16384 ~= Tan[23*Pi/128] ~= 0.633243016177569 */ \ t2 += (td*10375 + 8192) >> 14; \ /* 7405/8192 ~= Sin[23*Pi/64] ~= 0.903989293123443 */ \ td -= (t2*7405 + 4096) >> 13; \ /* 10375/16384 ~= Tan[23*Pi/128] ~= 0.633243016177569 */ \ t2 += (td*10375 + 8192) >> 14; \ /* 11539/16384 ~= Tan[25*Pi/128] ~= 0.704279460865044 */ \ tc += (t3*11539 + 8192) >> 14; \ /* 7713/8192 ~= Sin[25*Pi/64] ~= 0.941544065183021 */ \ t3 -= (tc*7713 + 4096) >> 13; \ /* 11539/16384 ~= Tan[25*Pi/128] ~= 0.704279460865044 */ \ tc += (t3*11539 + 8192) >> 14; \ /* 9281/16384 ~= Tan[21*Pi/128] ~= 0.566493002730344 */ \ ta += (t5*9281 + 8192) >> 14; \ /* 7027/8192 ~= Sin[21*Pi/64] ~= 0.857728610000272 */ \ t5 -= (ta*7027 + 4096) >> 13; \ /* 9281/16384 ~= Tan[21*Pi/128] ~= 0.566493002730344 */ \ ta += (t5*9281 + 8192) >> 14; \ /* 6393/8192 ~= Tan[27*Pi/128] ~= 0.780407659653944 */ \ t4 += (tb*6393 + 4096) >> 13; \ /* 3973/4096 ~= Sin[27*Pi/64] ~= 0.970031253194544 */ \ tb -= (t4*3973 + 2048) >> 12; \ /* 6393/8192 ~= Tan[27*Pi/128] ~= 0.780407659653944 */ \ t4 += (tb*6393 + 4096) >> 13; \ /* 7263/16384 ~= Tan[17*Pi/128] ~= 0.443269513890864 */ \ te += (t1*7263 + 8192) >> 14; \ /* 3035/4096 ~= Sin[17*Pi/64] ~= 0.740951125354959 */ \ t1 -= (te*3035 + 2048) >> 12; \ /* 14525/32768 ~= Tan[17*Pi/128] ~= 0.443269513890864 */ \ te += (t1*14525 + 16384) >> 15; \ /* 2417/32768 ~= Tan[3*Pi/128] ~= 0.0737644315224493 */ \ t8 -= (t7*2417 + 16384) >> 15; \ /* 601/4096 ~= Sin[3*Pi/64] ~= 0.146730474455362 */ \ t7 += (t8*601 + 2048) >> 12; \ /* 2417/32768 ~= Tan[3*Pi/128] ~= 0.0737644315224493 */ \ t8 -= (t7*2417 + 16384) >> 15; \ /* 7799/8192 ~= Tan[31*Pi/128] ~= 0.952079146700925 */ \ t0 += (tf*7799 + 4096) >> 13; \ /* 4091/4096 ~= Sin[31*Pi/64] ~= 0.998795456205172 */ \ tf -= (t0*4091 + 2048) >> 12; \ /* 7799/8192 ~= Tan[31*Pi/128] ~= 0.952079146700925 */ \ t0 += (tf*7799 + 4096) >> 13; \ /* TODO: Can we move these into another operation */ \ t1 = -t1; \ t3 = -t3; \ t5 = -t5; \ t9 = -t9; \ tb = -tb; \ td = -td; \ tf = -tf; \ t4 += ta; \ t4h = OD_DCT_RSHIFT(t4, 1); \ ta = t4h - ta; \ tb -= t5; \ tbh_ = OD_DCT_RSHIFT(tb, 1); \ t5 += tbh_; \ tc += t2; \ tch = OD_DCT_RSHIFT(tc, 1); \ t2 -= tch; \ t3 -= td; \ t3h_ = OD_DCT_RSHIFT(t3, 1); \ td += t3h_; \ t9 += t8; \ t9h_ = OD_DCT_RSHIFT(t9, 1); \ t8 -= t9h_; \ t6 -= t7; \ t6h = OD_DCT_RSHIFT(t6, 1); \ t7 += t6h; \ t1 += tf; \ t1h_ = OD_DCT_RSHIFT(t1, 1); \ tf -= t1h_; \ te -= t0; \ teh = OD_DCT_RSHIFT(te, 1); \ t0 += teh; \ ta += t9h_; \ t9 = ta - t9; \ t5 -= t6h; \ t6 += t5; \ td = teh - td; \ te = td - te; \ t2 = t1h_ - t2; \ t1 -= t2; \ t7 += t4h; \ t4 -= t7; \ t8 -= tbh_; \ tb += t8; \ t0 += tch; \ tc -= t0; \ tf -= t3h_; \ t3 += tf; \ /* TODO: Can we move this into another operation */ \ ta = -ta; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.820678790828660 */ \ td += (t2*6723 + 4096) >> 13; \ /* 16069/16384 ~= Sin[7*Pi/16] ~= 0.980785280403230 */ \ t2 -= (td*16069 + 8192) >> 14; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.820678790828660 */ \ td += (t2*6723 + 4096) >> 13; \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ t5 -= (ta*2485 + 4096) >> 13; \ /* 18205/32768 ~= Sin[3*Pi/16] ~= 0.555570233019602 */ \ ta += (t5*18205 + 16384) >> 15; \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ t5 -= (ta*2485 + 4096) >> 13; \ t2 += t5; \ t2h = OD_DCT_RSHIFT(t2, 1); \ t5 -= t2h; \ ta = td - ta; \ td -= OD_DCT_RSHIFT(ta, 1); \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ ta -= (t5*13573 + 8192) >> 14; \ /* 11585/32768 ~= Sin[Pi/4]/2 ~= 0.353553390593274 */ \ t5 += (ta*11585 + 16384) >> 15; \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ ta -= (t5*13573 + 8192) >> 14; \ /* 17515/32768 ~= Tan[5*Pi/32] ~= 0.534511135950792 */ \ t9 -= (t6*17515 + 16384) >> 15; \ /* 13623/16384 ~= Sin[5*Pi/16] ~= 0.831469612302545 */ \ t6 += (t9*13623 + 8192) >> 14; \ /* 17515/32768 ~= Tan[5*Pi/32]) ~= 0.534511135950792 */ \ t9 -= (t6*17515 + 16384) >> 15; \ /* 6723/8192 ~= Tan[7*Pi/32]) ~= 0.820678790828660 */ \ t1 -= (te*6723 + 4096) >> 13; \ /* 16069/16384 ~= Sin[7*Pi/16] ~= 0.980785280403230 */ \ te += (t1*16069 + 8192) >> 14; \ /* 6723/8192 ~= Tan[7*Pi/32]) ~= 0.820678790828660 */ \ t1 -= (te*6723 + 4096) >> 13; \ te += t6; \ teh = OD_DCT_RSHIFT(te, 1); \ t6 = teh - t6; \ t9 += t1; \ t1 -= OD_DCT_RSHIFT(t9, 1); \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ t9 -= (t6*19195 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ t6 -= (t9*11585 + 8192) >> 14; \ /* 7489/8192 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ t9 += (t6*7489 + 4096) >> 13; \ tb = tc - tb; \ tc = OD_DCT_RSHIFT(tb, 1) - tc; \ t3 += t4; \ t4 = OD_DCT_RSHIFT(t3, 1) - t4; \ /* TODO: Can we move this into another operation */ \ t3 = -t3; \ t8 += tf; \ tf = OD_DCT_RSHIFT(t8, 1) - tf; \ t0 += t7; \ t0h = OD_DCT_RSHIFT(t0, 1); \ t7 = t0h - t7; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ t3 += (tc*4161 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ tc -= (t3*15137 + 8192) >> 14; \ /* 14341/16384 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ t3 += (tc*14341 + 8192) >> 14; \ /* 14341/16384 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ t4 -= (tb*14341 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ tb += (t4*15137 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ t4 -= (tb*4161 + 8192) >> 14; \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ t8 += (t7*13573 + 8192) >> 14; \ /* 11585/32768 ~= Sin[Pi/4]/2 ~= 0.353553390593274 */ \ t7 -= (t8*11585 + 16384) >> 15; \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ t8 += (t7*13573 + 8192) >> 14; \ /* TODO: Can we move these into another operation */ \ t1 = -t1; \ t5 = -t5; \ t9 = -t9; \ tb = -tb; \ td = -td; \ } \ while (0) #define OD_FDCT_32(t0, tg, t8, to, t4, tk, tc, ts, t2, ti, ta, tq, t6, tm, \ te, tu, t1, th, t9, tp, t5, tl, td, tt, t3, tj, tb, tr, t7, tn, tf, tv) \ /* Embedded 32-point orthonormal Type-II fDCT. */ \ do { \ int tgh; \ int thh; \ int tih; \ int tkh; \ int tmh; \ int tnh; \ int toh; \ int tqh; \ int tsh; \ int tuh; \ int tvh; \ tv = t0 - tv; \ tvh = OD_DCT_RSHIFT(tv, 1); \ t0 -= tvh; \ tu += t1; \ tuh = OD_DCT_RSHIFT(tu, 1); \ t1 = tuh - t1; \ tt = t2 - tt; \ t2 -= OD_DCT_RSHIFT(tt, 1); \ ts += t3; \ tsh = OD_DCT_RSHIFT(ts, 1); \ t3 = tsh - t3; \ tr = t4 - tr; \ t4 -= OD_DCT_RSHIFT(tr, 1); \ tq += t5; \ tqh = OD_DCT_RSHIFT(tq, 1); \ t5 = tqh - t5; \ tp = t6 - tp; \ t6 -= OD_DCT_RSHIFT(tp, 1); \ to += t7; \ toh = OD_DCT_RSHIFT(to, 1); \ t7 = toh - t7; \ tn = t8 - tn; \ tnh = OD_DCT_RSHIFT(tn, 1); \ t8 -= tnh; \ tm += t9; \ tmh = OD_DCT_RSHIFT(tm, 1); \ t9 = tmh - t9; \ tl = ta - tl; \ ta -= OD_DCT_RSHIFT(tl, 1); \ tk += tb; \ tkh = OD_DCT_RSHIFT(tk, 1); \ tb = tkh - tb; \ tj = tc - tj; \ tc -= OD_DCT_RSHIFT(tj, 1); \ ti += td; \ tih = OD_DCT_RSHIFT(ti, 1); \ td = tih - td; \ th = te - th; \ thh = OD_DCT_RSHIFT(th, 1); \ te -= thh; \ tg += tf; \ tgh = OD_DCT_RSHIFT(tg, 1); \ tf = tgh - tf; \ OD_FDCT_16_ASYM(t0, tg, tgh, t8, to, toh, t4, tk, tkh, tc, ts, tsh, \ t2, ti, tih, ta, tq, tqh, t6, tm, tmh, te, tu, tuh); \ OD_FDST_16_ASYM(tv, tvh, tf, tn, tnh, t7, tr, tb, tj, t3, \ tt, td, tl, t5, tp, t9, th, thh, t1); \ } \ while (0) #define OD_IDCT_32(t0, tg, t8, to, t4, tk, tc, ts, t2, ti, ta, tq, t6, tm, \ te, tu, t1, th, t9, tp, t5, tl, td, tt, t3, tj, tb, tr, t7, tn, tf, tv) \ /* Embedded 32-point orthonormal Type-II iDCT. */ \ do { \ int t1h; \ int t3h; \ int t5h; \ int t7h; \ int t9h; \ int tbh; \ int tdh; \ int tfh; \ int thh; \ int tth; \ int tvh; \ OD_IDST_16_ASYM(tv, tvh, tn, tr, tj, tt, tth, tl, tp, th, thh, \ tu, tm, tq, ti, ts, tk, to, tg); \ OD_IDCT_16_ASYM(t0, t8, t4, tc, t2, ta, t6, te, \ t1, t1h, t9, t9h, t5, t5h, td, tdh, t3, t3h, tb, tbh, t7, t7h, tf, tfh); \ tu = t1h - tu; \ t1 -= tu; \ te += thh; \ th = te - th; \ tm = t9h - tm; \ t9 -= tm; \ t6 += OD_DCT_RSHIFT(tp, 1); \ tp = t6 - tp; \ tq = t5h - tq; \ t5 -= tq; \ ta += OD_DCT_RSHIFT(tl, 1); \ tl = ta - tl; \ ti = tdh - ti; \ td -= ti; \ t2 += tth; \ tt = t2 - tt; \ ts = t3h - ts; \ t3 -= ts; \ tc += OD_DCT_RSHIFT(tj, 1); \ tj = tc - tj; \ tk = tbh - tk; \ tb -= tk; \ t4 += OD_DCT_RSHIFT(tr, 1); \ tr = t4 - tr; \ to = t7h - to; \ t7 -= to; \ t8 += OD_DCT_RSHIFT(tn, 1); \ tn = t8 - tn; \ tg = tfh - tg; \ tf -= tg; \ t0 += tvh; \ tv = t0 - tv; \ } \ while (0) #if CONFIG_TX64X64 #define OD_FDCT_32_ASYM(t0, tg, tgh, t8, to, toh, t4, tk, tkh, tc, ts, tsh, \ t2, ti, tih, ta, tq, tqh, t6, tm, tmh, te, tu, tuh, t1, th, thh, \ t9, tp, tph, t5, tl, tlh, td, tt, tth, t3, tj, tjh, tb, tr, trh, \ t7, tn, tnh, tf, tv, tvh) \ /* Embedded 32-point asymmetric Type-II fDCT. */ \ do { \ t0 += tvh; \ tv = t0 - tv; \ t1 = tuh - t1; \ tu -= t1; \ t2 += tth; \ tt = t2 - tt; \ t3 = tsh - t3; \ ts -= t3; \ t4 += trh; \ tr = t4 - tr; \ t5 = tqh - t5; \ tq -= t5; \ t6 += tph; \ tp = t6 - tp; \ t7 = toh - t7; \ to -= t7; \ t8 += tnh; \ tn = t8 - tn; \ t9 = tmh - t9; \ tm -= t9; \ ta += tlh; \ tl = ta - tl; \ tb = tkh - tb; \ tk -= tb; \ tc += tjh; \ tj = tc - tj; \ td = tih - td; \ ti -= td; \ te += thh; \ th = te - th; \ tf = tgh - tf; \ tg -= tf; \ OD_FDCT_16(t0, tg, t8, to, t4, tk, tc, ts, \ t2, ti, ta, tq, t6, tm, te, tu); \ OD_FDST_16(tv, tf, tn, t7, tr, tb, tj, t3, \ tt, td, tl, t5, tp, t9, th, t1); \ } \ while (0) #define OD_IDCT_32_ASYM(t0, tg, t8, to, t4, tk, tc, ts, t2, ti, ta, tq, \ t6, tm, te, tu, t1, t1h, th, thh, t9, t9h, tp, tph, t5, t5h, tl, tlh, \ td, tdh, tt, tth, t3, t3h, tj, tjh, tb, tbh, tr, trh, t7, t7h, tn, tnh, \ tf, tfh, tv, tvh) \ /* Embedded 32-point asymmetric Type-II iDCT. */ \ do { \ OD_IDST_16(tv, tn, tr, tj, tt, tl, tp, th, \ tu, tm, tq, ti, ts, tk, to, tg); \ OD_IDCT_16(t0, t8, t4, tc, t2, ta, t6, te, \ t1, t9, t5, td, t3, tb, t7, tf); \ tv = t0 - tv; \ tvh = OD_DCT_RSHIFT(tv, 1); \ t0 -= tvh; \ t1 += tu; \ t1h = OD_DCT_RSHIFT(t1, 1); \ tu = t1h - tu; \ tt = t2 - tt; \ tth = OD_DCT_RSHIFT(tt, 1); \ t2 -= tth; \ t3 += ts; \ t3h = OD_DCT_RSHIFT(t3, 1); \ ts = t3h - ts; \ tr = t4 - tr; \ trh = OD_DCT_RSHIFT(tr, 1); \ t4 -= trh; \ t5 += tq; \ t5h = OD_DCT_RSHIFT(t5, 1); \ tq = t5h - tq; \ tp = t6 - tp; \ tph = OD_DCT_RSHIFT(tp, 1); \ t6 -= tph; \ t7 += to; \ t7h = OD_DCT_RSHIFT(t7, 1); \ to = t7h - to; \ tn = t8 - tn; \ tnh = OD_DCT_RSHIFT(tn, 1); \ t8 -= tnh; \ t9 += tm; \ t9h = OD_DCT_RSHIFT(t9, 1); \ tm = t9h - tm; \ tl = ta - tl; \ tlh = OD_DCT_RSHIFT(tl, 1); \ ta -= tlh; \ tb += tk; \ tbh = OD_DCT_RSHIFT(tb, 1); \ tk = tbh - tk; \ tj = tc - tj; \ tjh = OD_DCT_RSHIFT(tj, 1); \ tc -= tjh; \ td += ti; \ tdh = OD_DCT_RSHIFT(td, 1); \ ti = tdh - ti; \ th = te - th; \ thh = OD_DCT_RSHIFT(th, 1); \ te -= thh; \ tf += tg; \ tfh = OD_DCT_RSHIFT(tf, 1); \ tg = tfh - tg; \ } \ while (0) #define OD_FDST_32_ASYM(t0, tg, t8, to, t4, tk, tc, ts, t2, ti, ta, tq, t6, \ tm, te, tu, t1, th, t9, tp, t5, tl, td, tt, t3, tj, tb, tr, t7, tn, tf, tv) \ /* Embedded 32-point asymmetric Type-IV fDST. */ \ do { \ int t0h; \ int t1h; \ int t4h; \ int t5h; \ int tqh; \ int trh; \ int tuh; \ int tvh; \ \ tu = -tu; \ \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ OD_DCT_OVERFLOW_CHECK(tq, 13573, 8192, 271); \ t5 -= (tq*13573 + 8192) >> 14; \ /* 11585/32768 ~= Sin[Pi/4]/2 ~= 0.353553390593274 */ \ OD_DCT_OVERFLOW_CHECK(t5, 11585, 16384, 272); \ tq += (t5*11585 + 16384) >> 15; \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ OD_DCT_OVERFLOW_CHECK(tq, 13573, 8192, 273); \ t5 -= (tq*13573 + 8192) >> 14; \ /* 29957/32768 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ OD_DCT_OVERFLOW_CHECK(t6, 29957, 16384, 274); \ tp += (t6*29957 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ OD_DCT_OVERFLOW_CHECK(tp, 11585, 8192, 275); \ t6 -= (tp*11585 + 8192) >> 14; \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ OD_DCT_OVERFLOW_CHECK(t6, 19195, 16384, 276); \ tp -= (t6*19195 + 16384) >> 15; \ /* 29957/32768 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ OD_DCT_OVERFLOW_CHECK(t1, 29957, 16384, 277); \ tu += (t1*29957 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ OD_DCT_OVERFLOW_CHECK(tu, 11585, 8192, 278); \ t1 -= (tu*11585 + 8192) >> 14; \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ OD_DCT_OVERFLOW_CHECK(t1, 19195, 16384, 279); \ tu -= (t1*19195 + 16384) >> 15; \ /* 28681/32768 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ OD_DCT_OVERFLOW_CHECK(t2, 28681, 16384, 280); \ tt += (t2*28681 + 16384) >> 15; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ OD_DCT_OVERFLOW_CHECK(tt, 15137, 8192, 281); \ t2 -= (tt*15137 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ OD_DCT_OVERFLOW_CHECK(t2, 4161, 8192, 282); \ tt += (t2*4161 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ OD_DCT_OVERFLOW_CHECK(ts, 4161, 8192, 283); \ t3 += (ts*4161 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ OD_DCT_OVERFLOW_CHECK(t3, 15137, 8192, 284); \ ts -= (t3*15137 + 8192) >> 14; \ /* 14341/16384 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ OD_DCT_OVERFLOW_CHECK(ts, 14341, 8192, 285); \ t3 += (ts*14341 + 8192) >> 14; \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ OD_DCT_OVERFLOW_CHECK(tm, 19195, 16384, 286); \ t9 -= (tm*19195 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ OD_DCT_OVERFLOW_CHECK(t9, 11585, 8192, 287); \ tm -= (t9*11585 + 8192) >> 14; \ /* 7489/8192 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ OD_DCT_OVERFLOW_CHECK(tm, 7489, 4096, 288); \ t9 += (tm*7489 + 4096) >> 13; \ /* 3259/8192 ~= 2*Tan[Pi/16] ~= 0.397824734759316 */ \ OD_DCT_OVERFLOW_CHECK(tl, 3259, 4096, 289); \ ta += (tl*3259 + 4096) >> 13; \ /* 3135/16384 ~= Sin[Pi/8]/2 ~= 0.1913417161825449 */ \ OD_DCT_OVERFLOW_CHECK(ta, 3135, 8192, 290); \ tl -= (ta*3135 + 8192) >> 14; \ /* 3259/8192 ~= 2*Tan[Pi/16] ~= 0.397824734759316 */ \ OD_DCT_OVERFLOW_CHECK(tl, 3259, 4096, 291); \ ta += (tl*3259 + 4096) >> 13; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ OD_DCT_OVERFLOW_CHECK(tk, 4161, 8192, 292); \ tb += (tk*4161 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ OD_DCT_OVERFLOW_CHECK(tb, 15137, 8192, 293); \ tk -= (tb*15137 + 8192) >> 14; \ /* 14341/16384 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ OD_DCT_OVERFLOW_CHECK(tk, 14341, 8192, 294); \ tb += (tk*14341 + 8192) >> 14; \ /* 29957/32768 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ OD_DCT_OVERFLOW_CHECK(te, 29957, 16384, 295); \ th += (te*29957 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ OD_DCT_OVERFLOW_CHECK(th, 11585, 8192, 296); \ te -= (th*11585 + 8192) >> 14; \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ OD_DCT_OVERFLOW_CHECK(te, 19195, 16384, 297); \ th -= (te*19195 + 16384) >> 15; \ /* 28681/32768 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ OD_DCT_OVERFLOW_CHECK(tc, 28681, 16384, 298); \ tj += (tc*28681 + 16384) >> 15; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ OD_DCT_OVERFLOW_CHECK(tj, 15137, 8192, 299); \ tc -= (tj*15137 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ OD_DCT_OVERFLOW_CHECK(tc, 4161, 8192, 300); \ tj += (tc*4161 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ OD_DCT_OVERFLOW_CHECK(ti, 4161, 8192, 301); \ td += (ti*4161 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ OD_DCT_OVERFLOW_CHECK(td, 15137, 8192, 302); \ ti -= (td*15137 + 8192) >> 14; \ /* 14341/16384 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ OD_DCT_OVERFLOW_CHECK(ti, 14341, 8192, 303); \ td += (ti*14341 + 8192) >> 14; \ \ t1 = -t1; \ t2 = -t2; \ t3 = -t3; \ td = -td; \ tg = -tg; \ to = -to; \ ts = -ts; \ \ tr -= OD_DCT_RSHIFT(t5, 1); \ t5 += tr; \ tq -= OD_DCT_RSHIFT(t4, 1); /* pass */ \ t4 += tq; \ t6 -= OD_DCT_RSHIFT(t7, 1); \ t7 += t6; \ to -= OD_DCT_RSHIFT(tp, 1); /* pass */ \ tp += to; \ t1 += OD_DCT_RSHIFT(t0, 1); /* pass */ \ t0 -= t1; \ tv -= OD_DCT_RSHIFT(tu, 1); \ tu += tv; \ t3 -= OD_DCT_RSHIFT(tt, 1); \ tt += t3; \ t2 += OD_DCT_RSHIFT(ts, 1); \ ts -= t2; \ t9 -= OD_DCT_RSHIFT(t8, 1); /* pass */ \ t8 += t9; \ tn += OD_DCT_RSHIFT(tm, 1); \ tm -= tn; \ tb += OD_DCT_RSHIFT(ta, 1); \ ta -= tb; \ tl -= OD_DCT_RSHIFT(tk, 1); \ tk += tl; \ te -= OD_DCT_RSHIFT(tf, 1); /* pass */ \ tf += te; \ tg -= OD_DCT_RSHIFT(th, 1); \ th += tg; \ tc -= OD_DCT_RSHIFT(ti, 1); \ ti += tc; \ td += OD_DCT_RSHIFT(tj, 1); \ tj -= td; \ \ t4 = -t4; \ \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.8206787908286602 */ \ OD_DCT_OVERFLOW_CHECK(tr, 6723, 4096, 304); \ t4 += (tr*6723 + 4096) >> 13; \ /* 16069/16384 ~= Sin[7*Pi/16] ~= 0.9807852804032304 */ \ OD_DCT_OVERFLOW_CHECK(t4, 16069, 8192, 305); \ tr -= (t4*16069 + 8192) >> 14; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.8206787908286602 */ \ OD_DCT_OVERFLOW_CHECK(tr, 6723, 4096, 306); \ t4 += (tr*6723 + 4096) >> 13; \ /* 17515/32768 ~= Tan[5*Pi/32] ~= 0.5345111359507916 */ \ OD_DCT_OVERFLOW_CHECK(tq, 17515, 16384, 307); \ t5 += (tq*17515 + 16384) >> 15; \ /* 13623/16384 ~= Sin[5*Pi/16] ~= 0.8314696123025452 */ \ OD_DCT_OVERFLOW_CHECK(t5, 13623, 8192, 308); \ tq -= (t5*13623 + 8192) >> 14; \ /* 17515/32768 ~= Tan[5*Pi/32] ~= 0.5345111359507916 */ \ OD_DCT_OVERFLOW_CHECK(tq, 17515, 16384, 309); \ t5 += (tq*17515 + 16384) >> 15; \ /* 3227/32768 ~= Tan[Pi/32] ~= 0.09849140335716425 */ \ OD_DCT_OVERFLOW_CHECK(to, 3227, 16384, 310); \ t7 += (to*3227 + 16384) >> 15; \ /* 6393/32768 ~= Sin[Pi/16] ~= 0.19509032201612825 */ \ OD_DCT_OVERFLOW_CHECK(t7, 6393, 16384, 311); \ to -= (t7*6393 + 16384) >> 15; \ /* 3227/32768 ~= Tan[Pi/32] ~= 0.09849140335716425 */ \ OD_DCT_OVERFLOW_CHECK(to, 3227, 16384, 312); \ t7 += (to*3227 + 16384) >> 15; \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ OD_DCT_OVERFLOW_CHECK(tp, 2485, 4096, 313); \ t6 += (tp*2485 + 4096) >> 13; \ /* 18205/32768 ~= Sin[3*Pi/16] ~= 0.555570233019602 */ \ OD_DCT_OVERFLOW_CHECK(t6, 18205, 16384, 314); \ tp -= (t6*18205 + 16384) >> 15; \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ OD_DCT_OVERFLOW_CHECK(tp, 2485, 4096, 315); \ t6 += (tp*2485 + 4096) >> 13; \ \ t5 = -t5; \ \ tr += to; \ trh = OD_DCT_RSHIFT(tr, 1); \ to -= trh; \ t4 += t7; \ t4h = OD_DCT_RSHIFT(t4, 1); \ t7 -= t4h; \ t5 += tp; \ t5h = OD_DCT_RSHIFT(t5, 1); \ tp -= t5h; \ tq += t6; \ tqh = OD_DCT_RSHIFT(tq, 1); \ t6 -= tqh; \ t0 -= t3; \ t0h = OD_DCT_RSHIFT(t0, 1); \ t3 += t0h; \ tv -= ts; \ tvh = OD_DCT_RSHIFT(tv, 1); \ ts += tvh; \ tu += tt; \ tuh = OD_DCT_RSHIFT(tu, 1); \ tt -= tuh; \ t1 -= t2; \ t1h = OD_DCT_RSHIFT(t1, 1); \ t2 += t1h; \ t8 += tb; \ tb -= OD_DCT_RSHIFT(t8, 1); \ tn += tk; \ tk -= OD_DCT_RSHIFT(tn, 1); \ t9 += tl; \ tl -= OD_DCT_RSHIFT(t9, 1); \ tm -= ta; \ ta += OD_DCT_RSHIFT(tm, 1); \ tc -= tf; \ tf += OD_DCT_RSHIFT(tc, 1); \ tj += tg; \ tg -= OD_DCT_RSHIFT(tj, 1); \ td -= te; \ te += OD_DCT_RSHIFT(td, 1); \ ti += th; \ th -= OD_DCT_RSHIFT(ti, 1); \ \ t9 = -t9; \ tl = -tl; \ \ /* 805/16384 ~= Tan[Pi/64] ~= 0.04912684976946793 */ \ OD_DCT_OVERFLOW_CHECK(tn, 805, 8192, 316); \ t8 += (tn*805 + 8192) >> 14; \ /* 803/8192 ~= Sin[Pi/32] ~= 0.0980171403295606 */ \ OD_DCT_OVERFLOW_CHECK(t8, 803, 4096, 317); \ tn -= (t8*803 + 4096) >> 13; \ /* 805/16384 ~= Tan[Pi/64] ~= 0.04912684976946793 */ \ OD_DCT_OVERFLOW_CHECK(tn, 805, 8192, 318); \ t8 += (tn*805 + 8192) >> 14; \ /* 11725/32768 ~= Tan[7*Pi/64] ~= 0.3578057213145241 */ \ OD_DCT_OVERFLOW_CHECK(tb, 11725, 16384, 319); \ tk += (tb*11725 + 16384) >> 15; \ /* 5197/8192 ~= Sin[7*Pi/32] ~= 0.6343932841636455 */ \ OD_DCT_OVERFLOW_CHECK(tk, 5197, 4096, 320); \ tb -= (tk*5197 + 4096) >> 13; \ /* 11725/32768 ~= Tan[7*Pi/64] ~= 0.3578057213145241 */ \ OD_DCT_OVERFLOW_CHECK(tb, 11725, 16384, 321); \ tk += (tb*11725 + 16384) >> 15; \ /* 2455/4096 ~= Tan[11*Pi/64] ~= 0.5993769336819237 */ \ OD_DCT_OVERFLOW_CHECK(tl, 2455, 2048, 322); \ ta += (tl*2455 + 2048) >> 12; \ /* 14449/16384 ~= Sin[11*Pi/32] ~= 0.881921264348355 */ \ OD_DCT_OVERFLOW_CHECK(ta, 14449, 8192, 323); \ tl -= (ta*14449 + 8192) >> 14; \ /* 2455/4096 ~= Tan[11*Pi/64] ~= 0.5993769336819237 */ \ OD_DCT_OVERFLOW_CHECK(tl, 2455, 2048, 324); \ ta += (tl*2455 + 2048) >> 12; \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.14833598753834742 */ \ OD_DCT_OVERFLOW_CHECK(tm, 4861, 16384, 325); \ t9 += (tm*4861 + 16384) >> 15; \ /* 1189/4096 ~= Sin[3*Pi/32] ~= 0.29028467725446233 */ \ OD_DCT_OVERFLOW_CHECK(t9, 1189, 2048, 326); \ tm -= (t9*1189 + 2048) >> 12; \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.14833598753834742 */ \ OD_DCT_OVERFLOW_CHECK(tm, 4861, 16384, 327); \ t9 += (tm*4861 + 16384) >> 15; \ /* 805/16384 ~= Tan[Pi/64] ~= 0.04912684976946793 */ \ OD_DCT_OVERFLOW_CHECK(tg, 805, 8192, 328); \ tf += (tg*805 + 8192) >> 14; \ /* 803/8192 ~= Sin[Pi/32] ~= 0.0980171403295606 */ \ OD_DCT_OVERFLOW_CHECK(tf, 803, 4096, 329); \ tg -= (tf*803 + 4096) >> 13; \ /* 805/16384 ~= Tan[Pi/64] ~= 0.04912684976946793 */ \ OD_DCT_OVERFLOW_CHECK(tg, 805, 8192, 330); \ tf += (tg*805 + 8192) >> 14; \ /* 2931/8192 ~= Tan[7*Pi/64] ~= 0.3578057213145241 */ \ OD_DCT_OVERFLOW_CHECK(tj, 2931, 4096, 331); \ tc += (tj*2931 + 4096) >> 13; \ /* 5197/8192 ~= Sin[7*Pi/32] ~= 0.6343932841636455 */ \ OD_DCT_OVERFLOW_CHECK(tc, 5197, 4096, 332); \ tj -= (tc*5197 + 4096) >> 13; \ /* 2931/8192 ~= Tan[7*Pi/64] ~= 0.3578057213145241 */ \ OD_DCT_OVERFLOW_CHECK(tj, 2931, 4096, 333); \ tc += (tj*2931 + 4096) >> 13; \ /* 513/2048 ~= Tan[5*Pi/64] ~= 0.25048696019130545 */ \ OD_DCT_OVERFLOW_CHECK(ti, 513, 1024, 334); \ td += (ti*513 + 1024) >> 11; \ /* 7723/16384 ~= Sin[5*Pi/32] ~= 0.47139673682599764 */ \ OD_DCT_OVERFLOW_CHECK(td, 7723, 8192, 335); \ ti -= (td*7723 + 8192) >> 14; \ /* 513/2048 ~= Tan[5*Pi/64] ~= 0.25048696019130545 */ \ OD_DCT_OVERFLOW_CHECK(ti, 513, 1024, 336); \ td += (ti*513 + 1024) >> 11; \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.14833598753834742 */ \ OD_DCT_OVERFLOW_CHECK(th, 4861, 16384, 337); \ te += (th*4861 + 16384) >> 15; \ /* 1189/4096 ~= Sin[3*Pi/32] ~= 0.29028467725446233 */ \ OD_DCT_OVERFLOW_CHECK(te, 1189, 2048, 338); \ th -= (te*1189 + 2048) >> 12; \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.14833598753834742 */ \ OD_DCT_OVERFLOW_CHECK(th, 4861, 16384, 339); \ te += (th*4861 + 16384) >> 15; \ \ ta = -ta; \ tb = -tb; \ \ tt += t5h; \ t5 -= tt; \ t2 -= tqh; \ tq += t2; \ tp += t1h; \ t1 -= tp; \ t6 -= tuh; \ tu += t6; \ t7 += tvh; \ tv -= t7; \ to += t0h; \ t0 -= to; \ t3 -= t4h; \ t4 += t3; \ ts += trh; \ tr -= ts; \ tf -= OD_DCT_RSHIFT(tn, 1); \ tn += tf; \ tg -= OD_DCT_RSHIFT(t8, 1); \ t8 += tg; \ tk += OD_DCT_RSHIFT(tc, 1); \ tc -= tk; \ tb += OD_DCT_RSHIFT(tj, 1); \ tj -= tb; \ ta += OD_DCT_RSHIFT(ti, 1); \ ti -= ta; \ tl += OD_DCT_RSHIFT(td, 1); \ td -= tl; \ te -= OD_DCT_RSHIFT(tm, 1); \ tm += te; \ th -= OD_DCT_RSHIFT(t9, 1); \ t9 += th; \ ta -= t5; \ t5 += OD_DCT_RSHIFT(ta, 1); \ tq -= tl; \ tl += OD_DCT_RSHIFT(tq, 1); \ t2 -= ti; \ ti += OD_DCT_RSHIFT(t2, 1); \ td -= tt; \ tt += OD_DCT_RSHIFT(td, 1); \ tm += tp; \ tp -= OD_DCT_RSHIFT(tm, 1); \ t6 += t9; \ t9 -= OD_DCT_RSHIFT(t6, 1); \ te -= tu; \ tu += OD_DCT_RSHIFT(te, 1); \ t1 -= th; \ th += OD_DCT_RSHIFT(t1, 1); \ t0 -= tg; \ tg += OD_DCT_RSHIFT(t0, 1); \ tf += tv; \ tv -= OD_DCT_RSHIFT(tf, 1); \ t8 -= t7; \ t7 += OD_DCT_RSHIFT(t8, 1); \ to -= tn; \ tn += OD_DCT_RSHIFT(to, 1); \ t4 -= tk; \ tk += OD_DCT_RSHIFT(t4, 1); \ tb -= tr; \ tr += OD_DCT_RSHIFT(tb, 1); \ t3 -= tj; \ tj += OD_DCT_RSHIFT(t3, 1); \ tc -= ts; \ ts += OD_DCT_RSHIFT(tc, 1); \ \ tr = -tr; \ ts = -ts; \ tt = -tt; \ tu = -tu; \ \ /* 2847/4096 ~= (1/Sqrt[2] - Cos[63*Pi/128]/2)/Sin[63*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t0, 2847, 2048, 340); \ tv += (t0*2847 + 2048) >> 12; \ /* 5791/4096 ~= Sqrt[2]*Sin[63*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tv, 5791, 2048, 341); \ t0 -= (tv*5791 + 2048) >> 12; \ /* 5593/8192 ~= (1/Sqrt[2] - Cos[63*Pi/128])/Sin[63*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t0, 5593, 4096, 342); \ tv += (t0*5593 + 4096) >> 13; \ /* 4099/8192 ~= (1/Sqrt[2] - Cos[31*Pi/128]/2)/Sin[31*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tf, 4099, 4096, 343); \ tg -= (tf*4099 + 4096) >> 13; \ /* 1997/2048 ~= Sqrt[2]*Sin[31*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tg, 1997, 1024, 344); \ tf += (tg*1997 + 1024) >> 11; \ /* -815/32768 ~= (1/Sqrt[2] - Cos[31*Pi/128])/Sin[31*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tf, 815, 16384, 345); \ tg += (tf*815 + 16384) >> 15; \ /* 2527/4096 ~= (1/Sqrt[2] - Cos[17*Pi/128]/2)/Sin[17*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t8, 2527, 2048, 346); \ tn -= (t8*2527 + 2048) >> 12; \ /* 4695/8192 ~= Sqrt[2]*Sin[17*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tn, 4695, 4096, 347); \ t8 += (tn*4695 + 4096) >> 13; \ /* -4187/8192 ~= (1/Sqrt[2] - Cos[17*Pi/128])/Sin[17*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t8, 4187, 4096, 348); \ tn += (t8*4187 + 4096) >> 13; \ /* 5477/8192 ~= (1/Sqrt[2] - Cos[15*Pi/128]/2)/Sin[15*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(to, 5477, 4096, 349); \ t7 += (to*5477 + 4096) >> 13; \ /* 4169/8192 ~= Sqrt[2]*Sin[15*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t7, 4169, 4096, 350); \ to -= (t7*4169 + 4096) >> 13; \ /* -2571/4096 ~= (1/Sqrt[2] - Cos[15*Pi/128])/Sin[15*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(to, 2571, 2048, 351); \ t7 -= (to*2571 + 2048) >> 12; \ /* 5331/8192 ~= (1/Sqrt[2] - Cos[59*Pi/128]/2)/Sin[59*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t2, 5331, 4096, 352); \ tt += (t2*5331 + 4096) >> 13; \ /* 5749/4096 ~= Sqrt[2]*Sin[59*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tt, 5749, 2048, 353); \ t2 -= (tt*5749 + 2048) >> 12; \ /* 2413/4096 ~= (1/Sqrt[2] - Cos[59*Pi/128])/Sin[59*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t2, 2413, 2048, 354); \ tt += (t2*2413 + 2048) >> 12; \ /* 4167/8192 ~= (1/Sqrt[2] - Cos[27*Pi/128]/2)/Sin[27*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(td, 4167, 4096, 355); \ ti -= (td*4167 + 4096) >> 13; \ /* 891/1024 ~= Sqrt[2]*Sin[27*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(ti, 891, 512, 356); \ td += (ti*891 + 512) >> 10; \ /* -4327/32768 ~= (1/Sqrt[2] - Cos[27*Pi/128])/Sin[27*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(td, 4327, 16384, 357); \ ti += (td*4327 + 16384) >> 15; \ /* 2261/4096 ~= (1/Sqrt[2] - Cos[21*Pi/128]/2)/Sin[21*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(ta, 2261, 2048, 358); \ tl -= (ta*2261 + 2048) >> 12; \ /* 2855/4096 ~= Sqrt[2]*Sin[21*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tl, 2855, 2048, 359); \ ta += (tl*2855 + 2048) >> 12; \ /* -5417/16384 ~= (1/Sqrt[2] - Cos[21*Pi/128])/Sin[21*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(ta, 5417, 8192, 360); \ tl += (ta*5417 + 8192) >> 14; \ /* 3459/4096 ~= (1/Sqrt[2] - Cos[11*Pi/128]/2)/Sin[11*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tq, 3459, 2048, 361); \ t5 += (tq*3459 + 2048) >> 12; \ /* 1545/4096 ~= Sqrt[2]*Sin[11*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t5, 1545, 2048, 362); \ tq -= (t5*1545 + 2048) >> 12; \ /* -1971/2048 ~= (1/Sqrt[2] - Cos[11*Pi/128])/Sin[11*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tq, 1971, 1024, 363); \ t5 -= (tq*1971 + 1024) >> 11; \ /* 323/512 ~= (1/Sqrt[2] - Cos[57*Pi/128]/2)/Sin[57*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t3, 323, 256, 364); \ ts += (t3*323 + 256) >> 9; \ /* 5707/4096 ~= Sqrt[2]*Sin[57*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(ts, 5707, 2048, 365); \ t3 -= (ts*5707 + 2048) >> 12; \ /* 2229/4096 ~= (1/Sqrt[2] - Cos[57*Pi/128])/Sin[57*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t3, 2229, 2048, 366); \ ts += (t3*2229 + 2048) >> 12; \ /* 1061/2048 ~= (1/Sqrt[2] - Cos[25*Pi/128]/2)/Sin[25*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tc, 1061, 1024, 367); \ tj -= (tc*1061 + 1024) >> 11; \ /* 6671/8192 ~= Sqrt[2]*Sin[25*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tj, 6671, 4096, 368); \ tc += (tj*6671 + 4096) >> 13; \ /* -6287/32768 ~= (1/Sqrt[2] - Cos[25*Pi/128])/Sin[25*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tc, 6287, 16384, 369); \ tj += (tc*6287 + 16384) >> 15; \ /* 4359/8192 ~= (1/Sqrt[2] - Cos[23*Pi/128]/2)/Sin[23*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tb, 4359, 4096, 370); \ tk -= (tb*4359 + 4096) >> 13; \ /* 3099/4096 ~= Sqrt[2]*Sin[23*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tk, 3099, 2048, 371); \ tb += (tk*3099 + 2048) >> 12; \ /* -2109/8192 ~= (1/Sqrt[2] - Cos[23*Pi/128])/Sin[23*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tb, 2109, 4096, 372); \ tk += (tb*2109 + 4096) >> 13; \ /* 5017/8192 ~= (1/Sqrt[2] - Cos[55*Pi/128]/2)/Sin[55*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t4, 5017, 4096, 373); \ tr += (t4*5017 + 4096) >> 13; \ /* 1413/1024 ~= Sqrt[2]*Sin[55*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tr, 1413, 512, 374); \ t4 -= (tr*1413 + 512) >> 10; \ /* 8195/16384 ~= (1/Sqrt[2] - Cos[55*Pi/128])/Sin[55*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t4, 8195, 8192, 375); \ tr += (t4*8195 + 8192) >> 14; \ /* 2373/4096 ~= (1/Sqrt[2] - Cos[19*Pi/128]/2)/Sin[19*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tm, 2373, 2048, 376); \ t9 += (tm*2373 + 2048) >> 12; \ /* 5209/8192 ~= Sqrt[2]*Sin[19*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t9, 5209, 4096, 377); \ tm -= (t9*5209 + 4096) >> 13; \ /* -3391/8192 ~= (1/Sqrt[2] - Cos[19*Pi/128])/Sin[19*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tm, 3391, 4096, 378); \ t9 -= (tm*3391 + 4096) >> 13; \ /* 1517/2048 ~= (1/Sqrt[2] - Cos[13*Pi/128]/2)/Sin[13*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t6, 1517, 1024, 379); \ tp -= (t6*1517 + 1024) >> 11; \ /* 1817/4096 ~= Sqrt[2]*Sin[13*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tp, 1817, 2048, 380); \ t6 += (tp*1817 + 2048) >> 12; \ /* -6331/8192 ~= (1/Sqrt[2] - Cos[13*Pi/128])/Sin[13*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t6, 6331, 4096, 381); \ tp += (t6*6331 + 4096) >> 13; \ /* 515/1024 ~= (1/Sqrt[2] - Cos[29*Pi/128]/2)/Sin[29*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(te, 515, 512, 382); \ th -= (te*515 + 512) >> 10; \ /* 7567/8192 ~= Sqrt[2]*Sin[29*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(th, 7567, 4096, 383); \ te += (th*7567 + 4096) >> 13; \ /* -2513/32768 ~= (1/Sqrt[2] - Cos[29*Pi/128])/Sin[29*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(te, 2513, 16384, 384); \ th += (te*2513 + 16384) >> 15; \ /* 2753/4096 ~= (1/Sqrt[2] - Cos[61*Pi/128]/2)/Sin[61*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t1, 2753, 2048, 385); \ tu += (t1*2753 + 2048) >> 12; \ /* 5777/4096 ~= Sqrt[2]*Sin[61*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(tu, 5777, 2048, 386); \ t1 -= (tu*5777 + 2048) >> 12; \ /* 1301/2048 ~= (1/Sqrt[2] - Cos[61*Pi/128])/Sin[61*Pi/128] */ \ OD_DCT_OVERFLOW_CHECK(t1, 1301, 1024, 387); \ tu += (t1*1301 + 1024) >> 11; \ } \ while (0) #define OD_IDST_32_ASYM(t0, tg, t8, to, t4, tk, tc, ts, t2, ti, ta, tq, t6, \ tm, te, tu, t1, th, t9, tp, t5, tl, td, tt, t3, tj, tb, tr, t7, tn, tf, tv) \ /* Embedded 32-point asymmetric Type-IV iDST. */ \ do { \ int t0h; \ int t4h; \ int tbh; \ int tfh; \ int tgh; \ int tkh; \ int trh; \ int tvh; \ /* 1301/2048 ~= (1/Sqrt[2] - Cos[61*Pi/128])/Sin[61*Pi/128] */ \ tf -= (tg*1301 + 1024) >> 11; \ /* 5777/4096 ~= Sqrt[2]*Sin[61*Pi/128] */ \ tg += (tf*5777 + 2048) >> 12; \ /* 2753/4096 ~= (1/Sqrt[2] - Cos[61*Pi/128]/2)/Sin[61*Pi/128] */ \ tf -= (tg*2753 + 2048) >> 12; \ /* -2513/32768 ~= (1/Sqrt[2] - Cos[29*Pi/128])/Sin[29*Pi/128] */ \ th -= (te*2513 + 16384) >> 15; \ /* 7567/8192 ~= Sqrt[2]*Sin[29*Pi/128] */ \ te -= (th*7567 + 4096) >> 13; \ /* 515/1024 ~= (1/Sqrt[2] - Cos[29*Pi/128]/2)/Sin[29*Pi/128] */ \ th += (te*515 + 512) >> 10; \ /* -6331/8192 ~= (1/Sqrt[2] - Cos[13*Pi/128])/Sin[13*Pi/128] */ \ tj -= (tc*6331 + 4096) >> 13; \ /* 1817/4096 ~= Sqrt[2]*Sin[13*Pi/128] */ \ tc -= (tj*1817 + 2048) >> 12; \ /* 1517/2048 ~= (1/Sqrt[2] - Cos[13*Pi/128]/2)/Sin[13*Pi/128] */ \ tj += (tc*1517 + 1024) >> 11; \ /* -3391/8192 ~= (1/Sqrt[2] - Cos[19*Pi/128])/Sin[19*Pi/128] */ \ ti += (td*3391 + 4096) >> 13; \ /* 5209/8192 ~= Sqrt[2]*Sin[19*Pi/128] */ \ td += (ti*5209 + 4096) >> 13; \ /* 2373/4096 ~= (1/Sqrt[2] - Cos[19*Pi/128]/2)/Sin[19*Pi/128] */ \ ti -= (td*2373 + 2048) >> 12; \ /* 8195/16384 ~= (1/Sqrt[2] - Cos[55*Pi/128])/Sin[55*Pi/128] */ \ tr -= (t4*8195 + 8192) >> 14; \ /* 1413/1024 ~= Sqrt[2]*Sin[55*Pi/128] */ \ t4 += (tr*1413 + 512) >> 10; \ /* 5017/8192 ~= (1/Sqrt[2] - Cos[55*Pi/128]/2)/Sin[55*Pi/128] */ \ tr -= (t4*5017 + 4096) >> 13; \ /* -2109/8192 ~= (1/Sqrt[2] - Cos[23*Pi/128])/Sin[23*Pi/128] */ \ t5 -= (tq*2109 + 4096) >> 13; \ /* 3099/4096 ~= Sqrt[2]*Sin[23*Pi/128] */ \ tq -= (t5*3099 + 2048) >> 12; \ /* 4359/8192 ~= (1/Sqrt[2] - Cos[23*Pi/128]/2)/Sin[23*Pi/128] */ \ t5 += (tq*4359 + 4096) >> 13; \ /* -6287/32768 ~= (1/Sqrt[2] - Cos[25*Pi/128])/Sin[25*Pi/128] */ \ tp -= (t6*6287 + 16384) >> 15; \ /* 6671/8192 ~= Sqrt[2]*Sin[25*Pi/128] */ \ t6 -= (tp*6671 + 4096) >> 13; \ /* 1061/2048 ~= (1/Sqrt[2] - Cos[25*Pi/128]/2)/Sin[25*Pi/128] */ \ tp += (t6*1061 + 1024) >> 11; \ /* 2229/4096 ~= (1/Sqrt[2] - Cos[57*Pi/128])/Sin[57*Pi/128] */ \ t7 -= (to*2229 + 2048) >> 12; \ /* 5707/4096 ~= Sqrt[2]*Sin[57*Pi/128] */ \ to += (t7*5707 + 2048) >> 12; \ /* 323/512 ~= (1/Sqrt[2] - Cos[57*Pi/128]/2)/Sin[57*Pi/128] */ \ t7 -= (to*323 + 256) >> 9; \ /* -1971/2048 ~= (1/Sqrt[2] - Cos[11*Pi/128])/Sin[11*Pi/128] */ \ tk += (tb*1971 + 1024) >> 11; \ /* 1545/4096 ~= Sqrt[2]*Sin[11*Pi/128] */ \ tb += (tk*1545 + 2048) >> 12; \ /* 3459/4096 ~= (1/Sqrt[2] - Cos[11*Pi/128]/2)/Sin[11*Pi/128] */ \ tk -= (tb*3459 + 2048) >> 12; \ /* -5417/16384 ~= (1/Sqrt[2] - Cos[21*Pi/128])/Sin[21*Pi/128] */ \ tl -= (ta*5417 + 8192) >> 14; \ /* 2855/4096 ~= Sqrt[2]*Sin[21*Pi/128] */ \ ta -= (tl*2855 + 2048) >> 12; \ /* 2261/4096 ~= (1/Sqrt[2] - Cos[21*Pi/128]/2)/Sin[21*Pi/128] */ \ tl += (ta*2261 + 2048) >> 12; \ /* -4327/32768 ~= (1/Sqrt[2] - Cos[27*Pi/128])/Sin[27*Pi/128] */ \ t9 -= (tm*4327 + 16384) >> 15; \ /* 891/1024 ~= Sqrt[2]*Sin[27*Pi/128] */ \ tm -= (t9*891 + 512) >> 10; \ /* 4167/8192 ~= (1/Sqrt[2] - Cos[27*Pi/128]/2)/Sin[27*Pi/128] */ \ t9 += (tm*4167 + 4096) >> 13; \ /* 2413/4096 ~= (1/Sqrt[2] - Cos[59*Pi/128])/Sin[59*Pi/128] */ \ tn -= (t8*2413 + 2048) >> 12; \ /* 5749/4096 ~= Sqrt[2]*Sin[59*Pi/128] */ \ t8 += (tn*5749 + 2048) >> 12; \ /* 5331/8192 ~= (1/Sqrt[2] - Cos[59*Pi/128]/2)/Sin[59*Pi/128] */ \ tn -= (t8*5331 + 4096) >> 13; \ /* -2571/4096 ~= (1/Sqrt[2] - Cos[15*Pi/128])/Sin[15*Pi/128] */ \ ts += (t3*2571 + 2048) >> 12; \ /* 4169/8192 ~= Sqrt[2]*Sin[15*Pi/128] */ \ t3 += (ts*4169 + 4096) >> 13; \ /* 5477/8192 ~= (1/Sqrt[2] - Cos[15*Pi/128]/2)/Sin[15*Pi/128] */ \ ts -= (t3*5477 + 4096) >> 13; \ /* -4187/8192 ~= (1/Sqrt[2] - Cos[17*Pi/128])/Sin[17*Pi/128] */ \ tt -= (t2*4187 + 4096) >> 13; \ /* 4695/8192 ~= Sqrt[2]*Sin[17*Pi/128] */ \ t2 -= (tt*4695 + 4096) >> 13; \ /* 2527/4096 ~= (1/Sqrt[2] - Cos[17*Pi/128]/2)/Sin[17*Pi/128] */ \ tt += (t2*2527 + 2048) >> 12; \ /* -815/32768 ~= (1/Sqrt[2] - Cos[31*Pi/128])/Sin[31*Pi/128] */ \ t1 -= (tu*815 + 16384) >> 15; \ /* 1997/2048 ~= Sqrt[2]*Sin[31*Pi/128] */ \ tu -= (t1*1997 + 1024) >> 11; \ /* 4099/8192 ~= (1/Sqrt[2] - Cos[31*Pi/128]/2)/Sin[31*Pi/128] */ \ t1 += (tu*4099 + 4096) >> 13; \ /* 5593/8192 ~= (1/Sqrt[2] - Cos[63*Pi/128])/Sin[63*Pi/128] */ \ tv -= (t0*5593 + 4096) >> 13; \ /* 5791/4096 ~= Sqrt[2]*Sin[63*Pi/128] */ \ t0 += (tv*5791 + 2048) >> 12; \ /* 2847/4096 ~= (1/Sqrt[2] - Cos[63*Pi/128]/2)/Sin[63*Pi/128] */ \ tv -= (t0*2847 + 2048) >> 12; \ \ t7 = -t7; \ tf = -tf; \ tn = -tn; \ tr = -tr; \ \ t7 -= OD_DCT_RSHIFT(t6, 1); \ t6 += t7; \ tp -= OD_DCT_RSHIFT(to, 1); \ to += tp; \ tr -= OD_DCT_RSHIFT(tq, 1); \ tq += tr; \ t5 -= OD_DCT_RSHIFT(t4, 1); \ t4 += t5; \ tt -= OD_DCT_RSHIFT(t3, 1); \ t3 += tt; \ ts -= OD_DCT_RSHIFT(t2, 1); \ t2 += ts; \ tv += OD_DCT_RSHIFT(tu, 1); \ tu -= tv; \ t1 -= OD_DCT_RSHIFT(t0, 1); \ t0 += t1; \ th -= OD_DCT_RSHIFT(tg, 1); \ tg += th; \ tf -= OD_DCT_RSHIFT(te, 1); \ te += tf; \ ti += OD_DCT_RSHIFT(tc, 1); \ tc -= ti; \ tj += OD_DCT_RSHIFT(td, 1); \ td -= tj; \ tn -= OD_DCT_RSHIFT(tm, 1); \ tm += tn; \ t9 -= OD_DCT_RSHIFT(t8, 1); \ t8 += t9; \ tl -= OD_DCT_RSHIFT(tb, 1); \ tb += tl; \ tk -= OD_DCT_RSHIFT(ta, 1); \ ta += tk; \ \ ti -= th; \ th += OD_DCT_RSHIFT(ti, 1); \ td -= te; \ te += OD_DCT_RSHIFT(td, 1); \ tm += tl; \ tl -= OD_DCT_RSHIFT(tm, 1); \ t9 += ta; \ ta -= OD_DCT_RSHIFT(t9, 1); \ tp += tq; \ tq -= OD_DCT_RSHIFT(tp, 1); \ t6 += t5; \ t5 -= OD_DCT_RSHIFT(t6, 1); \ t2 -= t1; \ t1 += OD_DCT_RSHIFT(t2, 1); \ tt -= tu; \ tu += OD_DCT_RSHIFT(tt, 1); \ tr += t7; \ trh = OD_DCT_RSHIFT(tr, 1); \ t7 -= trh; \ t4 -= to; \ t4h = OD_DCT_RSHIFT(t4, 1); \ to += t4h; \ t0 += t3; \ t0h = OD_DCT_RSHIFT(t0, 1); \ t3 -= t0h; \ tv += ts; \ tvh = OD_DCT_RSHIFT(tv, 1); \ ts -= tvh; \ tf -= tc; \ tfh = OD_DCT_RSHIFT(tf, 1); \ tc += tfh; \ tg += tj; \ tgh = OD_DCT_RSHIFT(tg, 1); \ tj -= tgh; \ tb -= t8; \ tbh = OD_DCT_RSHIFT(tb, 1); \ t8 += tbh; \ tk += tn; \ tkh = OD_DCT_RSHIFT(tk, 1); \ tn -= tkh; \ \ ta = -ta; \ tq = -tq; \ \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.14833598753834742 */ \ te -= (th*4861 + 16384) >> 15; \ /* 1189/4096 ~= Sin[3*Pi/32] ~= 0.29028467725446233 */ \ th += (te*1189 + 2048) >> 12; \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.14833598753834742 */ \ te -= (th*4861 + 16384) >> 15; \ /* 513/2048 ~= Tan[5*Pi/64] ~= 0.25048696019130545 */ \ tm -= (t9*513 + 1024) >> 11; \ /* 7723/16384 ~= Sin[5*Pi/32] ~= 0.47139673682599764 */ \ t9 += (tm*7723 + 8192) >> 14; \ /* 513/2048 ~= Tan[5*Pi/64] ~= 0.25048696019130545 */ \ tm -= (t9*513 + 1024) >> 11; \ /* 2931/8192 ~= Tan[7*Pi/64] ~= 0.3578057213145241 */ \ t6 -= (tp*2931 + 4096) >> 13; \ /* 5197/8192 ~= Sin[7*Pi/32] ~= 0.6343932841636455 */ \ tp += (t6*5197 + 4096) >> 13; \ /* 2931/8192 ~= Tan[7*Pi/64] ~= 0.3578057213145241 */ \ t6 -= (tp*2931 + 4096) >> 13; \ /* 805/16384 ~= Tan[Pi/64] ~= 0.04912684976946793 */ \ tu -= (t1*805 + 8192) >> 14; \ /* 803/8192 ~= Sin[Pi/32] ~= 0.0980171403295606 */ \ t1 += (tu*803 + 4096) >> 13; \ /* 805/16384 ~= Tan[Pi/64] ~= 0.04912684976946793 */ \ tu -= (t1*805 + 8192) >> 14; \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.14833598753834742 */ \ ti -= (td*4861 + 16384) >> 15; \ /* 1189/4096 ~= Sin[3*Pi/32] ~= 0.29028467725446233 */ \ td += (ti*1189 + 2048) >> 12; \ /* 4861/32768 ~= Tan[3*Pi/64] ~= 0.14833598753834742 */ \ ti -= (td*4861 + 16384) >> 15; \ /* 2455/4096 ~= Tan[11*Pi/64] ~= 0.5993769336819237 */ \ ta -= (tl*2455 + 2048) >> 12; \ /* 14449/16384 ~= Sin[11*Pi/32] ~= 0.881921264348355 */ \ tl += (ta*14449 + 8192) >> 14; \ /* 2455/4096 ~= Tan[11*Pi/64] ~= 0.5993769336819237 */ \ ta -= (tl*2455 + 2048) >> 12; \ /* 11725/32768 ~= Tan[7*Pi/64] ~= 0.3578057213145241 */ \ t5 -= (tq*11725 + 16384) >> 15; \ /* 5197/8192 ~= Sin[7*Pi/32] ~= 0.6343932841636455 */ \ tq += (t5*5197 + 4096) >> 13; \ /* 11725/32768 ~= Tan[7*Pi/64] ~= 0.3578057213145241 */ \ t5 -= (tq*11725 + 16384) >> 15; \ /* 805/16384 ~= Tan[Pi/64] ~= 0.04912684976946793 */ \ t2 -= (tt*805 + 8192) >> 14; \ /* 803/8192 ~= Sin[Pi/32] ~= 0.0980171403295606 */ \ tt += (t2*803 + 4096) >> 13; \ /* 805/16384 ~= Tan[Pi/64] ~= 0.04912684976946793 */ \ t2 -= (tt*805 + 8192) >> 14; \ \ tl = -tl; \ ti = -ti; \ \ th += OD_DCT_RSHIFT(t9, 1); \ t9 -= th; \ te -= OD_DCT_RSHIFT(tm, 1); \ tm += te; \ t1 += OD_DCT_RSHIFT(tp, 1); \ tp -= t1; \ tu -= OD_DCT_RSHIFT(t6, 1); \ t6 += tu; \ ta -= OD_DCT_RSHIFT(td, 1); \ td += ta; \ tl += OD_DCT_RSHIFT(ti, 1); \ ti -= tl; \ t5 += OD_DCT_RSHIFT(tt, 1); \ tt -= t5; \ tq += OD_DCT_RSHIFT(t2, 1); \ t2 -= tq; \ \ t8 -= tgh; \ tg += t8; \ tn += tfh; \ tf -= tn; \ t7 -= tvh; \ tv += t7; \ to -= t0h; \ t0 += to; \ tc += tbh; \ tb -= tc; \ tj += tkh; \ tk -= tj; \ ts += t4h; \ t4 -= ts; \ t3 += trh; \ tr -= t3; \ \ tk = -tk; \ \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ tc -= (tj*2485 + 4096) >> 13; \ /* 18205/32768 ~= Sin[3*Pi/16] ~= 0.555570233019602 */ \ tj += (tc*18205 + 16384) >> 15; \ /* 2485/8192 ~= Tan[3*Pi/32] ~= 0.303346683607342 */ \ tc -= (tj*2485 + 4096) >> 13; \ /* 3227/32768 ~= Tan[Pi/32] ~= 0.09849140335716425 */ \ ts -= (t3*3227 + 16384) >> 15; \ /* 6393/32768 ~= Sin[Pi/16] ~= 0.19509032201612825 */ \ t3 += (ts*6393 + 16384) >> 15; \ /* 3227/32768 ~= Tan[Pi/32] ~= 0.09849140335716425 */ \ ts -= (t3*3227 + 16384) >> 15; \ /* 17515/32768 ~= Tan[5*Pi/32] ~= 0.5345111359507916 */ \ tk -= (tb*17515 + 16384) >> 15; \ /* 13623/16384 ~= Sin[5*Pi/16] ~= 0.8314696123025452 */ \ tb += (tk*13623 + 8192) >> 14; \ /* 17515/32768 ~= Tan[5*Pi/32] ~= 0.5345111359507916 */ \ tk -= (tb*17515 + 16384) >> 15; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.8206787908286602 */ \ t4 -= (tr*6723 + 4096) >> 13; \ /* 16069/16384 ~= Sin[7*Pi/16] ~= 0.9807852804032304 */ \ tr += (t4*16069 + 8192) >> 14; \ /* 6723/8192 ~= Tan[7*Pi/32] ~= 0.8206787908286602 */ \ t4 -= (tr*6723 + 4096) >> 13; \ \ t4 = -t4; \ \ tp += tm; \ tm -= OD_DCT_RSHIFT(tp, 1); \ t9 -= t6; \ t6 += OD_DCT_RSHIFT(t9, 1); \ th -= t1; \ t1 += OD_DCT_RSHIFT(th, 1); \ tu -= te; \ te += OD_DCT_RSHIFT(tu, 1); /* pass */ \ t5 -= tl; \ tl += OD_DCT_RSHIFT(t5, 1); \ ta += tq; \ tq -= OD_DCT_RSHIFT(ta, 1); \ td += tt; \ tt -= OD_DCT_RSHIFT(td, 1); \ t2 -= ti; \ ti += OD_DCT_RSHIFT(t2, 1); /* pass */ \ t7 += t8; \ t8 -= OD_DCT_RSHIFT(t7, 1); \ tn -= to; \ to += OD_DCT_RSHIFT(tn, 1); \ tf -= tv; \ tv += OD_DCT_RSHIFT(tf, 1); \ t0 += tg; \ tg -= OD_DCT_RSHIFT(t0, 1); /* pass */ \ tj -= t3; \ t3 += OD_DCT_RSHIFT(tj, 1); /* pass */ \ ts -= tc; \ tc += OD_DCT_RSHIFT(ts, 1); \ t4 -= tb; \ tb += OD_DCT_RSHIFT(t4, 1); /* pass */ \ tk -= tr; \ tr += OD_DCT_RSHIFT(tk, 1); \ \ t1 = -t1; \ t3 = -t3; \ t7 = -t7; \ t8 = -t8; \ tg = -tg; \ tm = -tm; \ to = -to; \ \ /* 14341/16384 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ tm -= (t9*14341 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ t9 += (tm*15137 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ tm -= (t9*4161 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ tp -= (t6*4161 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ t6 += (tp*15137 + 8192) >> 14; \ /* 28681/32768 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ tp -= (t6*28681 + 16384) >> 15; \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ th += (te*19195 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ te += (th*11585 + 8192) >> 14; \ /* 29957/32768 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ th -= (te*29957 + 16384) >> 15; \ /* 14341/16384 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ tq -= (t5*14341 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ t5 += (tq*15137 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ tq -= (t5*4161 + 8192) >> 14; \ /* 3259/8192 ~= 2*Tan[Pi/16] ~= 0.397824734759316 */ \ ta -= (tl*3259 + 4096) >> 13; \ /* 3135/16384 ~= Sin[Pi/8]/2 ~= 0.1913417161825449 */ \ tl += (ta*3135 + 8192) >> 14; \ /* 3259/8192 ~= 2*Tan[Pi/16] ~= 0.397824734759316 */ \ ta -= (tl*3259 + 4096) >> 13; \ /* 7489/8192 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ ti -= (td*7489 + 4096) >> 13; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ td += (ti*11585 + 8192) >> 14; \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ ti += (td*19195 + 16384) >> 15; \ /* 14341/16384 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ to -= (t7*14341 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ t7 += (to*15137 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ to -= (t7*4161 + 8192) >> 14; \ /* 4161/16384 ~= Tan[3*Pi/16] - Tan[Pi/8] ~= 0.253965075546204 */ \ tn -= (t8*4161 + 8192) >> 14; \ /* 15137/16384 ~= Sin[3*Pi/8] ~= 0.923879532511287 */ \ t8 += (tn*15137 + 8192) >> 14; \ /* 28681/32768 ~= Tan[3*Pi/16] + Tan[Pi/8]/2 ~= 0.875285419105846 */ \ tn -= (t8*28681 + 16384) >> 15; \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ tf += (tg*19195 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ tg += (tf*11585 + 8192) >> 14; \ /* 29957/32768 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ tf -= (tg*29957 + 16384) >> 15; \ /* -19195/32768 ~= Tan[Pi/8] - Tan[Pi/4] ~= -0.585786437626905 */ \ tj += (tc*19195 + 16384) >> 15; \ /* 11585/16384 ~= Sin[Pi/4] ~= 0.707106781186548 */ \ tc += (tj*11585 + 8192) >> 14; \ /* 29957/32768 ~= Tan[Pi/8] + Tan[Pi/4]/2 ~= 0.914213562373095 */ \ tj -= (tc*29957 + 16384) >> 15; \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ tk += (tb*13573 + 8192) >> 14; \ /* 11585/32768 ~= Sin[Pi/4]/2 ~= 0.353553390593274 */ \ tb -= (tk*11585 + 16384) >> 15; \ /* 13573/16384 ~= 2*Tan[Pi/8] ~= 0.828427124746190 */ \ tk += (tb*13573 + 8192) >> 14; \ \ tf = -tf; \ \ } \ while (0) #define OD_FDCT_64(u0, uw, ug, uM, u8, uE, uo, uU, u4, uA, uk, uQ, uc, uI, \ us, uY, u2, uy, ui, uO, ua, uG, uq, uW, u6, uC, um, uS, ue, uK, uu, u_, u1, \ ux, uh, uN, u9, uF, up, uV, u5, uB, ul, uR, ud, uJ, ut, uZ, u3, uz, uj, uP, \ ub, uH, ur, uX, u7, uD, un, uT, uf, uL, uv, u) \ /* Embedded 64-point orthonormal Type-II fDCT. */ \ do { \ int uwh; \ int uxh; \ int uyh; \ int uzh; \ int uAh; \ int uBh; \ int uCh; \ int uDh; \ int uEh; \ int uFh; \ int uGh; \ int uHh; \ int uIh; \ int uJh; \ int uKh; \ int uLh; \ int uMh; \ int uNh; \ int uOh; \ int uPh; \ int uQh; \ int uRh; \ int uSh; \ int uTh; \ int uUh; \ int uVh; \ int uWh; \ int uXh; \ int uYh; \ int uZh; \ int u_h; \ int uh_; \ u = u0 - u; \ uh_ = OD_DCT_RSHIFT(u, 1); \ u0 -= uh_; \ u_ += u1; \ u_h = OD_DCT_RSHIFT(u_, 1); \ u1 = u_h - u1; \ uZ = u2 - uZ; \ uZh = OD_DCT_RSHIFT(uZ, 1); \ u2 -= uZh; \ uY += u3; \ uYh = OD_DCT_RSHIFT(uY, 1); \ u3 = uYh - u3; \ uX = u4 - uX; \ uXh = OD_DCT_RSHIFT(uX, 1); \ u4 -= uXh; \ uW += u5; \ uWh = OD_DCT_RSHIFT(uW, 1); \ u5 = uWh - u5; \ uV = u6 - uV; \ uVh = OD_DCT_RSHIFT(uV, 1); \ u6 -= uVh; \ uU += u7; \ uUh = OD_DCT_RSHIFT(uU, 1); \ u7 = uUh - u7; \ uT = u8 - uT; \ uTh = OD_DCT_RSHIFT(uT, 1); \ u8 -= uTh; \ uS += u9; \ uSh = OD_DCT_RSHIFT(uS, 1); \ u9 = uSh - u9; \ uR = ua - uR; \ uRh = OD_DCT_RSHIFT(uR, 1); \ ua -= uRh; \ uQ += ub; \ uQh = OD_DCT_RSHIFT(uQ, 1); \ ub = uQh - ub; \ uP = uc - uP; \ uPh = OD_DCT_RSHIFT(uP, 1); \ uc -= uPh; \ uO += ud; \ uOh = OD_DCT_RSHIFT(uO, 1); \ ud = uOh - ud; \ uN = ue - uN; \ uNh = OD_DCT_RSHIFT(uN, 1); \ ue -= uNh; \ uM += uf; \ uMh = OD_DCT_RSHIFT(uM, 1); \ uf = uMh - uf; \ uL = ug - uL; \ uLh = OD_DCT_RSHIFT(uL, 1); \ ug -= uLh; \ uK += uh; \ uKh = OD_DCT_RSHIFT(uK, 1); \ uh = uKh - uh; \ uJ = ui - uJ; \ uJh = OD_DCT_RSHIFT(uJ, 1); \ ui -= uJh; \ uI += uj; \ uIh = OD_DCT_RSHIFT(uI, 1); \ uj = uIh - uj; \ uH = uk - uH; \ uHh = OD_DCT_RSHIFT(uH, 1); \ uk -= uHh; \ uG += ul; \ uGh = OD_DCT_RSHIFT(uG, 1); \ ul = uGh - ul; \ uF = um - uF; \ uFh = OD_DCT_RSHIFT(uF, 1); \ um -= uFh; \ uE += un; \ uEh = OD_DCT_RSHIFT(uE, 1); \ un = uEh - un; \ uD = uo - uD; \ uDh = OD_DCT_RSHIFT(uD, 1); \ uo -= uDh; \ uC += up; \ uCh = OD_DCT_RSHIFT(uC, 1); \ up = uCh - up; \ uB = uq - uB; \ uBh = OD_DCT_RSHIFT(uB, 1); \ uq -= uBh; \ uA += ur; \ uAh = OD_DCT_RSHIFT(uA, 1); \ ur = uAh - ur; \ uz = us - uz; \ uzh = OD_DCT_RSHIFT(uz, 1); \ us -= uzh; \ uy += ut; \ uyh = OD_DCT_RSHIFT(uy, 1); \ ut = uyh - ut; \ ux = uu - ux; \ uxh = OD_DCT_RSHIFT(ux, 1); \ uu -= uxh; \ uw += uv; \ uwh = OD_DCT_RSHIFT(uw, 1); \ uv = uwh - uv; \ OD_FDCT_32_ASYM(u0, uw, uwh, ug, uM, uMh, u8, uE, uEh, uo, uU, uUh, \ u4, uA, uAh, uk, uQ, uQh, uc, uI, uIh, us, uY, uYh, u2, uy, uyh, \ ui, uO, uOh, ua, uG, uGh, uq, uW, uWh, u6, uC, uCh, um, uS, uSh, \ ue, uK, uKh, uu, u_, u_h); \ OD_FDST_32_ASYM(u, uv, uL, uf, uT, un, uD, u7, uX, ur, uH, ub, uP, uj, \ uz, u3, uZ, ut, uJ, ud, uR, ul, uB, u5, uV, up, uF, u9, uN, uh, ux, u1); \ } \ while (0) #define OD_IDCT_64(u0, uw, ug, uM, u8, uE, uo, uU, u4, uA, uk, uQ, uc, uI, \ us, uY, u2, uy, ui, uO, ua, uG, uq, uW, u6, uC, um, uS, ue, uK, uu, u_, u1, \ ux, uh, uN, u9, uF, up, uV, u5, uB, ul, uR, ud, uJ, ut, uZ, u3, uz, uj, uP, \ ub, uH, ur, uX, u7, uD, un, uT, uf, uL, uv, u) \ /* Embedded 64-point orthonormal Type-II fDCT. */ \ do { \ int u1h; \ int u3h; \ int u5h; \ int u7h; \ int u9h; \ int ubh; \ int udh; \ int ufh; \ int uhh; \ int ujh; \ int ulh; \ int unh; \ int uph; \ int urh; \ int uth; \ int uvh; \ int uxh; \ int uzh; \ int uBh; \ int uDh; \ int uFh; \ int uHh; \ int uJh; \ int uLh; \ int uNh; \ int uPh; \ int uRh; \ int uTh; \ int uVh; \ int uXh; \ int uZh; \ int uh_; \ OD_IDST_32_ASYM(u, uL, uT, uD, uX, uH, uP, uz, uZ, uJ, uR, uB, uV, uF, \ uN, ux, u_, uK, uS, uC, uW, uG, uO, uy, uY, uI, uQ, uA, uU, uE, uM, uw); \ OD_IDCT_32_ASYM(u0, ug, u8, uo, u4, uk, uc, us, u2, ui, ua, uq, u6, um, \ ue, uu, u1, u1h, uh, uhh, u9, u9h, up, uph, u5, u5h, ul, ulh, ud, udh, \ ut, uth, u3, u3h, uj, ujh, ub, ubh, ur, urh, u7, u7h, un, unh, uf, ufh, \ uv, uvh); \ uh_ = OD_DCT_RSHIFT(u, 1); \ u0 += uh_; \ u = u0 - u; \ u_ = u1h - u_; \ u1 -= u_; \ uZh = OD_DCT_RSHIFT(uZ, 1); \ u2 += uZh; \ uZ = u2 - uZ; \ uY = u3h - uY; \ u3 -= uY; \ uXh = OD_DCT_RSHIFT(uX, 1); \ u4 += uXh; \ uX = u4 - uX; \ uW = u5h - uW; \ u5 -= uW; \ uVh = OD_DCT_RSHIFT(uV, 1); \ u6 += uVh; \ uV = u6 - uV; \ uU = u7h - uU; \ u7 -= uU; \ uTh = OD_DCT_RSHIFT(uT, 1); \ u8 += uTh; \ uT = u8 - uT; \ uS = u9h - uS; \ u9 -= uS; \ uRh = OD_DCT_RSHIFT(uR, 1); \ ua += uRh; \ uR = ua - uR; \ uQ = ubh - uQ; \ ub -= uQ; \ uPh = OD_DCT_RSHIFT(uP, 1); \ uc += uPh; \ uP = uc - uP; \ uO = udh - uO; \ ud -= uO; \ uNh = OD_DCT_RSHIFT(uN, 1); \ ue += uNh; \ uN = ue - uN; \ uM = ufh - uM; \ uf -= uM; \ uLh = OD_DCT_RSHIFT(uL, 1); \ ug += uLh; \ uL = ug - uL; \ uK = uhh - uK; \ uh -= uK; \ uJh = OD_DCT_RSHIFT(uJ, 1); \ ui += uJh; \ uJ = ui - uJ; \ uI = ujh - uI; \ uj -= uI; \ uHh = OD_DCT_RSHIFT(uH, 1); \ uk += uHh; \ uH = uk - uH; \ uG = ulh - uG; \ ul -= uG; \ uFh = OD_DCT_RSHIFT(uF, 1); \ um += uFh; \ uF = um - uF; \ uE = unh - uE; \ un -= uE; \ uDh = OD_DCT_RSHIFT(uD, 1); \ uo += uDh; \ uD = uo - uD; \ uC = uph - uC; \ up -= uC; \ uBh = OD_DCT_RSHIFT(uB, 1); \ uq += uBh; \ uB = uq - uB; \ uA = urh - uA; \ ur -= uA; \ uzh = OD_DCT_RSHIFT(uz, 1); \ us += uzh; \ uz = us - uz; \ uy = uth - uy; \ ut -= uy; \ uxh = OD_DCT_RSHIFT(ux, 1); \ uu += uxh; \ ux = uu - ux; \ uw = uvh - uw; \ uv -= uw; \ } while (0) #endif void od_bin_fdct4(od_coeff y[4], const od_coeff *x, int xstride) { int q0; int q1; int q2; int q3; q0 = x[0*xstride]; q2 = x[1*xstride]; q1 = x[2*xstride]; q3 = x[3*xstride]; OD_FDCT_4(q0, q2, q1, q3); y[0] = (od_coeff)q0; y[1] = (od_coeff)q1; y[2] = (od_coeff)q2; y[3] = (od_coeff)q3; } void od_bin_idct4(od_coeff *x, int xstride, const od_coeff y[4]) { int q0; int q1; int q2; int q3; q0 = y[0]; q2 = y[1]; q1 = y[2]; q3 = y[3]; OD_IDCT_4(q0, q2, q1, q3); x[0*xstride] = q0; x[1*xstride] = q1; x[2*xstride] = q2; x[3*xstride] = q3; } void od_bin_fdst4(od_coeff y[4], const od_coeff *x, int xstride) { int q0; int q1; int q2; int q3; q0 = x[3*xstride]; q2 = x[2*xstride]; q1 = x[1*xstride]; q3 = x[0*xstride]; OD_FDST_4(q0, q2, q1, q3); y[0] = (od_coeff)q3; y[1] = (od_coeff)q2; y[2] = (od_coeff)q1; y[3] = (od_coeff)q0; } void od_bin_idst4(od_coeff *x, int xstride, const od_coeff y[4]) { int q0; int q1; int q2; int q3; q0 = y[3]; q2 = y[2]; q1 = y[1]; q3 = y[0]; OD_IDST_4(q0, q2, q1, q3); x[0*xstride] = q3; x[1*xstride] = q2; x[2*xstride] = q1; x[3*xstride] = q0; } void od_bin_fdct8(od_coeff y[8], const od_coeff *x, int xstride) { int r0; int r1; int r2; int r3; int r4; int r5; int r6; int r7; r0 = x[0*xstride]; r4 = x[1*xstride]; r2 = x[2*xstride]; r6 = x[3*xstride]; r1 = x[4*xstride]; r5 = x[5*xstride]; r3 = x[6*xstride]; r7 = x[7*xstride]; OD_FDCT_8(r0, r4, r2, r6, r1, r5, r3, r7); y[0] = (od_coeff)r0; y[1] = (od_coeff)r1; y[2] = (od_coeff)r2; y[3] = (od_coeff)r3; y[4] = (od_coeff)r4; y[5] = (od_coeff)r5; y[6] = (od_coeff)r6; y[7] = (od_coeff)r7; } void od_bin_idct8(od_coeff *x, int xstride, const od_coeff y[8]) { int r0; int r1; int r2; int r3; int r4; int r5; int r6; int r7; r0 = y[0]; r4 = y[1]; r2 = y[2]; r6 = y[3]; r1 = y[4]; r5 = y[5]; r3 = y[6]; r7 = y[7]; OD_IDCT_8(r0, r4, r2, r6, r1, r5, r3, r7); x[0*xstride] = (od_coeff)r0; x[1*xstride] = (od_coeff)r1; x[2*xstride] = (od_coeff)r2; x[3*xstride] = (od_coeff)r3; x[4*xstride] = (od_coeff)r4; x[5*xstride] = (od_coeff)r5; x[6*xstride] = (od_coeff)r6; x[7*xstride] = (od_coeff)r7; } void od_bin_fdst8(od_coeff y[8], const od_coeff *x, int xstride) { int r0; int r1; int r2; int r3; int r4; int r5; int r6; int r7; r0 = x[0*xstride]; r4 = x[1*xstride]; r2 = x[2*xstride]; r6 = x[3*xstride]; r1 = x[4*xstride]; r5 = x[5*xstride]; r3 = x[6*xstride]; r7 = x[7*xstride]; OD_FDST_8(r0, r4, r2, r6, r1, r5, r3, r7); y[0] = (od_coeff)r0; y[1] = (od_coeff)r1; y[2] = (od_coeff)r2; y[3] = (od_coeff)r3; y[4] = (od_coeff)r4; y[5] = (od_coeff)r5; y[6] = (od_coeff)r6; y[7] = (od_coeff)r7; } void od_bin_idst8(od_coeff *x, int xstride, const od_coeff y[8]) { int r0; int r1; int r2; int r3; int r4; int r5; int r6; int r7; r0 = y[0]; r4 = y[1]; r2 = y[2]; r6 = y[3]; r1 = y[4]; r5 = y[5]; r3 = y[6]; r7 = y[7]; OD_IDST_8(r0, r4, r2, r6, r1, r5, r3, r7); x[0*xstride] = (od_coeff)r0; x[1*xstride] = (od_coeff)r1; x[2*xstride] = (od_coeff)r2; x[3*xstride] = (od_coeff)r3; x[4*xstride] = (od_coeff)r4; x[5*xstride] = (od_coeff)r5; x[6*xstride] = (od_coeff)r6; x[7*xstride] = (od_coeff)r7; } void od_bin_fdct16(od_coeff y[16], const od_coeff *x, int xstride) { int s0; int s1; int s2; int s3; int s4; int s5; int s6; int s7; int s8; int s9; int sa; int sb; int sc; int sd; int se; int sf; s0 = x[0*xstride]; s8 = x[1*xstride]; s4 = x[2*xstride]; sc = x[3*xstride]; s2 = x[4*xstride]; sa = x[5*xstride]; s6 = x[6*xstride]; se = x[7*xstride]; s1 = x[8*xstride]; s9 = x[9*xstride]; s5 = x[10*xstride]; sd = x[11*xstride]; s3 = x[12*xstride]; sb = x[13*xstride]; s7 = x[14*xstride]; sf = x[15*xstride]; OD_FDCT_16(s0, s8, s4, sc, s2, sa, s6, se, s1, s9, s5, sd, s3, sb, s7, sf); y[0] = (od_coeff)s0; y[1] = (od_coeff)s1; y[2] = (od_coeff)s2; y[3] = (od_coeff)s3; y[4] = (od_coeff)s4; y[5] = (od_coeff)s5; y[6] = (od_coeff)s6; y[7] = (od_coeff)s7; y[8] = (od_coeff)s8; y[9] = (od_coeff)s9; y[10] = (od_coeff)sa; y[11] = (od_coeff)sb; y[12] = (od_coeff)sc; y[13] = (od_coeff)sd; y[14] = (od_coeff)se; y[15] = (od_coeff)sf; } void od_bin_idct16(od_coeff *x, int xstride, const od_coeff y[16]) { int s0; int s1; int s2; int s3; int s4; int s5; int s6; int s7; int s8; int s9; int sa; int sb; int sc; int sd; int se; int sf; s0 = y[0]; s8 = y[1]; s4 = y[2]; sc = y[3]; s2 = y[4]; sa = y[5]; s6 = y[6]; se = y[7]; s1 = y[8]; s9 = y[9]; s5 = y[10]; sd = y[11]; s3 = y[12]; sb = y[13]; s7 = y[14]; sf = y[15]; OD_IDCT_16(s0, s8, s4, sc, s2, sa, s6, se, s1, s9, s5, sd, s3, sb, s7, sf); x[0*xstride] = (od_coeff)s0; x[1*xstride] = (od_coeff)s1; x[2*xstride] = (od_coeff)s2; x[3*xstride] = (od_coeff)s3; x[4*xstride] = (od_coeff)s4; x[5*xstride] = (od_coeff)s5; x[6*xstride] = (od_coeff)s6; x[7*xstride] = (od_coeff)s7; x[8*xstride] = (od_coeff)s8; x[9*xstride] = (od_coeff)s9; x[10*xstride] = (od_coeff)sa; x[11*xstride] = (od_coeff)sb; x[12*xstride] = (od_coeff)sc; x[13*xstride] = (od_coeff)sd; x[14*xstride] = (od_coeff)se; x[15*xstride] = (od_coeff)sf; } void od_bin_fdst16(od_coeff y[16], const od_coeff *x, int xstride) { int s0; int s1; int s2; int s3; int s4; int s5; int s6; int s7; int s8; int s9; int sa; int sb; int sc; int sd; int se; int sf; s0 = x[15*xstride]; s8 = x[14*xstride]; s4 = x[13*xstride]; sc = x[12*xstride]; s2 = x[11*xstride]; sa = x[10*xstride]; s6 = x[9*xstride]; se = x[8*xstride]; s1 = x[7*xstride]; s9 = x[6*xstride]; s5 = x[5*xstride]; sd = x[4*xstride]; s3 = x[3*xstride]; sb = x[2*xstride]; s7 = x[1*xstride]; sf = x[0*xstride]; OD_FDST_16(s0, s8, s4, sc, s2, sa, s6, se, s1, s9, s5, sd, s3, sb, s7, sf); y[0] = (od_coeff)sf; y[1] = (od_coeff)se; y[2] = (od_coeff)sd; y[3] = (od_coeff)sc; y[4] = (od_coeff)sb; y[5] = (od_coeff)sa; y[6] = (od_coeff)s9; y[7] = (od_coeff)s8; y[8] = (od_coeff)s7; y[9] = (od_coeff)s6; y[10] = (od_coeff)s5; y[11] = (od_coeff)s4; y[12] = (od_coeff)s3; y[13] = (od_coeff)s2; y[14] = (od_coeff)s1; y[15] = (od_coeff)s0; } void od_bin_idst16(od_coeff *x, int xstride, const od_coeff y[16]) { int s0; int s1; int s2; int s3; int s4; int s5; int s6; int s7; int s8; int s9; int sa; int sb; int sc; int sd; int se; int sf; s0 = y[15]; s8 = y[14]; s4 = y[13]; sc = y[12]; s2 = y[11]; sa = y[10]; s6 = y[9]; se = y[8]; s1 = y[7]; s9 = y[6]; s5 = y[5]; sd = y[4]; s3 = y[3]; sb = y[2]; s7 = y[1]; sf = y[0]; OD_IDST_16(s0, s8, s4, sc, s2, sa, s6, se, s1, s9, s5, sd, s3, sb, s7, sf); x[0*xstride] = (od_coeff)sf; x[1*xstride] = (od_coeff)se; x[2*xstride] = (od_coeff)sd; x[3*xstride] = (od_coeff)sc; x[4*xstride] = (od_coeff)sb; x[5*xstride] = (od_coeff)sa; x[6*xstride] = (od_coeff)s9; x[7*xstride] = (od_coeff)s8; x[8*xstride] = (od_coeff)s7; x[9*xstride] = (od_coeff)s6; x[10*xstride] = (od_coeff)s5; x[11*xstride] = (od_coeff)s4; x[12*xstride] = (od_coeff)s3; x[13*xstride] = (od_coeff)s2; x[14*xstride] = (od_coeff)s1; x[15*xstride] = (od_coeff)s0; } void od_bin_fdct32(od_coeff y[32], const od_coeff *x, int xstride) { /*215 adds, 38 shifts, 87 "muls".*/ int t0; int t1; int t2; int t3; int t4; int t5; int t6; int t7; int t8; int t9; int ta; int tb; int tc; int td; int te; int tf; int tg; int th; int ti; int tj; int tk; int tl; int tm; int tn; int to; int tp; int tq; int tr; int ts; int tt; int tu; int tv; t0 = x[0*xstride]; tg = x[1*xstride]; t8 = x[2*xstride]; to = x[3*xstride]; t4 = x[4*xstride]; tk = x[5*xstride]; tc = x[6*xstride]; ts = x[7*xstride]; t2 = x[8*xstride]; ti = x[9*xstride]; ta = x[10*xstride]; tq = x[11*xstride]; t6 = x[12*xstride]; tm = x[13*xstride]; te = x[14*xstride]; tu = x[15*xstride]; t1 = x[16*xstride]; th = x[17*xstride]; t9 = x[18*xstride]; tp = x[19*xstride]; t5 = x[20*xstride]; tl = x[21*xstride]; td = x[22*xstride]; tt = x[23*xstride]; t3 = x[24*xstride]; tj = x[25*xstride]; tb = x[26*xstride]; tr = x[27*xstride]; t7 = x[28*xstride]; tn = x[29*xstride]; tf = x[30*xstride]; tv = x[31*xstride]; OD_FDCT_32(t0, tg, t8, to, t4, tk, tc, ts, t2, ti, ta, tq, t6, tm, te, tu, t1, th, t9, tp, t5, tl, td, tt, t3, tj, tb, tr, t7, tn, tf, tv); y[0] = (od_coeff)t0; y[1] = (od_coeff)t1; y[2] = (od_coeff)t2; y[3] = (od_coeff)t3; y[4] = (od_coeff)t4; y[5] = (od_coeff)t5; y[6] = (od_coeff)t6; y[7] = (od_coeff)t7; y[8] = (od_coeff)t8; y[9] = (od_coeff)t9; y[10] = (od_coeff)ta; y[11] = (od_coeff)tb; y[12] = (od_coeff)tc; y[13] = (od_coeff)td; y[14] = (od_coeff)te; y[15] = (od_coeff)tf; y[16] = (od_coeff)tg; y[17] = (od_coeff)th; y[18] = (od_coeff)ti; y[19] = (od_coeff)tj; y[20] = (od_coeff)tk; y[21] = (od_coeff)tl; y[22] = (od_coeff)tm; y[23] = (od_coeff)tn; y[24] = (od_coeff)to; y[25] = (od_coeff)tp; y[26] = (od_coeff)tq; y[27] = (od_coeff)tr; y[28] = (od_coeff)ts; y[29] = (od_coeff)tt; y[30] = (od_coeff)tu; y[31] = (od_coeff)tv; } void od_bin_idct32(od_coeff *x, int xstride, const od_coeff y[32]) { int t0; int t1; int t2; int t3; int t4; int t5; int t6; int t7; int t8; int t9; int ta; int tb; int tc; int td; int te; int tf; int tg; int th; int ti; int tj; int tk; int tl; int tm; int tn; int to; int tp; int tq; int tr; int ts; int tt; int tu; int tv; t0 = y[0]; tg = y[1]; t8 = y[2]; to = y[3]; t4 = y[4]; tk = y[5]; tc = y[6]; ts = y[7]; t2 = y[8]; ti = y[9]; ta = y[10]; tq = y[11]; t6 = y[12]; tm = y[13]; te = y[14]; tu = y[15]; t1 = y[16]; th = y[17]; t9 = y[18]; tp = y[19]; t5 = y[20]; tl = y[21]; td = y[22]; tt = y[23]; t3 = y[24]; tj = y[25]; tb = y[26]; tr = y[27]; t7 = y[28]; tn = y[29]; tf = y[30]; tv = y[31]; OD_IDCT_32(t0, tg, t8, to, t4, tk, tc, ts, t2, ti, ta, tq, t6, tm, te, tu, t1, th, t9, tp, t5, tl, td, tt, t3, tj, tb, tr, t7, tn, tf, tv); x[0*xstride] = (od_coeff)t0; x[1*xstride] = (od_coeff)t1; x[2*xstride] = (od_coeff)t2; x[3*xstride] = (od_coeff)t3; x[4*xstride] = (od_coeff)t4; x[5*xstride] = (od_coeff)t5; x[6*xstride] = (od_coeff)t6; x[7*xstride] = (od_coeff)t7; x[8*xstride] = (od_coeff)t8; x[9*xstride] = (od_coeff)t9; x[10*xstride] = (od_coeff)ta; x[11*xstride] = (od_coeff)tb; x[12*xstride] = (od_coeff)tc; x[13*xstride] = (od_coeff)td; x[14*xstride] = (od_coeff)te; x[15*xstride] = (od_coeff)tf; x[16*xstride] = (od_coeff)tg; x[17*xstride] = (od_coeff)th; x[18*xstride] = (od_coeff)ti; x[19*xstride] = (od_coeff)tj; x[20*xstride] = (od_coeff)tk; x[21*xstride] = (od_coeff)tl; x[22*xstride] = (od_coeff)tm; x[23*xstride] = (od_coeff)tn; x[24*xstride] = (od_coeff)to; x[25*xstride] = (od_coeff)tp; x[26*xstride] = (od_coeff)tq; x[27*xstride] = (od_coeff)tr; x[28*xstride] = (od_coeff)ts; x[29*xstride] = (od_coeff)tt; x[30*xstride] = (od_coeff)tu; x[31*xstride] = (od_coeff)tv; } #if CONFIG_TX64X64 void od_bin_fdct64(od_coeff y[64], const od_coeff *x, int xstride) { int t0; int t1; int t2; int t3; int t4; int t5; int t6; int t7; int t8; int t9; int ta; int tb; int tc; int td; int te; int tf; int tg; int th; int ti; int tj; int tk; int tl; int tm; int tn; int to; int tp; int tq; int tr; int ts; int tt; int tu; int tv; int tw; int tx; int ty; int tz; int tA; int tB; int tC; int tD; int tE; int tF; int tG; int tH; int tI; int tJ; int tK; int tL; int tM; int tN; int tO; int tP; int tQ; int tR; int tS; int tT; int tU; int tV; int tW; int tX; int tY; int tZ; int t_; int t; t0 = x[0*xstride]; tw = x[1*xstride]; tg = x[2*xstride]; tM = x[3*xstride]; t8 = x[4*xstride]; tE = x[5*xstride]; to = x[6*xstride]; tU = x[7*xstride]; t4 = x[8*xstride]; tA = x[9*xstride]; tk = x[10*xstride]; tQ = x[11*xstride]; tc = x[12*xstride]; tI = x[13*xstride]; ts = x[14*xstride]; tY = x[15*xstride]; t2 = x[16*xstride]; ty = x[17*xstride]; ti = x[18*xstride]; tO = x[19*xstride]; ta = x[20*xstride]; tG = x[21*xstride]; tq = x[22*xstride]; tW = x[23*xstride]; t6 = x[24*xstride]; tC = x[25*xstride]; tm = x[26*xstride]; tS = x[27*xstride]; te = x[28*xstride]; tK = x[29*xstride]; tu = x[30*xstride]; t_ = x[31*xstride]; t1 = x[32*xstride]; tx = x[33*xstride]; th = x[34*xstride]; tN = x[35*xstride]; t9 = x[36*xstride]; tF = x[37*xstride]; tp = x[38*xstride]; tV = x[39*xstride]; t5 = x[40*xstride]; tB = x[41*xstride]; tl = x[42*xstride]; tR = x[43*xstride]; td = x[44*xstride]; tJ = x[45*xstride]; tt = x[46*xstride]; tZ = x[47*xstride]; t3 = x[48*xstride]; tz = x[49*xstride]; tj = x[50*xstride]; tP = x[51*xstride]; tb = x[52*xstride]; tH = x[53*xstride]; tr = x[54*xstride]; tX = x[55*xstride]; t7 = x[56*xstride]; tD = x[57*xstride]; tn = x[58*xstride]; tT = x[59*xstride]; tf = x[60*xstride]; tL = x[61*xstride]; tv = x[62*xstride]; t = x[63*xstride]; OD_FDCT_64(t0, tw, tg, tM, t8, tE, to, tU, t4, tA, tk, tQ, tc, tI, ts, tY, t2, ty, ti, tO, ta, tG, tq, tW, t6, tC, tm, tS, te, tK, tu, t_, t1, tx, th, tN, t9, tF, tp, tV, t5, tB, tl, tR, td, tJ, tt, tZ, t3, tz, tj, tP, tb, tH, tr, tX, t7, tD, tn, tT, tf, tL, tv, t); y[0] = (od_coeff)t0; y[1] = (od_coeff)t1; y[2] = (od_coeff)t2; y[3] = (od_coeff)t3; y[4] = (od_coeff)t4; y[5] = (od_coeff)t5; y[6] = (od_coeff)t6; y[7] = (od_coeff)t7; y[8] = (od_coeff)t8; y[9] = (od_coeff)t9; y[10] = (od_coeff)ta; y[11] = (od_coeff)tb; y[12] = (od_coeff)tc; y[13] = (od_coeff)td; y[14] = (od_coeff)te; y[15] = (od_coeff)tf; y[16] = (od_coeff)tg; y[17] = (od_coeff)th; y[18] = (od_coeff)ti; y[19] = (od_coeff)tj; y[20] = (od_coeff)tk; y[21] = (od_coeff)tl; y[22] = (od_coeff)tm; y[23] = (od_coeff)tn; y[24] = (od_coeff)to; y[25] = (od_coeff)tp; y[26] = (od_coeff)tq; y[27] = (od_coeff)tr; y[28] = (od_coeff)ts; y[29] = (od_coeff)tt; y[30] = (od_coeff)tu; y[31] = (od_coeff)tv; y[32] = (od_coeff)tw; y[33] = (od_coeff)tx; y[34] = (od_coeff)ty; y[35] = (od_coeff)tz; y[36] = (od_coeff)tA; y[37] = (od_coeff)tB; y[38] = (od_coeff)tC; y[39] = (od_coeff)tD; y[40] = (od_coeff)tE; y[41] = (od_coeff)tF; y[41] = (od_coeff)tF; y[42] = (od_coeff)tG; y[43] = (od_coeff)tH; y[44] = (od_coeff)tI; y[45] = (od_coeff)tJ; y[46] = (od_coeff)tK; y[47] = (od_coeff)tL; y[48] = (od_coeff)tM; y[49] = (od_coeff)tN; y[50] = (od_coeff)tO; y[51] = (od_coeff)tP; y[52] = (od_coeff)tQ; y[53] = (od_coeff)tR; y[54] = (od_coeff)tS; y[55] = (od_coeff)tT; y[56] = (od_coeff)tU; y[57] = (od_coeff)tV; y[58] = (od_coeff)tW; y[59] = (od_coeff)tX; y[60] = (od_coeff)tY; y[61] = (od_coeff)tZ; y[62] = (od_coeff)t_; y[63] = (od_coeff)t; } void od_bin_idct64(od_coeff *x, int xstride, const od_coeff y[64]) { int t0; int t1; int t2; int t3; int t4; int t5; int t6; int t7; int t8; int t9; int ta; int tb; int tc; int td; int te; int tf; int tg; int th; int ti; int tj; int tk; int tl; int tm; int tn; int to; int tp; int tq; int tr; int ts; int tt; int tu; int tv; int tw; int tx; int ty; int tz; int tA; int tB; int tC; int tD; int tE; int tF; int tG; int tH; int tI; int tJ; int tK; int tL; int tM; int tN; int tO; int tP; int tQ; int tR; int tS; int tT; int tU; int tV; int tW; int tX; int tY; int tZ; int t_; int t; t0 = y[0]; tw = y[1]; tg = y[2]; tM = y[3]; t8 = y[4]; tE = y[5]; to = y[6]; tU = y[7]; t4 = y[8]; tA = y[9]; tk = y[10]; tQ = y[11]; tc = y[12]; tI = y[13]; ts = y[14]; tY = y[15]; t2 = y[16]; ty = y[17]; ti = y[18]; tO = y[19]; ta = y[20]; tG = y[21]; tq = y[22]; tW = y[23]; t6 = y[24]; tC = y[25]; tm = y[26]; tS = y[27]; te = y[28]; tK = y[29]; tu = y[30]; t_ = y[31]; t1 = y[32]; tx = y[33]; th = y[34]; tN = y[35]; t9 = y[36]; tF = y[37]; tp = y[38]; tV = y[39]; t5 = y[40]; tB = y[41]; tl = y[42]; tR = y[43]; td = y[44]; tJ = y[45]; tt = y[46]; tZ = y[47]; t3 = y[48]; tz = y[49]; tj = y[50]; tP = y[51]; tb = y[52]; tH = y[53]; tr = y[54]; tX = y[55]; t7 = y[56]; tD = y[57]; tn = y[58]; tT = y[59]; tf = y[60]; tL = y[61]; tv = y[62]; t = y[63]; OD_IDCT_64(t0, tw, tg, tM, t8, tE, to, tU, t4, tA, tk, tQ, tc, tI, ts, tY, t2, ty, ti, tO, ta, tG, tq, tW, t6, tC, tm, tS, te, tK, tu, t_, t1, tx, th, tN, t9, tF, tp, tV, t5, tB, tl, tR, td, tJ, tt, tZ, t3, tz, tj, tP, tb, tH, tr, tX, t7, tD, tn, tT, tf, tL, tv, t); x[0*xstride] = (od_coeff)t0; x[1*xstride] = (od_coeff)t1; x[2*xstride] = (od_coeff)t2; x[3*xstride] = (od_coeff)t3; x[4*xstride] = (od_coeff)t4; x[5*xstride] = (od_coeff)t5; x[6*xstride] = (od_coeff)t6; x[7*xstride] = (od_coeff)t7; x[8*xstride] = (od_coeff)t8; x[9*xstride] = (od_coeff)t9; x[10*xstride] = (od_coeff)ta; x[11*xstride] = (od_coeff)tb; x[12*xstride] = (od_coeff)tc; x[13*xstride] = (od_coeff)td; x[14*xstride] = (od_coeff)te; x[15*xstride] = (od_coeff)tf; x[16*xstride] = (od_coeff)tg; x[17*xstride] = (od_coeff)th; x[18*xstride] = (od_coeff)ti; x[19*xstride] = (od_coeff)tj; x[20*xstride] = (od_coeff)tk; x[21*xstride] = (od_coeff)tl; x[22*xstride] = (od_coeff)tm; x[23*xstride] = (od_coeff)tn; x[24*xstride] = (od_coeff)to; x[25*xstride] = (od_coeff)tp; x[26*xstride] = (od_coeff)tq; x[27*xstride] = (od_coeff)tr; x[28*xstride] = (od_coeff)ts; x[29*xstride] = (od_coeff)tt; x[30*xstride] = (od_coeff)tu; x[31*xstride] = (od_coeff)tv; x[32*xstride] = (od_coeff)tw; x[33*xstride] = (od_coeff)tx; x[34*xstride] = (od_coeff)ty; x[35*xstride] = (od_coeff)tz; x[36*xstride] = (od_coeff)tA; x[37*xstride] = (od_coeff)tB; x[38*xstride] = (od_coeff)tC; x[39*xstride] = (od_coeff)tD; x[40*xstride] = (od_coeff)tE; x[41*xstride] = (od_coeff)tF; x[41*xstride] = (od_coeff)tF; x[42*xstride] = (od_coeff)tG; x[43*xstride] = (od_coeff)tH; x[44*xstride] = (od_coeff)tI; x[45*xstride] = (od_coeff)tJ; x[46*xstride] = (od_coeff)tK; x[47*xstride] = (od_coeff)tL; x[48*xstride] = (od_coeff)tM; x[49*xstride] = (od_coeff)tN; x[50*xstride] = (od_coeff)tO; x[51*xstride] = (od_coeff)tP; x[52*xstride] = (od_coeff)tQ; x[53*xstride] = (od_coeff)tR; x[54*xstride] = (od_coeff)tS; x[55*xstride] = (od_coeff)tT; x[56*xstride] = (od_coeff)tU; x[57*xstride] = (od_coeff)tV; x[58*xstride] = (od_coeff)tW; x[59*xstride] = (od_coeff)tX; x[60*xstride] = (od_coeff)tY; x[61*xstride] = (od_coeff)tZ; x[62*xstride] = (od_coeff)t_; x[63*xstride] = (od_coeff)t; } #endif void daala_fdct4(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[4]; od_coeff y[4]; for (i = 0; i < 4; i++) x[i] = (od_coeff)input[i]; od_bin_fdct4(y, x, 1); for (i = 0; i < 4; i++) output[i] = (tran_low_t)y[i]; } void daala_idct4(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[4]; od_coeff y[4]; for (i = 0; i < 4; i++) y[i] = input[i]; od_bin_idct4(x, 1, y); for (i = 0; i < 4; i++) output[i] = (tran_low_t)x[i]; } void daala_fdst4(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[4]; od_coeff y[4]; for (i = 0; i < 4; i++) x[i] = (od_coeff)input[i]; od_bin_fdst4(y, x, 1); for (i = 0; i < 4; i++) output[i] = (tran_low_t)y[i]; } void daala_idst4(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[4]; od_coeff y[4]; for (i = 0; i < 4; i++) y[i] = input[i]; od_bin_idst4(x, 1, y); for (i = 0; i < 4; i++) output[i] = (tran_low_t)x[i]; } void daala_idtx4(const tran_low_t *input, tran_low_t *output) { int i; for (i = 0; i < 4; i++) output[i] = input[i]; } void daala_fdct8(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[8]; od_coeff y[8]; for (i = 0; i < 8; i++) x[i] = (od_coeff)input[i]; od_bin_fdct8(y, x, 1); for (i = 0; i < 8; i++) output[i] = (tran_low_t)y[i]; } void daala_idct8(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[8]; od_coeff y[8]; for (i = 0; i < 8; i++) y[i] = (od_coeff)input[i]; od_bin_idct8(x, 1, y); for (i = 0; i < 8; i++) output[i] = (tran_low_t)x[i]; } void daala_fdst8(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[8]; od_coeff y[8]; for (i = 0; i < 8; i++) x[i] = (od_coeff)input[i]; od_bin_fdst8(y, x, 1); for (i = 0; i < 8; i++) output[i] = (tran_low_t)y[i]; } void daala_idst8(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[8]; od_coeff y[8]; for (i = 0; i < 8; i++) y[i] = (od_coeff)input[i]; od_bin_idst8(x, 1, y); for (i = 0; i < 8; i++) output[i] = (tran_low_t)x[i]; } void daala_idtx8(const tran_low_t *input, tran_low_t *output) { int i; for (i = 0; i < 8; i++) output[i] = input[i]; } void daala_fdct16(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[16]; od_coeff y[16]; for (i = 0; i < 16; i++) x[i] = (od_coeff)input[i]; od_bin_fdct16(y, x, 1); for (i = 0; i < 16; i++) output[i] = (tran_low_t)y[i]; } void daala_idct16(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[16]; od_coeff y[16]; for (i = 0; i < 16; i++) y[i] = (od_coeff)input[i]; od_bin_idct16(x, 1, y); for (i = 0; i < 16; i++) output[i] = (tran_low_t)x[i]; } void daala_fdst16(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[16]; od_coeff y[16]; for (i = 0; i < 16; i++) x[i] = (od_coeff)input[i]; od_bin_fdst16(y, x, 1); for (i = 0; i < 16; i++) output[i] = (tran_low_t)y[i]; } void daala_idst16(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[16]; od_coeff y[16]; for (i = 0; i < 16; i++) y[i] = (od_coeff)input[i]; od_bin_idst16(x, 1, y); for (i = 0; i < 16; i++) output[i] = (tran_low_t)x[i]; } void daala_idtx16(const tran_low_t *input, tran_low_t *output) { int i; for (i = 0; i < 16; i++) output[i] = input[i]; } void daala_fdct32(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[32]; od_coeff y[32]; for (i = 0; i < 32; i++) x[i] = (od_coeff)input[i]; od_bin_fdct32(y, x, 1); for (i = 0; i < 32; i++) output[i] = (tran_low_t)y[i]; } void daala_idct32(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[32]; od_coeff y[32]; for (i = 0; i < 32; i++) y[i] = (od_coeff)input[i]; od_bin_idct32(x, 1, y); for (i = 0; i < 32; i++) output[i] = (tran_low_t)x[i]; } /* Preserve the "half-right" transform behavior. */ void daala_fdst32(const tran_low_t *input, tran_low_t *output) { int i; tran_low_t inputhalf[16]; for (i = 0; i < 16; ++i) { output[16 + i] = input[i]; } for (i = 0; i < 16; ++i) { inputhalf[i] = input[i + 16]; } daala_fdct16(inputhalf, output); } /* Preserve the "half-right" transform behavior. */ void daala_idst32(const tran_low_t *input, tran_low_t *output) { int i; tran_low_t inputhalf[16]; for (i = 0; i < 16; ++i) { inputhalf[i] = input[i]; } for (i = 0; i < 16; ++i) { output[i] = input[16 + i]; } daala_idct16(inputhalf, output + 16); } void daala_idtx32(const tran_low_t *input, tran_low_t *output) { int i; for (i = 0; i < 32; i++) output[i] = input[i]; } #if CONFIG_TX64X64 void daala_fdct64(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[64]; od_coeff y[64]; for (i = 0; i < 64; i++) x[i] = (od_coeff)input[i]; od_bin_fdct64(y, x, 1); for (i = 0; i < 64; i++) output[i] = (tran_low_t)y[i]; } void daala_idct64(const tran_low_t *input, tran_low_t *output) { int i; od_coeff x[64]; od_coeff y[64]; for (i = 0; i < 64; i++) y[i] = (od_coeff)input[i]; od_bin_idct64(x, 1, y); for (i = 0; i < 64; i++) output[i] = (tran_low_t)x[i]; } /* Preserve the "half-right" transform behavior. */ void daala_fdst64(const tran_low_t *input, tran_low_t *output) { int i; tran_low_t inputhalf[32]; for (i = 0; i < 32; ++i) { output[32 + i] = input[i]; } for (i = 0; i < 32; ++i) { inputhalf[i] = input[i + 32]; } daala_fdct32(inputhalf, output); } /* Preserve the "half-right" transform behavior. */ void daala_idst64(const tran_low_t *input, tran_low_t *output) { int i; tran_low_t inputhalf[32]; for (i = 0; i < 32; ++i) { inputhalf[i] = input[i]; } for (i = 0; i < 32; ++i) { output[i] = input[32 + i]; } daala_idct32(inputhalf, output + 32); } void daala_idtx64(const tran_low_t *input, tran_low_t *output) { int i; for (i = 0; i < 64; i++) output[i] = input[i]; } #endif