1*600f14f4SXin Li /* libFLAC - Free Lossless Audio Codec library 2*600f14f4SXin Li * Copyright (C) 2000-2009 Josh Coalson 3*600f14f4SXin Li * Copyright (C) 2011-2023 Xiph.Org Foundation 4*600f14f4SXin Li * 5*600f14f4SXin Li * Redistribution and use in source and binary forms, with or without 6*600f14f4SXin Li * modification, are permitted provided that the following conditions 7*600f14f4SXin Li * are met: 8*600f14f4SXin Li * 9*600f14f4SXin Li * - Redistributions of source code must retain the above copyright 10*600f14f4SXin Li * notice, this list of conditions and the following disclaimer. 11*600f14f4SXin Li * 12*600f14f4SXin Li * - Redistributions in binary form must reproduce the above copyright 13*600f14f4SXin Li * notice, this list of conditions and the following disclaimer in the 14*600f14f4SXin Li * documentation and/or other materials provided with the distribution. 15*600f14f4SXin Li * 16*600f14f4SXin Li * - Neither the name of the Xiph.org Foundation nor the names of its 17*600f14f4SXin Li * contributors may be used to endorse or promote products derived from 18*600f14f4SXin Li * this software without specific prior written permission. 19*600f14f4SXin Li * 20*600f14f4SXin Li * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS 21*600f14f4SXin Li * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT 22*600f14f4SXin Li * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR 23*600f14f4SXin Li * A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR 24*600f14f4SXin Li * CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, 25*600f14f4SXin Li * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, 26*600f14f4SXin Li * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR 27*600f14f4SXin Li * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF 28*600f14f4SXin Li * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING 29*600f14f4SXin Li * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS 30*600f14f4SXin Li * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 31*600f14f4SXin Li */ 32*600f14f4SXin Li 33*600f14f4SXin Li #ifndef FLAC__PRIVATE__FIXED_H 34*600f14f4SXin Li #define FLAC__PRIVATE__FIXED_H 35*600f14f4SXin Li 36*600f14f4SXin Li #ifdef HAVE_CONFIG_H 37*600f14f4SXin Li #include <config.h> 38*600f14f4SXin Li #endif 39*600f14f4SXin Li 40*600f14f4SXin Li #include "private/cpu.h" 41*600f14f4SXin Li #include "private/float.h" 42*600f14f4SXin Li #include "FLAC/format.h" 43*600f14f4SXin Li 44*600f14f4SXin Li /* 45*600f14f4SXin Li * FLAC__fixed_compute_best_predictor() 46*600f14f4SXin Li * -------------------------------------------------------------------- 47*600f14f4SXin Li * Compute the best fixed predictor and the expected bits-per-sample 48*600f14f4SXin Li * of the residual signal for each order. The _wide() version uses 49*600f14f4SXin Li * 64-bit integers which is statistically necessary when bits-per- 50*600f14f4SXin Li * sample + log2(blocksize) > 30 51*600f14f4SXin Li * 52*600f14f4SXin Li * IN data[0,data_len-1] 53*600f14f4SXin Li * IN data_len 54*600f14f4SXin Li * OUT residual_bits_per_sample[0,FLAC__MAX_FIXED_ORDER] 55*600f14f4SXin Li */ 56*600f14f4SXin Li #ifndef FLAC__INTEGER_ONLY_LIBRARY 57*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor(const FLAC__int32 data[], uint32_t data_len, float residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 58*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor_wide(const FLAC__int32 data[], uint32_t data_len, float residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 59*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor_limit_residual(const FLAC__int32 data[], uint32_t data_len, float residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 60*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor_limit_residual_33bit(const FLAC__int64 data[], uint32_t data_len, float residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 61*600f14f4SXin Li # ifndef FLAC__NO_ASM 62*600f14f4SXin Li # if (defined FLAC__CPU_IA32 || defined FLAC__CPU_X86_64) && FLAC__HAS_X86INTRIN 63*600f14f4SXin Li # ifdef FLAC__SSE2_SUPPORTED 64*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor_intrin_sse2(const FLAC__int32 data[], uint32_t data_len, float residual_bits_per_sample[FLAC__MAX_FIXED_ORDER + 1]); 65*600f14f4SXin Li # endif 66*600f14f4SXin Li # ifdef FLAC__SSSE3_SUPPORTED 67*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor_intrin_ssse3(const FLAC__int32 data[], uint32_t data_len, float residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 68*600f14f4SXin Li # endif 69*600f14f4SXin Li # ifdef FLAC__SSE4_2_SUPPORTED 70*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor_limit_residual_intrin_sse42(const FLAC__int32 data[], uint32_t data_len, float residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 71*600f14f4SXin Li # endif 72*600f14f4SXin Li # ifdef FLAC__AVX2_SUPPORTED 73*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor_wide_intrin_avx2(const FLAC__int32 data[], uint32_t data_len, float residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 74*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor_limit_residual_intrin_avx2(const FLAC__int32 data[], uint32_t data_len, float residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 75*600f14f4SXin Li # endif 76*600f14f4SXin Li # endif 77*600f14f4SXin Li # endif 78*600f14f4SXin Li #else 79*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor(const FLAC__int32 data[], uint32_t data_len, FLAC__fixedpoint residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 80*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor_wide(const FLAC__int32 data[], uint32_t data_len, FLAC__fixedpoint residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 81*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor_limit_residual(const FLAC__int32 data[], uint32_t data_len, FLAC__fixedpoint residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 82*600f14f4SXin Li uint32_t FLAC__fixed_compute_best_predictor_limit_residual_33bit(const FLAC__int64 data[], uint32_t data_len, FLAC__fixedpoint residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]); 83*600f14f4SXin Li #endif 84*600f14f4SXin Li 85*600f14f4SXin Li /* 86*600f14f4SXin Li * FLAC__fixed_compute_residual() 87*600f14f4SXin Li * -------------------------------------------------------------------- 88*600f14f4SXin Li * Compute the residual signal obtained from sutracting the predicted 89*600f14f4SXin Li * signal from the original. 90*600f14f4SXin Li * 91*600f14f4SXin Li * IN data[-order,data_len-1] original signal (NOTE THE INDICES!) 92*600f14f4SXin Li * IN data_len length of original signal 93*600f14f4SXin Li * IN order <= FLAC__MAX_FIXED_ORDER fixed-predictor order 94*600f14f4SXin Li * OUT residual[0,data_len-1] residual signal 95*600f14f4SXin Li */ 96*600f14f4SXin Li void FLAC__fixed_compute_residual(const FLAC__int32 data[], uint32_t data_len, uint32_t order, FLAC__int32 residual[]); 97*600f14f4SXin Li void FLAC__fixed_compute_residual_wide(const FLAC__int32 data[], uint32_t data_len, uint32_t order, FLAC__int32 residual[]); 98*600f14f4SXin Li void FLAC__fixed_compute_residual_wide_33bit(const FLAC__int64 data[], uint32_t data_len, uint32_t order, FLAC__int32 residual[]); 99*600f14f4SXin Li 100*600f14f4SXin Li /* 101*600f14f4SXin Li * FLAC__fixed_restore_signal() 102*600f14f4SXin Li * -------------------------------------------------------------------- 103*600f14f4SXin Li * Restore the original signal by summing the residual and the 104*600f14f4SXin Li * predictor. 105*600f14f4SXin Li * 106*600f14f4SXin Li * IN residual[0,data_len-1] residual signal 107*600f14f4SXin Li * IN data_len length of original signal 108*600f14f4SXin Li * IN order <= FLAC__MAX_FIXED_ORDER fixed-predictor order 109*600f14f4SXin Li * *** IMPORTANT: the caller must pass in the historical samples: 110*600f14f4SXin Li * IN data[-order,-1] previously-reconstructed historical samples 111*600f14f4SXin Li * OUT data[0,data_len-1] original signal 112*600f14f4SXin Li */ 113*600f14f4SXin Li void FLAC__fixed_restore_signal(const FLAC__int32 residual[], uint32_t data_len, uint32_t order, FLAC__int32 data[]); 114*600f14f4SXin Li void FLAC__fixed_restore_signal_wide(const FLAC__int32 residual[], uint32_t data_len, uint32_t order, FLAC__int32 data[]); 115*600f14f4SXin Li void FLAC__fixed_restore_signal_wide_33bit(const FLAC__int32 residual[], uint32_t data_len, uint32_t order, FLAC__int64 data[]); 116*600f14f4SXin Li 117*600f14f4SXin Li #endif 118