x86/g722dsp: add ff_g722_apply_qmf_sse2
[ffmpeg.git] / libavcodec / g722dsp.c
index 5d3aedb..f148053 100644 (file)
@@ -1,20 +1,20 @@
 /*
  * Copyright (c) 2015 Peter Meerwald <pmeerw@pmeerw.net>
  *
 /*
  * Copyright (c) 2015 Peter Meerwald <pmeerw@pmeerw.net>
  *
- * This file is part of Libav.
+ * This file is part of FFmpeg.
  *
  *
- * Libav is free software; you can redistribute it and/or
+ * 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.
  *
  * 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.
  *
- * Libav is distributed in the hope that it will be useful,
+ * 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
  * 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 Libav; if not, write to the Free Software
+ * License along with FFmpeg; if not, write to the Free Software
  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  */
 
  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  */
 
  * quadrature mirror filter (QMF) coefficients (ITU-T G.722 Table 11) inlined
  * in code below: 3, -11, 12, 32, -210, 951, 3876, -805, 362, -156, 53, -11
  */
  * quadrature mirror filter (QMF) coefficients (ITU-T G.722 Table 11) inlined
  * in code below: 3, -11, 12, 32, -210, 951, 3876, -805, 362, -156, 53, -11
  */
-static const int16_t qmf_coeffs[12] = {
-    3, -11, 12, 32, -210, 951, 3876, -805, 362, -156, 53, -11,
-};
 
 static void g722_apply_qmf(const int16_t *prev_samples, int xout[2])
 {
 
 static void g722_apply_qmf(const int16_t *prev_samples, int xout[2])
 {
-    int i;
-
-    xout[0] = 0;
-    xout[1] = 0;
-    for (i = 0; i < 12; i++) {
-        MAC16(xout[1], prev_samples[2*i  ], qmf_coeffs[i   ]);
-        MAC16(xout[0], prev_samples[2*i+1], qmf_coeffs[11-i]);
-    }
+    xout[1] = MUL16(*prev_samples++, 3);
+    xout[0] = MUL16(*prev_samples++, -11);
+
+    MAC16(xout[1], *prev_samples++, -11);
+    MAC16(xout[0], *prev_samples++, 53);
+
+    MAC16(xout[1], *prev_samples++, 12);
+    MAC16(xout[0], *prev_samples++, -156);
+
+    MAC16(xout[1], *prev_samples++, 32);
+    MAC16(xout[0], *prev_samples++, 362);
+
+    MAC16(xout[1], *prev_samples++, -210);
+    MAC16(xout[0], *prev_samples++, -805);
+
+    MAC16(xout[1], *prev_samples++, 951);
+    MAC16(xout[0], *prev_samples++, 3876);
+
+    MAC16(xout[1], *prev_samples++, 3876);
+    MAC16(xout[0], *prev_samples++, 951);
+
+    MAC16(xout[1], *prev_samples++, -805);
+    MAC16(xout[0], *prev_samples++, -210);
+
+    MAC16(xout[1], *prev_samples++, 362);
+    MAC16(xout[0], *prev_samples++, 32);
+
+    MAC16(xout[1], *prev_samples++, -156);
+    MAC16(xout[0], *prev_samples++, 12);
+
+    MAC16(xout[1], *prev_samples++, 53);
+    MAC16(xout[0], *prev_samples++, -11);
+
+    MAC16(xout[1], *prev_samples++, -11);
+    MAC16(xout[0], *prev_samples++, 3);
 }
 
 av_cold void ff_g722dsp_init(G722DSPContext *c)
 {
     c->apply_qmf = g722_apply_qmf;
 }
 
 av_cold void ff_g722dsp_init(G722DSPContext *c)
 {
     c->apply_qmf = g722_apply_qmf;
+
+    if (ARCH_ARM)
+        ff_g722dsp_init_arm(c);
+    if (ARCH_X86)
+        ff_g722dsp_init_x86(c);
 }
 }