|
|
|
/*
|
|
|
|
* RV40 decoder motion compensation functions x86-optimised
|
|
|
|
* Copyright (c) 2008 Konstantin Shishkov
|
|
|
|
*
|
|
|
|
* This file is part of FFmpeg.
|
|
|
|
*
|
|
|
|
* FFmpeg is free software; you can redistribute it and/or
|
|
|
|
* modify it under the terms of the GNU Lesser General Public
|
|
|
|
* License as published by the Free Software Foundation; either
|
|
|
|
* version 2.1 of the License, or (at your option) any later version.
|
|
|
|
*
|
|
|
|
* FFmpeg is distributed in the hope that it will be useful,
|
|
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
|
|
* Lesser General Public License for more details.
|
|
|
|
*
|
|
|
|
* You should have received a copy of the GNU Lesser General Public
|
|
|
|
* License along with FFmpeg; if not, write to the Free Software
|
|
|
|
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
|
|
|
|
*/
|
|
|
|
|
|
|
|
/**
|
|
|
|
* @file
|
|
|
|
* RV40 decoder motion compensation functions x86-optimised
|
|
|
|
* 2,0 and 0,2 have h264 equivalents.
|
|
|
|
* 3,3 is bugged in the rv40 format and maps to _xy2 version
|
|
|
|
*/
|
|
|
|
|
|
|
|
#include "libavcodec/rv34dsp.h"
|
|
|
|
#include "libavutil/attributes.h"
|
|
|
|
#include "libavutil/mem.h"
|
|
|
|
#include "libavutil/x86/cpu.h"
|
|
|
|
#include "hpeldsp.h"
|
|
|
|
|
|
|
|
#define DEFINE_FN(op, size, insn) \
|
|
|
|
static void op##_rv40_qpel##size##_mc33_##insn(uint8_t *dst, const uint8_t *src, \
|
|
|
|
ptrdiff_t stride) \
|
|
|
|
{ \
|
|
|
|
ff_##op##_pixels##size##_xy2_##insn(dst, src, stride, size); \
|
|
|
|
}
|
|
|
|
|
|
|
|
#if HAVE_YASM
|
|
|
|
void ff_put_rv40_chroma_mc8_mmx (uint8_t *dst, uint8_t *src,
|
|
|
|
int stride, int h, int x, int y);
|
|
|
|
void ff_avg_rv40_chroma_mc8_mmxext(uint8_t *dst, uint8_t *src,
|
|
|
|
int stride, int h, int x, int y);
|
|
|
|
void ff_avg_rv40_chroma_mc8_3dnow(uint8_t *dst, uint8_t *src,
|
|
|
|
int stride, int h, int x, int y);
|
|
|
|
|
|
|
|
void ff_put_rv40_chroma_mc4_mmx (uint8_t *dst, uint8_t *src,
|
|
|
|
int stride, int h, int x, int y);
|
|
|
|
void ff_avg_rv40_chroma_mc4_mmxext(uint8_t *dst, uint8_t *src,
|
|
|
|
int stride, int h, int x, int y);
|
|
|
|
void ff_avg_rv40_chroma_mc4_3dnow(uint8_t *dst, uint8_t *src,
|
|
|
|
int stride, int h, int x, int y);
|
|
|
|
|
rv40: x86 SIMD for biweight
Provide MMX, SSE2 and SSSE3 versions, with a fast-path when the weights are
multiples of 512 (which is often the case when the values round up nicely).
*_TIMER report for the 16x16 and 8x8 cases:
C:
9015 decicycles in 16, 524257 runs, 31 skips
2656 decicycles in 8, 524271 runs, 17 skips
MMX:
4156 decicycles in 16, 262090 runs, 54 skips
1206 decicycles in 8, 262131 runs, 13 skips
MMX on fast-path:
2760 decicycles in 16, 524222 runs, 66 skips
995 decicycles in 8, 524252 runs, 36 skips
SSE2:
2163 decicycles in 16, 262131 runs, 13 skips
832 decicycles in 8, 262137 runs, 7 skips
SSE2 with fast path:
1783 decicycles in 16, 524276 runs, 12 skips
711 decicycles in 8, 524283 runs, 5 skips
SSSE3:
2117 decicycles in 16, 262136 runs, 8 skips
814 decicycles in 8, 262143 runs, 1 skips
SSSE3 with fast path:
1315 decicycles in 16, 524285 runs, 3 skips
578 decicycles in 8, 524286 runs, 2 skips
This means around a 4% speedup for some sequences.
Signed-off-by: Diego Biurrun <diego@biurrun.de>
13 years ago
|
|
|
#define DECLARE_WEIGHT(opt) \
|
|
|
|
void ff_rv40_weight_func_rnd_16_##opt(uint8_t *dst, uint8_t *src1, uint8_t *src2, \
|
|
|
|
int w1, int w2, ptrdiff_t stride); \
|
|
|
|
void ff_rv40_weight_func_rnd_8_##opt (uint8_t *dst, uint8_t *src1, uint8_t *src2, \
|
|
|
|
int w1, int w2, ptrdiff_t stride); \
|
|
|
|
void ff_rv40_weight_func_nornd_16_##opt(uint8_t *dst, uint8_t *src1, uint8_t *src2, \
|
|
|
|
int w1, int w2, ptrdiff_t stride); \
|
|
|
|
void ff_rv40_weight_func_nornd_8_##opt (uint8_t *dst, uint8_t *src1, uint8_t *src2, \
|
|
|
|
int w1, int w2, ptrdiff_t stride);
|
|
|
|
DECLARE_WEIGHT(mmxext)
|
rv40: x86 SIMD for biweight
Provide MMX, SSE2 and SSSE3 versions, with a fast-path when the weights are
multiples of 512 (which is often the case when the values round up nicely).
*_TIMER report for the 16x16 and 8x8 cases:
C:
9015 decicycles in 16, 524257 runs, 31 skips
2656 decicycles in 8, 524271 runs, 17 skips
MMX:
4156 decicycles in 16, 262090 runs, 54 skips
1206 decicycles in 8, 262131 runs, 13 skips
MMX on fast-path:
2760 decicycles in 16, 524222 runs, 66 skips
995 decicycles in 8, 524252 runs, 36 skips
SSE2:
2163 decicycles in 16, 262131 runs, 13 skips
832 decicycles in 8, 262137 runs, 7 skips
SSE2 with fast path:
1783 decicycles in 16, 524276 runs, 12 skips
711 decicycles in 8, 524283 runs, 5 skips
SSSE3:
2117 decicycles in 16, 262136 runs, 8 skips
814 decicycles in 8, 262143 runs, 1 skips
SSSE3 with fast path:
1315 decicycles in 16, 524285 runs, 3 skips
578 decicycles in 8, 524286 runs, 2 skips
This means around a 4% speedup for some sequences.
Signed-off-by: Diego Biurrun <diego@biurrun.de>
13 years ago
|
|
|
DECLARE_WEIGHT(sse2)
|
|
|
|
DECLARE_WEIGHT(ssse3)
|
|
|
|
|
|
|
|
/** @{ */
|
|
|
|
/**
|
|
|
|
* Define one qpel function.
|
|
|
|
* LOOPSIZE must be already set to the number of pixels processed per
|
|
|
|
* iteration in the inner loop of the called functions.
|
|
|
|
* COFF(x) must be already defined so as to provide the offset into any
|
|
|
|
* array of coeffs used by the called function for the qpel position x.
|
|
|
|
*/
|
|
|
|
#define QPEL_FUNC_DECL(OP, SIZE, PH, PV, OPT) \
|
|
|
|
static void OP ## rv40_qpel ##SIZE ##_mc ##PH ##PV ##OPT(uint8_t *dst, \
|
|
|
|
const uint8_t *src, \
|
|
|
|
ptrdiff_t stride) \
|
|
|
|
{ \
|
|
|
|
int i; \
|
|
|
|
if (PH && PV) { \
|
|
|
|
LOCAL_ALIGNED(16, uint8_t, tmp, [SIZE * (SIZE + 5)]); \
|
|
|
|
uint8_t *tmpptr = tmp + SIZE * 2; \
|
|
|
|
src -= stride * 2; \
|
|
|
|
\
|
|
|
|
for (i = 0; i < SIZE; i += LOOPSIZE) \
|
|
|
|
ff_put_rv40_qpel_h ##OPT(tmp + i, SIZE, src + i, stride, \
|
|
|
|
SIZE + 5, HCOFF(PH)); \
|
|
|
|
for (i = 0; i < SIZE; i += LOOPSIZE) \
|
|
|
|
ff_ ##OP ##rv40_qpel_v ##OPT(dst + i, stride, tmpptr + i, \
|
|
|
|
SIZE, SIZE, VCOFF(PV)); \
|
|
|
|
} else if (PV) { \
|
|
|
|
for (i = 0; i < SIZE; i += LOOPSIZE) \
|
|
|
|
ff_ ##OP ##rv40_qpel_v ## OPT(dst + i, stride, src + i, \
|
|
|
|
stride, SIZE, VCOFF(PV)); \
|
|
|
|
} else { \
|
|
|
|
for (i = 0; i < SIZE; i += LOOPSIZE) \
|
|
|
|
ff_ ##OP ##rv40_qpel_h ## OPT(dst + i, stride, src + i, \
|
|
|
|
stride, SIZE, HCOFF(PH)); \
|
|
|
|
} \
|
|
|
|
}
|
|
|
|
|
|
|
|
/** Declare functions for sizes 8 and 16 and given operations
|
|
|
|
* and qpel position. */
|
|
|
|
#define QPEL_FUNCS_DECL(OP, PH, PV, OPT) \
|
|
|
|
QPEL_FUNC_DECL(OP, 8, PH, PV, OPT) \
|
|
|
|
QPEL_FUNC_DECL(OP, 16, PH, PV, OPT)
|
|
|
|
|
|
|
|
/** Declare all functions for all sizes and qpel positions */
|
|
|
|
#define QPEL_MC_DECL(OP, OPT) \
|
|
|
|
void ff_ ##OP ##rv40_qpel_h ##OPT(uint8_t *dst, ptrdiff_t dstStride, \
|
|
|
|
const uint8_t *src, \
|
|
|
|
ptrdiff_t srcStride, \
|
|
|
|
int len, int m); \
|
|
|
|
void ff_ ##OP ##rv40_qpel_v ##OPT(uint8_t *dst, ptrdiff_t dstStride, \
|
|
|
|
const uint8_t *src, \
|
|
|
|
ptrdiff_t srcStride, \
|
|
|
|
int len, int m); \
|
|
|
|
QPEL_FUNCS_DECL(OP, 0, 1, OPT) \
|
|
|
|
QPEL_FUNCS_DECL(OP, 0, 3, OPT) \
|
|
|
|
QPEL_FUNCS_DECL(OP, 1, 0, OPT) \
|
|
|
|
QPEL_FUNCS_DECL(OP, 1, 1, OPT) \
|
|
|
|
QPEL_FUNCS_DECL(OP, 1, 2, OPT) \
|
|
|
|
QPEL_FUNCS_DECL(OP, 1, 3, OPT) \
|
|
|
|
QPEL_FUNCS_DECL(OP, 2, 1, OPT) \
|
|
|
|
QPEL_FUNCS_DECL(OP, 2, 2, OPT) \
|
|
|
|
QPEL_FUNCS_DECL(OP, 2, 3, OPT) \
|
|
|
|
QPEL_FUNCS_DECL(OP, 3, 0, OPT) \
|
|
|
|
QPEL_FUNCS_DECL(OP, 3, 1, OPT) \
|
|
|
|
QPEL_FUNCS_DECL(OP, 3, 2, OPT)
|
|
|
|
/** @} */
|
|
|
|
|
|
|
|
#define LOOPSIZE 8
|
|
|
|
#define HCOFF(x) (32 * ((x) - 1))
|
|
|
|
#define VCOFF(x) (32 * ((x) - 1))
|
|
|
|
QPEL_MC_DECL(put_, _ssse3)
|
|
|
|
QPEL_MC_DECL(avg_, _ssse3)
|
|
|
|
|
|
|
|
#undef LOOPSIZE
|
|
|
|
#undef HCOFF
|
|
|
|
#undef VCOFF
|
|
|
|
#define LOOPSIZE 8
|
|
|
|
#define HCOFF(x) (64 * ((x) - 1))
|
|
|
|
#define VCOFF(x) (64 * ((x) - 1))
|
|
|
|
QPEL_MC_DECL(put_, _sse2)
|
|
|
|
QPEL_MC_DECL(avg_, _sse2)
|
|
|
|
|
|
|
|
#if ARCH_X86_32
|
|
|
|
#undef LOOPSIZE
|
|
|
|
#undef HCOFF
|
|
|
|
#undef VCOFF
|
|
|
|
#define LOOPSIZE 4
|
|
|
|
#define HCOFF(x) (64 * ((x) - 1))
|
|
|
|
#define VCOFF(x) (64 * ((x) - 1))
|
|
|
|
|
|
|
|
QPEL_MC_DECL(put_, _mmx)
|
|
|
|
|
|
|
|
#define ff_put_rv40_qpel_h_mmxext ff_put_rv40_qpel_h_mmx
|
|
|
|
#define ff_put_rv40_qpel_v_mmxext ff_put_rv40_qpel_v_mmx
|
|
|
|
QPEL_MC_DECL(avg_, _mmxext)
|
|
|
|
|
|
|
|
#define ff_put_rv40_qpel_h_3dnow ff_put_rv40_qpel_h_mmx
|
|
|
|
#define ff_put_rv40_qpel_v_3dnow ff_put_rv40_qpel_v_mmx
|
|
|
|
QPEL_MC_DECL(avg_, _3dnow)
|
|
|
|
#endif
|
|
|
|
|
|
|
|
/** @{ */
|
|
|
|
/** Set one function */
|
|
|
|
#define QPEL_FUNC_SET(OP, SIZE, PH, PV, OPT) \
|
|
|
|
c-> OP ## pixels_tab[2 - SIZE / 8][4 * PV + PH] = OP ## rv40_qpel ##SIZE ## _mc ##PH ##PV ##OPT;
|
|
|
|
|
|
|
|
/** Set functions put and avg for sizes 8 and 16 and a given qpel position */
|
|
|
|
#define QPEL_FUNCS_SET(OP, PH, PV, OPT) \
|
|
|
|
QPEL_FUNC_SET(OP, 8, PH, PV, OPT) \
|
|
|
|
QPEL_FUNC_SET(OP, 16, PH, PV, OPT)
|
|
|
|
|
|
|
|
/** Set all functions for all sizes and qpel positions */
|
|
|
|
#define QPEL_MC_SET(OP, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 0, 1, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 0, 3, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 1, 0, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 1, 1, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 1, 2, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 1, 3, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 2, 1, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 2, 2, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 2, 3, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 3, 0, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 3, 1, OPT) \
|
|
|
|
QPEL_FUNCS_SET (OP, 3, 2, OPT)
|
|
|
|
/** @} */
|
|
|
|
|
|
|
|
DEFINE_FN(put, 8, ssse3)
|
|
|
|
|
|
|
|
DEFINE_FN(put, 16, sse2)
|
|
|
|
DEFINE_FN(put, 16, ssse3)
|
|
|
|
|
|
|
|
DEFINE_FN(avg, 8, mmxext)
|
|
|
|
DEFINE_FN(avg, 8, ssse3)
|
|
|
|
|
|
|
|
DEFINE_FN(avg, 16, sse2)
|
|
|
|
DEFINE_FN(avg, 16, ssse3)
|
|
|
|
#endif /* HAVE_YASM */
|
|
|
|
|
|
|
|
#if HAVE_MMX_INLINE
|
|
|
|
DEFINE_FN(put, 8, mmx)
|
|
|
|
DEFINE_FN(avg, 8, mmx)
|
|
|
|
DEFINE_FN(put, 16, mmx)
|
|
|
|
DEFINE_FN(avg, 16, mmx)
|
|
|
|
#endif
|
|
|
|
|
|
|
|
av_cold void ff_rv40dsp_init_x86(RV34DSPContext *c)
|
|
|
|
{
|
|
|
|
av_unused int cpu_flags = av_get_cpu_flags();
|
|
|
|
|
|
|
|
#if HAVE_MMX_INLINE
|
|
|
|
if (INLINE_MMX(cpu_flags)) {
|
|
|
|
c->put_pixels_tab[0][15] = put_rv40_qpel16_mc33_mmx;
|
|
|
|
c->put_pixels_tab[1][15] = put_rv40_qpel8_mc33_mmx;
|
|
|
|
c->avg_pixels_tab[0][15] = avg_rv40_qpel16_mc33_mmx;
|
|
|
|
c->avg_pixels_tab[1][15] = avg_rv40_qpel8_mc33_mmx;
|
|
|
|
}
|
|
|
|
#endif /* HAVE_MMX_INLINE */
|
|
|
|
|
|
|
|
#if HAVE_YASM
|
|
|
|
if (EXTERNAL_MMX(cpu_flags)) {
|
|
|
|
c->put_chroma_pixels_tab[0] = ff_put_rv40_chroma_mc8_mmx;
|
|
|
|
c->put_chroma_pixels_tab[1] = ff_put_rv40_chroma_mc4_mmx;
|
|
|
|
#if ARCH_X86_32
|
|
|
|
QPEL_MC_SET(put_, _mmx)
|
|
|
|
#endif
|
|
|
|
}
|
|
|
|
if (EXTERNAL_AMD3DNOW(cpu_flags)) {
|
|
|
|
c->avg_chroma_pixels_tab[0] = ff_avg_rv40_chroma_mc8_3dnow;
|
|
|
|
c->avg_chroma_pixels_tab[1] = ff_avg_rv40_chroma_mc4_3dnow;
|
|
|
|
#if ARCH_X86_32
|
|
|
|
QPEL_MC_SET(avg_, _3dnow)
|
|
|
|
#endif
|
|
|
|
}
|
|
|
|
if (EXTERNAL_MMXEXT(cpu_flags)) {
|
|
|
|
c->avg_pixels_tab[1][15] = avg_rv40_qpel8_mc33_mmxext;
|
|
|
|
c->avg_chroma_pixels_tab[0] = ff_avg_rv40_chroma_mc8_mmxext;
|
|
|
|
c->avg_chroma_pixels_tab[1] = ff_avg_rv40_chroma_mc4_mmxext;
|
|
|
|
c->rv40_weight_pixels_tab[0][0] = ff_rv40_weight_func_rnd_16_mmxext;
|
|
|
|
c->rv40_weight_pixels_tab[0][1] = ff_rv40_weight_func_rnd_8_mmxext;
|
|
|
|
c->rv40_weight_pixels_tab[1][0] = ff_rv40_weight_func_nornd_16_mmxext;
|
|
|
|
c->rv40_weight_pixels_tab[1][1] = ff_rv40_weight_func_nornd_8_mmxext;
|
|
|
|
#if ARCH_X86_32
|
|
|
|
QPEL_MC_SET(avg_, _mmxext)
|
|
|
|
#endif
|
|
|
|
}
|
|
|
|
if (EXTERNAL_SSE2(cpu_flags)) {
|
|
|
|
c->put_pixels_tab[0][15] = put_rv40_qpel16_mc33_sse2;
|
|
|
|
c->avg_pixels_tab[0][15] = avg_rv40_qpel16_mc33_sse2;
|
|
|
|
c->rv40_weight_pixels_tab[0][0] = ff_rv40_weight_func_rnd_16_sse2;
|
|
|
|
c->rv40_weight_pixels_tab[0][1] = ff_rv40_weight_func_rnd_8_sse2;
|
|
|
|
c->rv40_weight_pixels_tab[1][0] = ff_rv40_weight_func_nornd_16_sse2;
|
|
|
|
c->rv40_weight_pixels_tab[1][1] = ff_rv40_weight_func_nornd_8_sse2;
|
|
|
|
QPEL_MC_SET(put_, _sse2)
|
|
|
|
QPEL_MC_SET(avg_, _sse2)
|
rv40: x86 SIMD for biweight
Provide MMX, SSE2 and SSSE3 versions, with a fast-path when the weights are
multiples of 512 (which is often the case when the values round up nicely).
*_TIMER report for the 16x16 and 8x8 cases:
C:
9015 decicycles in 16, 524257 runs, 31 skips
2656 decicycles in 8, 524271 runs, 17 skips
MMX:
4156 decicycles in 16, 262090 runs, 54 skips
1206 decicycles in 8, 262131 runs, 13 skips
MMX on fast-path:
2760 decicycles in 16, 524222 runs, 66 skips
995 decicycles in 8, 524252 runs, 36 skips
SSE2:
2163 decicycles in 16, 262131 runs, 13 skips
832 decicycles in 8, 262137 runs, 7 skips
SSE2 with fast path:
1783 decicycles in 16, 524276 runs, 12 skips
711 decicycles in 8, 524283 runs, 5 skips
SSSE3:
2117 decicycles in 16, 262136 runs, 8 skips
814 decicycles in 8, 262143 runs, 1 skips
SSSE3 with fast path:
1315 decicycles in 16, 524285 runs, 3 skips
578 decicycles in 8, 524286 runs, 2 skips
This means around a 4% speedup for some sequences.
Signed-off-by: Diego Biurrun <diego@biurrun.de>
13 years ago
|
|
|
}
|
|
|
|
if (EXTERNAL_SSSE3(cpu_flags)) {
|
|
|
|
c->put_pixels_tab[0][15] = put_rv40_qpel16_mc33_ssse3;
|
|
|
|
c->put_pixels_tab[1][15] = put_rv40_qpel8_mc33_ssse3;
|
|
|
|
c->avg_pixels_tab[0][15] = avg_rv40_qpel16_mc33_ssse3;
|
|
|
|
c->avg_pixels_tab[1][15] = avg_rv40_qpel8_mc33_ssse3;
|
|
|
|
c->rv40_weight_pixels_tab[0][0] = ff_rv40_weight_func_rnd_16_ssse3;
|
|
|
|
c->rv40_weight_pixels_tab[0][1] = ff_rv40_weight_func_rnd_8_ssse3;
|
|
|
|
c->rv40_weight_pixels_tab[1][0] = ff_rv40_weight_func_nornd_16_ssse3;
|
|
|
|
c->rv40_weight_pixels_tab[1][1] = ff_rv40_weight_func_nornd_8_ssse3;
|
|
|
|
QPEL_MC_SET(put_, _ssse3)
|
|
|
|
QPEL_MC_SET(avg_, _ssse3)
|
rv40: x86 SIMD for biweight
Provide MMX, SSE2 and SSSE3 versions, with a fast-path when the weights are
multiples of 512 (which is often the case when the values round up nicely).
*_TIMER report for the 16x16 and 8x8 cases:
C:
9015 decicycles in 16, 524257 runs, 31 skips
2656 decicycles in 8, 524271 runs, 17 skips
MMX:
4156 decicycles in 16, 262090 runs, 54 skips
1206 decicycles in 8, 262131 runs, 13 skips
MMX on fast-path:
2760 decicycles in 16, 524222 runs, 66 skips
995 decicycles in 8, 524252 runs, 36 skips
SSE2:
2163 decicycles in 16, 262131 runs, 13 skips
832 decicycles in 8, 262137 runs, 7 skips
SSE2 with fast path:
1783 decicycles in 16, 524276 runs, 12 skips
711 decicycles in 8, 524283 runs, 5 skips
SSSE3:
2117 decicycles in 16, 262136 runs, 8 skips
814 decicycles in 8, 262143 runs, 1 skips
SSSE3 with fast path:
1315 decicycles in 16, 524285 runs, 3 skips
578 decicycles in 8, 524286 runs, 2 skips
This means around a 4% speedup for some sequences.
Signed-off-by: Diego Biurrun <diego@biurrun.de>
13 years ago
|
|
|
}
|
|
|
|
#endif /* HAVE_YASM */
|
|
|
|
}
|