/****************************************************************************** * * Copyright (C) 2015 The Android Open Source Project * * 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. * ***************************************************************************** * Originally developed and contributed by Ittiam Systems Pvt. Ltd, Bangalore */ /** ******************************************************************************* * @file * ideint_cac_ssse3.c * * @brief * This file include the definitions of the combing artifact check function * of the de-interlacer and some variant of that. * * @author * Ittiam * * @par List of Functions: * cac_4x8() * ideint_cac() * * @remarks * In the de-interlacer workspace, cac is not a seperate assembly module as * it comes along with the de_int_decision() function. But in C-Model, to * keep the things cleaner, it was made to be a separate function during * cac experiments long after the assembly was written by Mudit. * ******************************************************************************* */ /*****************************************************************************/ /* File Includes */ /*****************************************************************************/ /* System include files */ #include #include #include #include #include /* User include files */ #include "icv_datatypes.h" #include "icv_macros.h" #include "icv.h" #include "icv_variance.h" #include "icv_sad.h" #include "ideint.h" #include "ideint_defs.h" #include "ideint_structs.h" #include "ideint_cac.h" /** ******************************************************************************* * * @brief * Combing artifact check function for 8x8 block * * @par Description * Determines CAC for 8x8 block by calling 8x4 CAC function * * @param[in] pu1_top * Top field * * @param[in] pu1_bot * Bottom field * * @param[in] top_strd * Top field Stride * * @param[in] bot_strd * Bottom field stride * * @returns * combing artifact flag (1 = detected, 0 = not detected) * * @remarks * ******************************************************************************* */ WORD32 ideint_cac_8x8_ssse3(UWORD8 *pu1_top, UWORD8 *pu1_bot, WORD32 top_strd, WORD32 bot_strd) { WORD32 ca; /* combing artifact result */ WORD32 i; WORD32 adj[2] = {0}; WORD32 alt[2] = {0}; WORD32 sum_1, sum_2, sum_3, sum_4; WORD32 sum_diff, diff_sum; __m128i top[4]; __m128i bot[4]; __m128i sum_t[4]; __m128i sum_b[4]; __m128i zero; zero = _mm_setzero_si128(); for(i = 0; i < 4; i++) { /* Load top */ top[i] = (__m128i)_mm_loadl_epi64((__m128i *) (pu1_top)); pu1_top += top_strd; /* Load bottom */ bot[i] = (__m128i)_mm_loadl_epi64((__m128i *) (pu1_bot)); pu1_bot += bot_strd; /* Unpack */ top[i] = _mm_unpacklo_epi8(top[i], zero); bot[i] = _mm_unpacklo_epi8(bot[i], zero); /* Compute row sums */ sum_t[i] = _mm_sad_epu8(top[i], zero); sum_b[i] = _mm_sad_epu8(bot[i], zero); } /* Compute row based alt and adj */ for(i = 0; i < 4; i += 2) { sum_1 = _mm_cvtsi128_si32(sum_t[i + 0]); sum_2 = _mm_cvtsi128_si32(sum_b[i + 0]); sum_diff = ABS_DIF(sum_1, sum_2); if(sum_diff >= RSUM_CSUM_THRESH) adj[0] += sum_diff; sum_3 = _mm_cvtsi128_si32(sum_t[i + 1]); sum_4 = _mm_cvtsi128_si32(sum_b[i + 1]); sum_diff = ABS_DIF(sum_3, sum_4); if(sum_diff >= RSUM_CSUM_THRESH) adj[0] += sum_diff; alt[0] += ABS_DIF(sum_1, sum_3); alt[0] += ABS_DIF(sum_2, sum_4); sum_1 = _mm_cvtsi128_si32(_mm_srli_si128(sum_t[i + 0], 8)); sum_2 = _mm_cvtsi128_si32(_mm_srli_si128(sum_b[i + 0], 8)); sum_diff = ABS_DIF(sum_1, sum_2); if(sum_diff >= RSUM_CSUM_THRESH) adj[1] += sum_diff; sum_3 = _mm_cvtsi128_si32(_mm_srli_si128(sum_t[i + 1], 8)); sum_4 = _mm_cvtsi128_si32(_mm_srli_si128(sum_b[i + 1], 8)); sum_diff = ABS_DIF(sum_3, sum_4); if(sum_diff >= RSUM_CSUM_THRESH) adj[1] += sum_diff; alt[1] += ABS_DIF(sum_1, sum_3); alt[1] += ABS_DIF(sum_2, sum_4); } /* Compute column based adj */ { __m128i avg1, avg2; __m128i top_avg, bot_avg; __m128i min, max, diff, thresh; __m128i mask; avg1 = _mm_avg_epu8(top[0], top[1]); avg2 = _mm_avg_epu8(top[2], top[3]); top_avg = _mm_avg_epu8(avg1, avg2); avg1 = _mm_avg_epu8(bot[0], bot[1]); avg2 = _mm_avg_epu8(bot[2], bot[3]); bot_avg = _mm_avg_epu8(avg1, avg2); min = _mm_min_epu8(top_avg, bot_avg); max = _mm_max_epu8(top_avg, bot_avg); diff = _mm_sub_epi16(max, min); thresh = _mm_set1_epi16((RSUM_CSUM_THRESH >> 2) - 1); mask = _mm_cmpgt_epi16(diff, thresh); diff = _mm_and_si128(diff, mask); diff_sum = _mm_extract_epi16(diff, 0); diff_sum += _mm_extract_epi16(diff, 1); diff_sum += _mm_extract_epi16(diff, 2); diff_sum += _mm_extract_epi16(diff, 3); adj[0] += diff_sum << 2; diff_sum = _mm_extract_epi16(diff, 4); diff_sum += _mm_extract_epi16(diff, 5); diff_sum += _mm_extract_epi16(diff, 6); diff_sum += _mm_extract_epi16(diff, 7); adj[1] += diff_sum << 2; } /* Compute column based alt */ { __m128i avg1, avg2; __m128i even_avg, odd_avg, diff; avg1 = _mm_avg_epu8(top[0], bot[0]); avg2 = _mm_avg_epu8(top[2], bot[2]); even_avg = _mm_avg_epu8(avg1, avg2); avg1 = _mm_avg_epu8(top[1], bot[1]); avg2 = _mm_avg_epu8(top[3], bot[3]); odd_avg = _mm_avg_epu8(avg1, avg2); diff = _mm_sad_epu8(even_avg, odd_avg); diff_sum = _mm_cvtsi128_si32(diff); alt[0] += diff_sum << 2; diff_sum = _mm_cvtsi128_si32(_mm_srli_si128(diff, 8)); alt[1] += diff_sum << 2; } alt[0] += (alt[0] >> SAD_BIAS_MULT_SHIFT) + (SAD_BIAS_ADDITIVE >> 1); alt[1] += (alt[1] >> SAD_BIAS_MULT_SHIFT) + (SAD_BIAS_ADDITIVE >> 1); ca = (alt[0] < adj[0]); ca |= (alt[1] < adj[1]); return ca; }