Source code
Revision control
Copy as Markdown
Other Tools
/*
* Copyright (c) 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.
*/
#include <stddef.h>
#include <stdint.h>
#include "config/aom_config.h"
#include "config/aom_dsp_rtcd.h"
#include "aom_dsp/aom_dsp_common.h"
// The 2 unused parameters are place holders for PIC enabled build.
// These definitions are for functions defined in subpel_variance.asm
#define DECL(w, opt) \
int aom_sub_pixel_variance##w##xh_##opt( \
const uint8_t *src, ptrdiff_t src_stride, int x_offset, int y_offset, \
const uint8_t *dst, ptrdiff_t dst_stride, int height, unsigned int *sse, \
void *unused0, void *unused)
#define DECLS(opt) \
DECL(4, opt); \
DECL(8, opt); \
DECL(16, opt)
DECLS(ssse3);
#undef DECLS
#undef DECL
#define FN(w, h, wf, wlog2, hlog2, opt, cast_prod, cast) \
unsigned int aom_sub_pixel_variance##w##x##h##_##opt( \
const uint8_t *src, int src_stride, int x_offset, int y_offset, \
const uint8_t *dst, int dst_stride, unsigned int *sse_ptr) { \
/*Avoid overflow in helper by capping height.*/ \
const int hf = AOMMIN(h, 64); \
unsigned int sse = 0; \
int se = 0; \
for (int i = 0; i < (w / wf); ++i) { \
const uint8_t *src_ptr = src; \
const uint8_t *dst_ptr = dst; \
for (int j = 0; j < (h / hf); ++j) { \
unsigned int sse2; \
const int se2 = aom_sub_pixel_variance##wf##xh_##opt( \
src_ptr, src_stride, x_offset, y_offset, dst_ptr, dst_stride, hf, \
&sse2, NULL, NULL); \
dst_ptr += hf * dst_stride; \
src_ptr += hf * src_stride; \
se += se2; \
sse += sse2; \
} \
src += wf; \
dst += wf; \
} \
*sse_ptr = sse; \
return sse - (unsigned int)(cast_prod(cast se * se) >> (wlog2 + hlog2)); \
}
#if !CONFIG_REALTIME_ONLY
#define FNS(opt) \
FN(128, 128, 16, 7, 7, opt, (int64_t), (int64_t)) \
FN(128, 64, 16, 7, 6, opt, (int64_t), (int64_t)) \
FN(64, 128, 16, 6, 7, opt, (int64_t), (int64_t)) \
FN(64, 64, 16, 6, 6, opt, (int64_t), (int64_t)) \
FN(64, 32, 16, 6, 5, opt, (int64_t), (int64_t)) \
FN(32, 64, 16, 5, 6, opt, (int64_t), (int64_t)) \
FN(32, 32, 16, 5, 5, opt, (int64_t), (int64_t)) \
FN(32, 16, 16, 5, 4, opt, (int64_t), (int64_t)) \
FN(16, 32, 16, 4, 5, opt, (int64_t), (int64_t)) \
FN(16, 16, 16, 4, 4, opt, (uint32_t), (int64_t)) \
FN(16, 8, 16, 4, 3, opt, (int32_t), (int32_t)) \
FN(8, 16, 8, 3, 4, opt, (int32_t), (int32_t)) \
FN(8, 8, 8, 3, 3, opt, (int32_t), (int32_t)) \
FN(8, 4, 8, 3, 2, opt, (int32_t), (int32_t)) \
FN(4, 8, 4, 2, 3, opt, (int32_t), (int32_t)) \
FN(4, 4, 4, 2, 2, opt, (int32_t), (int32_t)) \
FN(4, 16, 4, 2, 4, opt, (int32_t), (int32_t)) \
FN(16, 4, 16, 4, 2, opt, (int32_t), (int32_t)) \
FN(8, 32, 8, 3, 5, opt, (uint32_t), (int64_t)) \
FN(32, 8, 16, 5, 3, opt, (uint32_t), (int64_t)) \
FN(16, 64, 16, 4, 6, opt, (int64_t), (int64_t)) \
FN(64, 16, 16, 6, 4, opt, (int64_t), (int64_t))
#else
#define FNS(opt) \
FN(128, 128, 16, 7, 7, opt, (int64_t), (int64_t)) \
FN(128, 64, 16, 7, 6, opt, (int64_t), (int64_t)) \
FN(64, 128, 16, 6, 7, opt, (int64_t), (int64_t)) \
FN(64, 64, 16, 6, 6, opt, (int64_t), (int64_t)) \
FN(64, 32, 16, 6, 5, opt, (int64_t), (int64_t)) \
FN(32, 64, 16, 5, 6, opt, (int64_t), (int64_t)) \
FN(32, 32, 16, 5, 5, opt, (int64_t), (int64_t)) \
FN(32, 16, 16, 5, 4, opt, (int64_t), (int64_t)) \
FN(16, 32, 16, 4, 5, opt, (int64_t), (int64_t)) \
FN(16, 16, 16, 4, 4, opt, (uint32_t), (int64_t)) \
FN(16, 8, 16, 4, 3, opt, (int32_t), (int32_t)) \
FN(8, 16, 8, 3, 4, opt, (int32_t), (int32_t)) \
FN(8, 8, 8, 3, 3, opt, (int32_t), (int32_t)) \
FN(8, 4, 8, 3, 2, opt, (int32_t), (int32_t)) \
FN(4, 8, 4, 2, 3, opt, (int32_t), (int32_t)) \
FN(4, 4, 4, 2, 2, opt, (int32_t), (int32_t))
#endif
FNS(ssse3)
#undef FNS
#undef FN
// The 2 unused parameters are place holders for PIC enabled build.
#define DECL(w, opt) \
int aom_sub_pixel_avg_variance##w##xh_##opt( \
const uint8_t *src, ptrdiff_t src_stride, int x_offset, int y_offset, \
const uint8_t *dst, ptrdiff_t dst_stride, const uint8_t *sec, \
ptrdiff_t sec_stride, int height, unsigned int *sse, void *unused0, \
void *unused)
#define DECLS(opt) \
DECL(4, opt); \
DECL(8, opt); \
DECL(16, opt)
DECLS(ssse3);
#undef DECL
#undef DECLS
#define FN(w, h, wf, wlog2, hlog2, opt, cast_prod, cast) \
unsigned int aom_sub_pixel_avg_variance##w##x##h##_##opt( \
const uint8_t *src, int src_stride, int x_offset, int y_offset, \
const uint8_t *dst, int dst_stride, unsigned int *sse_ptr, \
const uint8_t *sec) { \
/*Avoid overflow in helper by capping height.*/ \
const int hf = AOMMIN(h, 64); \
unsigned int sse = 0; \
int se = 0; \
for (int i = 0; i < (w / wf); ++i) { \
const uint8_t *src_ptr = src; \
const uint8_t *dst_ptr = dst; \
const uint8_t *sec_ptr = sec; \
for (int j = 0; j < (h / hf); ++j) { \
unsigned int sse2; \
const int se2 = aom_sub_pixel_avg_variance##wf##xh_##opt( \
src_ptr, src_stride, x_offset, y_offset, dst_ptr, dst_stride, \
sec_ptr, w, hf, &sse2, NULL, NULL); \
dst_ptr += hf * dst_stride; \
src_ptr += hf * src_stride; \
sec_ptr += hf * w; \
se += se2; \
sse += sse2; \
} \
src += wf; \
dst += wf; \
sec += wf; \
} \
*sse_ptr = sse; \
return sse - (unsigned int)(cast_prod(cast se * se) >> (wlog2 + hlog2)); \
}
#if !CONFIG_REALTIME_ONLY
#define FNS(opt) \
FN(128, 128, 16, 7, 7, opt, (int64_t), (int64_t)) \
FN(128, 64, 16, 7, 6, opt, (int64_t), (int64_t)) \
FN(64, 128, 16, 6, 7, opt, (int64_t), (int64_t)) \
FN(64, 64, 16, 6, 6, opt, (int64_t), (int64_t)) \
FN(64, 32, 16, 6, 5, opt, (int64_t), (int64_t)) \
FN(32, 64, 16, 5, 6, opt, (int64_t), (int64_t)) \
FN(32, 32, 16, 5, 5, opt, (int64_t), (int64_t)) \
FN(32, 16, 16, 5, 4, opt, (int64_t), (int64_t)) \
FN(16, 32, 16, 4, 5, opt, (int64_t), (int64_t)) \
FN(16, 16, 16, 4, 4, opt, (uint32_t), (int64_t)) \
FN(16, 8, 16, 4, 3, opt, (uint32_t), (int32_t)) \
FN(8, 16, 8, 3, 4, opt, (uint32_t), (int32_t)) \
FN(8, 8, 8, 3, 3, opt, (uint32_t), (int32_t)) \
FN(8, 4, 8, 3, 2, opt, (uint32_t), (int32_t)) \
FN(4, 8, 4, 2, 3, opt, (uint32_t), (int32_t)) \
FN(4, 4, 4, 2, 2, opt, (uint32_t), (int32_t)) \
FN(4, 16, 4, 2, 4, opt, (int32_t), (int32_t)) \
FN(16, 4, 16, 4, 2, opt, (int32_t), (int32_t)) \
FN(8, 32, 8, 3, 5, opt, (uint32_t), (int64_t)) \
FN(32, 8, 16, 5, 3, opt, (uint32_t), (int64_t)) \
FN(16, 64, 16, 4, 6, opt, (int64_t), (int64_t)) \
FN(64, 16, 16, 6, 4, opt, (int64_t), (int64_t))
#else
#define FNS(opt) \
FN(128, 128, 16, 7, 7, opt, (int64_t), (int64_t)) \
FN(128, 64, 16, 7, 6, opt, (int64_t), (int64_t)) \
FN(64, 128, 16, 6, 7, opt, (int64_t), (int64_t)) \
FN(64, 64, 16, 6, 6, opt, (int64_t), (int64_t)) \
FN(64, 32, 16, 6, 5, opt, (int64_t), (int64_t)) \
FN(32, 64, 16, 5, 6, opt, (int64_t), (int64_t)) \
FN(32, 32, 16, 5, 5, opt, (int64_t), (int64_t)) \
FN(32, 16, 16, 5, 4, opt, (int64_t), (int64_t)) \
FN(16, 32, 16, 4, 5, opt, (int64_t), (int64_t)) \
FN(16, 16, 16, 4, 4, opt, (uint32_t), (int64_t)) \
FN(16, 8, 16, 4, 3, opt, (uint32_t), (int32_t)) \
FN(8, 16, 8, 3, 4, opt, (uint32_t), (int32_t)) \
FN(8, 8, 8, 3, 3, opt, (uint32_t), (int32_t)) \
FN(8, 4, 8, 3, 2, opt, (uint32_t), (int32_t)) \
FN(4, 8, 4, 2, 3, opt, (uint32_t), (int32_t)) \
FN(4, 4, 4, 2, 2, opt, (uint32_t), (int32_t))
#endif
FNS(ssse3)
#undef FNS
#undef FN