/* * Copyright (c) 2001-2016, Alliance for Open Media. All rights reserved * * This source code is subject to the terms of the BSD 2 Clause License and * the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License * was not distributed with this source code in the LICENSE file, you can * obtain it at www.aomedia.org/license/software. If the Alliance for Open * Media Patent License 1.0 was not distributed with this source code in the * PATENTS file, you can obtain it at www.aomedia.org/license/patent. */ /* clang-format off */ #ifdef HAVE_CONFIG_H # include "config.h" #endif #include "enums.h" #include "odintrin.h" #include "partition.h" #include "zigzag.h" OD_EXTERN const index_pair *OD_ZIGZAG4[4] = { OD_ZIGZAG4_DCT_DCT, OD_ZIGZAG4_ADST_DCT, OD_ZIGZAG4_DCT_ADST, OD_ZIGZAG4_ADST_ADST }; OD_EXTERN const index_pair *OD_ZIGZAG8[4] = { OD_ZIGZAG8_DCT_DCT, OD_ZIGZAG8_ADST_DCT, OD_ZIGZAG8_DCT_ADST, OD_ZIGZAG8_ADST_ADST }; OD_EXTERN const index_pair *OD_ZIGZAG16[4] = { OD_ZIGZAG16_DCT_DCT, OD_ZIGZAG16_ADST_DCT, OD_ZIGZAG16_DCT_ADST, OD_ZIGZAG16_ADST_ADST }; OD_EXTERN const index_pair *OD_ZIGZAG32[4] = { OD_ZIGZAG32_DCT_DCT, OD_ZIGZAG32_DCT_DCT, OD_ZIGZAG32_DCT_DCT, OD_ZIGZAG32_DCT_DCT }; /* The tables below specify how coefficient blocks are translated to and from PVQ partition coding scan order for 4x4, 8x8 and 16x16 */ static const int OD_LAYOUT32_OFFSETS[4] = { 0, 128, 256, 768 }; const band_layout OD_LAYOUT32 = { OD_ZIGZAG32, 32, 3, OD_LAYOUT32_OFFSETS }; static const int OD_LAYOUT16_OFFSETS[4] = { 0, 32, 64, 192 }; const band_layout OD_LAYOUT16 = { OD_ZIGZAG16, 16, 3, OD_LAYOUT16_OFFSETS }; const int OD_LAYOUT8_OFFSETS[4] = { 0, 8, 16, 48 }; const band_layout OD_LAYOUT8 = { OD_ZIGZAG8, 8, 3, OD_LAYOUT8_OFFSETS }; static const int OD_LAYOUT4_OFFSETS[2] = { 0, 15 }; const band_layout OD_LAYOUT4 = { OD_ZIGZAG4, 4, 1, OD_LAYOUT4_OFFSETS }; /* First element is the number of bands, followed by the list all the band boundaries. */ static const int OD_BAND_OFFSETS4[] = {1, 1, 16}; static const int OD_BAND_OFFSETS8[] = {4, 1, 16, 24, 32, 64}; static const int OD_BAND_OFFSETS16[] = {7, 1, 16, 24, 32, 64, 96, 128, 256}; static const int OD_BAND_OFFSETS32[] = {10, 1, 16, 24, 32, 64, 96, 128, 256, 384, 512, 1024}; static const int OD_BAND_OFFSETS64[] = {13, 1, 16, 24, 32, 64, 96, 128, 256, 384, 512, 1024, 1536, 2048, 4096}; const int *const OD_BAND_OFFSETS[OD_TXSIZES + 1] = { OD_BAND_OFFSETS4, OD_BAND_OFFSETS8, OD_BAND_OFFSETS16, OD_BAND_OFFSETS32, OD_BAND_OFFSETS64 }; /** Perform a single stage of conversion from a coefficient block in * raster order into coding scan order * * @param [in] layout scan order specification * @param [out] dst destination vector * @param [in] src source coefficient block * @param [int] int source vector row stride */ static void od_band_from_raster(const band_layout *layout, tran_low_t *dst, const tran_low_t *src, int stride, TX_TYPE tx_type) { int i; int len; len = layout->band_offsets[layout->nb_bands]; for (i = 0; i < len; i++) { dst[i] = src[layout->dst_table[tx_type][i][1]*stride + layout->dst_table[tx_type][i][0]]; } } /** Perform a single stage of conversion from a vector in coding scan order back into a coefficient block in raster order * * @param [in] layout scan order specification * @param [out] dst destination coefficient block * @param [in] src source vector * @param [int] stride destination vector row stride */ static void od_raster_from_band(const band_layout *layout, tran_low_t *dst, int stride, TX_TYPE tx_type, const tran_low_t *src) { int i; int len; len = layout->band_offsets[layout->nb_bands]; for (i = 0; i < len; i++) { dst[layout->dst_table[tx_type][i][1]*stride + layout->dst_table[tx_type][i][0]] = src[i]; } } static const band_layout *const OD_LAYOUTS[] = {&OD_LAYOUT4, &OD_LAYOUT8, &OD_LAYOUT16, &OD_LAYOUT32}; /** Converts a coefficient block in raster order into a vector in * coding scan order with the PVQ partitions laid out one after * another. This works in stages; the 4x4 conversion is applied to * the coefficients nearest DC, then the 8x8 applied to the 8x8 block * nearest DC that was not already coded by 4x4, then 16x16 following * the same pattern. * * @param [out] dst destination vector * @param [in] n block size (along one side) * @param [in] ty_type transfrom type * @param [in] src source coefficient block * @param [in] stride source vector row stride */ void od_raster_to_coding_order(tran_low_t *dst, int n, TX_TYPE ty_type, const tran_low_t *src, int stride) { int bs; /* dst + 1 because DC is not included for 4x4 blocks. */ od_band_from_raster(OD_LAYOUTS[0], dst + 1, src, stride, ty_type); for (bs = 1; bs < OD_TXSIZES; bs++) { int size; int offset; /* Length of block size > 4. */ size = 1 << (OD_LOG_BSIZE0 + bs); /* Offset is the size of the previous block squared. */ offset = 1 << 2*(OD_LOG_BSIZE0 - 1 + bs); if (n >= size) { /* 3 16x16 bands come after 3 8x8 bands, which come after 2 4x4 bands. */ od_band_from_raster(OD_LAYOUTS[bs], dst + offset, src, stride, ty_type); } } dst[0] = src[0]; } /** Converts a vector in coding scan order witht he PVQ partitions * laid out one after another into a coefficient block in raster * order. This works in stages in the reverse order of raster->scan * order; the 16x16 conversion is applied to the coefficients that * don't appear in an 8x8 block, then the 8x8 applied to the 8x8 block * sans the 4x4 block it contains, then 4x4 is converted sans DC. * * @param [out] dst destination coefficient block * @param [in] stride destination vector row stride * @param [in] src source vector * @param [in] n block size (along one side) */ void od_coding_order_to_raster(tran_low_t *dst, int stride, TX_TYPE ty_type, const tran_low_t *src, int n) { int bs; /* src + 1 because DC is not included for 4x4 blocks. */ od_raster_from_band(OD_LAYOUTS[0], dst, stride, ty_type, src + 1); for (bs = 1; bs < OD_TXSIZES; bs++) { int size; int offset; /* Length of block size > 4 */ size = 1 << (OD_LOG_BSIZE0 + bs); /* Offset is the size of the previous block squared. */ offset = 1 << 2*(OD_LOG_BSIZE0 - 1 + bs); if (n >= size) { /* 3 16x16 bands come after 3 8x8 bands, which come after 2 4x4 bands. */ od_raster_from_band(OD_LAYOUTS[bs], dst, stride, ty_type, src + offset); } } dst[0] = src[0]; } /** Perform a single stage of conversion from a coefficient block in * raster order into coding scan order * * @param [in] layout scan order specification * @param [out] dst destination vector * @param [in] src source coefficient block * @param [int] int source vector row stride */ static void od_band_from_raster_16(const band_layout *layout, int16_t *dst, const int16_t *src, int stride) { int i; int len; len = layout->band_offsets[layout->nb_bands]; for (i = 0; i < len; i++) { dst[i] = src[layout->dst_table[DCT_DCT][i][1]*stride + layout->dst_table[DCT_DCT][i][0]]; } } /** Converts a coefficient block in raster order into a vector in * coding scan order with the PVQ partitions laid out one after * another. This works in stages; the 4x4 conversion is applied to * the coefficients nearest DC, then the 8x8 applied to the 8x8 block * nearest DC that was not already coded by 4x4, then 16x16 following * the same pattern. * * @param [out] dst destination vector * @param [in] n block size (along one side) * @param [in] src source coefficient block * @param [in] stride source vector row stride */ void od_raster_to_coding_order_16(int16_t *dst, int n, const int16_t *src, int stride) { int bs; /* dst + 1 because DC is not included for 4x4 blocks. */ od_band_from_raster_16(OD_LAYOUTS[0], dst + 1, src, stride); for (bs = 1; bs < OD_TXSIZES; bs++) { int size; int offset; /* Length of block size > 4. */ size = 1 << (OD_LOG_BSIZE0 + bs); /* Offset is the size of the previous block squared. */ offset = 1 << 2*(OD_LOG_BSIZE0 - 1 + bs); if (n >= size) { /* 3 16x16 bands come after 3 8x8 bands, which come after 2 4x4 bands. */ od_band_from_raster_16(OD_LAYOUTS[bs], dst + offset, src, stride); } } dst[0] = src[0]; }