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