git.sesse.net Git - ffmpeg/blob - libavcodec/lossless_audiodsp.c

   1 /*
   2  * Monkey's Audio lossless audio decoder
   3  * Copyright (c) 2007 Benjamin Zores <ben@geexbox.org>
   4  *  based upon libdemac from Dave Chapman.
   5  *
   6  * This file is part of FFmpeg.
   7  *
   8  * FFmpeg is free software; you can redistribute it and/or
   9  * modify it under the terms of the GNU Lesser General Public
  10  * License as published by the Free Software Foundation; either
  11  * version 2.1 of the License, or (at your option) any later version.
  12  *
  13  * FFmpeg is distributed in the hope that it will be useful,
  14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  16  * Lesser General Public License for more details.
  17  *
  18  * You should have received a copy of the GNU Lesser General Public
  19  * License along with FFmpeg; if not, write to the Free Software
  20  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  21  */
  22
  23 #include "avcodec.h"
  24 #include "lossless_audiodsp.h"
  25
  26 static int32_t scalarproduct_and_madd_int16_c(int16_t *v1, const int16_t *v2,
  27                                               const int16_t *v3,
  28                                               int order, int mul)
  29 {
  30     int res = 0;
  31
  32     while (order--) {
  33         res   += *v1 * *v2++;
  34         *v1++ += mul * *v3++;
  35     }
  36     return res;
  37 }
  38
  39 static int32_t scalarproduct_and_madd_int32_c(int32_t *v1, const int32_t *v2,
  40                                               const int32_t *v3,
  41                                               int order, int mul)
  42 {
  43     int res = 0;
  44
  45     while (order--) {
  46         res   += *v1 * *v2++;
  47         *v1++ += mul * *v3++;
  48     }
  49     return res;
  50 }
  51
  52 av_cold void ff_llauddsp_init(LLAudDSPContext *c)
  53 {
  54     c->scalarproduct_and_madd_int16 = scalarproduct_and_madd_int16_c;
  55     c->scalarproduct_and_madd_int32 = scalarproduct_and_madd_int32_c;
  56
  57     if (ARCH_ARM)
  58         ff_llauddsp_init_arm(c);
  59     if (ARCH_PPC)
  60         ff_llauddsp_init_ppc(c);
  61     if (ARCH_X86)
  62         ff_llauddsp_init_x86(c);
  63 }