michael@0: /*********************************************************************** michael@0: Copyright (c) 2006-2011, Skype Limited. All rights reserved. michael@0: Redistribution and use in source and binary forms, with or without michael@0: modification, are permitted provided that the following conditions michael@0: are met: michael@0: - Redistributions of source code must retain the above copyright notice, michael@0: this list of conditions and the following disclaimer. michael@0: - Redistributions in binary form must reproduce the above copyright michael@0: notice, this list of conditions and the following disclaimer in the michael@0: documentation and/or other materials provided with the distribution. michael@0: - Neither the name of Internet Society, IETF or IETF Trust, nor the michael@0: names of specific contributors, may be used to endorse or promote michael@0: products derived from this software without specific prior written michael@0: permission. michael@0: THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" michael@0: AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE michael@0: IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE michael@0: ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE michael@0: LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR michael@0: CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF michael@0: SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS michael@0: INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN michael@0: CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) michael@0: ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE michael@0: POSSIBILITY OF SUCH DAMAGE. michael@0: ***********************************************************************/ michael@0: michael@0: #ifdef HAVE_CONFIG_H michael@0: #include "config.h" michael@0: #endif michael@0: michael@0: #include "main.h" michael@0: michael@0: /*#define silk_enc_map(a) ((a) > 0 ? 1 : 0)*/ michael@0: /*#define silk_dec_map(a) ((a) > 0 ? 1 : -1)*/ michael@0: /* shifting avoids if-statement */ michael@0: #define silk_enc_map(a) ( silk_RSHIFT( (a), 15 ) + 1 ) michael@0: #define silk_dec_map(a) ( silk_LSHIFT( (a), 1 ) - 1 ) michael@0: michael@0: /* Encodes signs of excitation */ michael@0: void silk_encode_signs( michael@0: ec_enc *psRangeEnc, /* I/O Compressor data structure */ michael@0: const opus_int8 pulses[], /* I pulse signal */ michael@0: opus_int length, /* I length of input */ michael@0: const opus_int signalType, /* I Signal type */ michael@0: const opus_int quantOffsetType, /* I Quantization offset type */ michael@0: const opus_int sum_pulses[ MAX_NB_SHELL_BLOCKS ] /* I Sum of absolute pulses per block */ michael@0: ) michael@0: { michael@0: opus_int i, j, p; michael@0: opus_uint8 icdf[ 2 ]; michael@0: const opus_int8 *q_ptr; michael@0: const opus_uint8 *icdf_ptr; michael@0: michael@0: icdf[ 1 ] = 0; michael@0: q_ptr = pulses; michael@0: i = silk_SMULBB( 7, silk_ADD_LSHIFT( quantOffsetType, signalType, 1 ) ); michael@0: icdf_ptr = &silk_sign_iCDF[ i ]; michael@0: length = silk_RSHIFT( length + SHELL_CODEC_FRAME_LENGTH/2, LOG2_SHELL_CODEC_FRAME_LENGTH ); michael@0: for( i = 0; i < length; i++ ) { michael@0: p = sum_pulses[ i ]; michael@0: if( p > 0 ) { michael@0: icdf[ 0 ] = icdf_ptr[ silk_min( p & 0x1F, 6 ) ]; michael@0: for( j = 0; j < SHELL_CODEC_FRAME_LENGTH; j++ ) { michael@0: if( q_ptr[ j ] != 0 ) { michael@0: ec_enc_icdf( psRangeEnc, silk_enc_map( q_ptr[ j ]), icdf, 8 ); michael@0: } michael@0: } michael@0: } michael@0: q_ptr += SHELL_CODEC_FRAME_LENGTH; michael@0: } michael@0: } michael@0: michael@0: /* Decodes signs of excitation */ michael@0: void silk_decode_signs( michael@0: ec_dec *psRangeDec, /* I/O Compressor data structure */ michael@0: opus_int pulses[], /* I/O pulse signal */ michael@0: opus_int length, /* I length of input */ michael@0: const opus_int signalType, /* I Signal type */ michael@0: const opus_int quantOffsetType, /* I Quantization offset type */ michael@0: const opus_int sum_pulses[ MAX_NB_SHELL_BLOCKS ] /* I Sum of absolute pulses per block */ michael@0: ) michael@0: { michael@0: opus_int i, j, p; michael@0: opus_uint8 icdf[ 2 ]; michael@0: opus_int *q_ptr; michael@0: const opus_uint8 *icdf_ptr; michael@0: michael@0: icdf[ 1 ] = 0; michael@0: q_ptr = pulses; michael@0: i = silk_SMULBB( 7, silk_ADD_LSHIFT( quantOffsetType, signalType, 1 ) ); michael@0: icdf_ptr = &silk_sign_iCDF[ i ]; michael@0: length = silk_RSHIFT( length + SHELL_CODEC_FRAME_LENGTH/2, LOG2_SHELL_CODEC_FRAME_LENGTH ); michael@0: for( i = 0; i < length; i++ ) { michael@0: p = sum_pulses[ i ]; michael@0: if( p > 0 ) { michael@0: icdf[ 0 ] = icdf_ptr[ silk_min( p & 0x1F, 6 ) ]; michael@0: for( j = 0; j < SHELL_CODEC_FRAME_LENGTH; j++ ) { michael@0: if( q_ptr[ j ] > 0 ) { michael@0: /* attach sign */ michael@0: #if 0 michael@0: /* conditional implementation */ michael@0: if( ec_dec_icdf( psRangeDec, icdf, 8 ) == 0 ) { michael@0: q_ptr[ j ] = -q_ptr[ j ]; michael@0: } michael@0: #else michael@0: /* implementation with shift, subtraction, multiplication */ michael@0: q_ptr[ j ] *= silk_dec_map( ec_dec_icdf( psRangeDec, icdf, 8 ) ); michael@0: #endif michael@0: } michael@0: } michael@0: } michael@0: q_ptr += SHELL_CODEC_FRAME_LENGTH; michael@0: } michael@0: }