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: #ifndef AVCODEC_FFT_INTERNAL_H
yading@10: #define AVCODEC_FFT_INTERNAL_H
yading@10: 
yading@10: #if CONFIG_FFT_FLOAT
yading@10: 
yading@10: #define FIX15(v) (v)
yading@10: #define sqrthalf (float)M_SQRT1_2
yading@10: 
yading@10: #define BF(x, y, a, b) do {                     \
yading@10:         x = a - b;                              \
yading@10:         y = a + b;                              \
yading@10:     } while (0)
yading@10: 
yading@10: #define CMUL(dre, dim, are, aim, bre, bim) do { \
yading@10:         (dre) = (are) * (bre) - (aim) * (bim);  \
yading@10:         (dim) = (are) * (bim) + (aim) * (bre);  \
yading@10:     } while (0)
yading@10: 
yading@10: #else
yading@10: 
yading@10: #include "fft.h"
yading@10: #include "mathops.h"
yading@10: 
yading@10: void ff_mdct_calcw_c(FFTContext *s, FFTDouble *output, const FFTSample *input);
yading@10: 
yading@10: #define SCALE_FLOAT(a, bits) lrint((a) * (double)(1 << (bits)))
yading@10: #define FIX15(a) av_clip(SCALE_FLOAT(a, 15), -32767, 32767)
yading@10: 
yading@10: #define sqrthalf ((int16_t)((1<<15)*M_SQRT1_2))
yading@10: 
yading@10: #define BF(x, y, a, b) do {                     \
yading@10:         x = (a - b) >> 1;                       \
yading@10:         y = (a + b) >> 1;                       \
yading@10:     } while (0)
yading@10: 
yading@10: #define CMULS(dre, dim, are, aim, bre, bim, sh) do {            \
yading@10:         (dre) = (MUL16(are, bre) - MUL16(aim, bim)) >> sh;      \
yading@10:         (dim) = (MUL16(are, bim) + MUL16(aim, bre)) >> sh;      \
yading@10:     } while (0)
yading@10: 
yading@10: #define CMUL(dre, dim, are, aim, bre, bim)      \
yading@10:     CMULS(dre, dim, are, aim, bre, bim, 15)
yading@10: 
yading@10: #define CMULL(dre, dim, are, aim, bre, bim)     \
yading@10:     CMULS(dre, dim, are, aim, bre, bim, 0)
yading@10: 
yading@10: #endif /* CONFIG_FFT_FLOAT */
yading@10: 
yading@10: #define ff_imdct_calc_c FFT_NAME(ff_imdct_calc_c)
yading@10: #define ff_imdct_half_c FFT_NAME(ff_imdct_half_c)
yading@10: #define ff_mdct_calc_c  FFT_NAME(ff_mdct_calc_c)
yading@10: 
yading@10: void ff_imdct_calc_c(FFTContext *s, FFTSample *output, const FFTSample *input);
yading@10: void ff_imdct_half_c(FFTContext *s, FFTSample *output, const FFTSample *input);
yading@10: void ff_mdct_calc_c(FFTContext *s, FFTSample *output, const FFTSample *input);
yading@10: 
yading@10: #endif /* AVCODEC_FFT_INTERNAL_H */