247 lines
8.6 KiB
C
247 lines
8.6 KiB
C
/*
|
|
* Copyright (c) 2004 Romain Dolbeau <romain@dolbeau.org>
|
|
*
|
|
* 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
|
|
*/
|
|
|
|
#include "libavutil/mem.h"
|
|
#include "libavutil/ppc/util_altivec.h"
|
|
|
|
/* this code assume that stride % 16 == 0 */
|
|
|
|
#define CHROMA_MC8_ALTIVEC_CORE(BIAS1, BIAS2) \
|
|
vsrc2ssH = (vec_s16)VEC_MERGEH(zero_u8v,(vec_u8)vsrc2uc);\
|
|
vsrc3ssH = (vec_s16)VEC_MERGEH(zero_u8v,(vec_u8)vsrc3uc);\
|
|
\
|
|
psum = vec_mladd(vA, vsrc0ssH, BIAS1);\
|
|
psum = vec_mladd(vB, vsrc1ssH, psum);\
|
|
psum = vec_mladd(vC, vsrc2ssH, psum);\
|
|
psum = vec_mladd(vD, vsrc3ssH, psum);\
|
|
psum = BIAS2(psum);\
|
|
psum = vec_sr(psum, v6us);\
|
|
\
|
|
vdst = vec_ld(0, dst);\
|
|
ppsum = (vec_u8)vec_pack(psum, psum);\
|
|
vfdst = vec_perm(vdst, ppsum, fperm);\
|
|
\
|
|
OP_U8_ALTIVEC(fsum, vfdst, vdst);\
|
|
\
|
|
vec_st(fsum, 0, dst);\
|
|
\
|
|
vsrc0ssH = vsrc2ssH;\
|
|
vsrc1ssH = vsrc3ssH;\
|
|
\
|
|
dst += stride;\
|
|
src += stride;
|
|
|
|
#define CHROMA_MC8_ALTIVEC_CORE_SIMPLE \
|
|
\
|
|
vsrc0ssH = (vec_s16)VEC_MERGEH(zero_u8v,(vec_u8)vsrc0uc);\
|
|
vsrc1ssH = (vec_s16)VEC_MERGEH(zero_u8v,(vec_u8)vsrc1uc);\
|
|
\
|
|
psum = vec_mladd(vA, vsrc0ssH, v32ss);\
|
|
psum = vec_mladd(vE, vsrc1ssH, psum);\
|
|
psum = vec_sr(psum, v6us);\
|
|
\
|
|
vdst = vec_ld(0, dst);\
|
|
ppsum = (vec_u8)vec_pack(psum, psum);\
|
|
vfdst = vec_perm(vdst, ppsum, fperm);\
|
|
\
|
|
OP_U8_ALTIVEC(fsum, vfdst, vdst);\
|
|
\
|
|
vec_st(fsum, 0, dst);\
|
|
\
|
|
dst += stride;\
|
|
src += stride;
|
|
|
|
#define noop(a) a
|
|
#define add28(a) vec_add(v28ss, a)
|
|
|
|
#if HAVE_BIGENDIAN
|
|
#define GET_VSRC1(vs0, off, b, perm0, s){ \
|
|
vec_u8 vsrcCuc, vsrcDuc; \
|
|
vsrcCuc = vec_ld(off, s); \
|
|
if (loadSecond){ \
|
|
vsrcDuc = vec_ld(off + b, s); \
|
|
} else \
|
|
vsrcDuc = vsrcCuc; \
|
|
\
|
|
vs0 = vec_perm(vsrcCuc, vsrcDuc, perm0); \
|
|
}
|
|
#define GET_VSRC(vs0, vs1, off, b, perm0, perm1, s){ \
|
|
vec_u8 vsrcCuc, vsrcDuc; \
|
|
vsrcCuc = vec_ld(off, s); \
|
|
if (loadSecond){ \
|
|
vsrcDuc = vec_ld(off + b, s); \
|
|
} else \
|
|
vsrcDuc = vsrcCuc; \
|
|
\
|
|
vs0 = vec_perm(vsrcCuc, vsrcDuc, perm0); \
|
|
if (reallyBadAlign){ \
|
|
vs1 = vsrcDuc; \
|
|
} else \
|
|
vs1 = vec_perm(vsrcCuc, vsrcDuc, perm1); \
|
|
}
|
|
|
|
#else
|
|
|
|
#define GET_VSRC1(vs0, off, b, perm0, s){ \
|
|
vs0 = vec_vsx_ld(off, s); \
|
|
}
|
|
#define GET_VSRC(vs0, vs1, off, b, perm0, perm1, s){ \
|
|
vs0 = vec_vsx_ld(off, s); \
|
|
vs1 = vec_vsx_ld(off + 1, s); \
|
|
}
|
|
#endif /* HAVE_BIGENDIAN */
|
|
|
|
#ifdef PREFIX_h264_chroma_mc8_altivec
|
|
static void PREFIX_h264_chroma_mc8_altivec(uint8_t * dst, uint8_t * src,
|
|
ptrdiff_t stride, int h,
|
|
int x, int y)
|
|
{
|
|
DECLARE_ALIGNED(16, signed int, ABCD)[4] =
|
|
{((8 - x) * (8 - y)),
|
|
(( x) * (8 - y)),
|
|
((8 - x) * ( y)),
|
|
(( x) * ( y))};
|
|
register int i;
|
|
vec_u8 fperm;
|
|
LOAD_ZERO;
|
|
const vec_s32 vABCD = vec_ld(0, ABCD);
|
|
const vec_s16 vA = VEC_SPLAT16(vABCD, 1);
|
|
const vec_s16 vB = VEC_SPLAT16(vABCD, 3);
|
|
const vec_s16 vC = VEC_SPLAT16(vABCD, 5);
|
|
const vec_s16 vD = VEC_SPLAT16(vABCD, 7);
|
|
const vec_s16 v32ss = vec_sl(vec_splat_s16(1),vec_splat_u16(5));
|
|
const vec_u16 v6us = vec_splat_u16(6);
|
|
|
|
vec_u8 vsrcperm0, vsrcperm1;
|
|
vec_u8 vsrc0uc, vsrc1uc;
|
|
vec_s16 vsrc0ssH, vsrc1ssH;
|
|
vec_u8 vsrc2uc, vsrc3uc;
|
|
vec_s16 vsrc2ssH, vsrc3ssH, psum;
|
|
vec_u8 vdst, ppsum, vfdst, fsum;
|
|
#if HAVE_BIGENDIAN
|
|
register int loadSecond = (((unsigned long)src) % 16) <= 7 ? 0 : 1;
|
|
register int reallyBadAlign = (((unsigned long)src) % 16) == 15 ? 1 : 0;
|
|
vsrcperm0 = vec_lvsl(0, src);
|
|
vsrcperm1 = vec_lvsl(1, src);
|
|
#endif
|
|
|
|
if (((unsigned long)dst) % 16 == 0) {
|
|
fperm = (vec_u8){0x10, 0x11, 0x12, 0x13,
|
|
0x14, 0x15, 0x16, 0x17,
|
|
0x08, 0x09, 0x0A, 0x0B,
|
|
0x0C, 0x0D, 0x0E, 0x0F};
|
|
} else {
|
|
fperm = (vec_u8){0x00, 0x01, 0x02, 0x03,
|
|
0x04, 0x05, 0x06, 0x07,
|
|
0x18, 0x19, 0x1A, 0x1B,
|
|
0x1C, 0x1D, 0x1E, 0x1F};
|
|
}
|
|
|
|
GET_VSRC(vsrc0uc, vsrc1uc, 0, 16, vsrcperm0, vsrcperm1, src);
|
|
|
|
vsrc0ssH = (vec_s16)VEC_MERGEH(zero_u8v,(vec_u8)vsrc0uc);
|
|
vsrc1ssH = (vec_s16)VEC_MERGEH(zero_u8v,(vec_u8)vsrc1uc);
|
|
|
|
if (ABCD[3]) {
|
|
for (i = 0 ; i < h ; i++) {
|
|
GET_VSRC(vsrc2uc, vsrc3uc, stride, 16, vsrcperm0, vsrcperm1, src);
|
|
CHROMA_MC8_ALTIVEC_CORE(v32ss, noop);
|
|
}
|
|
} else {
|
|
const vec_s16 vE = vec_add(vB, vC);
|
|
if (ABCD[2]) { // x == 0 B == 0
|
|
for (i = 0 ; i < h ; i++) {
|
|
GET_VSRC1(vsrc1uc, stride, 15, vsrcperm0, src);
|
|
CHROMA_MC8_ALTIVEC_CORE_SIMPLE;
|
|
vsrc0uc = vsrc1uc;
|
|
}
|
|
} else { // y == 0 C == 0
|
|
for (i = 0 ; i < h ; i++) {
|
|
GET_VSRC(vsrc0uc, vsrc1uc, 0, 15, vsrcperm0, vsrcperm1, src);
|
|
CHROMA_MC8_ALTIVEC_CORE_SIMPLE;
|
|
}
|
|
}
|
|
}
|
|
}
|
|
#endif
|
|
|
|
/* this code assume that stride % 16 == 0 */
|
|
#ifdef PREFIX_no_rnd_vc1_chroma_mc8_altivec
|
|
static void PREFIX_no_rnd_vc1_chroma_mc8_altivec(uint8_t *dst, uint8_t *src,
|
|
ptrdiff_t stride, int h,
|
|
int x, int y)
|
|
{
|
|
DECLARE_ALIGNED(16, signed int, ABCD)[4] =
|
|
{((8 - x) * (8 - y)),
|
|
(( x) * (8 - y)),
|
|
((8 - x) * ( y)),
|
|
(( x) * ( y))};
|
|
register int i;
|
|
vec_u8 fperm;
|
|
LOAD_ZERO;
|
|
const vec_s32 vABCD = vec_ld(0, ABCD);
|
|
const vec_s16 vA = VEC_SPLAT16(vABCD, 1);
|
|
const vec_s16 vB = VEC_SPLAT16(vABCD, 3);
|
|
const vec_s16 vC = VEC_SPLAT16(vABCD, 5);
|
|
const vec_s16 vD = VEC_SPLAT16(vABCD, 7);
|
|
const vec_s16 v28ss = vec_sub(vec_sl(vec_splat_s16(1),vec_splat_u16(5)),vec_splat_s16(4));
|
|
const vec_u16 v6us = vec_splat_u16(6);
|
|
|
|
vec_u8 vsrcperm0, vsrcperm1;
|
|
vec_u8 vsrc0uc, vsrc1uc;
|
|
vec_s16 vsrc0ssH, vsrc1ssH;
|
|
vec_u8 vsrc2uc, vsrc3uc;
|
|
vec_s16 vsrc2ssH, vsrc3ssH, psum;
|
|
vec_u8 vdst, ppsum, vfdst, fsum;
|
|
#if HAVE_BIGENDIAN
|
|
register int loadSecond = (((unsigned long)src) % 16) <= 7 ? 0 : 1;
|
|
register int reallyBadAlign = (((unsigned long)src) % 16) == 15 ? 1 : 0;
|
|
vsrcperm0 = vec_lvsl(0, src);
|
|
vsrcperm1 = vec_lvsl(1, src);
|
|
#endif
|
|
|
|
if (((unsigned long)dst) % 16 == 0) {
|
|
fperm = (vec_u8){0x10, 0x11, 0x12, 0x13,
|
|
0x14, 0x15, 0x16, 0x17,
|
|
0x08, 0x09, 0x0A, 0x0B,
|
|
0x0C, 0x0D, 0x0E, 0x0F};
|
|
} else {
|
|
fperm = (vec_u8){0x00, 0x01, 0x02, 0x03,
|
|
0x04, 0x05, 0x06, 0x07,
|
|
0x18, 0x19, 0x1A, 0x1B,
|
|
0x1C, 0x1D, 0x1E, 0x1F};
|
|
}
|
|
|
|
GET_VSRC(vsrc0uc, vsrc1uc, 0, 16, vsrcperm0, vsrcperm1, src);
|
|
|
|
vsrc0ssH = (vec_s16)VEC_MERGEH(zero_u8v, (vec_u8)vsrc0uc);
|
|
vsrc1ssH = (vec_s16)VEC_MERGEH(zero_u8v, (vec_u8)vsrc1uc);
|
|
|
|
for (i = 0 ; i < h ; i++) {
|
|
GET_VSRC(vsrc2uc, vsrc3uc, stride, 16, vsrcperm0, vsrcperm1, src);
|
|
CHROMA_MC8_ALTIVEC_CORE(vec_splat_s16(0), add28);
|
|
}
|
|
}
|
|
#endif
|
|
|
|
#undef noop
|
|
#undef add28
|
|
#undef CHROMA_MC8_ALTIVEC_CORE
|