cd8bbb90fcdedd257914114a2ca20b2ac62ebd8b
[ffmpeg.git] / libavcodec / ppc / vorbisdsp_altivec.c
1 /*
2  * Copyright (c) 2006 Luca Barbato <lu_zero@gentoo.org>
3  *
4  * This file is part of Libav.
5  *
6  * Libav 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  * Libav 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 Libav; 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/cpu.h"
26 #include "libavutil/ppc/types_altivec.h"
27 #include "libavutil/ppc/util_altivec.h"
28 #include "libavcodec/vorbisdsp.h"
29
30 #if HAVE_ALTIVEC
31 static void vorbis_inverse_coupling_altivec(float *mag, float *ang,
32                                             int 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
55 void ff_vorbisdsp_init_ppc(VorbisDSPContext* c)
56 {
57 #if HAVE_ALTIVEC
58     if (av_get_cpu_flags() & AV_CPU_FLAG_ALTIVEC) {
59       c->vorbis_inverse_coupling = vorbis_inverse_coupling_altivec;
60     }
61 #endif /* HAVE_ALTIVEC */
62 }