457 lines
15 KiB
C
457 lines
15 KiB
C
/******************************************************************************
|
|
*
|
|
* 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_iquant_itrans_recon.c
|
|
*
|
|
* @brief
|
|
* Contains function definitions for inverse quantization, inverse
|
|
* transform and reconstruction
|
|
*
|
|
* @author
|
|
* 100470
|
|
*
|
|
* @par List of Functions:
|
|
* - ihevc_iquant_itrans_recon_4x4_ttype1()
|
|
* - ihevc_iquant_itrans_recon_4x4()
|
|
*
|
|
* @remarks
|
|
* None
|
|
*
|
|
*******************************************************************************
|
|
*/
|
|
#include <stdio.h>
|
|
#include <string.h>
|
|
#include "ihevc_typedefs.h"
|
|
#include "ihevc_macros.h"
|
|
#include "ihevc_platform_macros.h"
|
|
#include "ihevc_defs.h"
|
|
#include "ihevc_trans_tables.h"
|
|
#include "ihevc_iquant_itrans_recon.h"
|
|
#include "ihevc_func_selector.h"
|
|
#include "ihevc_trans_macros.h"
|
|
|
|
/* All the functions here are replicated from ihevc_itrans.c and modified to */
|
|
/* include reconstruction */
|
|
|
|
/**
|
|
*******************************************************************************
|
|
*
|
|
* @brief
|
|
* This function performs inverse quantization, inverse transform
|
|
* type1(DST) and reconstruction for 4x4 input block
|
|
*
|
|
* @par Description:
|
|
* Performs inverse quantization , inverse transform type 1 and adds
|
|
* prediction data and clips output to 8 bit
|
|
*
|
|
* @param[in] pi2_src
|
|
* Input 4x4 coefficients
|
|
*
|
|
* @param[in] pi2_tmp
|
|
* Temporary 4x4 buffer for storing inverse
|
|
* transform 1st stage output
|
|
*
|
|
* @param[in] pu1_pred
|
|
* Prediction 4x4 block
|
|
*
|
|
* @param[in] pi2_dequant_coeff
|
|
* Dequant Coeffs
|
|
*
|
|
* @param[out] pu1_dst
|
|
* Output 4x4 block
|
|
*
|
|
* @param[in] qp_div
|
|
* Quantization parameter / 6
|
|
*
|
|
* @param[in] qp_rem
|
|
* Quantization parameter % 6
|
|
*
|
|
* @param[in] src_strd
|
|
* Input stride
|
|
*
|
|
* @param[in] pred_strd
|
|
* Prediction stride
|
|
*
|
|
* @param[in] dst_strd
|
|
* Output Stride
|
|
*
|
|
* @param[in] zero_cols
|
|
* Zero columns in pi2_src
|
|
*
|
|
* @param[in] zero_rows
|
|
* Zero Rows in pi2_src
|
|
*
|
|
* @returns Void
|
|
*
|
|
* @remarks
|
|
* None
|
|
*
|
|
*******************************************************************************
|
|
*/
|
|
|
|
void ihevc_iquant_itrans_recon_4x4_ttype1(WORD16 *pi2_src,
|
|
WORD16 *pi2_tmp,
|
|
UWORD8 *pu1_pred,
|
|
WORD16 *pi2_dequant_coeff,
|
|
UWORD8 *pu1_dst,
|
|
WORD32 qp_div, /* qpscaled / 6 */
|
|
WORD32 qp_rem, /* qpscaled % 6 */
|
|
WORD32 src_strd,
|
|
WORD32 pred_strd,
|
|
WORD32 dst_strd,
|
|
WORD32 zero_cols,
|
|
WORD32 zero_rows)
|
|
{
|
|
UNUSED(zero_rows);
|
|
/* Inverse Quant and Inverse Transform and Reconstruction */
|
|
{
|
|
WORD32 i, c[4];
|
|
WORD32 add;
|
|
WORD32 shift;
|
|
WORD16 *pi2_tmp_orig;
|
|
WORD32 shift_iq;
|
|
WORD32 trans_size;
|
|
/* Inverse Quantization constants */
|
|
{
|
|
WORD32 log2_trans_size, bit_depth;
|
|
|
|
log2_trans_size = 2;
|
|
bit_depth = 8 + 0;
|
|
shift_iq = bit_depth + log2_trans_size - 5;
|
|
}
|
|
|
|
trans_size = TRANS_SIZE_4;
|
|
pi2_tmp_orig = pi2_tmp;
|
|
|
|
/* Inverse Transform 1st stage */
|
|
shift = IT_SHIFT_STAGE_1;
|
|
add = 1 << (shift - 1);
|
|
|
|
for(i = 0; i < trans_size; i++)
|
|
{
|
|
/* Checking for Zero Cols */
|
|
if((zero_cols & 1) == 1)
|
|
{
|
|
memset(pi2_tmp, 0, trans_size * sizeof(WORD16));
|
|
}
|
|
else
|
|
{
|
|
WORD32 iq_tmp_1, iq_tmp_2, iq_tmp_3;
|
|
// Intermediate Variables
|
|
IQUANT_4x4(iq_tmp_1,
|
|
pi2_src[0 * src_strd],
|
|
pi2_dequant_coeff[0 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
IQUANT_4x4(iq_tmp_2,
|
|
pi2_src[2 * src_strd],
|
|
pi2_dequant_coeff[2 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
c[0] = iq_tmp_1 + iq_tmp_2;
|
|
|
|
IQUANT_4x4(iq_tmp_1,
|
|
pi2_src[2 * src_strd],
|
|
pi2_dequant_coeff[2 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
IQUANT_4x4(iq_tmp_2,
|
|
pi2_src[3 * src_strd],
|
|
pi2_dequant_coeff[3 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
c[1] = iq_tmp_1 + iq_tmp_2;
|
|
|
|
IQUANT_4x4(iq_tmp_1,
|
|
pi2_src[0 * src_strd],
|
|
pi2_dequant_coeff[0 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
IQUANT_4x4(iq_tmp_2,
|
|
pi2_src[3 * src_strd],
|
|
pi2_dequant_coeff[3 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
c[2] = iq_tmp_1 - iq_tmp_2;
|
|
|
|
IQUANT_4x4(iq_tmp_1,
|
|
pi2_src[1 * src_strd],
|
|
pi2_dequant_coeff[1 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
c[3] = 74 * iq_tmp_1;
|
|
|
|
pi2_tmp[0] =
|
|
CLIP_S16((29 * c[0] + 55 * c[1] + c[3] + add) >> shift);
|
|
pi2_tmp[1] =
|
|
CLIP_S16((55 * c[2] - 29 * c[1] + c[3] + add) >> shift);
|
|
|
|
IQUANT_4x4(iq_tmp_1,
|
|
pi2_src[0 * src_strd],
|
|
pi2_dequant_coeff[0 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
IQUANT_4x4(iq_tmp_2,
|
|
pi2_src[2 * src_strd],
|
|
pi2_dequant_coeff[2 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
IQUANT_4x4(iq_tmp_3,
|
|
pi2_src[3 * src_strd],
|
|
pi2_dequant_coeff[3 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
|
|
pi2_tmp[2] =
|
|
CLIP_S16((74 * (iq_tmp_1 - iq_tmp_2 + iq_tmp_3) + add) >> shift);
|
|
pi2_tmp[3] =
|
|
CLIP_S16((55 * c[0] + 29 * c[2] - c[3] + add) >> shift);
|
|
}
|
|
pi2_src++;
|
|
pi2_dequant_coeff++;
|
|
pi2_tmp += trans_size;
|
|
zero_cols = zero_cols >> 1;
|
|
}
|
|
|
|
pi2_tmp = pi2_tmp_orig;
|
|
|
|
/* Inverse Transform 2nd stage */
|
|
shift = IT_SHIFT_STAGE_2;
|
|
add = 1 << (shift - 1);
|
|
|
|
for(i = 0; i < trans_size; i++)
|
|
{
|
|
WORD32 itrans_out;
|
|
|
|
// Intermediate Variables
|
|
c[0] = pi2_tmp[0] + pi2_tmp[2 * trans_size];
|
|
c[1] = pi2_tmp[2 * trans_size] + pi2_tmp[3 * trans_size];
|
|
c[2] = pi2_tmp[0] - pi2_tmp[3 * trans_size];
|
|
c[3] = 74 * pi2_tmp[trans_size];
|
|
|
|
itrans_out =
|
|
CLIP_S16((29 * c[0] + 55 * c[1] + c[3] + add) >> shift);
|
|
pu1_dst[0] = CLIP_U8((itrans_out + pu1_pred[0]));
|
|
|
|
itrans_out =
|
|
CLIP_S16((55 * c[2] - 29 * c[1] + c[3] + add) >> shift);
|
|
pu1_dst[1] = CLIP_U8((itrans_out + pu1_pred[1]));
|
|
|
|
itrans_out =
|
|
CLIP_S16((74 * (pi2_tmp[0] - pi2_tmp[2 * trans_size] + pi2_tmp[3 * trans_size]) + add) >> shift);
|
|
pu1_dst[2] = CLIP_U8((itrans_out + pu1_pred[2]));
|
|
|
|
itrans_out =
|
|
CLIP_S16((55 * c[0] + 29 * c[2] - c[3] + add) >> shift);
|
|
pu1_dst[3] = CLIP_U8((itrans_out + pu1_pred[3]));
|
|
pi2_tmp++;
|
|
pu1_pred += pred_strd;
|
|
pu1_dst += dst_strd;
|
|
}
|
|
}
|
|
}
|
|
|
|
/**
|
|
*******************************************************************************
|
|
*
|
|
* @brief
|
|
* This function performs inverse quantization, inverse transform and
|
|
* reconstruction for 4x4 input block
|
|
*
|
|
* @par Description:
|
|
* Performs inverse quantization , inverse transform and adds the
|
|
* prediction data and clips output to 8 bit
|
|
*
|
|
* @param[in] pi2_src
|
|
* Input 4x4 coefficients
|
|
*
|
|
* @param[in] pi2_tmp
|
|
* Temporary 4x4 buffer for storing inverse
|
|
* transform 1st stage output
|
|
*
|
|
* @param[in] pu1_pred
|
|
* Prediction 4x4 block
|
|
*
|
|
* @param[in] pi2_dequant_coeff
|
|
* Dequant Coeffs
|
|
*
|
|
* @param[out] pu1_dst
|
|
* Output 4x4 block
|
|
*
|
|
* @param[in] qp_div
|
|
* Quantization parameter / 6
|
|
*
|
|
* @param[in] qp_rem
|
|
* Quantization parameter % 6
|
|
*
|
|
* @param[in] src_strd
|
|
* Input stride
|
|
*
|
|
* @param[in] pred_strd
|
|
* Prediction stride
|
|
*
|
|
* @param[in] dst_strd
|
|
* Output Stride
|
|
*
|
|
* @param[in] zero_cols
|
|
* Zero columns in pi2_src
|
|
*
|
|
* @param[in] zero_rows
|
|
* Zero Rows in pi2_src
|
|
*
|
|
* @returns Void
|
|
*
|
|
* @remarks
|
|
* None
|
|
*
|
|
*******************************************************************************
|
|
*/
|
|
|
|
void ihevc_iquant_itrans_recon_4x4(WORD16 *pi2_src,
|
|
WORD16 *pi2_tmp,
|
|
UWORD8 *pu1_pred,
|
|
WORD16 *pi2_dequant_coeff,
|
|
UWORD8 *pu1_dst,
|
|
WORD32 qp_div, /* qpscaled / 6 */
|
|
WORD32 qp_rem, /* qpscaled % 6 */
|
|
WORD32 src_strd,
|
|
WORD32 pred_strd,
|
|
WORD32 dst_strd,
|
|
WORD32 zero_cols,
|
|
WORD32 zero_rows)
|
|
{
|
|
UNUSED(zero_rows);
|
|
/* Inverse Transform */
|
|
{
|
|
WORD32 j;
|
|
WORD32 e[2], o[2];
|
|
WORD32 add;
|
|
WORD32 shift;
|
|
WORD16 *pi2_tmp_orig;
|
|
WORD32 shift_iq;
|
|
WORD32 trans_size;
|
|
/* Inverse Quantization constants */
|
|
{
|
|
WORD32 log2_trans_size, bit_depth;
|
|
|
|
log2_trans_size = 2;
|
|
bit_depth = 8 + 0;
|
|
shift_iq = bit_depth + log2_trans_size - 5;
|
|
}
|
|
|
|
trans_size = TRANS_SIZE_4;
|
|
pi2_tmp_orig = pi2_tmp;
|
|
|
|
/* Inverse Transform 1st stage */
|
|
shift = IT_SHIFT_STAGE_1;
|
|
add = 1 << (shift - 1);
|
|
|
|
for(j = 0; j < trans_size; j++)
|
|
{
|
|
/* Checking for Zero Cols */
|
|
if((zero_cols & 1) == 1)
|
|
{
|
|
memset(pi2_tmp, 0, trans_size * sizeof(WORD16));
|
|
}
|
|
else
|
|
{
|
|
WORD32 iq_tmp_1, iq_tmp_2;
|
|
/* Utilizing symmetry properties to the maximum to minimize the number of multiplications */
|
|
IQUANT_4x4(iq_tmp_1,
|
|
pi2_src[1 * src_strd],
|
|
pi2_dequant_coeff[1 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
IQUANT_4x4(iq_tmp_2,
|
|
pi2_src[3 * src_strd],
|
|
pi2_dequant_coeff[3 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
|
|
o[0] = g_ai2_ihevc_trans_4[1][0] * iq_tmp_1
|
|
+ g_ai2_ihevc_trans_4[3][0] * iq_tmp_2;
|
|
o[1] = g_ai2_ihevc_trans_4[1][1] * iq_tmp_1
|
|
+ g_ai2_ihevc_trans_4[3][1] * iq_tmp_2;
|
|
|
|
IQUANT_4x4(iq_tmp_1,
|
|
pi2_src[0 * src_strd],
|
|
pi2_dequant_coeff[0 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
IQUANT_4x4(iq_tmp_2,
|
|
pi2_src[2 * src_strd],
|
|
pi2_dequant_coeff[2 * trans_size] * g_ihevc_iquant_scales[qp_rem],
|
|
shift_iq, qp_div);
|
|
|
|
e[0] = g_ai2_ihevc_trans_4[0][0] * iq_tmp_1
|
|
+ g_ai2_ihevc_trans_4[2][0] * iq_tmp_2;
|
|
e[1] = g_ai2_ihevc_trans_4[0][1] * iq_tmp_1
|
|
+ g_ai2_ihevc_trans_4[2][1] * iq_tmp_2;
|
|
|
|
pi2_tmp[0] =
|
|
CLIP_S16(((e[0] + o[0] + add) >> shift));
|
|
pi2_tmp[1] =
|
|
CLIP_S16(((e[1] + o[1] + add) >> shift));
|
|
pi2_tmp[2] =
|
|
CLIP_S16(((e[1] - o[1] + add) >> shift));
|
|
pi2_tmp[3] =
|
|
CLIP_S16(((e[0] - o[0] + add) >> shift));
|
|
}
|
|
pi2_src++;
|
|
pi2_dequant_coeff++;
|
|
pi2_tmp += trans_size;
|
|
zero_cols = zero_cols >> 1;
|
|
}
|
|
|
|
pi2_tmp = pi2_tmp_orig;
|
|
|
|
/* Inverse Transform 2nd stage */
|
|
shift = IT_SHIFT_STAGE_2;
|
|
add = 1 << (shift - 1);
|
|
|
|
for(j = 0; j < trans_size; j++)
|
|
{
|
|
WORD32 itrans_out;
|
|
|
|
/* Utilizing symmetry properties to the maximum to minimize the number of multiplications */
|
|
o[0] = g_ai2_ihevc_trans_4[1][0] * pi2_tmp[trans_size]
|
|
+ g_ai2_ihevc_trans_4[3][0]
|
|
* pi2_tmp[3 * trans_size];
|
|
o[1] = g_ai2_ihevc_trans_4[1][1] * pi2_tmp[trans_size]
|
|
+ g_ai2_ihevc_trans_4[3][1]
|
|
* pi2_tmp[3 * trans_size];
|
|
e[0] = g_ai2_ihevc_trans_4[0][0] * pi2_tmp[0]
|
|
+ g_ai2_ihevc_trans_4[2][0]
|
|
* pi2_tmp[2 * trans_size];
|
|
e[1] = g_ai2_ihevc_trans_4[0][1] * pi2_tmp[0]
|
|
+ g_ai2_ihevc_trans_4[2][1]
|
|
* pi2_tmp[2 * trans_size];
|
|
|
|
itrans_out =
|
|
CLIP_S16(((e[0] + o[0] + add) >> shift));
|
|
pu1_dst[0] = CLIP_U8((itrans_out + pu1_pred[0]));
|
|
|
|
itrans_out =
|
|
CLIP_S16(((e[1] + o[1] + add) >> shift));
|
|
pu1_dst[1] = CLIP_U8((itrans_out + pu1_pred[1]));
|
|
|
|
itrans_out =
|
|
CLIP_S16(((e[1] - o[1] + add) >> shift));
|
|
pu1_dst[2] = CLIP_U8((itrans_out + pu1_pred[2]));
|
|
|
|
itrans_out =
|
|
CLIP_S16(((e[0] - o[0] + add) >> shift));
|
|
pu1_dst[3] = CLIP_U8((itrans_out + pu1_pred[3]));
|
|
|
|
pi2_tmp++;
|
|
pu1_pred += pred_strd;
|
|
pu1_dst += dst_strd;
|
|
|
|
}
|
|
}
|
|
}
|