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     unsigned res = 0;
  31
  32     do {
  33         res   += *v1 * *v2++;
  34         *v1++ += mul * *v3++;
  35         res   += *v1 * *v2++;
  36         *v1++ += mul * *v3++;
  37     } while (order-=2);
  38     return res;
  39 }
  40
  41 static int32_t scalarproduct_and_madd_int32_c(int16_t *v1, const int32_t *v2,
  42                                               const int16_t *v3,
  43                                               int order, int mul)
  44 {
  45     int res = 0;
  46
  47     do {
  48         res   += *v1 * (uint32_t)*v2++;
  49         *v1++ += mul * *v3++;
  50         res   += *v1 * (uint32_t)*v2++;
  51         *v1++ += mul * *v3++;
  52     } while (order-=2);
  53     return res;
  54 }
  55
  56 av_cold void ff_llauddsp_init(LLAudDSPContext *c)
  57 {
  58     c->scalarproduct_and_madd_int16 = scalarproduct_and_madd_int16_c;
  59     c->scalarproduct_and_madd_int32 = scalarproduct_and_madd_int32_c;
  60
  61     if (ARCH_ARM)
  62         ff_llauddsp_init_arm(c);
  63     if (ARCH_PPC)
  64         ff_llauddsp_init_ppc(c);
  65     if (ARCH_X86)
  66         ff_llauddsp_init_x86(c);
  67 }