liboil/dct/imdct32_f32.c

   1 /*
   2  * LIBOIL - Library of Optimized Inner Loops
   3  * Copyright (c) 2003,2004 David A. Schleef <ds@schleef.org>
   4  * All rights reserved.
   5  *
   6  * Redistribution and use in source and binary forms, with or without
   7  * modification, are permitted provided that the following conditions
   8  * are met:
   9  * 1. Redistributions of source code must retain the above copyright
  10  *    notice, this list of conditions and the following disclaimer.
  11  * 2. Redistributions in binary form must reproduce the above copyright
  12  *    notice, this list of conditions and the following disclaimer in the
  13  *    documentation and/or other materials provided with the distribution.
  14  *
  15  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
  16  * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
  17  * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  18  * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT,
  19  * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
  20  * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
  21  * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  22  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
  23  * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING
  24  * IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  25  * POSSIBILITY OF SUCH DAMAGE.
  26  */
  27
  28 #ifdef HAVE_CONFIG_H
  29 #include "config.h"
  30 #endif
  31
  32 #include <liboil/liboilfunction.h>
  33 #include <liboil/dct/dct.h>
  34 #include <math.h>
  35
  36
  37 /**
  38  * oil_imdct32_f32:
  39  * @d_32:
  40  * @s_32:
  41  *
  42  * Performs a Inverse Modified Discrete Cosine Transform (IMDCT)
  43  * on the source array @s_32 and places the result in @d_32.
  44  */
  45 OIL_DEFINE_CLASS (imdct32_f32, "float *d_32, float *s_32");
  46
  47 static void imdct32_f32_ref (float *dest, const float *src)
  48 {
  49         double x;
  50         int i,j;
  51         double coeff;
  52
  53         for(i=0;i<32;i++){
  54                 x = 0;
  55                 for(j=0;j<32;j++){
  56                         coeff = cos((M_PI/32)*i*(j+0.5));
  57                         x += coeff * src[j];
  58                 }
  59                 dest[i] = x;
  60         }
  61 }
  62
  63 OIL_DEFINE_IMPL_REF (imdct32_f32_ref, imdct32_f32);
  64