1/*
2 * Copyright (c) 2006 Luca Barbato <lu_zero@gentoo.org>
3 *
4 * This file is part of FFmpeg.
5 *
6 * FFmpeg is free software; you can redistribute it and/or
7 * modify it under the terms of the GNU Lesser General Public
8 * License as published by the Free Software Foundation; either
9 * version 2.1 of the License, or (at your option) any later version.
10 *
11 * FFmpeg is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14 * Lesser General Public License for more details.
15 *
16 * You should have received a copy of the GNU Lesser General Public
17 * License along with FFmpeg; if not, write to the Free Software
18 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19 */
20
21#include "config.h"
22#if HAVE_ALTIVEC_H
23#include <altivec.h>
24#endif
25#include "libavutil/attributes.h"
26#include "libavutil/cpu.h"
27#include "libavutil/ppc/cpu.h"
28#include "libavcodec/vorbisdsp.h"
29
30#if HAVE_ALTIVEC
31static void vorbis_inverse_coupling_altivec(float *mag, float *ang,
32                                            intptr_t blocksize)
33{
34    int i;
35    vector float m, a;
36    vector bool int t0, t1;
37    const vector unsigned int v_31 = //XXX
38        vec_add(vec_add(vec_splat_u32(15),vec_splat_u32(15)),vec_splat_u32(1));
39    for (i = 0; i < blocksize; i += 4) {
40        m = vec_ld(0, mag+i);
41        a = vec_ld(0, ang+i);
42        t0 = vec_cmple(m, (vector float)vec_splat_u32(0));
43        t1 = vec_cmple(a, (vector float)vec_splat_u32(0));
44        a = vec_xor(a, (vector float) vec_sl((vector unsigned int)t0, v_31));
45        t0 = (vector bool int)vec_and(a, t1);
46        t1 = (vector bool int)vec_andc(a, t1);
47        a = vec_sub(m, (vector float)t1);
48        m = vec_add(m, (vector float)t0);
49        vec_stl(a, 0, ang+i);
50        vec_stl(m, 0, mag+i);
51    }
52}
53#endif /* HAVE_ALTIVEC */
54
55av_cold void ff_vorbisdsp_init_ppc(VorbisDSPContext *c)
56{
57#if HAVE_ALTIVEC
58    if (!PPC_ALTIVEC(av_get_cpu_flags()))
59        return;
60
61    c->vorbis_inverse_coupling = vorbis_inverse_coupling_altivec;
62#endif /* HAVE_ALTIVEC */
63}
64