michael@0: /*
michael@0:  *  Copyright (c) 2010 The WebM project authors. All Rights Reserved.
michael@0:  *
michael@0:  *  Use of this source code is governed by a BSD-style license
michael@0:  *  that can be found in the LICENSE file in the root of the source
michael@0:  *  tree. An additional intellectual property rights grant can be found
michael@0:  *  in the file PATENTS.  All contributing project authors may
michael@0:  *  be found in the AUTHORS file in the root of the source tree.
michael@0:  */
michael@0: 
michael@0: #ifndef VP9_COMMON_VP9_IDCT_H_
michael@0: #define VP9_COMMON_VP9_IDCT_H_
michael@0: 
michael@0: #include <assert.h>
michael@0: 
michael@0: #include "./vpx_config.h"
michael@0: #include "vpx/vpx_integer.h"
michael@0: #include "vp9/common/vp9_common.h"
michael@0: #include "vp9/common/vp9_enums.h"
michael@0: 
michael@0: 
michael@0: // Constants and Macros used by all idct/dct functions
michael@0: #define DCT_CONST_BITS 14
michael@0: #define DCT_CONST_ROUNDING  (1 << (DCT_CONST_BITS - 1))
michael@0: 
michael@0: #define UNIT_QUANT_SHIFT 2
michael@0: #define UNIT_QUANT_FACTOR (1 << UNIT_QUANT_SHIFT)
michael@0: 
michael@0: #define pair_set_epi16(a, b) \
michael@0:   _mm_set_epi16(b, a, b, a, b, a, b, a)
michael@0: 
michael@0: #define pair_set_epi32(a, b) \
michael@0:   _mm_set_epi32(b, a, b, a)
michael@0: 
michael@0: // Constants:
michael@0: //  for (int i = 1; i< 32; ++i)
michael@0: //    printf("static const int cospi_%d_64 = %.0f;\n", i,
michael@0: //           round(16384 * cos(i*M_PI/64)));
michael@0: // Note: sin(k*Pi/64) = cos((32-k)*Pi/64)
michael@0: static const int cospi_1_64  = 16364;
michael@0: static const int cospi_2_64  = 16305;
michael@0: static const int cospi_3_64  = 16207;
michael@0: static const int cospi_4_64  = 16069;
michael@0: static const int cospi_5_64  = 15893;
michael@0: static const int cospi_6_64  = 15679;
michael@0: static const int cospi_7_64  = 15426;
michael@0: static const int cospi_8_64  = 15137;
michael@0: static const int cospi_9_64  = 14811;
michael@0: static const int cospi_10_64 = 14449;
michael@0: static const int cospi_11_64 = 14053;
michael@0: static const int cospi_12_64 = 13623;
michael@0: static const int cospi_13_64 = 13160;
michael@0: static const int cospi_14_64 = 12665;
michael@0: static const int cospi_15_64 = 12140;
michael@0: static const int cospi_16_64 = 11585;
michael@0: static const int cospi_17_64 = 11003;
michael@0: static const int cospi_18_64 = 10394;
michael@0: static const int cospi_19_64 = 9760;
michael@0: static const int cospi_20_64 = 9102;
michael@0: static const int cospi_21_64 = 8423;
michael@0: static const int cospi_22_64 = 7723;
michael@0: static const int cospi_23_64 = 7005;
michael@0: static const int cospi_24_64 = 6270;
michael@0: static const int cospi_25_64 = 5520;
michael@0: static const int cospi_26_64 = 4756;
michael@0: static const int cospi_27_64 = 3981;
michael@0: static const int cospi_28_64 = 3196;
michael@0: static const int cospi_29_64 = 2404;
michael@0: static const int cospi_30_64 = 1606;
michael@0: static const int cospi_31_64 = 804;
michael@0: 
michael@0: //  16384 * sqrt(2) * sin(kPi/9) * 2 / 3
michael@0: static const int sinpi_1_9 = 5283;
michael@0: static const int sinpi_2_9 = 9929;
michael@0: static const int sinpi_3_9 = 13377;
michael@0: static const int sinpi_4_9 = 15212;
michael@0: 
michael@0: static INLINE int dct_const_round_shift(int input) {
michael@0:   int rv = ROUND_POWER_OF_TWO(input, DCT_CONST_BITS);
michael@0:   return (int16_t)rv;
michael@0: }
michael@0: 
michael@0: typedef void (*transform_1d)(const int16_t*, int16_t*);
michael@0: 
michael@0: typedef struct {
michael@0:   transform_1d cols, rows;  // vertical and horizontal
michael@0: } transform_2d;
michael@0: 
michael@0: void vp9_iwht4x4_add(const int16_t *input, uint8_t *dest, int stride, int eob);
michael@0: 
michael@0: void vp9_idct4x4_add(const int16_t *input, uint8_t *dest, int stride, int eob);
michael@0: void vp9_idct8x8_add(const int16_t *input, uint8_t *dest, int stride, int eob);
michael@0: void vp9_idct16x16_add(const int16_t *input, uint8_t *dest, int stride, int
michael@0:                        eob);
michael@0: void vp9_idct32x32_add(const int16_t *input, uint8_t *dest, int stride,
michael@0:                        int eob);
michael@0: 
michael@0: void vp9_iht4x4_add(TX_TYPE tx_type, const int16_t *input, uint8_t *dest,
michael@0:                     int stride, int eob);
michael@0: void vp9_iht8x8_add(TX_TYPE tx_type, const int16_t *input, uint8_t *dest,
michael@0:                     int stride, int eob);
michael@0: void vp9_iht16x16_add(TX_TYPE tx_type, const int16_t *input, uint8_t *dest,
michael@0:                       int stride, int eob);
michael@0: 
michael@0: 
michael@0: #endif  // VP9_COMMON_VP9_IDCT_H_