1  /* { dg-do compile } */
       2  /* { dg-options "-O3 -Wuninitialized --param vect-max-version-for-alias-checks=20" } */
       3  
       4  #include <stdint.h>
       5  
       6  #define A1  2896 /* (1/sqrt(2))<<12 */
       7  #define A2  2217
       8  #define A3  3784
       9  #define A4 -5352
      10  
      11  #define IDCT_TRANSFORM(dest,s0,s1,s2,s3,s4,s5,s6,s7,d0,d1,d2,d3,d4,d5,d6,d7,munge,src) {\
      12      const int a0 = (src)[s0] + (src)[s4]; \
      13      const int a1 = (src)[s0] - (src)[s4]; \
      14      const int a2 = (src)[s2] + (src)[s6]; \
      15      const int a3 = (A1*((src)[s2] - (src)[s6])) >> 11; \
      16      const int a4 = (src)[s5] + (src)[s3]; \
      17      const int a5 = (src)[s5] - (src)[s3]; \
      18      const int a6 = (src)[s1] + (src)[s7]; \
      19      const int a7 = (src)[s1] - (src)[s7]; \
      20      const int b0 = a4 + a6; \
      21      const int b1 = (A3*(a5 + a7)) >> 11; \
      22      const int b2 = ((A4*a5) >> 11) - b0 + b1; \
      23      const int b3 = (A1*(a6 - a4) >> 11) - b2; \
      24      const int b4 = ((A2*a7) >> 11) + b3 - b1; \
      25      (dest)[d0] = munge(a0+a2   +b0); \
      26      (dest)[d1] = munge(a1+a3-a2+b2); \
      27      (dest)[d2] = munge(a1-a3+a2+b3); \
      28      (dest)[d3] = munge(a0-a2   -b4); \
      29      (dest)[d4] = munge(a0-a2   +b4); \
      30      (dest)[d5] = munge(a1-a3+a2-b3); \
      31      (dest)[d6] = munge(a1+a3-a2-b2); \
      32      (dest)[d7] = munge(a0+a2   -b0); \
      33  }
      34  
      35  #define MUNGE_NONE(x) (x)
      36  #define IDCT_COL(dest,src) IDCT_TRANSFORM(dest,0,8,16,24,32,40,48,56,0,8,16,24,32,40,48,56,MUNGE_NONE,src)
      37  
      38  #define MUNGE_ROW(x) (((x) + 0x7F)>>8)
      39  #define IDCT_ROW(dest,src) IDCT_TRANSFORM(dest,0,1,2,3,4,5,6,7,0,1,2,3,4,5,6,7,MUNGE_ROW,src)
      40  
      41  static inline void bink_idct_col(int *dest, const int32_t *src)
      42  {
      43      if ((src[8]|src[16]|src[24]|src[32]|src[40]|src[48]|src[56])==0) {
      44          dest[0]  =
      45          dest[8]  =
      46          dest[16] =
      47          dest[24] =
      48          dest[32] =
      49          dest[40] =
      50          dest[48] =
      51          dest[56] = src[0];
      52      } else {
      53          IDCT_COL(dest, src);
      54      }
      55  }
      56  
      57  int bink_idct_put_c(uint8_t *dest, int linesize, int32_t *block)
      58  {
      59      int i;
      60      int temp[64];
      61      for (i = 0; i < 8; i++)
      62          bink_idct_col(&temp[i], &block[i]);
      63      for (i = 0; i < 8; i++) {
      64          IDCT_ROW( (&dest[i*linesize]), (&temp[8*i]) );
      65      }
      66  
      67      return 0;
      68  }
      69