Merge commit 'eecd29b3fd7fee221580798346d6582b75c7ade4'
authorMichael Niedermayer <michaelni@gmx.at>
Sun, 31 May 2015 21:23:58 +0000 (23:23 +0200)
committerMichael Niedermayer <michaelni@gmx.at>
Sun, 31 May 2015 21:42:36 +0000 (23:42 +0200)
* commit 'eecd29b3fd7fee221580798346d6582b75c7ade4':
  ppc: vsx: Implement float_dsp

Conflicts:
libavutil/ppc/float_dsp_init.c

Merged-by: Michael Niedermayer <michaelni@gmx.at>
1  2 
arch.mak
libavutil/ppc/float_dsp_init.c
libavutil/ppc/float_dsp_vsx.c
libavutil/ppc/float_dsp_vsx.h

diff --cc arch.mak
+++ b/arch.mak
@@@ -4,13 -4,8 +4,14 @@@ OBJS-$(HAVE_ARMV8)   += $(ARMV8-OBJS
  OBJS-$(HAVE_VFP)     += $(VFP-OBJS)     $(VFP-OBJS-yes)
  OBJS-$(HAVE_NEON)    += $(NEON-OBJS)    $(NEON-OBJS-yes)
  
 +OBJS-$(HAVE_MIPSFPU)   += $(MIPSFPU-OBJS)    $(MIPSFPU-OBJS-yes)
 +OBJS-$(HAVE_MIPSDSPR1) += $(MIPSDSPR1-OBJS)  $(MIPSDSPR1-OBJS-yes)
 +OBJS-$(HAVE_MIPSDSPR2) += $(MIPSDSPR2-OBJS)  $(MIPSDSPR2-OBJS-yes)
 +OBJS-$(HAVE_MSA)       += $(MSA-OBJS)        $(MSA-OBJS-yes)
 +OBJS-$(HAVE_LOONGSON3) += $(LOONGSON3-OBJS)  $(LOONGSON3-OBJS-yes)
 +
  OBJS-$(HAVE_ALTIVEC) += $(ALTIVEC-OBJS) $(ALTIVEC-OBJS-yes)
+ OBJS-$(HAVE_VSX)     += $(VSX-OBJS) $(VSX-OBJS-yes)
  
  OBJS-$(HAVE_MMX)     += $(MMX-OBJS)     $(MMX-OBJS-yes)
  OBJS-$(HAVE_YASM)    += $(YASM-OBJS)    $(YASM-OBJS-yes)
  
  av_cold void ff_float_dsp_init_ppc(AVFloatDSPContext *fdsp, int bit_exact)
  {
-     if (!PPC_ALTIVEC(av_get_cpu_flags()))
-         return;
 -    if (HAVE_BIGENDIAN && PPC_ALTIVEC(av_get_cpu_flags())) {
++    if (PPC_ALTIVEC(av_get_cpu_flags())) {
+         fdsp->vector_fmul = ff_vector_fmul_altivec;
+         fdsp->vector_fmul_add = ff_vector_fmul_add_altivec;
+         fdsp->vector_fmul_reverse = ff_vector_fmul_reverse_altivec;
  
-     fdsp->vector_fmul = ff_vector_fmul_altivec;
-     fdsp->vector_fmul_add = ff_vector_fmul_add_altivec;
-     fdsp->vector_fmul_reverse = ff_vector_fmul_reverse_altivec;
+         if (!bit_exact) {
+             fdsp->vector_fmul_window = ff_vector_fmul_window_altivec;
+         }
+     }
+     if (PPC_VSX(av_get_cpu_flags())) {
+         fdsp->vector_fmul = ff_vector_fmul_vsx;
+         fdsp->vector_fmul_add = ff_vector_fmul_add_vsx;
+         fdsp->vector_fmul_reverse = ff_vector_fmul_reverse_vsx;
  
-     if (!bit_exact) {
-         fdsp->vector_fmul_window = ff_vector_fmul_window_altivec;
+         if (!bit_exact) {
+             fdsp->vector_fmul_window = ff_vector_fmul_window_vsx;
+         }
      }
  }
index 0000000,4135db1..a91e7c1
mode 000000,100644..100644
--- /dev/null
@@@ -1,0 -1,117 +1,117 @@@
 - * This file is part of Libav.
+ /*
+  * Copyright (c) 2015 Luca Barbato <lu_zero@gentoo.org>
+  *
 - * Libav is free software; you can redistribute it and/or
++ * This file is part of FFmpeg.
+  *
 - * Libav is distributed in the hope that it will be useful,
++ * 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.
+  *
 - * License along with Libav; if not, write to the Free Software
++ * 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 "util_altivec.h"
+ #include "float_dsp_vsx.h"
+ void ff_vector_fmul_vsx(float *dst,
+                         const float *src0, const float *src1,
+                         int len)
+ {
+     int i;
+     vec_f d0, d1, zero = (vec_f)vec_splat_u32(0);
+     for (i = 0; i < len - 7; i += 8) {
+         d0 = vec_vsx_ld( 0, src0 + i);
+         d1 = vec_vsx_ld(16, src0 + i);
+         d0 = vec_madd(d0, vec_vsx_ld( 0, src1 + i), zero);
+         d1 = vec_madd(d1, vec_vsx_ld(16, src1 + i), zero);
+         vec_vsx_st(d0,  0, dst + i);
+         vec_vsx_st(d1, 16, dst + i);
+     }
+ }
+ void ff_vector_fmul_window_vsx(float *dst, const float *src0,
+                                const float *src1, const float *win,
+                                int len)
+ {
+     vec_f zero, t0, t1, s0, s1, wi, wj;
+     const vec_u8 reverse = vcprm(3, 2, 1, 0);
+     int i, j;
+     dst  += len;
+     win  += len;
+     src0 += len;
+     zero = (vec_f)vec_splat_u32(0);
+     for (i = -len * 4, j = len * 4 - 16; i < 0; i += 16, j -= 16) {
+         s0 = vec_vsx_ld(i, src0);
+         s1 = vec_vsx_ld(j, src1);
+         wi = vec_vsx_ld(i, win);
+         wj = vec_vsx_ld(j, win);
+         s1 = vec_perm(s1, s1, reverse);
+         wj = vec_perm(wj, wj, reverse);
+         t0 = vec_madd(s0, wj, zero);
+         t0 = vec_nmsub(s1, wi, t0);
+         t1 = vec_madd(s0, wi, zero);
+         t1 = vec_madd(s1, wj, t1);
+         t1 = vec_perm(t1, t1, reverse);
+         vec_vsx_st(t0, i, dst);
+         vec_vsx_st(t1, j, dst);
+     }
+ }
+ void ff_vector_fmul_add_vsx(float *dst, const float *src0,
+                             const float *src1, const float *src2,
+                             int len)
+ {
+     int i;
+     vec_f d, s0, s1, s2;
+     for (i = 0; i < len - 3; i += 4) {
+         s0 = vec_vsx_ld(0, src0 + i);
+         s1 = vec_vsx_ld(0, src1 + i);
+         s2 = vec_vsx_ld(0, src2 + i);
+         d = vec_madd(s0, s1, s2);
+         vec_vsx_st(d, 0, dst + i);
+     }
+ }
+ void ff_vector_fmul_reverse_vsx(float *dst, const float *src0,
+                                 const float *src1, int len)
+ {
+     int i;
+     vec_f d, s0, s1, h0, l0, s2, s3;
+     vec_f zero = (vec_f)vec_splat_u32(0);
+     src1 += len - 4;
+     for (i = 0; i < len - 7; i += 8) {
+         s1 = vec_vsx_ld(0, src1 - i);              // [a,b,c,d]
+         s0 = vec_vsx_ld(0, src0 + i);
+         l0 = vec_mergel(s1, s1);               // [c,c,d,d]
+         s3 = vec_vsx_ld(-16, src1 - i);
+         h0 = vec_mergeh(s1, s1);               // [a,a,b,b]
+         s2 = vec_vsx_ld(16, src0 + i);
+         s1 = vec_mergeh(vec_mergel(l0, h0),    // [d,b,d,b]
+                         vec_mergeh(l0, h0));   // [c,a,c,a]
+         // [d,c,b,a]
+         l0 = vec_mergel(s3, s3);
+         d = vec_madd(s0, s1, zero);
+         h0 = vec_mergeh(s3, s3);
+         vec_vsx_st(d, 0, dst + i);
+         s3 = vec_mergeh(vec_mergel(l0, h0),
+                         vec_mergeh(l0, h0));
+         d = vec_madd(s2, s3, zero);
+         vec_vsx_st(d, 16, dst + i);
+     }
+ }
index 0000000,32860dc..5d3cc5e
mode 000000,100644..100644
--- /dev/null
@@@ -1,0 -1,38 +1,38 @@@
 - * This file is part of Libav.
+ /*
+  * Copyright (c) 2015 Luca Barbato <lu_zero@gentoo.org>
+  *
 - * Libav is free software; you can redistribute it and/or
++ * This file is part of FFmpeg.
+  *
 - * Libav is distributed in the hope that it will be useful,
++ * 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.
+  *
 - * License along with Libav; if not, write to the Free Software
++ * 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
+  */
+ #ifndef AVUTIL_PPC_FLOAT_DSP_VSX_H
+ #define AVUTIL_PPC_FLOAT_DSP_VSX_H
+ void ff_vector_fmul_vsx(float *dst, const float *src0,
+                         const float *src1, int len);
+ void ff_vector_fmul_window_vsx(float *dst, const float *src0,
+                                const float *src1, const float *win,
+                                int len);
+ void ff_vector_fmul_add_vsx(float *dst, const float *src0,
+                             const float *src1, const float *src2,
+                             int len);
+ void ff_vector_fmul_reverse_vsx(float *dst, const float *src0,
+                                 const float *src1, int len);
+ #endif /* AVUTIL_PPC_FLOAT_DSP_VSX_H */