blob: 079784de7764a9d2121f78ff0242a7f1018c500e [file] [log] [blame]
/******************************************************************************
*
* Copyright (C) 2012 Ittiam Systems Pvt Ltd, Bangalore
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at:
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
******************************************************************************/
/**
*******************************************************************************
* @file
* ihevc_trans_macros.h
*
* @brief
* Macros used in the forward transform and inverse transform functions
*
* @author
* Ittiam
*
* @remarks
* None
*
*******************************************************************************
*/
#ifndef IHEVC_TRANS_MACROS_H_
#define IHEVC_TRANS_MACROS_H_
#define QUANT(out, inp, quant_coeff, qp_div, log2_trans_size, q_add) \
{ \
LWORD64 tmp; \
WORD32 sign; \
WORD32 bit_depth,transform_shift; \
WORD32 q_bits, quant_multiplier; \
\
/* q_bits and q_add calculation*/ \
/* To be moved outside in neon. To be computer once per transform call */ \
bit_depth = 8; \
transform_shift = MAX_TR_DYNAMIC_RANGE - bit_depth - log2_trans_size; \
quant_multiplier = 4 ; /* because quant_coeff are multiplied by 16. Instead of multiplying, we can reduce the division factor q_bits by 4 */ \
q_bits = QUANT_SHIFT + qp_div + transform_shift + SCALING_Q_SHIFT - quant_multiplier ; \
\
sign = (inp)<0 ? -1:1; \
\
tmp = (LWORD64)(abs(inp)); \
tmp = tmp * (quant_coeff); \
tmp = tmp + (((LWORD64)q_add) << (q_bits - QUANT_ROUND_FACTOR_Q)); \
tmp = tmp >> q_bits; \
\
tmp = tmp * sign; \
out = (WORD16) CLIP_S16(tmp); \
} \
#define QUANT_HBD(out, inp, quant_coeff, qp_div, log2_trans_size, q_add, bit_depth) \
{ \
LWORD64 tmp; \
WORD32 sign; \
WORD32 transform_shift; \
WORD32 q_bits, quant_multiplier; \
\
/* q_bits and q_add calculation*/ \
/* To be moved outside in neon. To be computer once per transform call */ \
\
transform_shift = MAX_TR_DYNAMIC_RANGE - bit_depth - log2_trans_size; \
quant_multiplier = 4 ; /* because quant_coeff are multiplied by 16. Instead of multiplying, we can reduce the division factor q_bits by 4 */ \
q_bits = QUANT_SHIFT + qp_div + transform_shift + SCALING_Q_SHIFT - quant_multiplier ; \
\
sign = (inp)<0 ? -1:1; \
\
tmp = (LWORD64)(abs(inp)); \
tmp = tmp * (quant_coeff); \
tmp = tmp + (((LWORD64)q_add) << (q_bits - QUANT_ROUND_FACTOR_Q)); \
tmp = tmp >> q_bits; \
\
tmp = tmp * sign; \
out = (WORD16) CLIP_S16(tmp); \
}
/* added by 100028 */
#define QUANT_NO_WEIGHTMAT(out, inp, quant_coeff, qp_div, log2_trans_size, q_add) \
{ \
WORD32 tmp; \
WORD32 sign; \
WORD32 bit_depth,transform_shift; \
WORD32 q_bits, quant_multiplier; \
\
/* q_bits and q_add calculation*/ \
/* To be moved outside in neon. To be computer once per transform call */ \
bit_depth = 8; \
transform_shift = MAX_TR_DYNAMIC_RANGE - bit_depth - log2_trans_size; \
quant_multiplier = 4 ; /* because quant_coeff are multiplied by 16. Instead of multiplying, we can reduce the division factor q_bits by 4 */ \
q_bits = QUANT_SHIFT + qp_div + transform_shift + SCALING_Q_SHIFT - quant_multiplier - FLAT_RESCALE_MAT_Q_SHIFT /* 2048 */; \
\
sign = (inp)<0 ? -1:1; \
\
tmp = (WORD32)(abs(inp)); \
tmp = tmp * (quant_coeff); \
tmp = tmp + (((WORD32)q_add) << (q_bits - QUANT_ROUND_FACTOR_Q)); \
tmp = tmp >> q_bits; \
\
tmp = tmp * sign; \
out = (WORD16) CLIP_S16(tmp); \
}
#define QUANT_NO_WEIGHTMAT_HBD(out, inp, quant_coeff, qp_div, log2_trans_size, q_add, bit_depth) \
{ \
WORD32 tmp; \
WORD32 sign; \
WORD32 transform_shift; \
WORD32 q_bits, quant_multiplier; \
\
/* q_bits and q_add calculation*/ \
/* To be moved outside in neon. To be computer once per transform call */ \
\
transform_shift = MAX_TR_DYNAMIC_RANGE - bit_depth - log2_trans_size; \
quant_multiplier = 4 ; /* because quant_coeff are multiplied by 16. Instead of multiplying, we can reduce the division factor q_bits by 4 */ \
q_bits = QUANT_SHIFT + qp_div + transform_shift + SCALING_Q_SHIFT - quant_multiplier - FLAT_RESCALE_MAT_Q_SHIFT /* 2048 */; \
\
sign = (inp)<0 ? -1:1; \
\
tmp = (WORD32)(abs(inp)); \
tmp = tmp * (quant_coeff); \
tmp = tmp + (((WORD32)q_add) << (q_bits - QUANT_ROUND_FACTOR_Q)); \
tmp = tmp >> q_bits; \
\
tmp = tmp * sign; \
out = (WORD16) CLIP_S16(tmp); \
}
/* Reference Inverse Quantization: "pi2_src"(Coefficients) will be clipped to 15 or 14 bits when (qp_div > shift_iq). Spec doesn't have any clip mentioned */
/* Inverse quantization other than 4x4 */
/* No clipping is needed for "pi2_src"(coefficients) */
#define IQUANT(res, coeff /*pi2_src[index*src_strd]*/, dequant_coeff /*pi2_dequant_coeff[index*trans_size] * g_ihevc_iquant_scales[qp_rem] */, shift_iq, qp_div) \
{ \
WORD32 tmp, add_iq; \
\
add_iq = SHL_NEG(1 , (shift_iq - qp_div - 1)); /* To be moved outside in neon. To be computed once per transform call */ \
\
tmp = coeff * dequant_coeff ; \
tmp = tmp + add_iq; \
tmp = SHR_NEG(tmp,(shift_iq - qp_div)); \
\
res = CLIP_S16(tmp); \
}
/* 4x4 inverse quantization */
/* Options : */
/* 1. Clip "pi2_src"(coefficients) to 10 bits if "(qp_div >= shift_iq)" or 16 bits if "(qp_div < shift_iq)"*/
/* 2. Increasing precision of "pi2_src"(coefficients) to 64 bits */
#define IQUANT_4x4(res, coeff /*pi2_src[index*src_strd]*/, dequant_coeff /*pi2_dequant_coeff[index*trans_size] * g_ihevc_iquant_scales[qp_rem] */, shift_iq, qp_div) \
{ \
WORD32 clip_coeff, tmp; \
WORD32 coeff_min,coeff_max; \
WORD32 coeff_bit_range; \
WORD32 add_iq; \
add_iq = SHL_NEG(1 , (shift_iq - qp_div - 1)); /* To be moved outside in neon. To be computed once per transform call */ \
\
coeff_bit_range = 16; \
if(qp_div > shift_iq) \
coeff_bit_range = 10; \
\
coeff_min = -(1<<(coeff_bit_range-1)); \
coeff_max = (1<<(coeff_bit_range-1)) - 1; \
\
clip_coeff = CLIP3(coeff,coeff_min,coeff_max); \
\
tmp = clip_coeff * dequant_coeff ; \
tmp = tmp + add_iq; \
tmp = SHR_NEG(tmp,(shift_iq - qp_div)); \
\
res = CLIP_S16(tmp); \
}
#endif /* IHEVC_TRANS_MACROS_H_ */