yading@10: /*
yading@10:  * Copyright (c) 2012 Mans Rullgard <mans@mansr.com>
yading@10:  *
yading@10:  * This file is part of FFmpeg.
yading@10:  *
yading@10:  * FFmpeg is free software; you can redistribute it and/or
yading@10:  * modify it under the terms of the GNU Lesser General Public
yading@10:  * License as published by the Free Software Foundation; either
yading@10:  * version 2.1 of the License, or (at your option) any later version.
yading@10:  *
yading@10:  * FFmpeg is distributed in the hope that it will be useful,
yading@10:  * but WITHOUT ANY WARRANTY; without even the implied warranty of
yading@10:  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
yading@10:  * Lesser General Public License for more details.
yading@10:  *
yading@10:  * You should have received a copy of the GNU Lesser General Public
yading@10:  * License along with FFmpeg; if not, write to the Free Software
yading@10:  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
yading@10:  */
yading@10: 
yading@10: #include <stdint.h>
yading@10: #include "libavutil/avutil.h"
yading@10: 
yading@10: #undef FUNC
yading@10: #undef FSUF
yading@10: #undef sample
yading@10: #undef sample_type
yading@10: #undef OUT
yading@10: #undef S
yading@10: 
yading@10: #if SAMPLE_SIZE == 32
yading@10: #   define sample_type  int32_t
yading@10: #else
yading@10: #   define sample_type  int16_t
yading@10: #endif
yading@10: 
yading@10: #if PLANAR
yading@10: #   define FSUF   AV_JOIN(SAMPLE_SIZE, p)
yading@10: #   define sample sample_type *
yading@10: #   define OUT(n) n
yading@10: #   define S(s, c, i) (s[c][i])
yading@10: #else
yading@10: #   define FSUF   SAMPLE_SIZE
yading@10: #   define sample sample_type
yading@10: #   define OUT(n) n[0]
yading@10: #   define S(s, c, i) (*s++)
yading@10: #endif
yading@10: 
yading@10: #define FUNC(n) AV_JOIN(n ## _, FSUF)
yading@10: 
yading@10: static void FUNC(flac_decorrelate_indep_c)(uint8_t **out, int32_t **in,
yading@10:                                            int channels, int len, int shift)
yading@10: {
yading@10:     sample *samples = (sample *) OUT(out);
yading@10:     int i, j;
yading@10: 
yading@10:     for (j = 0; j < len; j++)
yading@10:         for (i = 0; i < channels; i++)
yading@10:             S(samples, i, j) = in[i][j] << shift;
yading@10: }
yading@10: 
yading@10: static void FUNC(flac_decorrelate_ls_c)(uint8_t **out, int32_t **in,
yading@10:                                         int channels, int len, int shift)
yading@10: {
yading@10:     sample *samples = (sample *) OUT(out);
yading@10:     int i;
yading@10: 
yading@10:     for (i = 0; i < len; i++) {
yading@10:         int a = in[0][i];
yading@10:         int b = in[1][i];
yading@10:         S(samples, 0, i) =  a      << shift;
yading@10:         S(samples, 1, i) = (a - b) << shift;
yading@10:     }
yading@10: }
yading@10: 
yading@10: static void FUNC(flac_decorrelate_rs_c)(uint8_t **out, int32_t **in,
yading@10:                                         int channels, int len, int shift)
yading@10: {
yading@10:     sample *samples = (sample *) OUT(out);
yading@10:     int i;
yading@10: 
yading@10:     for (i = 0; i < len; i++) {
yading@10:         int a = in[0][i];
yading@10:         int b = in[1][i];
yading@10:         S(samples, 0, i) = (a + b) << shift;
yading@10:         S(samples, 1, i) =  b      << shift;
yading@10:     }
yading@10: }
yading@10: 
yading@10: static void FUNC(flac_decorrelate_ms_c)(uint8_t **out, int32_t **in,
yading@10:                                         int channels, int len, int shift)
yading@10: {
yading@10:     sample *samples = (sample *) OUT(out);
yading@10:     int i;
yading@10: 
yading@10:     for (i = 0; i < len; i++) {
yading@10:         int a = in[0][i];
yading@10:         int b = in[1][i];
yading@10:         a -= b >> 1;
yading@10:         S(samples, 0, i) = (a + b) << shift;
yading@10:         S(samples, 1, i) =  a      << shift;
yading@10:     }
yading@10: }