blob: 634a9edd59fcfe28ecf0eb5524471fd4acf96f25 [file] [log] [blame]
/*
* 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];
}