/****************************************************************************** * * 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_itrans_recon_atom_intr.c * * @brief * Contains function definitions for inverse quantization, inverse * transform and reconstruction * * @author * 100470 * 100592 (edited by) * * @par List of Functions: * - ihevc_itrans_recon_4x4_ttype1_ssse3() * - ihevc_itrans_recon_4x4_ssse3() * - ihevc_itrans_recon_8x8_ssse3() * * @remarks * None * ******************************************************************************* */ #include <stdio.h> #include <string.h> #include "ihevc_typedefs.h" #include "ihevc_platform_macros.h" #include "ihevc_macros.h" #include "ihevc_defs.h" #include "ihevc_func_selector.h" #include "ihevc_trans_tables.h" #include "ihevc_iquant_itrans_recon.h" #include "ihevc_trans_macros.h" #include <immintrin.h> #include <emmintrin.h> #include <tmmintrin.h> /** ******************************************************************************* * * @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 * * @returns Void * * @remarks * None * ******************************************************************************* */ void ihevc_itrans_recon_4x4_ttype1_ssse3(WORD16 *pi2_src, WORD16 *pi2_tmp, UWORD8 *pu1_pred, UWORD8 *pu1_dst, WORD32 src_strd, WORD32 pred_strd, WORD32 dst_strd, WORD32 zero_cols, WORD32 zero_rows) { __m128i m_temp_reg_0; __m128i m_temp_reg_1; __m128i m_temp_reg_2; __m128i m_temp_reg_3; __m128i m_temp_reg_4; __m128i m_temp_reg_10; __m128i m_temp_reg_11; __m128i m_temp_reg_12; __m128i m_temp_reg_13; __m128i m_temp_reg_14; __m128i m_temp_reg_20; __m128i m_temp_reg_21; __m128i m_temp_reg_22; __m128i m_temp_reg_23; __m128i m_temp_reg_24; __m128i m_temp_reg_25; __m128i m_temp_reg_30; __m128i m_temp_reg_31; __m128i m_temp_reg_32; __m128i m_temp_reg_33; __m128i m_temp_reg_34; __m128i m_temp_reg_35; __m128i m_temp_reg_36; __m128i m_rdng_factor; __m128i m_count; __m128i m_ge_zero16b_flag_row0; __m128i m_ge_zero16b_flag_row1; __m128i m_ge_zero16b_flag_row2; __m128i m_ge_zero16b_flag_row3; __m128i m_zero = _mm_setzero_si128(); WORD32 i4_shift = IT_SHIFT_STAGE_1; UNUSED(zero_cols); UNUSED(zero_rows); UNUSED(pi2_tmp); m_temp_reg_0 = _mm_loadl_epi64((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_1 = _mm_loadl_epi64((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_2 = _mm_loadl_epi64((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_3 = _mm_loadl_epi64((__m128i *)pi2_src); m_ge_zero16b_flag_row0 = _mm_cmpgt_epi16(m_zero, m_temp_reg_0); m_ge_zero16b_flag_row1 = _mm_cmpgt_epi16(m_zero, m_temp_reg_1); m_ge_zero16b_flag_row2 = _mm_cmpgt_epi16(m_zero, m_temp_reg_2); m_ge_zero16b_flag_row3 = _mm_cmpgt_epi16(m_zero, m_temp_reg_3); m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_0, m_ge_zero16b_flag_row0); m_temp_reg_1 = _mm_unpacklo_epi16(m_temp_reg_1, m_ge_zero16b_flag_row1); m_temp_reg_2 = _mm_unpacklo_epi16(m_temp_reg_2, m_ge_zero16b_flag_row2); m_temp_reg_3 = _mm_unpacklo_epi16(m_temp_reg_3, m_ge_zero16b_flag_row3); /*m_temp_reg_0 = _mm_cvtepi16_epi32(m_temp_reg_0); m_temp_reg_2 = _mm_cvtepi16_epi32(m_temp_reg_2); m_temp_reg_1 = _mm_cvtepi16_epi32(m_temp_reg_1); m_temp_reg_3 = _mm_cvtepi16_epi32(m_temp_reg_3);*/ /* c[4] in m_temp_reg_14 */ /* c[4] = src[0] - src[2] + src[3] */ { m_temp_reg_14 = _mm_sub_epi32(m_temp_reg_0, m_temp_reg_2); } /* c[3] in m_temp_reg_13 */ { m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_1, 6); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_1, 3); m_temp_reg_22 = _mm_slli_epi32(m_temp_reg_1, 1); m_temp_reg_23 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_21); m_temp_reg_13 = _mm_add_epi32(m_temp_reg_22, m_temp_reg_23); //m_temp_reg_13 = _mm_mullo_epi32(m_temp_reg_1, m_coeff3); } /* c[0] in m_temp_reg_10 */ { m_temp_reg_10 = _mm_add_epi32(m_temp_reg_0, m_temp_reg_2); } /* c[1] in m_temp_reg_11 */ { m_temp_reg_11 = _mm_add_epi32(m_temp_reg_2, m_temp_reg_3); } /* c[2] in m_temp_reg_12 */ { m_temp_reg_12 = _mm_sub_epi32(m_temp_reg_0, m_temp_reg_3); } /* c[4] in m_temp_reg_14 */ /* c[4] = src[0] - src[2] + src[3] */ { m_temp_reg_14 = _mm_add_epi32(m_temp_reg_14, m_temp_reg_3); } /* Stage 1 outputs stored in m_temp_reg_20-23 */ { m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_10, 5); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_10, 1); m_temp_reg_23 = _mm_sub_epi32(m_temp_reg_20, m_temp_reg_10); m_temp_reg_30 = _mm_sub_epi32(m_temp_reg_23, m_temp_reg_21); //m_temp_reg_30 = _mm_mullo_epi32(m_temp_reg_10, m_coeff1);//29*c0 m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_11, 6); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_11, 3); m_temp_reg_23 = _mm_sub_epi32(m_temp_reg_20, m_temp_reg_11); m_temp_reg_31 = _mm_sub_epi32(m_temp_reg_23, m_temp_reg_21); //m_temp_reg_31 = _mm_mullo_epi32(m_temp_reg_11, m_coeff2);//55*c1 m_rdng_factor = _mm_cvtsi32_si128((1 << (i4_shift - 1))); m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_11, 5); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_11, 1); m_temp_reg_23 = _mm_sub_epi32(m_temp_reg_20, m_temp_reg_11); m_temp_reg_32 = _mm_sub_epi32(m_temp_reg_23, m_temp_reg_21); //m_temp_reg_32 = _mm_mullo_epi32(m_temp_reg_11, m_coeff1);//29*c1 m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_12, 6); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_12, 3); m_temp_reg_23 = _mm_sub_epi32(m_temp_reg_20, m_temp_reg_12); m_temp_reg_33 = _mm_sub_epi32(m_temp_reg_23, m_temp_reg_21); //m_temp_reg_33 = _mm_mullo_epi32(m_temp_reg_12, m_coeff2);//55*c2 m_rdng_factor = _mm_unpacklo_epi32(m_rdng_factor, m_rdng_factor); m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_10, 6); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_10, 3); m_temp_reg_23 = _mm_sub_epi32(m_temp_reg_20, m_temp_reg_10); m_temp_reg_34 = _mm_sub_epi32(m_temp_reg_23, m_temp_reg_21); //m_temp_reg_34 = _mm_mullo_epi32(m_temp_reg_10, m_coeff2);//55*c0 m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_12, 5); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_12, 1); m_temp_reg_23 = _mm_sub_epi32(m_temp_reg_20, m_temp_reg_12); m_temp_reg_35 = _mm_sub_epi32(m_temp_reg_23, m_temp_reg_21); //m_temp_reg_35 = _mm_mullo_epi32(m_temp_reg_12, m_coeff1);//29*c2 m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_14, 6); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_14, 3); m_temp_reg_22 = _mm_slli_epi32(m_temp_reg_14, 1); m_temp_reg_23 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_21); m_temp_reg_36 = _mm_add_epi32(m_temp_reg_22, m_temp_reg_23); //m_temp_reg_36 = _mm_mullo_epi32(m_temp_reg_14, m_coeff3);//74*c4 m_rdng_factor = _mm_unpacklo_epi64(m_rdng_factor, m_rdng_factor); m_count = _mm_cvtsi32_si128(i4_shift); m_temp_reg_20 = _mm_add_epi32(m_temp_reg_30, m_temp_reg_31); m_temp_reg_4 = _mm_add_epi32(m_rdng_factor, m_temp_reg_13); m_temp_reg_20 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_4); m_temp_reg_21 = _mm_sub_epi32(m_temp_reg_33, m_temp_reg_32); m_temp_reg_21 = _mm_add_epi32(m_temp_reg_21, m_temp_reg_4); m_temp_reg_23 = _mm_add_epi32(m_temp_reg_34, m_temp_reg_35); m_temp_reg_4 = _mm_sub_epi32(m_rdng_factor, m_temp_reg_13); m_temp_reg_23 = _mm_add_epi32(m_temp_reg_23, m_temp_reg_4); m_temp_reg_22 = _mm_add_epi32(m_temp_reg_36, m_rdng_factor); m_temp_reg_20 = _mm_sra_epi32(m_temp_reg_20, m_count); m_temp_reg_21 = _mm_sra_epi32(m_temp_reg_21, m_count); m_temp_reg_23 = _mm_sra_epi32(m_temp_reg_23, m_count); m_temp_reg_22 = _mm_sra_epi32(m_temp_reg_22, m_count); m_temp_reg_20 = _mm_packs_epi32(m_temp_reg_20, m_temp_reg_21); m_temp_reg_21 = _mm_packs_epi32(m_temp_reg_22, m_temp_reg_23); m_temp_reg_22 = _mm_srli_si128(m_temp_reg_20, 8); m_temp_reg_23 = _mm_srli_si128(m_temp_reg_21, 8); m_temp_reg_24 = _mm_unpacklo_epi16(m_temp_reg_20, m_temp_reg_22); m_temp_reg_25 = _mm_unpacklo_epi16(m_temp_reg_21, m_temp_reg_23); m_temp_reg_30 = _mm_unpacklo_epi32(m_temp_reg_24, m_temp_reg_25); m_temp_reg_31 = _mm_unpackhi_epi32(m_temp_reg_24, m_temp_reg_25); } /* Stage 2 */ { i4_shift = IT_SHIFT_STAGE_2; /*m_temp_reg_22 = _mm_srli_si128(m_temp_reg_20, 8); m_temp_reg_20 = _mm_cvtepi16_epi32(m_temp_reg_20); m_temp_reg_23 = _mm_srli_si128(m_temp_reg_21, 8); m_temp_reg_21 = _mm_cvtepi16_epi32(m_temp_reg_21); m_temp_reg_22 = _mm_cvtepi16_epi32(m_temp_reg_22); m_temp_reg_23 = _mm_cvtepi16_epi32(m_temp_reg_23);*/ m_ge_zero16b_flag_row0 = _mm_cmpgt_epi16(m_zero, m_temp_reg_30); m_ge_zero16b_flag_row1 = _mm_cmpgt_epi16(m_zero, m_temp_reg_31); m_temp_reg_20 = _mm_unpacklo_epi16(m_temp_reg_30, m_ge_zero16b_flag_row0); m_temp_reg_21 = _mm_unpacklo_epi16(m_temp_reg_31, m_ge_zero16b_flag_row1); m_temp_reg_22 = _mm_unpackhi_epi16(m_temp_reg_30, m_ge_zero16b_flag_row0); m_temp_reg_23 = _mm_unpackhi_epi16(m_temp_reg_31, m_ge_zero16b_flag_row1); /* c[4] stored in m_temp_reg_4 */ { m_temp_reg_4 = _mm_sub_epi32(m_temp_reg_20, m_temp_reg_21); } /* c[3] stored in m_temp_reg_3 */ { m_temp_reg_10 = _mm_slli_epi32(m_temp_reg_22, 6); m_temp_reg_11 = _mm_slli_epi32(m_temp_reg_22, 3); m_temp_reg_12 = _mm_slli_epi32(m_temp_reg_22, 1); m_temp_reg_13 = _mm_add_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_12, m_temp_reg_13); //m_temp_reg_3 = _mm_mullo_epi32(m_temp_reg_22, m_coeff3); } /* c[0] stored in m_temp_reg_0 */ { m_temp_reg_0 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_21); } /* c[1] stored in m_temp_reg_1 */ { m_temp_reg_1 = _mm_add_epi32(m_temp_reg_23, m_temp_reg_21); } /* c[2] stored in m_temp_reg_2 */ { m_temp_reg_2 = _mm_sub_epi32(m_temp_reg_20, m_temp_reg_23); } /* c[4] stored in m_temp_reg_4 */ { m_temp_reg_4 = _mm_add_epi32(m_temp_reg_4, m_temp_reg_23); } /* Stage 2 output generation */ { m_temp_reg_10 = _mm_slli_epi32(m_temp_reg_0, 5); m_temp_reg_11 = _mm_slli_epi32(m_temp_reg_0, 1); m_temp_reg_13 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_0); m_temp_reg_30 = _mm_sub_epi32(m_temp_reg_13, m_temp_reg_11); //m_temp_reg_30 = _mm_mullo_epi32(m_temp_reg_0, m_coeff1);//29*c0 m_temp_reg_10 = _mm_slli_epi32(m_temp_reg_1, 6); m_temp_reg_11 = _mm_slli_epi32(m_temp_reg_1, 3); m_temp_reg_13 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_1); m_temp_reg_31 = _mm_sub_epi32(m_temp_reg_13, m_temp_reg_11); //m_temp_reg_31 = _mm_mullo_epi32(m_temp_reg_1, m_coeff2);//55*c1 m_rdng_factor = _mm_cvtsi32_si128((1 << (i4_shift - 1))); m_temp_reg_10 = _mm_slli_epi32(m_temp_reg_1, 5); m_temp_reg_11 = _mm_slli_epi32(m_temp_reg_1, 1); m_temp_reg_13 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_1); m_temp_reg_32 = _mm_sub_epi32(m_temp_reg_13, m_temp_reg_11); //m_temp_reg_32 = _mm_mullo_epi32(m_temp_reg_1, m_coeff1);//29*c1 m_temp_reg_10 = _mm_slli_epi32(m_temp_reg_2, 6); m_temp_reg_11 = _mm_slli_epi32(m_temp_reg_2, 3); m_temp_reg_13 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_2); m_temp_reg_33 = _mm_sub_epi32(m_temp_reg_13, m_temp_reg_11); //m_temp_reg_33 = _mm_mullo_epi32(m_temp_reg_2, m_coeff2);//55*c2 m_rdng_factor = _mm_unpacklo_epi32(m_rdng_factor, m_rdng_factor); m_temp_reg_10 = _mm_slli_epi32(m_temp_reg_0, 6); m_temp_reg_11 = _mm_slli_epi32(m_temp_reg_0, 3); m_temp_reg_13 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_0); m_temp_reg_34 = _mm_sub_epi32(m_temp_reg_13, m_temp_reg_11); //m_temp_reg_34 = _mm_mullo_epi32(m_temp_reg_0, m_coeff2);//55*c0 m_temp_reg_10 = _mm_slli_epi32(m_temp_reg_2, 5); m_temp_reg_11 = _mm_slli_epi32(m_temp_reg_2, 1); m_temp_reg_13 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_2); m_temp_reg_35 = _mm_sub_epi32(m_temp_reg_13, m_temp_reg_11); //m_temp_reg_35 = _mm_mullo_epi32(m_temp_reg_2, m_coeff1);//29*c2 m_temp_reg_10 = _mm_slli_epi32(m_temp_reg_4, 6); m_temp_reg_11 = _mm_slli_epi32(m_temp_reg_4, 3); m_temp_reg_12 = _mm_slli_epi32(m_temp_reg_4, 1); m_temp_reg_13 = _mm_add_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_36 = _mm_add_epi32(m_temp_reg_12, m_temp_reg_13); //m_temp_reg_36 = _mm_mullo_epi32(m_temp_reg_4, m_coeff3);//74*c4 m_rdng_factor = _mm_unpacklo_epi64(m_rdng_factor, m_rdng_factor); m_count = _mm_cvtsi32_si128(i4_shift); m_temp_reg_4 = _mm_add_epi32(m_rdng_factor, m_temp_reg_3); m_temp_reg_20 = _mm_add_epi32(m_temp_reg_30, m_temp_reg_31); m_temp_reg_20 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_4); m_temp_reg_21 = _mm_sub_epi32(m_temp_reg_33, m_temp_reg_32); m_temp_reg_21 = _mm_add_epi32(m_temp_reg_21, m_temp_reg_4); m_temp_reg_4 = _mm_sub_epi32(m_rdng_factor, m_temp_reg_3); m_temp_reg_23 = _mm_add_epi32(m_temp_reg_34, m_temp_reg_35); m_temp_reg_23 = _mm_add_epi32(m_temp_reg_23, m_temp_reg_4); m_temp_reg_22 = _mm_add_epi32(m_temp_reg_36, m_rdng_factor); m_temp_reg_20 = _mm_sra_epi32(m_temp_reg_20, m_count); m_temp_reg_21 = _mm_sra_epi32(m_temp_reg_21, m_count); m_temp_reg_23 = _mm_sra_epi32(m_temp_reg_23, m_count); m_temp_reg_22 = _mm_sra_epi32(m_temp_reg_22, m_count); m_temp_reg_20 = _mm_packs_epi32(m_temp_reg_20, m_temp_reg_21); m_temp_reg_21 = _mm_packs_epi32(m_temp_reg_22, m_temp_reg_23); m_temp_reg_22 = _mm_srli_si128(m_temp_reg_20, 8); m_temp_reg_23 = _mm_srli_si128(m_temp_reg_21, 8); m_temp_reg_24 = _mm_unpacklo_epi16(m_temp_reg_20, m_temp_reg_22); m_temp_reg_25 = _mm_unpacklo_epi16(m_temp_reg_21, m_temp_reg_23); m_temp_reg_20 = _mm_unpacklo_epi32(m_temp_reg_24, m_temp_reg_25); m_temp_reg_21 = _mm_unpackhi_epi32(m_temp_reg_24, m_temp_reg_25); } /* Recon and store */ { WORD32 *pi4_dst = (WORD32 *)pu1_dst; m_temp_reg_0 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_1 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_2 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_3 = _mm_loadl_epi64((__m128i *)pu1_pred); m_temp_reg_0 = _mm_unpacklo_epi8(m_temp_reg_0, m_zero); m_temp_reg_1 = _mm_unpacklo_epi8(m_temp_reg_1, m_zero); m_temp_reg_2 = _mm_unpacklo_epi8(m_temp_reg_2, m_zero); m_temp_reg_3 = _mm_unpacklo_epi8(m_temp_reg_3, m_zero); /*m_temp_reg_0 = _mm_cvtepu8_epi16(m_temp_reg_0); m_temp_reg_1 = _mm_cvtepu8_epi16(m_temp_reg_1); m_temp_reg_2 = _mm_cvtepu8_epi16(m_temp_reg_2); m_temp_reg_3 = _mm_cvtepu8_epi16(m_temp_reg_3);*/ m_temp_reg_0 = _mm_unpacklo_epi64(m_temp_reg_0, m_temp_reg_1); m_temp_reg_1 = _mm_unpacklo_epi64(m_temp_reg_2, m_temp_reg_3); m_temp_reg_20 = _mm_add_epi16(m_temp_reg_20, m_temp_reg_0); m_temp_reg_21 = _mm_add_epi16(m_temp_reg_21, m_temp_reg_1); m_temp_reg_0 = _mm_packus_epi16(m_temp_reg_20, m_temp_reg_21); *pi4_dst = _mm_cvtsi128_si32(m_temp_reg_0); m_temp_reg_1 = _mm_srli_si128(m_temp_reg_0, 4); m_temp_reg_2 = _mm_srli_si128(m_temp_reg_0, 8); m_temp_reg_3 = _mm_srli_si128(m_temp_reg_0, 12); pu1_dst += dst_strd; pi4_dst = (WORD32 *)(pu1_dst); *pi4_dst = _mm_cvtsi128_si32(m_temp_reg_1); pu1_dst += dst_strd; pi4_dst = (WORD32 *)(pu1_dst); *pi4_dst = _mm_cvtsi128_si32(m_temp_reg_2); pu1_dst += dst_strd; pi4_dst = (WORD32 *)(pu1_dst); *pi4_dst = _mm_cvtsi128_si32(m_temp_reg_3); } } } /** ******************************************************************************* * * @brief * This function performs inverse quantization, inverse transform * (DCT) and reconstruction for 4x4 input block * * @par Description: * Performs inverse quantization , inverse transform 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 * * @returns Void * * @remarks * None * ******************************************************************************* */ void ihevc_itrans_recon_4x4_ssse3(WORD16 *pi2_src, WORD16 *pi2_tmp, UWORD8 *pu1_pred, UWORD8 *pu1_dst, WORD32 src_strd, WORD32 pred_strd, WORD32 dst_strd, WORD32 zero_cols, WORD32 zero_rows) { __m128i m_temp_reg_0; __m128i m_temp_reg_1; __m128i m_temp_reg_2; __m128i m_temp_reg_3; __m128i m_temp_reg_4; __m128i m_temp_reg_10; __m128i m_temp_reg_11; __m128i m_temp_reg_12; __m128i m_temp_reg_13; __m128i m_temp_reg_14; __m128i m_temp_reg_15; __m128i m_temp_reg_20; __m128i m_temp_reg_21; __m128i m_temp_reg_22; __m128i m_temp_reg_23; __m128i m_temp_reg_24; __m128i m_temp_reg_25; __m128i m_temp_reg_30; __m128i m_temp_reg_31; __m128i m_temp_reg_33; __m128i m_temp_reg_34; __m128i m_rdng_factor; __m128i m_count; __m128i m_ge_zero16b_flag_row0; __m128i m_ge_zero16b_flag_row1; __m128i m_ge_zero16b_flag_row2; __m128i m_ge_zero16b_flag_row3; __m128i m_zero = _mm_setzero_si128(); WORD32 i4_shift = IT_SHIFT_STAGE_1; UNUSED(zero_rows); UNUSED(zero_cols); UNUSED(pi2_tmp); m_temp_reg_0 = _mm_loadl_epi64((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_1 = _mm_loadl_epi64((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_2 = _mm_loadl_epi64((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_3 = _mm_loadl_epi64((__m128i *)pi2_src); m_ge_zero16b_flag_row0 = _mm_cmpgt_epi16(m_zero, m_temp_reg_0); m_ge_zero16b_flag_row1 = _mm_cmpgt_epi16(m_zero, m_temp_reg_1); m_ge_zero16b_flag_row2 = _mm_cmpgt_epi16(m_zero, m_temp_reg_2); m_ge_zero16b_flag_row3 = _mm_cmpgt_epi16(m_zero, m_temp_reg_3); m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_0, m_ge_zero16b_flag_row0); m_temp_reg_1 = _mm_unpacklo_epi16(m_temp_reg_1, m_ge_zero16b_flag_row1); m_temp_reg_2 = _mm_unpacklo_epi16(m_temp_reg_2, m_ge_zero16b_flag_row2); m_temp_reg_3 = _mm_unpacklo_epi16(m_temp_reg_3, m_ge_zero16b_flag_row3); /*m_temp_reg_0 = _mm_cvtepi16_epi32(m_temp_reg_0); m_temp_reg_2 = _mm_cvtepi16_epi32(m_temp_reg_2); m_temp_reg_1 = _mm_cvtepi16_epi32(m_temp_reg_1); m_temp_reg_3 = _mm_cvtepi16_epi32(m_temp_reg_3);*/ /* e */ { m_temp_reg_10 = _mm_slli_epi32(m_temp_reg_0, 6); m_temp_reg_11 = _mm_slli_epi32(m_temp_reg_2, 6); } /* o */ { m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_1, 5); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_1, 2); m_temp_reg_12 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_21); //m_temp_reg_12 = _mm_mullo_epi32(m_temp_reg_1, m_coeff1);//src[1]*36 m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_3, 6); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_3, 4); m_temp_reg_22 = _mm_slli_epi32(m_temp_reg_3, 1); m_temp_reg_23 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_3); m_temp_reg_24 = _mm_add_epi32(m_temp_reg_21, m_temp_reg_22); m_temp_reg_13 = _mm_add_epi32(m_temp_reg_23, m_temp_reg_24); //m_temp_reg_13 = _mm_mullo_epi32(m_temp_reg_3, m_coeff3);//src[3]*83 m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_1, 6); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_1, 4); m_temp_reg_22 = _mm_slli_epi32(m_temp_reg_1, 1); m_temp_reg_23 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_1); m_temp_reg_24 = _mm_add_epi32(m_temp_reg_21, m_temp_reg_22); m_temp_reg_14 = _mm_add_epi32(m_temp_reg_23, m_temp_reg_24); //m_temp_reg_14 = _mm_mullo_epi32(m_temp_reg_1, m_coeff3);//src[1]*83 m_temp_reg_20 = _mm_slli_epi32(m_temp_reg_3, 5); m_temp_reg_21 = _mm_slli_epi32(m_temp_reg_3, 2); m_temp_reg_15 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_21); //m_temp_reg_15 = _mm_mullo_epi32(m_temp_reg_3, m_coeff1);//src[3]*36 } m_rdng_factor = _mm_cvtsi32_si128((1 << (i4_shift - 1))); /* e1 stored in m_temp_reg_31 */ { m_temp_reg_31 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_11); } m_rdng_factor = _mm_unpacklo_epi32(m_rdng_factor, m_rdng_factor); /* e0 stored in m_temp_reg_30 */ { m_temp_reg_30 = _mm_add_epi32(m_temp_reg_10, m_temp_reg_11); } m_rdng_factor = _mm_unpacklo_epi64(m_rdng_factor, m_rdng_factor); m_count = _mm_cvtsi32_si128(i4_shift); /* o1 stored in m_temp_reg_33 */ { m_temp_reg_33 = _mm_sub_epi32(m_temp_reg_12, m_temp_reg_13); } /* e1 + add */ { m_temp_reg_31 = _mm_add_epi32(m_temp_reg_31, m_rdng_factor); } /* e0 + add */ { m_temp_reg_30 = _mm_add_epi32(m_temp_reg_30, m_rdng_factor); } /* o0 stored in m_temp_reg_34 */ { m_temp_reg_34 = _mm_add_epi32(m_temp_reg_14, m_temp_reg_15); } /* Stage 1 outputs */ { m_temp_reg_21 = _mm_add_epi32(m_temp_reg_31, m_temp_reg_33); m_temp_reg_22 = _mm_sub_epi32(m_temp_reg_31, m_temp_reg_33); m_temp_reg_20 = _mm_add_epi32(m_temp_reg_30, m_temp_reg_34); m_temp_reg_23 = _mm_sub_epi32(m_temp_reg_30, m_temp_reg_34); m_temp_reg_21 = _mm_sra_epi32(m_temp_reg_21, m_count); m_temp_reg_20 = _mm_sra_epi32(m_temp_reg_20, m_count); m_temp_reg_22 = _mm_sra_epi32(m_temp_reg_22, m_count); m_temp_reg_23 = _mm_sra_epi32(m_temp_reg_23, m_count); m_temp_reg_20 = _mm_packs_epi32(m_temp_reg_20, m_temp_reg_21); m_temp_reg_21 = _mm_packs_epi32(m_temp_reg_22, m_temp_reg_23); m_temp_reg_22 = _mm_srli_si128(m_temp_reg_20, 8); m_temp_reg_23 = _mm_srli_si128(m_temp_reg_21, 8); m_temp_reg_24 = _mm_unpacklo_epi16(m_temp_reg_20, m_temp_reg_22); m_temp_reg_25 = _mm_unpacklo_epi16(m_temp_reg_21, m_temp_reg_23); m_temp_reg_30 = _mm_unpacklo_epi32(m_temp_reg_24, m_temp_reg_25); m_temp_reg_31 = _mm_unpackhi_epi32(m_temp_reg_24, m_temp_reg_25); } /* Stage 2 */ { i4_shift = IT_SHIFT_STAGE_2; /*m_temp_reg_22 = _mm_srli_si128(m_temp_reg_20, 8); m_temp_reg_23 = _mm_srli_si128(m_temp_reg_21, 8);*/ m_ge_zero16b_flag_row0 = _mm_cmpgt_epi16(m_zero, m_temp_reg_30); m_ge_zero16b_flag_row1 = _mm_cmpgt_epi16(m_zero, m_temp_reg_31); m_temp_reg_20 = _mm_unpacklo_epi16(m_temp_reg_30, m_ge_zero16b_flag_row0); m_temp_reg_21 = _mm_unpacklo_epi16(m_temp_reg_31, m_ge_zero16b_flag_row1); m_temp_reg_22 = _mm_unpackhi_epi16(m_temp_reg_30, m_ge_zero16b_flag_row0); m_temp_reg_23 = _mm_unpackhi_epi16(m_temp_reg_31, m_ge_zero16b_flag_row1); /*m_temp_reg_20 = _mm_cvtepi16_epi32(m_temp_reg_20); m_temp_reg_21 = _mm_cvtepi16_epi32(m_temp_reg_21); m_temp_reg_22 = _mm_cvtepi16_epi32(m_temp_reg_22); m_temp_reg_23 = _mm_cvtepi16_epi32(m_temp_reg_23);*/ /* e */ { m_temp_reg_10 = _mm_slli_epi32(m_temp_reg_20, 6); } /* o */ /*{ m_temp_reg_12 = _mm_mullo_epi32(m_temp_reg_22, m_coeff1);//src[1]*36 m_temp_reg_14 = _mm_mullo_epi32(m_temp_reg_22, m_coeff3);//src[1]*83 m_temp_reg_13 = _mm_mullo_epi32(m_temp_reg_23, m_coeff3);//src[3]*83 m_temp_reg_15 = _mm_mullo_epi32(m_temp_reg_23, m_coeff1);//src[3]*36 }*/ { m_temp_reg_0 = _mm_slli_epi32(m_temp_reg_22, 5); m_temp_reg_1 = _mm_slli_epi32(m_temp_reg_22, 2); m_temp_reg_12 = _mm_add_epi32(m_temp_reg_0, m_temp_reg_1); //m_temp_reg_12 = _mm_mullo_epi32(m_temp_reg_1, m_coeff1);//src[1]*36 m_temp_reg_0 = _mm_slli_epi32(m_temp_reg_23, 6); m_temp_reg_1 = _mm_slli_epi32(m_temp_reg_23, 4); m_temp_reg_2 = _mm_slli_epi32(m_temp_reg_23, 1); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_0, m_temp_reg_23); m_temp_reg_4 = _mm_add_epi32(m_temp_reg_1, m_temp_reg_2); m_temp_reg_13 = _mm_add_epi32(m_temp_reg_3, m_temp_reg_4); //m_temp_reg_13 = _mm_mullo_epi32(m_temp_reg_3, m_coeff3);//src[3]*83 m_temp_reg_0 = _mm_slli_epi32(m_temp_reg_22, 6); m_temp_reg_1 = _mm_slli_epi32(m_temp_reg_22, 4); m_temp_reg_2 = _mm_slli_epi32(m_temp_reg_22, 1); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_0, m_temp_reg_22); m_temp_reg_4 = _mm_add_epi32(m_temp_reg_1, m_temp_reg_2); m_temp_reg_14 = _mm_add_epi32(m_temp_reg_3, m_temp_reg_4); //m_temp_reg_14 = _mm_mullo_epi32(m_temp_reg_1, m_coeff3);//src[1]*83 m_temp_reg_0 = _mm_slli_epi32(m_temp_reg_23, 5); m_temp_reg_1 = _mm_slli_epi32(m_temp_reg_23, 2); m_temp_reg_15 = _mm_add_epi32(m_temp_reg_0, m_temp_reg_1); //m_temp_reg_15 = _mm_mullo_epi32(m_temp_reg_3, m_coeff1);//src[3]*36 } /* e */ { m_temp_reg_11 = _mm_slli_epi32(m_temp_reg_21, 6); } m_rdng_factor = _mm_cvtsi32_si128((1 << (i4_shift - 1))); /* e1 stored in m_temp_reg_31 */ { m_temp_reg_31 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_11); } m_rdng_factor = _mm_unpacklo_epi32(m_rdng_factor, m_rdng_factor); /* e0 stored in m_temp_reg_30 */ { m_temp_reg_30 = _mm_add_epi32(m_temp_reg_10, m_temp_reg_11); } m_rdng_factor = _mm_unpacklo_epi64(m_rdng_factor, m_rdng_factor); m_count = _mm_cvtsi32_si128(i4_shift); /* o1 stored in m_temp_reg_33 */ { m_temp_reg_33 = _mm_sub_epi32(m_temp_reg_12, m_temp_reg_13); } /* e1 + add */ { m_temp_reg_31 = _mm_add_epi32(m_temp_reg_31, m_rdng_factor); } /* e0 + add */ { m_temp_reg_30 = _mm_add_epi32(m_temp_reg_30, m_rdng_factor); } /* o0 stored in m_temp_reg_34 */ { m_temp_reg_34 = _mm_add_epi32(m_temp_reg_14, m_temp_reg_15); } /* Stage 2 outputs */ { m_temp_reg_21 = _mm_add_epi32(m_temp_reg_31, m_temp_reg_33); m_temp_reg_22 = _mm_sub_epi32(m_temp_reg_31, m_temp_reg_33); m_temp_reg_20 = _mm_add_epi32(m_temp_reg_30, m_temp_reg_34); m_temp_reg_23 = _mm_sub_epi32(m_temp_reg_30, m_temp_reg_34); m_temp_reg_21 = _mm_sra_epi32(m_temp_reg_21, m_count); m_temp_reg_22 = _mm_sra_epi32(m_temp_reg_22, m_count); m_temp_reg_20 = _mm_sra_epi32(m_temp_reg_20, m_count); m_temp_reg_23 = _mm_sra_epi32(m_temp_reg_23, m_count); m_temp_reg_20 = _mm_packs_epi32(m_temp_reg_20, m_temp_reg_21); m_temp_reg_21 = _mm_packs_epi32(m_temp_reg_22, m_temp_reg_23); m_temp_reg_22 = _mm_srli_si128(m_temp_reg_20, 8); m_temp_reg_23 = _mm_srli_si128(m_temp_reg_21, 8); m_temp_reg_24 = _mm_unpacklo_epi16(m_temp_reg_20, m_temp_reg_22); m_temp_reg_25 = _mm_unpacklo_epi16(m_temp_reg_21, m_temp_reg_23); m_temp_reg_20 = _mm_unpacklo_epi32(m_temp_reg_24, m_temp_reg_25); m_temp_reg_21 = _mm_unpackhi_epi32(m_temp_reg_24, m_temp_reg_25); } /* Recon and store */ { UWORD32 *pu4_dst = (UWORD32 *)pu1_dst; m_temp_reg_0 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_1 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_2 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_3 = _mm_loadl_epi64((__m128i *)pu1_pred); //m_temp_reg_0 = _mm_cvtepu8_epi16(m_temp_reg_0); //m_temp_reg_1 = _mm_cvtepu8_epi16(m_temp_reg_1); //m_temp_reg_2 = _mm_cvtepu8_epi16(m_temp_reg_2); //m_temp_reg_3 = _mm_cvtepu8_epi16(m_temp_reg_3); m_temp_reg_0 = _mm_unpacklo_epi8(m_temp_reg_0, m_zero); m_temp_reg_1 = _mm_unpacklo_epi8(m_temp_reg_1, m_zero); m_temp_reg_2 = _mm_unpacklo_epi8(m_temp_reg_2, m_zero); m_temp_reg_3 = _mm_unpacklo_epi8(m_temp_reg_3, m_zero); m_temp_reg_0 = _mm_unpacklo_epi64(m_temp_reg_0, m_temp_reg_1); m_temp_reg_1 = _mm_unpacklo_epi64(m_temp_reg_2, m_temp_reg_3); m_temp_reg_20 = _mm_add_epi16(m_temp_reg_20, m_temp_reg_0); m_temp_reg_21 = _mm_add_epi16(m_temp_reg_21, m_temp_reg_1); m_temp_reg_0 = _mm_packus_epi16(m_temp_reg_20, m_temp_reg_21); *pu4_dst = _mm_cvtsi128_si32(m_temp_reg_0); m_temp_reg_1 = _mm_srli_si128(m_temp_reg_0, 4); m_temp_reg_2 = _mm_srli_si128(m_temp_reg_0, 8); m_temp_reg_3 = _mm_srli_si128(m_temp_reg_0, 12); pu1_dst += dst_strd; pu4_dst = (UWORD32 *)(pu1_dst); *pu4_dst = _mm_cvtsi128_si32(m_temp_reg_1); pu1_dst += dst_strd; pu4_dst = (UWORD32 *)(pu1_dst); *pu4_dst = _mm_cvtsi128_si32(m_temp_reg_2); pu1_dst += dst_strd; pu4_dst = (UWORD32 *)(pu1_dst); *pu4_dst = _mm_cvtsi128_si32(m_temp_reg_3); } } } /** ******************************************************************************* * * @brief * This function performs inverse quantization, inverse transform and * reconstruction for 8c8 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 8x8 coefficients * * @param[in] pi2_tmp * Temporary 8x8 buffer for storing inverse * transform 1st stage output * * @param[in] pu1_pred * Prediction 8x8 block * * @param[in] pi2_dequant_coeff * Dequant Coeffs * * @param[out] pu1_dst * Output 8x8 block * * @param[in] src_strd * Input stride * * @param[in] qp_div * Quantization parameter / 6 * * @param[in] qp_rem * Quantization parameter % 6 * * @param[in] pred_strd * Prediction stride * * @param[in] dst_strd * Output Stride * * @param[in] zero_cols * Zero columns in pi2_src * * @returns Void * * @remarks * None * ******************************************************************************* */ void ihevc_itrans_recon_8x8_ssse3(WORD16 *pi2_src, WORD16 *pi2_tmp, UWORD8 *pu1_pred, UWORD8 *pu1_dst, WORD32 src_strd, WORD32 pred_strd, WORD32 dst_strd, WORD32 zero_cols, WORD32 zero_rows) { __m128i m_temp_reg_0; __m128i m_temp_reg_1; __m128i m_temp_reg_2; __m128i m_temp_reg_3; __m128i m_temp_reg_5; __m128i m_temp_reg_6; __m128i m_temp_reg_7; __m128i m_temp_reg_4; __m128i m_temp_reg_10; __m128i m_temp_reg_11; __m128i m_temp_reg_12; __m128i m_temp_reg_13; __m128i m_temp_reg_14; __m128i m_temp_reg_15; __m128i m_temp_reg_16; __m128i m_temp_reg_17; __m128i m_temp_reg_20; __m128i m_temp_reg_21; __m128i m_temp_reg_22; __m128i m_temp_reg_23; __m128i m_temp_reg_24; __m128i m_temp_reg_25; __m128i m_temp_reg_26; __m128i m_temp_reg_27; __m128i m_temp_reg_30; __m128i m_temp_reg_31; __m128i m_temp_reg_32; __m128i m_temp_reg_33; __m128i m_temp_reg_34; __m128i m_temp_reg_35; __m128i m_temp_reg_36; __m128i m_temp_reg_37; __m128i m_temp_reg_40; __m128i m_temp_reg_41; __m128i m_temp_reg_42; __m128i m_temp_reg_43; __m128i m_temp_reg_44; __m128i m_temp_reg_45; __m128i m_temp_reg_46; __m128i m_temp_reg_47; __m128i m_temp_reg_50; __m128i m_temp_reg_51; __m128i m_temp_reg_52; __m128i m_temp_reg_53; __m128i m_temp_reg_54; __m128i m_temp_reg_55; __m128i m_temp_reg_56; __m128i m_temp_reg_57; __m128i m_temp_reg_60; __m128i m_temp_reg_61; __m128i m_temp_reg_62; __m128i m_temp_reg_63; __m128i m_temp_reg_64; __m128i m_temp_reg_65; __m128i m_temp_reg_66; __m128i m_temp_reg_67; __m128i m_temp_reg_70; __m128i m_temp_reg_71; __m128i m_temp_reg_72; __m128i m_temp_reg_73; __m128i m_temp_reg_74; __m128i m_temp_reg_75; __m128i m_temp_reg_76; __m128i m_temp_reg_77; __m128i m_coeff1, m_coeff2, m_coeff3, m_coeff4; WORD32 check_row_stage_1; /* Lokesh */ WORD32 check_row_stage_2; /* Lokesh */ __m128i m_rdng_factor; //__m128i m_count; WORD32 i4_shift = IT_SHIFT_STAGE_1; UNUSED(zero_rows); UNUSED(zero_cols); UNUSED(pi2_tmp); check_row_stage_1 = ((zero_rows & 0xF0) != 0xF0) ? 1 : 0; check_row_stage_2 = ((zero_cols & 0xF0) != 0xF0) ? 1 : 0; m_temp_reg_70 = _mm_load_si128((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_71 = _mm_load_si128((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_72 = _mm_load_si128((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_73 = _mm_load_si128((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_74 = _mm_load_si128((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_75 = _mm_load_si128((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_76 = _mm_load_si128((__m128i *)pi2_src); pi2_src += src_strd; m_temp_reg_77 = _mm_load_si128((__m128i *)pi2_src); if(!check_row_stage_2) { if(!check_row_stage_1) { /* ee0 is present in the registers m_temp_reg_10 and m_temp_reg_11 */ /* ee1 is present in the registers m_temp_reg_12 and m_temp_reg_13 */ { //Interleaving 0,4 row in 0 , 1 Rishab /*coef2 for m_temp_reg_12 and m_temp_reg_13 , coef1 for m_temp_reg_10 and m_temp_reg_11*/ m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[3][0]); m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[0][0]); m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_70, m_temp_reg_74); m_temp_reg_10 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_12 = _mm_madd_epi16(m_temp_reg_0, m_coeff2); } /* eo0 is present in the registers m_temp_reg_14 and m_temp_reg_15 */ /* eo1 is present in the registers m_temp_reg_16 and m_temp_reg_17 */ /* as upper 8 bytes are zeros so m_temp_reg_15 and m_temp_reg_17 are not used*/ { m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[1][0]); //sub 2B*36-6B*83 ,2T*36-6T*83 m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[2][0]); //add 2B*83+6B*36 ,2T*83+6T*36 /* Combining instructions to eliminate them based on zero_rows : Lokesh */ //Interleaving 2,6 row in 4, 5 Rishab m_temp_reg_4 = _mm_unpacklo_epi16(m_temp_reg_72, m_temp_reg_76); m_temp_reg_16 = _mm_madd_epi16(m_temp_reg_4, m_coeff1); m_temp_reg_14 = _mm_madd_epi16(m_temp_reg_4, m_coeff2); /* Loading coeff for computing o0, o1, o2 and o3 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[2][0]); m_coeff4 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[3][0]); m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[0][0]); m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[1][0]); /* e */ /* e0 stored in m_temp_reg_40 and m_temp_reg_41 */ /* e1 stored in m_temp_reg_42 and m_temp_reg_43 */ /* e3 stored in m_temp_reg_46 and m_temp_reg_47 */ /* e2 stored in m_temp_reg_44 and m_temp_reg_45 */ m_temp_reg_42 = _mm_add_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_44 = _mm_sub_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_40 = _mm_add_epi32(m_temp_reg_10, m_temp_reg_14); m_temp_reg_46 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_14); } /* o */ { /* o0 stored in m_temp_reg_30 and m_temp_reg_31 */ { m_temp_reg_60 = _mm_unpacklo_epi16(m_temp_reg_71, m_temp_reg_73); //o0:1B*89+3B*75,5B*50+7B*18 m_temp_reg_30 = _mm_madd_epi16(m_temp_reg_60, m_coeff1); m_rdng_factor = _mm_cvtsi32_si128((1 << (i4_shift - 1))); m_rdng_factor = _mm_shuffle_epi32(m_rdng_factor, 0x0000); /* Column 0 of destination computed here */ /* It is stored in m_temp_reg_50 */ /* Column 7 of destination computed here */ /* It is stored in m_temp_reg_57 */ /* Upper 8 bytes of both registers are zero due to zero_cols*/ m_temp_reg_62 = _mm_add_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_63 = _mm_setzero_si128(); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); //o1:1B*75-3B*18,5B*89+7B*50 m_temp_reg_32 = _mm_madd_epi16(m_temp_reg_60, m_coeff3); m_temp_reg_50 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_57 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_63); /* Loading coeff for computing o2 in the next block */ m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[4][0]); m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[5][0]); /* o1 stored in m_temp_reg_32 and m_temp_reg_33 */ /* Column 1 of destination computed here */ /* It is stored in m_temp_reg_51 */ /* Column 6 of destination computed here */ /* It is stored in m_temp_reg_56 */ m_temp_reg_62 = _mm_add_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); //o2:1B*50-3B*89,5B*18+7B*75 m_temp_reg_34 = _mm_madd_epi16(m_temp_reg_60, m_coeff1); m_temp_reg_51 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_56 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_63); /* o2 stored in m_temp_reg_34 and m_temp_reg_35 */ /* Loading coeff for computing o3 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[6][0]); m_coeff4 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[7][0]); /* Column 2 of destination computed here */ /* It is stored in m_temp_reg_52 */ /* Column 5 of destination computed here */ /* It is stored in m_temp_reg_55 */ m_temp_reg_62 = _mm_add_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); //o3:1B*18-3B*50,5B*75-7B*89 m_temp_reg_36 = _mm_madd_epi16(m_temp_reg_60, m_coeff3); m_temp_reg_52 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_55 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_63); /* o3 stored in m_temp_reg_36 and m_temp_reg_37 */ /* Column 3 of destination computed here */ /* It is stored in m_temp_reg_53 */ /* Column 4 of destination computed here */ /* It is stored in m_temp_reg_54 */ m_temp_reg_62 = _mm_add_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); m_temp_reg_53 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_54 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_63); } } /* Transpose of the destination 8x8 matrix done here */ /* and ultimately stored in registers m_temp_reg_50 to m_temp_reg_57 */ /* respectively */ { m_temp_reg_10 = _mm_unpacklo_epi16(m_temp_reg_50, m_temp_reg_51); m_temp_reg_11 = _mm_unpacklo_epi16(m_temp_reg_52, m_temp_reg_53); //m_temp_reg_14 = _mm_unpackhi_epi16(m_temp_reg_50, m_temp_reg_51); //m_temp_reg_15 = _mm_unpackhi_epi16(m_temp_reg_52, m_temp_reg_53); m_temp_reg_0 = _mm_unpacklo_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_1 = _mm_unpackhi_epi32(m_temp_reg_10, m_temp_reg_11); //m_temp_reg_2 = _mm_unpacklo_epi32(m_temp_reg_14, m_temp_reg_15); //m_temp_reg_3 = _mm_unpackhi_epi32(m_temp_reg_14, m_temp_reg_15); m_temp_reg_12 = _mm_unpacklo_epi16(m_temp_reg_54, m_temp_reg_55); m_temp_reg_13 = _mm_unpacklo_epi16(m_temp_reg_56, m_temp_reg_57); //m_temp_reg_16 = _mm_unpackhi_epi16(m_temp_reg_54, m_temp_reg_55); //m_temp_reg_17 = _mm_unpackhi_epi16(m_temp_reg_56, m_temp_reg_57); m_temp_reg_4 = _mm_unpacklo_epi32(m_temp_reg_12, m_temp_reg_13); m_temp_reg_5 = _mm_unpackhi_epi32(m_temp_reg_12, m_temp_reg_13); //m_temp_reg_6 = _mm_unpacklo_epi32(m_temp_reg_16, m_temp_reg_17); //m_temp_reg_7 = _mm_unpackhi_epi32(m_temp_reg_16, m_temp_reg_17); m_temp_reg_50 = _mm_unpacklo_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_51 = _mm_unpackhi_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_52 = _mm_unpacklo_epi64(m_temp_reg_1, m_temp_reg_5); m_temp_reg_53 = _mm_unpackhi_epi64(m_temp_reg_1, m_temp_reg_5); /*m_temp_reg_54 = _mm_unpacklo_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_55 = _mm_unpackhi_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_56 = _mm_unpacklo_epi64(m_temp_reg_3, m_temp_reg_7); m_temp_reg_57 = _mm_unpackhi_epi64(m_temp_reg_3, m_temp_reg_7); */ m_temp_reg_54 = _mm_setzero_si128(); m_temp_reg_55 = _mm_setzero_si128(); m_temp_reg_56 = _mm_setzero_si128(); m_temp_reg_57 = _mm_setzero_si128(); } } else { /* ee0 is present in the registers m_temp_reg_10 and m_temp_reg_11 */ /* ee1 is present in the registers m_temp_reg_12 and m_temp_reg_13 */ { //Interleaving 0,4 row in 0 , 1 Rishab /*coef2 for m_temp_reg_12 and m_temp_reg_13 , coef1 for m_temp_reg_10 and m_temp_reg_11*/ m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[3][0]); m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[0][0]); m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_70, m_temp_reg_74); m_temp_reg_10 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_12 = _mm_madd_epi16(m_temp_reg_0, m_coeff2); } /* eo0 is present in the registers m_temp_reg_14 and m_temp_reg_15 */ /* eo1 is present in the registers m_temp_reg_16 and m_temp_reg_17 */ /* as upper 8 bytes are zeros so m_temp_reg_15 and m_temp_reg_17 are not used*/ { m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[1][0]); //sub 2B*36-6B*83 ,2T*36-6T*83 m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[2][0]); //add 2B*83+6B*36 ,2T*83+6T*36 /* Combining instructions to eliminate them based on zero_rows : Lokesh */ //Interleaving 2,6 row in 4, 5 Rishab m_temp_reg_4 = _mm_unpacklo_epi16(m_temp_reg_72, m_temp_reg_76); m_temp_reg_16 = _mm_madd_epi16(m_temp_reg_4, m_coeff1); m_temp_reg_14 = _mm_madd_epi16(m_temp_reg_4, m_coeff2); /* Loading coeff for computing o0, o1, o2 and o3 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[2][0]); m_coeff4 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[3][0]); m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[0][0]); m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[1][0]); /* e */ /* e0 stored in m_temp_reg_40 and m_temp_reg_41 */ /* e1 stored in m_temp_reg_42 and m_temp_reg_43 */ /* e3 stored in m_temp_reg_46 and m_temp_reg_47 */ /* e2 stored in m_temp_reg_44 and m_temp_reg_45 */ m_temp_reg_42 = _mm_add_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_44 = _mm_sub_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_40 = _mm_add_epi32(m_temp_reg_10, m_temp_reg_14); m_temp_reg_46 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_14); } /* o */ { /* o0 stored in m_temp_reg_30 and m_temp_reg_31 */ { m_temp_reg_60 = _mm_unpacklo_epi16(m_temp_reg_71, m_temp_reg_73); m_temp_reg_64 = _mm_unpacklo_epi16(m_temp_reg_75, m_temp_reg_77); //o0:1B*89+3B*75,5B*50+7B*18 m_temp_reg_20 = _mm_madd_epi16(m_temp_reg_60, m_coeff1); m_temp_reg_24 = _mm_madd_epi16(m_temp_reg_64, m_coeff2); m_rdng_factor = _mm_cvtsi32_si128((1 << (i4_shift - 1))); m_rdng_factor = _mm_shuffle_epi32(m_rdng_factor, 0x0000); m_temp_reg_30 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_24); /* Column 0 of destination computed here */ /* It is stored in m_temp_reg_50 */ /* Column 7 of destination computed here */ /* It is stored in m_temp_reg_57 */ /* Upper 8 bytes of both registers are zero due to zero_cols*/ m_temp_reg_62 = _mm_add_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_63 = _mm_setzero_si128(); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); //o1:1B*75-3B*18,5B*89+7B*50 m_temp_reg_22 = _mm_madd_epi16(m_temp_reg_60, m_coeff3); m_temp_reg_26 = _mm_madd_epi16(m_temp_reg_64, m_coeff4); m_temp_reg_50 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_57 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_63); /* Loading coeff for computing o2 in the next block */ m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[4][0]); m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[5][0]); /* o1 stored in m_temp_reg_32 and m_temp_reg_33 */ m_temp_reg_32 = _mm_sub_epi32(m_temp_reg_22, m_temp_reg_26); /* Column 1 of destination computed here */ /* It is stored in m_temp_reg_51 */ /* Column 6 of destination computed here */ /* It is stored in m_temp_reg_56 */ m_temp_reg_62 = _mm_add_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); //o2:1B*50-3B*89,5B*18+7B*75 m_temp_reg_20 = _mm_madd_epi16(m_temp_reg_60, m_coeff1); m_temp_reg_24 = _mm_madd_epi16(m_temp_reg_64, m_coeff2); m_temp_reg_51 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_56 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_63); /* o2 stored in m_temp_reg_34 and m_temp_reg_35 */ /* Loading coeff for computing o3 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[6][0]); m_coeff4 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[7][0]); m_temp_reg_34 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_24); /* Column 2 of destination computed here */ /* It is stored in m_temp_reg_52 */ /* Column 5 of destination computed here */ /* It is stored in m_temp_reg_55 */ m_temp_reg_62 = _mm_add_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); //o3:1B*18-3B*50,5B*75-7B*89 m_temp_reg_22 = _mm_madd_epi16(m_temp_reg_60, m_coeff3); m_temp_reg_26 = _mm_madd_epi16(m_temp_reg_64, m_coeff4); m_temp_reg_52 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_55 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_63); /* o3 stored in m_temp_reg_36 and m_temp_reg_37 */ m_temp_reg_36 = _mm_add_epi32(m_temp_reg_22, m_temp_reg_26); /* Column 3 of destination computed here */ /* It is stored in m_temp_reg_53 */ /* Column 4 of destination computed here */ /* It is stored in m_temp_reg_54 */ m_temp_reg_62 = _mm_add_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); m_temp_reg_53 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_54 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_63); } } /* Transpose of the destination 8x8 matrix done here */ /* and ultimately stored in registers m_temp_reg_50 to m_temp_reg_57 */ /* respectively */ { m_temp_reg_10 = _mm_unpacklo_epi16(m_temp_reg_50, m_temp_reg_51); m_temp_reg_11 = _mm_unpacklo_epi16(m_temp_reg_52, m_temp_reg_53); //m_temp_reg_14 = _mm_unpackhi_epi16(m_temp_reg_50, m_temp_reg_51); //m_temp_reg_15 = _mm_unpackhi_epi16(m_temp_reg_52, m_temp_reg_53); m_temp_reg_0 = _mm_unpacklo_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_1 = _mm_unpackhi_epi32(m_temp_reg_10, m_temp_reg_11); //m_temp_reg_2 = _mm_unpacklo_epi32(m_temp_reg_14, m_temp_reg_15); //m_temp_reg_3 = _mm_unpackhi_epi32(m_temp_reg_14, m_temp_reg_15); m_temp_reg_12 = _mm_unpacklo_epi16(m_temp_reg_54, m_temp_reg_55); m_temp_reg_13 = _mm_unpacklo_epi16(m_temp_reg_56, m_temp_reg_57); //m_temp_reg_16 = _mm_unpackhi_epi16(m_temp_reg_54, m_temp_reg_55); //m_temp_reg_17 = _mm_unpackhi_epi16(m_temp_reg_56, m_temp_reg_57); m_temp_reg_4 = _mm_unpacklo_epi32(m_temp_reg_12, m_temp_reg_13); m_temp_reg_5 = _mm_unpackhi_epi32(m_temp_reg_12, m_temp_reg_13); //m_temp_reg_6 = _mm_unpacklo_epi32(m_temp_reg_16, m_temp_reg_17); //m_temp_reg_7 = _mm_unpackhi_epi32(m_temp_reg_16, m_temp_reg_17); m_temp_reg_50 = _mm_unpacklo_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_51 = _mm_unpackhi_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_52 = _mm_unpacklo_epi64(m_temp_reg_1, m_temp_reg_5); m_temp_reg_53 = _mm_unpackhi_epi64(m_temp_reg_1, m_temp_reg_5); /*m_temp_reg_54 = _mm_unpacklo_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_55 = _mm_unpackhi_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_56 = _mm_unpacklo_epi64(m_temp_reg_3, m_temp_reg_7); m_temp_reg_57 = _mm_unpackhi_epi64(m_temp_reg_3, m_temp_reg_7); */ m_temp_reg_54 = _mm_setzero_si128(); m_temp_reg_55 = _mm_setzero_si128(); m_temp_reg_56 = _mm_setzero_si128(); m_temp_reg_57 = _mm_setzero_si128(); } } /* Stage 2 */ i4_shift = IT_SHIFT_STAGE_2; { /* ee0 is present in the registers m_temp_reg_10 and m_temp_reg_11 */ /* ee1 is present in the registers m_temp_reg_12 and m_temp_reg_13 */ { m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[0][0]); //add m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[3][0]); //sub m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_50, m_temp_reg_54); m_temp_reg_1 = _mm_unpackhi_epi16(m_temp_reg_50, m_temp_reg_54); m_temp_reg_10 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_12 = _mm_madd_epi16(m_temp_reg_0, m_coeff2); m_temp_reg_11 = _mm_madd_epi16(m_temp_reg_1, m_coeff1); m_temp_reg_13 = _mm_madd_epi16(m_temp_reg_1, m_coeff2); m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[1][0]); m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[2][0]); } /* eo0 is present in the registers m_temp_reg_14 and m_temp_reg_15 */ /* eo1 is present in the registers m_temp_reg_16 and m_temp_reg_17 */ { m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_52, m_temp_reg_56); m_temp_reg_1 = _mm_unpackhi_epi16(m_temp_reg_52, m_temp_reg_56); m_temp_reg_16 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_14 = _mm_madd_epi16(m_temp_reg_0, m_coeff2); m_temp_reg_17 = _mm_madd_epi16(m_temp_reg_1, m_coeff1); m_temp_reg_15 = _mm_madd_epi16(m_temp_reg_1, m_coeff2); /* Loading coeff for computing o0 in the next block */ m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[0][0]); m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_51, m_temp_reg_53); m_temp_reg_1 = _mm_unpackhi_epi16(m_temp_reg_51, m_temp_reg_53); /* e */ /* e0 stored in m_temp_reg_40 and m_temp_reg_41 */ /* e1 stored in m_temp_reg_42 and m_temp_reg_43 */ /* e3 stored in m_temp_reg_46 and m_temp_reg_47 */ /* e2 stored in m_temp_reg_44 and m_temp_reg_45 */ m_temp_reg_42 = _mm_add_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_44 = _mm_sub_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_40 = _mm_add_epi32(m_temp_reg_10, m_temp_reg_14); m_temp_reg_46 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_14); m_temp_reg_43 = _mm_add_epi32(m_temp_reg_13, m_temp_reg_17); m_temp_reg_45 = _mm_sub_epi32(m_temp_reg_13, m_temp_reg_17); m_temp_reg_41 = _mm_add_epi32(m_temp_reg_11, m_temp_reg_15); m_temp_reg_47 = _mm_sub_epi32(m_temp_reg_11, m_temp_reg_15); } /* o */ { //m_temp_reg_4 = _mm_unpacklo_epi16(m_temp_reg_55,m_temp_reg_57); //m_temp_reg_5 = _mm_unpackhi_epi16(m_temp_reg_55,m_temp_reg_57); /* o0 stored in m_temp_reg_30 and m_temp_reg_31 */ { //o0:1B*89+3B*75,1T*89+3T*75 m_temp_reg_30 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_31 = _mm_madd_epi16(m_temp_reg_1, m_coeff1); m_rdng_factor = _mm_cvtsi32_si128((1 << (i4_shift - 1))); m_rdng_factor = _mm_shuffle_epi32(m_rdng_factor, 0x0000); /* Loading coeff for computing o1 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[2][0]); /* Column 0 of destination computed here */ /* It is stored in m_temp_reg_50 */ /* Column 7 of destination computed here */ /* It is stored in m_temp_reg_57 */ m_temp_reg_2 = _mm_add_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_6 = _mm_sub_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_41, m_temp_reg_31); m_temp_reg_7 = _mm_sub_epi32(m_temp_reg_41, m_temp_reg_31); m_temp_reg_2 = _mm_add_epi32(m_temp_reg_2, m_rdng_factor); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_3, m_rdng_factor); m_temp_reg_6 = _mm_add_epi32(m_temp_reg_6, m_rdng_factor); m_temp_reg_7 = _mm_add_epi32(m_temp_reg_7, m_rdng_factor); //o1:1B*75-3B*18,1T*75-3T*18 m_temp_reg_32 = _mm_madd_epi16(m_temp_reg_0, m_coeff3); m_temp_reg_33 = _mm_madd_epi16(m_temp_reg_1, m_coeff3); m_temp_reg_2 = _mm_srai_epi32(m_temp_reg_2, i4_shift); m_temp_reg_3 = _mm_srai_epi32(m_temp_reg_3, i4_shift); m_temp_reg_6 = _mm_srai_epi32(m_temp_reg_6, i4_shift); m_temp_reg_7 = _mm_srai_epi32(m_temp_reg_7, i4_shift); m_temp_reg_50 = _mm_packs_epi32(m_temp_reg_2, m_temp_reg_3); m_temp_reg_57 = _mm_packs_epi32(m_temp_reg_6, m_temp_reg_7); /* o1 stored in m_temp_reg_32 and m_temp_reg_33 */ /* Loading coeff for computing o2 in the next block */ m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[4][0]); /* Column 1 of destination computed here */ /* It is stored in m_temp_reg_51 */ /* Column 6 of destination computed here */ /* It is stored in m_temp_reg_56 */ m_temp_reg_2 = _mm_add_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_6 = _mm_sub_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_43, m_temp_reg_33); m_temp_reg_7 = _mm_sub_epi32(m_temp_reg_43, m_temp_reg_33); m_temp_reg_2 = _mm_add_epi32(m_temp_reg_2, m_rdng_factor); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_3, m_rdng_factor); m_temp_reg_6 = _mm_add_epi32(m_temp_reg_6, m_rdng_factor); m_temp_reg_7 = _mm_add_epi32(m_temp_reg_7, m_rdng_factor); //o2:1B*50-3B*89,5T*18+7T*75. m_temp_reg_34 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_35 = _mm_madd_epi16(m_temp_reg_1, m_coeff1); m_temp_reg_2 = _mm_srai_epi32(m_temp_reg_2, i4_shift); m_temp_reg_3 = _mm_srai_epi32(m_temp_reg_3, i4_shift); m_temp_reg_6 = _mm_srai_epi32(m_temp_reg_6, i4_shift); m_temp_reg_7 = _mm_srai_epi32(m_temp_reg_7, i4_shift); m_temp_reg_51 = _mm_packs_epi32(m_temp_reg_2, m_temp_reg_3); m_temp_reg_56 = _mm_packs_epi32(m_temp_reg_6, m_temp_reg_7); /* o2 stored in m_temp_reg_34 and m_temp_reg_35 */ /* Loading coeff for computing o3 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[6][0]); /* Column 2 of destination computed here */ /* It is stored in m_temp_reg_52 */ /* Column 5 of destination computed here */ /* It is stored in m_temp_reg_55 */ m_temp_reg_2 = _mm_add_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_6 = _mm_sub_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_45, m_temp_reg_35); m_temp_reg_7 = _mm_sub_epi32(m_temp_reg_45, m_temp_reg_35); m_temp_reg_2 = _mm_add_epi32(m_temp_reg_2, m_rdng_factor); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_3, m_rdng_factor); m_temp_reg_6 = _mm_add_epi32(m_temp_reg_6, m_rdng_factor); m_temp_reg_7 = _mm_add_epi32(m_temp_reg_7, m_rdng_factor); //o3:1B*18-3B*50,1T*18-3T*50 m_temp_reg_36 = _mm_madd_epi16(m_temp_reg_0, m_coeff3); m_temp_reg_37 = _mm_madd_epi16(m_temp_reg_1, m_coeff3); m_temp_reg_2 = _mm_srai_epi32(m_temp_reg_2, i4_shift); m_temp_reg_3 = _mm_srai_epi32(m_temp_reg_3, i4_shift); m_temp_reg_6 = _mm_srai_epi32(m_temp_reg_6, i4_shift); m_temp_reg_7 = _mm_srai_epi32(m_temp_reg_7, i4_shift); m_temp_reg_52 = _mm_packs_epi32(m_temp_reg_2, m_temp_reg_3); m_temp_reg_55 = _mm_packs_epi32(m_temp_reg_6, m_temp_reg_7); /* o3 stored in m_temp_reg_36 and m_temp_reg_37 */ /* Column 3 of destination computed here */ /* It is stored in m_temp_reg_53 */ /* Column 4 of destination computed here */ /* It is stored in m_temp_reg_54 */ m_temp_reg_20 = _mm_add_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_22 = _mm_sub_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_21 = _mm_add_epi32(m_temp_reg_47, m_temp_reg_37); m_temp_reg_23 = _mm_sub_epi32(m_temp_reg_47, m_temp_reg_37); m_temp_reg_20 = _mm_add_epi32(m_temp_reg_20, m_rdng_factor); m_temp_reg_21 = _mm_add_epi32(m_temp_reg_21, m_rdng_factor); m_temp_reg_22 = _mm_add_epi32(m_temp_reg_22, m_rdng_factor); m_temp_reg_23 = _mm_add_epi32(m_temp_reg_23, m_rdng_factor); m_temp_reg_20 = _mm_srai_epi32(m_temp_reg_20, i4_shift); m_temp_reg_21 = _mm_srai_epi32(m_temp_reg_21, i4_shift); m_temp_reg_22 = _mm_srai_epi32(m_temp_reg_22, i4_shift); m_temp_reg_23 = _mm_srai_epi32(m_temp_reg_23, i4_shift); m_temp_reg_53 = _mm_packs_epi32(m_temp_reg_20, m_temp_reg_21); m_temp_reg_54 = _mm_packs_epi32(m_temp_reg_22, m_temp_reg_23); } } /* Transpose of the destination 8x8 matrix done here */ /* and ultimately stored in registers m_temp_reg_50 to m_temp_reg_57 */ /* respectively */ { m_temp_reg_10 = _mm_unpacklo_epi16(m_temp_reg_50, m_temp_reg_51); m_temp_reg_11 = _mm_unpacklo_epi16(m_temp_reg_52, m_temp_reg_53); m_temp_reg_14 = _mm_unpackhi_epi16(m_temp_reg_50, m_temp_reg_51); m_temp_reg_15 = _mm_unpackhi_epi16(m_temp_reg_52, m_temp_reg_53); m_temp_reg_0 = _mm_unpacklo_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_1 = _mm_unpackhi_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_2 = _mm_unpacklo_epi32(m_temp_reg_14, m_temp_reg_15); m_temp_reg_3 = _mm_unpackhi_epi32(m_temp_reg_14, m_temp_reg_15); m_temp_reg_12 = _mm_unpacklo_epi16(m_temp_reg_54, m_temp_reg_55); m_temp_reg_13 = _mm_unpacklo_epi16(m_temp_reg_56, m_temp_reg_57); m_temp_reg_16 = _mm_unpackhi_epi16(m_temp_reg_54, m_temp_reg_55); m_temp_reg_17 = _mm_unpackhi_epi16(m_temp_reg_56, m_temp_reg_57); m_temp_reg_4 = _mm_unpacklo_epi32(m_temp_reg_12, m_temp_reg_13); m_temp_reg_5 = _mm_unpackhi_epi32(m_temp_reg_12, m_temp_reg_13); m_temp_reg_6 = _mm_unpacklo_epi32(m_temp_reg_16, m_temp_reg_17); m_temp_reg_7 = _mm_unpackhi_epi32(m_temp_reg_16, m_temp_reg_17); m_temp_reg_10 = _mm_unpacklo_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_11 = _mm_unpackhi_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_12 = _mm_unpacklo_epi64(m_temp_reg_1, m_temp_reg_5); m_temp_reg_13 = _mm_unpackhi_epi64(m_temp_reg_1, m_temp_reg_5); m_temp_reg_14 = _mm_unpacklo_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_15 = _mm_unpackhi_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_16 = _mm_unpacklo_epi64(m_temp_reg_3, m_temp_reg_7); m_temp_reg_17 = _mm_unpackhi_epi64(m_temp_reg_3, m_temp_reg_7); } /* Recon and store */ { m_temp_reg_0 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_1 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_2 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_3 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_4 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_5 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_6 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_7 = _mm_loadl_epi64((__m128i *)pu1_pred); m_temp_reg_50 = _mm_setzero_si128(); m_temp_reg_0 = _mm_unpacklo_epi8(m_temp_reg_0, m_temp_reg_50); m_temp_reg_1 = _mm_unpacklo_epi8(m_temp_reg_1, m_temp_reg_50); m_temp_reg_2 = _mm_unpacklo_epi8(m_temp_reg_2, m_temp_reg_50); m_temp_reg_3 = _mm_unpacklo_epi8(m_temp_reg_3, m_temp_reg_50); m_temp_reg_4 = _mm_unpacklo_epi8(m_temp_reg_4, m_temp_reg_50); m_temp_reg_5 = _mm_unpacklo_epi8(m_temp_reg_5, m_temp_reg_50); m_temp_reg_6 = _mm_unpacklo_epi8(m_temp_reg_6, m_temp_reg_50); m_temp_reg_7 = _mm_unpacklo_epi8(m_temp_reg_7, m_temp_reg_50); m_temp_reg_50 = _mm_add_epi16(m_temp_reg_10, m_temp_reg_0); m_temp_reg_51 = _mm_add_epi16(m_temp_reg_11, m_temp_reg_1); m_temp_reg_52 = _mm_add_epi16(m_temp_reg_12, m_temp_reg_2); m_temp_reg_53 = _mm_add_epi16(m_temp_reg_13, m_temp_reg_3); m_temp_reg_54 = _mm_add_epi16(m_temp_reg_14, m_temp_reg_4); m_temp_reg_55 = _mm_add_epi16(m_temp_reg_15, m_temp_reg_5); m_temp_reg_56 = _mm_add_epi16(m_temp_reg_16, m_temp_reg_6); m_temp_reg_57 = _mm_add_epi16(m_temp_reg_17, m_temp_reg_7); m_temp_reg_50 = _mm_packus_epi16(m_temp_reg_50, m_temp_reg_50); m_temp_reg_51 = _mm_packus_epi16(m_temp_reg_51, m_temp_reg_51); m_temp_reg_52 = _mm_packus_epi16(m_temp_reg_52, m_temp_reg_52); m_temp_reg_53 = _mm_packus_epi16(m_temp_reg_53, m_temp_reg_53); m_temp_reg_54 = _mm_packus_epi16(m_temp_reg_54, m_temp_reg_54); m_temp_reg_55 = _mm_packus_epi16(m_temp_reg_55, m_temp_reg_55); m_temp_reg_56 = _mm_packus_epi16(m_temp_reg_56, m_temp_reg_56); m_temp_reg_57 = _mm_packus_epi16(m_temp_reg_57, m_temp_reg_57); _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_50); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_51); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_52); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_53); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_54); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_55); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_56); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_57); pu1_dst += dst_strd; } } } else { /* ee0 is present in the registers m_temp_reg_10 and m_temp_reg_11 */ /* ee1 is present in the registers m_temp_reg_12 and m_temp_reg_13 */ if(!check_row_stage_1) { /* ee0 is present in the registers m_temp_reg_10 and m_temp_reg_11 */ /* ee1 is present in the registers m_temp_reg_12 and m_temp_reg_13 */ { //Interleaving 0,4 row in 0 , 1 Rishab /*coef2 for m_temp_reg_12 and m_temp_reg_13 , coef1 for m_temp_reg_10 and m_temp_reg_11*/ m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[3][0]); m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[0][0]); m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_70, m_temp_reg_74); m_temp_reg_1 = _mm_unpackhi_epi16(m_temp_reg_70, m_temp_reg_74); m_temp_reg_10 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_12 = _mm_madd_epi16(m_temp_reg_0, m_coeff2); m_temp_reg_11 = _mm_madd_epi16(m_temp_reg_1, m_coeff1); m_temp_reg_13 = _mm_madd_epi16(m_temp_reg_1, m_coeff2); } /* eo0 is present in the registers m_temp_reg_14 and m_temp_reg_15 */ /* eo1 is present in the registers m_temp_reg_16 and m_temp_reg_17 */ { m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[1][0]); //sub 2B*36-6B*83 ,2T*36-6T*83 m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[2][0]); //add 2B*83+6B*36 ,2T*83+6T*36 /* Combining instructions to eliminate them based on zero_rows : Lokesh */ //Interleaving 2,6 row in 4, 5 Rishab m_temp_reg_4 = _mm_unpacklo_epi16(m_temp_reg_72, m_temp_reg_76); m_temp_reg_5 = _mm_unpackhi_epi16(m_temp_reg_72, m_temp_reg_76); m_temp_reg_16 = _mm_madd_epi16(m_temp_reg_4, m_coeff1); m_temp_reg_14 = _mm_madd_epi16(m_temp_reg_4, m_coeff2); m_temp_reg_17 = _mm_madd_epi16(m_temp_reg_5, m_coeff1); m_temp_reg_15 = _mm_madd_epi16(m_temp_reg_5, m_coeff2); /* Loading coeff for computing o0, o1, o2 and o3 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[2][0]); //m_coeff4 = _mm_loadu_si128((__m128i *) &g_ai2_ihevc_trans_intr_odd_8[3][0]); m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[0][0]); //m_coeff2 = _mm_loadu_si128((__m128i *) &g_ai2_ihevc_trans_intr_odd_8[1][0]); } /* e */ { /* e0 stored in m_temp_reg_40 and m_temp_reg_41 */ /* e1 stored in m_temp_reg_42 and m_temp_reg_43 */ /* e3 stored in m_temp_reg_46 and m_temp_reg_47 */ /* e2 stored in m_temp_reg_44 and m_temp_reg_45 */ m_temp_reg_42 = _mm_add_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_44 = _mm_sub_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_40 = _mm_add_epi32(m_temp_reg_10, m_temp_reg_14); m_temp_reg_46 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_14); m_temp_reg_43 = _mm_add_epi32(m_temp_reg_13, m_temp_reg_17); m_temp_reg_45 = _mm_sub_epi32(m_temp_reg_13, m_temp_reg_17); m_temp_reg_41 = _mm_add_epi32(m_temp_reg_11, m_temp_reg_15); m_temp_reg_47 = _mm_sub_epi32(m_temp_reg_11, m_temp_reg_15); } /* o */ { /* o0 stored in m_temp_reg_30 and m_temp_reg_31 */ { m_temp_reg_60 = _mm_unpacklo_epi16(m_temp_reg_71, m_temp_reg_73); m_temp_reg_61 = _mm_unpackhi_epi16(m_temp_reg_71, m_temp_reg_73); // m_temp_reg_64 = _mm_unpacklo_epi16(m_temp_reg_75,m_temp_reg_77); // m_temp_reg_65 = _mm_unpackhi_epi16(m_temp_reg_75,m_temp_reg_77); //o0:1B*89+3B*75,1T*89+3T*75 m_temp_reg_30 = _mm_madd_epi16(m_temp_reg_60, m_coeff1); m_temp_reg_31 = _mm_madd_epi16(m_temp_reg_61, m_coeff1); //m_temp_reg_24 = _mm_madd_epi16(m_temp_reg_64, m_coeff2); //m_temp_reg_25 = _mm_madd_epi16(m_temp_reg_65, m_coeff2); m_rdng_factor = _mm_cvtsi32_si128((1 << (i4_shift - 1))); m_rdng_factor = _mm_shuffle_epi32(m_rdng_factor, 0x0000); //m_temp_reg_30 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_24); //m_temp_reg_31 = _mm_add_epi32(m_temp_reg_21, m_temp_reg_25); } /* Column 0 of destination computed here */ /* It is stored in m_temp_reg_50 */ /* Column 7 of destination computed here */ /* It is stored in m_temp_reg_57 */ { m_temp_reg_62 = _mm_add_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_41, m_temp_reg_31); m_temp_reg_67 = _mm_sub_epi32(m_temp_reg_41, m_temp_reg_31); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_63, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_67 = _mm_add_epi32(m_temp_reg_67, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_63 = _mm_srai_epi32(m_temp_reg_63, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); m_temp_reg_67 = _mm_srai_epi32(m_temp_reg_67, i4_shift); //o1:1B*75-3B*18,1T*75-3T*18,5B*89+7B*50,5T*89+7T*50 m_temp_reg_32 = _mm_madd_epi16(m_temp_reg_60, m_coeff3); //m_temp_reg_26 = _mm_madd_epi16(m_temp_reg_64, m_coeff4); m_temp_reg_33 = _mm_madd_epi16(m_temp_reg_61, m_coeff3); //m_temp_reg_27 = _mm_madd_epi16(m_temp_reg_65, m_coeff4); m_temp_reg_50 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_57 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_67); /* Loading coeff for computing o2 in the next block */ m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[4][0]); /* o1 stored in m_temp_reg_32 and m_temp_reg_33 */ //m_temp_reg_32 = _mm_sub_epi32(m_temp_reg_22, m_temp_reg_26); //m_temp_reg_33 = _mm_sub_epi32(m_temp_reg_23, m_temp_reg_27); } /* Column 1 of destination computed here */ /* It is stored in m_temp_reg_51 */ /* Column 6 of destination computed here */ /* It is stored in m_temp_reg_56 */ { m_temp_reg_62 = _mm_add_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_43, m_temp_reg_33); m_temp_reg_67 = _mm_sub_epi32(m_temp_reg_43, m_temp_reg_33); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_63, m_rdng_factor); m_temp_reg_67 = _mm_add_epi32(m_temp_reg_67, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_63 = _mm_srai_epi32(m_temp_reg_63, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); m_temp_reg_67 = _mm_srai_epi32(m_temp_reg_67, i4_shift); //o2:1B*50-3B*89,1T*50-3T*89 m_temp_reg_34 = _mm_madd_epi16(m_temp_reg_60, m_coeff1); m_temp_reg_35 = _mm_madd_epi16(m_temp_reg_61, m_coeff1); m_temp_reg_51 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_56 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_67); /* o2 stored in m_temp_reg_34 and m_temp_reg_35 */ /* Loading coeff for computing o3 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[6][0]); //m_temp_reg_34 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_24); //m_temp_reg_35 = _mm_add_epi32(m_temp_reg_21, m_temp_reg_25); } /* Column 2 of destination computed here */ /* It is stored in m_temp_reg_52 */ /* Column 5 of destination computed here */ /* It is stored in m_temp_reg_55 */ { m_temp_reg_62 = _mm_add_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_45, m_temp_reg_35); m_temp_reg_67 = _mm_sub_epi32(m_temp_reg_45, m_temp_reg_35); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_63, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_67 = _mm_add_epi32(m_temp_reg_67, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_63 = _mm_srai_epi32(m_temp_reg_63, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); m_temp_reg_67 = _mm_srai_epi32(m_temp_reg_67, i4_shift); //o3:1B*18-3B*50,1T*18-3T*50 m_temp_reg_36 = _mm_madd_epi16(m_temp_reg_60, m_coeff3); m_temp_reg_37 = _mm_madd_epi16(m_temp_reg_61, m_coeff3); m_temp_reg_52 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_55 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_67); /* o3 stored in m_temp_reg_36 and m_temp_reg_37 */ } /* Column 3 of destination computed here */ /* It is stored in m_temp_reg_53 */ /* Column 4 of destination computed here */ /* It is stored in m_temp_reg_54 */ { m_temp_reg_62 = _mm_add_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_47, m_temp_reg_37); m_temp_reg_67 = _mm_sub_epi32(m_temp_reg_47, m_temp_reg_37); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_63, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_67 = _mm_add_epi32(m_temp_reg_67, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_63 = _mm_srai_epi32(m_temp_reg_63, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); m_temp_reg_67 = _mm_srai_epi32(m_temp_reg_67, i4_shift); m_temp_reg_53 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_54 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_67); } } /* Transpose of the destination 8x8 matrix done here */ /* and ultimately stored in registers m_temp_reg_50 to m_temp_reg_57 */ /* respectively */ { m_temp_reg_10 = _mm_unpacklo_epi16(m_temp_reg_50, m_temp_reg_51); m_temp_reg_11 = _mm_unpacklo_epi16(m_temp_reg_52, m_temp_reg_53); m_temp_reg_14 = _mm_unpackhi_epi16(m_temp_reg_50, m_temp_reg_51); m_temp_reg_15 = _mm_unpackhi_epi16(m_temp_reg_52, m_temp_reg_53); m_temp_reg_0 = _mm_unpacklo_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_1 = _mm_unpackhi_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_2 = _mm_unpacklo_epi32(m_temp_reg_14, m_temp_reg_15); m_temp_reg_3 = _mm_unpackhi_epi32(m_temp_reg_14, m_temp_reg_15); m_temp_reg_12 = _mm_unpacklo_epi16(m_temp_reg_54, m_temp_reg_55); m_temp_reg_13 = _mm_unpacklo_epi16(m_temp_reg_56, m_temp_reg_57); m_temp_reg_16 = _mm_unpackhi_epi16(m_temp_reg_54, m_temp_reg_55); m_temp_reg_17 = _mm_unpackhi_epi16(m_temp_reg_56, m_temp_reg_57); m_temp_reg_4 = _mm_unpacklo_epi32(m_temp_reg_12, m_temp_reg_13); m_temp_reg_5 = _mm_unpackhi_epi32(m_temp_reg_12, m_temp_reg_13); m_temp_reg_6 = _mm_unpacklo_epi32(m_temp_reg_16, m_temp_reg_17); m_temp_reg_7 = _mm_unpackhi_epi32(m_temp_reg_16, m_temp_reg_17); m_temp_reg_50 = _mm_unpacklo_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_51 = _mm_unpackhi_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_52 = _mm_unpacklo_epi64(m_temp_reg_1, m_temp_reg_5); m_temp_reg_53 = _mm_unpackhi_epi64(m_temp_reg_1, m_temp_reg_5); m_temp_reg_54 = _mm_unpacklo_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_55 = _mm_unpackhi_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_56 = _mm_unpacklo_epi64(m_temp_reg_3, m_temp_reg_7); m_temp_reg_57 = _mm_unpackhi_epi64(m_temp_reg_3, m_temp_reg_7); } } else { /* ee0 is present in the registers m_temp_reg_10 and m_temp_reg_11 */ /* ee1 is present in the registers m_temp_reg_12 and m_temp_reg_13 */ { //Interleaving 0,4 row in 0 , 1 Rishab /*coef2 for m_temp_reg_12 and m_temp_reg_13 , coef1 for m_temp_reg_10 and m_temp_reg_11*/ m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[3][0]); m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[0][0]); m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_70, m_temp_reg_74); m_temp_reg_1 = _mm_unpackhi_epi16(m_temp_reg_70, m_temp_reg_74); m_temp_reg_10 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_12 = _mm_madd_epi16(m_temp_reg_0, m_coeff2); m_temp_reg_11 = _mm_madd_epi16(m_temp_reg_1, m_coeff1); m_temp_reg_13 = _mm_madd_epi16(m_temp_reg_1, m_coeff2); } /* eo0 is present in the registers m_temp_reg_14 and m_temp_reg_15 */ /* eo1 is present in the registers m_temp_reg_16 and m_temp_reg_17 */ { m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[1][0]); //sub 2B*36-6B*83 ,2T*36-6T*83 m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[2][0]); //add 2B*83+6B*36 ,2T*83+6T*36 /* Combining instructions to eliminate them based on zero_rows : Lokesh */ //Interleaving 2,6 row in 4, 5 Rishab m_temp_reg_4 = _mm_unpacklo_epi16(m_temp_reg_72, m_temp_reg_76); m_temp_reg_5 = _mm_unpackhi_epi16(m_temp_reg_72, m_temp_reg_76); m_temp_reg_16 = _mm_madd_epi16(m_temp_reg_4, m_coeff1); m_temp_reg_14 = _mm_madd_epi16(m_temp_reg_4, m_coeff2); m_temp_reg_17 = _mm_madd_epi16(m_temp_reg_5, m_coeff1); m_temp_reg_15 = _mm_madd_epi16(m_temp_reg_5, m_coeff2); /* Loading coeff for computing o0, o1, o2 and o3 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[2][0]); m_coeff4 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[3][0]); m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[0][0]); m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[1][0]); } /* e */ { /* e0 stored in m_temp_reg_40 and m_temp_reg_41 */ /* e1 stored in m_temp_reg_42 and m_temp_reg_43 */ /* e3 stored in m_temp_reg_46 and m_temp_reg_47 */ /* e2 stored in m_temp_reg_44 and m_temp_reg_45 */ m_temp_reg_42 = _mm_add_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_44 = _mm_sub_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_40 = _mm_add_epi32(m_temp_reg_10, m_temp_reg_14); m_temp_reg_46 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_14); m_temp_reg_43 = _mm_add_epi32(m_temp_reg_13, m_temp_reg_17); m_temp_reg_45 = _mm_sub_epi32(m_temp_reg_13, m_temp_reg_17); m_temp_reg_41 = _mm_add_epi32(m_temp_reg_11, m_temp_reg_15); m_temp_reg_47 = _mm_sub_epi32(m_temp_reg_11, m_temp_reg_15); } /* o */ { /* o0 stored in m_temp_reg_30 and m_temp_reg_31 */ { m_temp_reg_60 = _mm_unpacklo_epi16(m_temp_reg_71, m_temp_reg_73); m_temp_reg_61 = _mm_unpackhi_epi16(m_temp_reg_71, m_temp_reg_73); m_temp_reg_64 = _mm_unpacklo_epi16(m_temp_reg_75, m_temp_reg_77); m_temp_reg_65 = _mm_unpackhi_epi16(m_temp_reg_75, m_temp_reg_77); //o0:1B*89+3B*75,1T*89+3T*75,5B*50+7B*18,5T*50+7T*18 m_temp_reg_20 = _mm_madd_epi16(m_temp_reg_60, m_coeff1); m_temp_reg_21 = _mm_madd_epi16(m_temp_reg_61, m_coeff1); m_temp_reg_24 = _mm_madd_epi16(m_temp_reg_64, m_coeff2); m_temp_reg_25 = _mm_madd_epi16(m_temp_reg_65, m_coeff2); m_rdng_factor = _mm_cvtsi32_si128((1 << (i4_shift - 1))); m_rdng_factor = _mm_shuffle_epi32(m_rdng_factor, 0x0000); m_temp_reg_30 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_24); m_temp_reg_31 = _mm_add_epi32(m_temp_reg_21, m_temp_reg_25); } /* Column 0 of destination computed here */ /* It is stored in m_temp_reg_50 */ /* Column 7 of destination computed here */ /* It is stored in m_temp_reg_57 */ { m_temp_reg_62 = _mm_add_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_41, m_temp_reg_31); m_temp_reg_67 = _mm_sub_epi32(m_temp_reg_41, m_temp_reg_31); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_63, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_67 = _mm_add_epi32(m_temp_reg_67, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_63 = _mm_srai_epi32(m_temp_reg_63, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); m_temp_reg_67 = _mm_srai_epi32(m_temp_reg_67, i4_shift); //o1:1B*75-3B*18,1T*75-3T*18,5B*89+7B*50,5T*89+7T*50 m_temp_reg_22 = _mm_madd_epi16(m_temp_reg_60, m_coeff3); m_temp_reg_26 = _mm_madd_epi16(m_temp_reg_64, m_coeff4); m_temp_reg_23 = _mm_madd_epi16(m_temp_reg_61, m_coeff3); m_temp_reg_27 = _mm_madd_epi16(m_temp_reg_65, m_coeff4); m_temp_reg_50 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_57 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_67); /* Loading coeff for computing o2 in the next block */ m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[4][0]); m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[5][0]); /* o1 stored in m_temp_reg_32 and m_temp_reg_33 */ m_temp_reg_32 = _mm_sub_epi32(m_temp_reg_22, m_temp_reg_26); m_temp_reg_33 = _mm_sub_epi32(m_temp_reg_23, m_temp_reg_27); } /* Column 1 of destination computed here */ /* It is stored in m_temp_reg_51 */ /* Column 6 of destination computed here */ /* It is stored in m_temp_reg_56 */ { m_temp_reg_62 = _mm_add_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_43, m_temp_reg_33); m_temp_reg_67 = _mm_sub_epi32(m_temp_reg_43, m_temp_reg_33); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_63, m_rdng_factor); m_temp_reg_67 = _mm_add_epi32(m_temp_reg_67, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_63 = _mm_srai_epi32(m_temp_reg_63, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); m_temp_reg_67 = _mm_srai_epi32(m_temp_reg_67, i4_shift); //o2:1B*50-3B*89,1T*50-3T*89,5B*18+7B*75,5T*18+7T*75 m_temp_reg_20 = _mm_madd_epi16(m_temp_reg_60, m_coeff1); m_temp_reg_24 = _mm_madd_epi16(m_temp_reg_64, m_coeff2); m_temp_reg_21 = _mm_madd_epi16(m_temp_reg_61, m_coeff1); m_temp_reg_25 = _mm_madd_epi16(m_temp_reg_65, m_coeff2); m_temp_reg_51 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_56 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_67); /* o2 stored in m_temp_reg_34 and m_temp_reg_35 */ /* Loading coeff for computing o3 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[6][0]); m_coeff4 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[7][0]); m_temp_reg_34 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_24); m_temp_reg_35 = _mm_add_epi32(m_temp_reg_21, m_temp_reg_25); } /* Column 2 of destination computed here */ /* It is stored in m_temp_reg_52 */ /* Column 5 of destination computed here */ /* It is stored in m_temp_reg_55 */ { m_temp_reg_62 = _mm_add_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_45, m_temp_reg_35); m_temp_reg_67 = _mm_sub_epi32(m_temp_reg_45, m_temp_reg_35); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_63, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_67 = _mm_add_epi32(m_temp_reg_67, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_63 = _mm_srai_epi32(m_temp_reg_63, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); m_temp_reg_67 = _mm_srai_epi32(m_temp_reg_67, i4_shift); //o3:1B*18-3B*50,1T*18-3T*50,5B*75-7B*89,5T*75-7T*89 m_temp_reg_22 = _mm_madd_epi16(m_temp_reg_60, m_coeff3); m_temp_reg_26 = _mm_madd_epi16(m_temp_reg_64, m_coeff4); m_temp_reg_23 = _mm_madd_epi16(m_temp_reg_61, m_coeff3); m_temp_reg_27 = _mm_madd_epi16(m_temp_reg_65, m_coeff4); m_temp_reg_52 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_55 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_67); /* o3 stored in m_temp_reg_36 and m_temp_reg_37 */ m_temp_reg_36 = _mm_add_epi32(m_temp_reg_22, m_temp_reg_26); m_temp_reg_37 = _mm_add_epi32(m_temp_reg_23, m_temp_reg_27); } /* Column 3 of destination computed here */ /* It is stored in m_temp_reg_53 */ /* Column 4 of destination computed here */ /* It is stored in m_temp_reg_54 */ { m_temp_reg_62 = _mm_add_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_66 = _mm_sub_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_47, m_temp_reg_37); m_temp_reg_67 = _mm_sub_epi32(m_temp_reg_47, m_temp_reg_37); m_temp_reg_62 = _mm_add_epi32(m_temp_reg_62, m_rdng_factor); m_temp_reg_63 = _mm_add_epi32(m_temp_reg_63, m_rdng_factor); m_temp_reg_66 = _mm_add_epi32(m_temp_reg_66, m_rdng_factor); m_temp_reg_67 = _mm_add_epi32(m_temp_reg_67, m_rdng_factor); m_temp_reg_62 = _mm_srai_epi32(m_temp_reg_62, i4_shift); m_temp_reg_63 = _mm_srai_epi32(m_temp_reg_63, i4_shift); m_temp_reg_66 = _mm_srai_epi32(m_temp_reg_66, i4_shift); m_temp_reg_67 = _mm_srai_epi32(m_temp_reg_67, i4_shift); m_temp_reg_53 = _mm_packs_epi32(m_temp_reg_62, m_temp_reg_63); m_temp_reg_54 = _mm_packs_epi32(m_temp_reg_66, m_temp_reg_67); } } /* Transpose of the destination 8x8 matrix done here */ /* and ultimately stored in registers m_temp_reg_50 to m_temp_reg_57 */ /* respectively */ { m_temp_reg_10 = _mm_unpacklo_epi16(m_temp_reg_50, m_temp_reg_51); m_temp_reg_11 = _mm_unpacklo_epi16(m_temp_reg_52, m_temp_reg_53); m_temp_reg_14 = _mm_unpackhi_epi16(m_temp_reg_50, m_temp_reg_51); m_temp_reg_15 = _mm_unpackhi_epi16(m_temp_reg_52, m_temp_reg_53); m_temp_reg_0 = _mm_unpacklo_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_1 = _mm_unpackhi_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_2 = _mm_unpacklo_epi32(m_temp_reg_14, m_temp_reg_15); m_temp_reg_3 = _mm_unpackhi_epi32(m_temp_reg_14, m_temp_reg_15); m_temp_reg_12 = _mm_unpacklo_epi16(m_temp_reg_54, m_temp_reg_55); m_temp_reg_13 = _mm_unpacklo_epi16(m_temp_reg_56, m_temp_reg_57); m_temp_reg_16 = _mm_unpackhi_epi16(m_temp_reg_54, m_temp_reg_55); m_temp_reg_17 = _mm_unpackhi_epi16(m_temp_reg_56, m_temp_reg_57); m_temp_reg_4 = _mm_unpacklo_epi32(m_temp_reg_12, m_temp_reg_13); m_temp_reg_5 = _mm_unpackhi_epi32(m_temp_reg_12, m_temp_reg_13); m_temp_reg_6 = _mm_unpacklo_epi32(m_temp_reg_16, m_temp_reg_17); m_temp_reg_7 = _mm_unpackhi_epi32(m_temp_reg_16, m_temp_reg_17); m_temp_reg_50 = _mm_unpacklo_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_51 = _mm_unpackhi_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_52 = _mm_unpacklo_epi64(m_temp_reg_1, m_temp_reg_5); m_temp_reg_53 = _mm_unpackhi_epi64(m_temp_reg_1, m_temp_reg_5); m_temp_reg_54 = _mm_unpacklo_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_55 = _mm_unpackhi_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_56 = _mm_unpacklo_epi64(m_temp_reg_3, m_temp_reg_7); m_temp_reg_57 = _mm_unpackhi_epi64(m_temp_reg_3, m_temp_reg_7); } } /* Stage 2 */ i4_shift = IT_SHIFT_STAGE_2; { /* ee0 is present in the registers m_temp_reg_10 and m_temp_reg_11 */ /* ee1 is present in the registers m_temp_reg_12 and m_temp_reg_13 */ { m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[0][0]); //add m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[3][0]); //sub m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_50, m_temp_reg_54); m_temp_reg_1 = _mm_unpackhi_epi16(m_temp_reg_50, m_temp_reg_54); m_temp_reg_10 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_12 = _mm_madd_epi16(m_temp_reg_0, m_coeff2); m_temp_reg_11 = _mm_madd_epi16(m_temp_reg_1, m_coeff1); m_temp_reg_13 = _mm_madd_epi16(m_temp_reg_1, m_coeff2); m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[1][0]); m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_even_8[2][0]); } /* eo0 is present in the registers m_temp_reg_14 and m_temp_reg_15 */ /* eo1 is present in the registers m_temp_reg_16 and m_temp_reg_17 */ { //m_temp_reg_66 = _mm_mullo_epi32(m_temp_reg_2, m_coeff1); //m_temp_reg_64 = _mm_mullo_epi32(m_temp_reg_0, m_coeff2); //m_temp_reg_62 = _mm_mullo_epi32(m_temp_reg_2, m_coeff2); //m_temp_reg_60 = _mm_mullo_epi32(m_temp_reg_0, m_coeff1); m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_52, m_temp_reg_56); m_temp_reg_1 = _mm_unpackhi_epi16(m_temp_reg_52, m_temp_reg_56); m_temp_reg_16 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_14 = _mm_madd_epi16(m_temp_reg_0, m_coeff2); m_temp_reg_17 = _mm_madd_epi16(m_temp_reg_1, m_coeff1); m_temp_reg_15 = _mm_madd_epi16(m_temp_reg_1, m_coeff2); //m_temp_reg_16 = _mm_sub_epi32(m_temp_reg_64, m_temp_reg_66); //m_temp_reg_14 = _mm_add_epi32(m_temp_reg_60, m_temp_reg_62); /* Loading coeff for computing o0 in the next block */ m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[0][0]); m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[1][0]); m_temp_reg_0 = _mm_unpacklo_epi16(m_temp_reg_51, m_temp_reg_53); m_temp_reg_1 = _mm_unpackhi_epi16(m_temp_reg_51, m_temp_reg_53); /*m_temp_reg_3 = _mm_srli_si128(m_temp_reg_53, 8); m_temp_reg_1 = _mm_cvtepi16_epi32(m_temp_reg_1); m_temp_reg_3 = _mm_cvtepi16_epi32(m_temp_reg_3); */ } /* e */ { /* e0 stored in m_temp_reg_40 and m_temp_reg_41 */ /* e1 stored in m_temp_reg_42 and m_temp_reg_43 */ /* e3 stored in m_temp_reg_46 and m_temp_reg_47 */ /* e2 stored in m_temp_reg_44 and m_temp_reg_45 */ m_temp_reg_42 = _mm_add_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_44 = _mm_sub_epi32(m_temp_reg_12, m_temp_reg_16); m_temp_reg_40 = _mm_add_epi32(m_temp_reg_10, m_temp_reg_14); m_temp_reg_46 = _mm_sub_epi32(m_temp_reg_10, m_temp_reg_14); m_temp_reg_43 = _mm_add_epi32(m_temp_reg_13, m_temp_reg_17); m_temp_reg_45 = _mm_sub_epi32(m_temp_reg_13, m_temp_reg_17); m_temp_reg_41 = _mm_add_epi32(m_temp_reg_11, m_temp_reg_15); m_temp_reg_47 = _mm_sub_epi32(m_temp_reg_11, m_temp_reg_15); } /* o */ { /* m_temp_reg_4 = _mm_cvtepi16_epi32(m_temp_reg_55); m_temp_reg_5 = _mm_srli_si128(m_temp_reg_55, 8); m_temp_reg_6 = _mm_cvtepi16_epi32(m_temp_reg_57); m_temp_reg_7 = _mm_srli_si128(m_temp_reg_57, 8); m_temp_reg_5 = _mm_cvtepi16_epi32(m_temp_reg_5); m_temp_reg_7 = _mm_cvtepi16_epi32(m_temp_reg_7); */ m_temp_reg_4 = _mm_unpacklo_epi16(m_temp_reg_55, m_temp_reg_57); m_temp_reg_5 = _mm_unpackhi_epi16(m_temp_reg_55, m_temp_reg_57); /* o0 stored in m_temp_reg_30 and m_temp_reg_31 */ { //o0:1B*89+3B*75,1T*89+3T*75,5B*50+7B*18,5T*50+7T*18 m_temp_reg_20 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_21 = _mm_madd_epi16(m_temp_reg_1, m_coeff1); m_temp_reg_24 = _mm_madd_epi16(m_temp_reg_4, m_coeff2); m_temp_reg_25 = _mm_madd_epi16(m_temp_reg_5, m_coeff2); m_rdng_factor = _mm_cvtsi32_si128((1 << (i4_shift - 1))); m_rdng_factor = _mm_shuffle_epi32(m_rdng_factor, 0x0000); /* Loading coeff for computing o1 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[2][0]); m_coeff4 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[3][0]); m_temp_reg_30 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_24); m_temp_reg_31 = _mm_add_epi32(m_temp_reg_21, m_temp_reg_25); } /* Column 0 of destination computed here */ /* It is stored in m_temp_reg_50 */ /* Column 7 of destination computed here */ /* It is stored in m_temp_reg_57 */ { m_temp_reg_2 = _mm_add_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_6 = _mm_sub_epi32(m_temp_reg_40, m_temp_reg_30); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_41, m_temp_reg_31); m_temp_reg_7 = _mm_sub_epi32(m_temp_reg_41, m_temp_reg_31); m_temp_reg_2 = _mm_add_epi32(m_temp_reg_2, m_rdng_factor); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_3, m_rdng_factor); m_temp_reg_6 = _mm_add_epi32(m_temp_reg_6, m_rdng_factor); m_temp_reg_7 = _mm_add_epi32(m_temp_reg_7, m_rdng_factor); m_temp_reg_2 = _mm_srai_epi32(m_temp_reg_2, i4_shift); m_temp_reg_3 = _mm_srai_epi32(m_temp_reg_3, i4_shift); m_temp_reg_6 = _mm_srai_epi32(m_temp_reg_6, i4_shift); m_temp_reg_7 = _mm_srai_epi32(m_temp_reg_7, i4_shift); //o1:1B*75-3B*18,1T*75-3T*18,5B*89+7B*50,5T*89+7T*50 m_temp_reg_22 = _mm_madd_epi16(m_temp_reg_0, m_coeff3); m_temp_reg_26 = _mm_madd_epi16(m_temp_reg_4, m_coeff4); m_temp_reg_23 = _mm_madd_epi16(m_temp_reg_1, m_coeff3); m_temp_reg_27 = _mm_madd_epi16(m_temp_reg_5, m_coeff4); m_temp_reg_50 = _mm_packs_epi32(m_temp_reg_2, m_temp_reg_3); m_temp_reg_57 = _mm_packs_epi32(m_temp_reg_6, m_temp_reg_7); /* o1 stored in m_temp_reg_32 and m_temp_reg_33 */ /* Loading coeff for computing o2 in the next block */ m_coeff1 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[4][0]); m_coeff2 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[5][0]); m_temp_reg_32 = _mm_sub_epi32(m_temp_reg_22, m_temp_reg_26); m_temp_reg_33 = _mm_sub_epi32(m_temp_reg_23, m_temp_reg_27); } /* Column 1 of destination computed here */ /* It is stored in m_temp_reg_51 */ /* Column 6 of destination computed here */ /* It is stored in m_temp_reg_56 */ { m_temp_reg_2 = _mm_add_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_6 = _mm_sub_epi32(m_temp_reg_42, m_temp_reg_32); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_43, m_temp_reg_33); m_temp_reg_7 = _mm_sub_epi32(m_temp_reg_43, m_temp_reg_33); m_temp_reg_2 = _mm_add_epi32(m_temp_reg_2, m_rdng_factor); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_3, m_rdng_factor); m_temp_reg_6 = _mm_add_epi32(m_temp_reg_6, m_rdng_factor); m_temp_reg_7 = _mm_add_epi32(m_temp_reg_7, m_rdng_factor); m_temp_reg_2 = _mm_srai_epi32(m_temp_reg_2, i4_shift); m_temp_reg_3 = _mm_srai_epi32(m_temp_reg_3, i4_shift); m_temp_reg_6 = _mm_srai_epi32(m_temp_reg_6, i4_shift); m_temp_reg_7 = _mm_srai_epi32(m_temp_reg_7, i4_shift); //o2:1B*50-3B*89,1T*50-3T*89,5B*18+7B*75,5T*18+7T*75 m_temp_reg_20 = _mm_madd_epi16(m_temp_reg_0, m_coeff1); m_temp_reg_24 = _mm_madd_epi16(m_temp_reg_4, m_coeff2); m_temp_reg_21 = _mm_madd_epi16(m_temp_reg_1, m_coeff1); m_temp_reg_25 = _mm_madd_epi16(m_temp_reg_5, m_coeff2); m_temp_reg_51 = _mm_packs_epi32(m_temp_reg_2, m_temp_reg_3); m_temp_reg_56 = _mm_packs_epi32(m_temp_reg_6, m_temp_reg_7); /* o2 stored in m_temp_reg_34 and m_temp_reg_35 */ /* Loading coeff for computing o3 in the next block */ m_coeff3 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[6][0]); m_coeff4 = _mm_load_si128((__m128i *)&g_ai2_ihevc_trans_intr_odd_8[7][0]); m_temp_reg_34 = _mm_add_epi32(m_temp_reg_20, m_temp_reg_24); m_temp_reg_35 = _mm_add_epi32(m_temp_reg_21, m_temp_reg_25); } /* Column 2 of destination computed here */ /* It is stored in m_temp_reg_52 */ /* Column 5 of destination computed here */ /* It is stored in m_temp_reg_55 */ { m_temp_reg_2 = _mm_add_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_6 = _mm_sub_epi32(m_temp_reg_44, m_temp_reg_34); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_45, m_temp_reg_35); m_temp_reg_7 = _mm_sub_epi32(m_temp_reg_45, m_temp_reg_35); m_temp_reg_2 = _mm_add_epi32(m_temp_reg_2, m_rdng_factor); m_temp_reg_3 = _mm_add_epi32(m_temp_reg_3, m_rdng_factor); m_temp_reg_6 = _mm_add_epi32(m_temp_reg_6, m_rdng_factor); m_temp_reg_7 = _mm_add_epi32(m_temp_reg_7, m_rdng_factor); m_temp_reg_2 = _mm_srai_epi32(m_temp_reg_2, i4_shift); m_temp_reg_3 = _mm_srai_epi32(m_temp_reg_3, i4_shift); m_temp_reg_6 = _mm_srai_epi32(m_temp_reg_6, i4_shift); m_temp_reg_7 = _mm_srai_epi32(m_temp_reg_7, i4_shift); //o3:1B*18-3B*50,1T*18-3T*50,5B*75-7B*89,5T*75-7T*89 m_temp_reg_22 = _mm_madd_epi16(m_temp_reg_0, m_coeff3); m_temp_reg_26 = _mm_madd_epi16(m_temp_reg_4, m_coeff4); m_temp_reg_23 = _mm_madd_epi16(m_temp_reg_1, m_coeff3); m_temp_reg_27 = _mm_madd_epi16(m_temp_reg_5, m_coeff4); m_temp_reg_52 = _mm_packs_epi32(m_temp_reg_2, m_temp_reg_3); m_temp_reg_55 = _mm_packs_epi32(m_temp_reg_6, m_temp_reg_7); /* o3 stored in m_temp_reg_36 and m_temp_reg_37 */ m_temp_reg_36 = _mm_add_epi32(m_temp_reg_22, m_temp_reg_26); m_temp_reg_37 = _mm_add_epi32(m_temp_reg_23, m_temp_reg_27); } /* Column 3 of destination computed here */ /* It is stored in m_temp_reg_53 */ /* Column 4 of destination computed here */ /* It is stored in m_temp_reg_54 */ { m_temp_reg_20 = _mm_add_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_22 = _mm_sub_epi32(m_temp_reg_46, m_temp_reg_36); m_temp_reg_21 = _mm_add_epi32(m_temp_reg_47, m_temp_reg_37); m_temp_reg_23 = _mm_sub_epi32(m_temp_reg_47, m_temp_reg_37); m_temp_reg_20 = _mm_add_epi32(m_temp_reg_20, m_rdng_factor); m_temp_reg_21 = _mm_add_epi32(m_temp_reg_21, m_rdng_factor); m_temp_reg_22 = _mm_add_epi32(m_temp_reg_22, m_rdng_factor); m_temp_reg_23 = _mm_add_epi32(m_temp_reg_23, m_rdng_factor); m_temp_reg_20 = _mm_srai_epi32(m_temp_reg_20, i4_shift); m_temp_reg_21 = _mm_srai_epi32(m_temp_reg_21, i4_shift); m_temp_reg_22 = _mm_srai_epi32(m_temp_reg_22, i4_shift); m_temp_reg_23 = _mm_srai_epi32(m_temp_reg_23, i4_shift); m_temp_reg_53 = _mm_packs_epi32(m_temp_reg_20, m_temp_reg_21); m_temp_reg_54 = _mm_packs_epi32(m_temp_reg_22, m_temp_reg_23); } } /* Transpose of the destination 8x8 matrix done here */ /* and ultimately stored in registers m_temp_reg_50 to m_temp_reg_57 */ /* respectively */ { m_temp_reg_10 = _mm_unpacklo_epi16(m_temp_reg_50, m_temp_reg_51); m_temp_reg_11 = _mm_unpacklo_epi16(m_temp_reg_52, m_temp_reg_53); m_temp_reg_14 = _mm_unpackhi_epi16(m_temp_reg_50, m_temp_reg_51); m_temp_reg_15 = _mm_unpackhi_epi16(m_temp_reg_52, m_temp_reg_53); m_temp_reg_0 = _mm_unpacklo_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_1 = _mm_unpackhi_epi32(m_temp_reg_10, m_temp_reg_11); m_temp_reg_2 = _mm_unpacklo_epi32(m_temp_reg_14, m_temp_reg_15); m_temp_reg_3 = _mm_unpackhi_epi32(m_temp_reg_14, m_temp_reg_15); m_temp_reg_12 = _mm_unpacklo_epi16(m_temp_reg_54, m_temp_reg_55); m_temp_reg_13 = _mm_unpacklo_epi16(m_temp_reg_56, m_temp_reg_57); m_temp_reg_16 = _mm_unpackhi_epi16(m_temp_reg_54, m_temp_reg_55); m_temp_reg_17 = _mm_unpackhi_epi16(m_temp_reg_56, m_temp_reg_57); m_temp_reg_4 = _mm_unpacklo_epi32(m_temp_reg_12, m_temp_reg_13); m_temp_reg_5 = _mm_unpackhi_epi32(m_temp_reg_12, m_temp_reg_13); m_temp_reg_6 = _mm_unpacklo_epi32(m_temp_reg_16, m_temp_reg_17); m_temp_reg_7 = _mm_unpackhi_epi32(m_temp_reg_16, m_temp_reg_17); m_temp_reg_10 = _mm_unpacklo_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_11 = _mm_unpackhi_epi64(m_temp_reg_0, m_temp_reg_4); m_temp_reg_12 = _mm_unpacklo_epi64(m_temp_reg_1, m_temp_reg_5); m_temp_reg_13 = _mm_unpackhi_epi64(m_temp_reg_1, m_temp_reg_5); m_temp_reg_14 = _mm_unpacklo_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_15 = _mm_unpackhi_epi64(m_temp_reg_2, m_temp_reg_6); m_temp_reg_16 = _mm_unpacklo_epi64(m_temp_reg_3, m_temp_reg_7); m_temp_reg_17 = _mm_unpackhi_epi64(m_temp_reg_3, m_temp_reg_7); } /* Recon and store */ { m_temp_reg_0 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_1 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_2 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_3 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_4 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_5 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_6 = _mm_loadl_epi64((__m128i *)pu1_pred); pu1_pred += pred_strd; m_temp_reg_7 = _mm_loadl_epi64((__m128i *)pu1_pred); m_temp_reg_50 = _mm_setzero_si128(); m_temp_reg_0 = _mm_unpacklo_epi8(m_temp_reg_0, m_temp_reg_50); m_temp_reg_1 = _mm_unpacklo_epi8(m_temp_reg_1, m_temp_reg_50); m_temp_reg_2 = _mm_unpacklo_epi8(m_temp_reg_2, m_temp_reg_50); m_temp_reg_3 = _mm_unpacklo_epi8(m_temp_reg_3, m_temp_reg_50); m_temp_reg_4 = _mm_unpacklo_epi8(m_temp_reg_4, m_temp_reg_50); m_temp_reg_5 = _mm_unpacklo_epi8(m_temp_reg_5, m_temp_reg_50); m_temp_reg_6 = _mm_unpacklo_epi8(m_temp_reg_6, m_temp_reg_50); m_temp_reg_7 = _mm_unpacklo_epi8(m_temp_reg_7, m_temp_reg_50); m_temp_reg_50 = _mm_add_epi16(m_temp_reg_10, m_temp_reg_0); m_temp_reg_51 = _mm_add_epi16(m_temp_reg_11, m_temp_reg_1); m_temp_reg_52 = _mm_add_epi16(m_temp_reg_12, m_temp_reg_2); m_temp_reg_53 = _mm_add_epi16(m_temp_reg_13, m_temp_reg_3); m_temp_reg_54 = _mm_add_epi16(m_temp_reg_14, m_temp_reg_4); m_temp_reg_55 = _mm_add_epi16(m_temp_reg_15, m_temp_reg_5); m_temp_reg_56 = _mm_add_epi16(m_temp_reg_16, m_temp_reg_6); m_temp_reg_57 = _mm_add_epi16(m_temp_reg_17, m_temp_reg_7); m_temp_reg_50 = _mm_packus_epi16(m_temp_reg_50, m_temp_reg_50); m_temp_reg_51 = _mm_packus_epi16(m_temp_reg_51, m_temp_reg_51); m_temp_reg_52 = _mm_packus_epi16(m_temp_reg_52, m_temp_reg_52); m_temp_reg_53 = _mm_packus_epi16(m_temp_reg_53, m_temp_reg_53); m_temp_reg_54 = _mm_packus_epi16(m_temp_reg_54, m_temp_reg_54); m_temp_reg_55 = _mm_packus_epi16(m_temp_reg_55, m_temp_reg_55); m_temp_reg_56 = _mm_packus_epi16(m_temp_reg_56, m_temp_reg_56); m_temp_reg_57 = _mm_packus_epi16(m_temp_reg_57, m_temp_reg_57); _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_50); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_51); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_52); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_53); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_54); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_55); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_56); pu1_dst += dst_strd; _mm_storel_epi64((__m128i *)pu1_dst, m_temp_reg_57); pu1_dst += dst_strd; } } } }