cannam@154: /* Copyright (c) 2007-2008 CSIRO cannam@154: Copyright (c) 2007-2009 Xiph.Org Foundation cannam@154: Written by Jean-Marc Valin */ cannam@154: /* cannam@154: Redistribution and use in source and binary forms, with or without cannam@154: modification, are permitted provided that the following conditions cannam@154: are met: cannam@154: cannam@154: - Redistributions of source code must retain the above copyright cannam@154: notice, this list of conditions and the following disclaimer. cannam@154: cannam@154: - Redistributions in binary form must reproduce the above copyright cannam@154: notice, this list of conditions and the following disclaimer in the cannam@154: documentation and/or other materials provided with the distribution. cannam@154: cannam@154: THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS cannam@154: ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT cannam@154: LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR cannam@154: A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER cannam@154: OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, cannam@154: EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, cannam@154: PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR cannam@154: PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF cannam@154: LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING cannam@154: NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS cannam@154: SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. cannam@154: */ cannam@154: cannam@154: #ifdef HAVE_CONFIG_H cannam@154: #include "config.h" cannam@154: #endif cannam@154: cannam@154: #include "quant_bands.h" cannam@154: #include "laplace.h" cannam@154: #include cannam@154: #include "os_support.h" cannam@154: #include "arch.h" cannam@154: #include "mathops.h" cannam@154: #include "stack_alloc.h" cannam@154: #include "rate.h" cannam@154: cannam@154: #ifdef FIXED_POINT cannam@154: /* Mean energy in each band quantized in Q4 */ cannam@154: const signed char eMeans[25] = { cannam@154: 103,100, 92, 85, 81, cannam@154: 77, 72, 70, 78, 75, cannam@154: 73, 71, 78, 74, 69, cannam@154: 72, 70, 74, 76, 71, cannam@154: 60, 60, 60, 60, 60 cannam@154: }; cannam@154: #else cannam@154: /* Mean energy in each band quantized in Q4 and converted back to float */ cannam@154: const opus_val16 eMeans[25] = { cannam@154: 6.437500f, 6.250000f, 5.750000f, 5.312500f, 5.062500f, cannam@154: 4.812500f, 4.500000f, 4.375000f, 4.875000f, 4.687500f, cannam@154: 4.562500f, 4.437500f, 4.875000f, 4.625000f, 4.312500f, cannam@154: 4.500000f, 4.375000f, 4.625000f, 4.750000f, 4.437500f, cannam@154: 3.750000f, 3.750000f, 3.750000f, 3.750000f, 3.750000f cannam@154: }; cannam@154: #endif cannam@154: /* prediction coefficients: 0.9, 0.8, 0.65, 0.5 */ cannam@154: #ifdef FIXED_POINT cannam@154: static const opus_val16 pred_coef[4] = {29440, 26112, 21248, 16384}; cannam@154: static const opus_val16 beta_coef[4] = {30147, 22282, 12124, 6554}; cannam@154: static const opus_val16 beta_intra = 4915; cannam@154: #else cannam@154: static const opus_val16 pred_coef[4] = {29440/32768., 26112/32768., 21248/32768., 16384/32768.}; cannam@154: static const opus_val16 beta_coef[4] = {30147/32768., 22282/32768., 12124/32768., 6554/32768.}; cannam@154: static const opus_val16 beta_intra = 4915/32768.; cannam@154: #endif cannam@154: cannam@154: /*Parameters of the Laplace-like probability models used for the coarse energy. cannam@154: There is one pair of parameters for each frame size, prediction type cannam@154: (inter/intra), and band number. cannam@154: The first number of each pair is the probability of 0, and the second is the cannam@154: decay rate, both in Q8 precision.*/ cannam@154: static const unsigned char e_prob_model[4][2][42] = { cannam@154: /*120 sample frames.*/ cannam@154: { cannam@154: /*Inter*/ cannam@154: { cannam@154: 72, 127, 65, 129, 66, 128, 65, 128, 64, 128, 62, 128, 64, 128, cannam@154: 64, 128, 92, 78, 92, 79, 92, 78, 90, 79, 116, 41, 115, 40, cannam@154: 114, 40, 132, 26, 132, 26, 145, 17, 161, 12, 176, 10, 177, 11 cannam@154: }, cannam@154: /*Intra*/ cannam@154: { cannam@154: 24, 179, 48, 138, 54, 135, 54, 132, 53, 134, 56, 133, 55, 132, cannam@154: 55, 132, 61, 114, 70, 96, 74, 88, 75, 88, 87, 74, 89, 66, cannam@154: 91, 67, 100, 59, 108, 50, 120, 40, 122, 37, 97, 43, 78, 50 cannam@154: } cannam@154: }, cannam@154: /*240 sample frames.*/ cannam@154: { cannam@154: /*Inter*/ cannam@154: { cannam@154: 83, 78, 84, 81, 88, 75, 86, 74, 87, 71, 90, 73, 93, 74, cannam@154: 93, 74, 109, 40, 114, 36, 117, 34, 117, 34, 143, 17, 145, 18, cannam@154: 146, 19, 162, 12, 165, 10, 178, 7, 189, 6, 190, 8, 177, 9 cannam@154: }, cannam@154: /*Intra*/ cannam@154: { cannam@154: 23, 178, 54, 115, 63, 102, 66, 98, 69, 99, 74, 89, 71, 91, cannam@154: 73, 91, 78, 89, 86, 80, 92, 66, 93, 64, 102, 59, 103, 60, cannam@154: 104, 60, 117, 52, 123, 44, 138, 35, 133, 31, 97, 38, 77, 45 cannam@154: } cannam@154: }, cannam@154: /*480 sample frames.*/ cannam@154: { cannam@154: /*Inter*/ cannam@154: { cannam@154: 61, 90, 93, 60, 105, 42, 107, 41, 110, 45, 116, 38, 113, 38, cannam@154: 112, 38, 124, 26, 132, 27, 136, 19, 140, 20, 155, 14, 159, 16, cannam@154: 158, 18, 170, 13, 177, 10, 187, 8, 192, 6, 175, 9, 159, 10 cannam@154: }, cannam@154: /*Intra*/ cannam@154: { cannam@154: 21, 178, 59, 110, 71, 86, 75, 85, 84, 83, 91, 66, 88, 73, cannam@154: 87, 72, 92, 75, 98, 72, 105, 58, 107, 54, 115, 52, 114, 55, cannam@154: 112, 56, 129, 51, 132, 40, 150, 33, 140, 29, 98, 35, 77, 42 cannam@154: } cannam@154: }, cannam@154: /*960 sample frames.*/ cannam@154: { cannam@154: /*Inter*/ cannam@154: { cannam@154: 42, 121, 96, 66, 108, 43, 111, 40, 117, 44, 123, 32, 120, 36, cannam@154: 119, 33, 127, 33, 134, 34, 139, 21, 147, 23, 152, 20, 158, 25, cannam@154: 154, 26, 166, 21, 173, 16, 184, 13, 184, 10, 150, 13, 139, 15 cannam@154: }, cannam@154: /*Intra*/ cannam@154: { cannam@154: 22, 178, 63, 114, 74, 82, 84, 83, 92, 82, 103, 62, 96, 72, cannam@154: 96, 67, 101, 73, 107, 72, 113, 55, 118, 52, 125, 52, 118, 52, cannam@154: 117, 55, 135, 49, 137, 39, 157, 32, 145, 29, 97, 33, 77, 40 cannam@154: } cannam@154: } cannam@154: }; cannam@154: cannam@154: static const unsigned char small_energy_icdf[3]={2,1,0}; cannam@154: cannam@154: static opus_val32 loss_distortion(const opus_val16 *eBands, opus_val16 *oldEBands, int start, int end, int len, int C) cannam@154: { cannam@154: int c, i; cannam@154: opus_val32 dist = 0; cannam@154: c=0; do { cannam@154: for (i=start;inbEBands]; cannam@154: oldE = MAX16(-QCONST16(9.f,DB_SHIFT), oldEBands[i+c*m->nbEBands]); cannam@154: #ifdef FIXED_POINT cannam@154: f = SHL32(EXTEND32(x),7) - PSHR32(MULT16_16(coef,oldE), 8) - prev[c]; cannam@154: /* Rounding to nearest integer here is really important! */ cannam@154: qi = (f+QCONST32(.5f,DB_SHIFT+7))>>(DB_SHIFT+7); cannam@154: decay_bound = EXTRACT16(MAX32(-QCONST16(28.f,DB_SHIFT), cannam@154: SUB32((opus_val32)oldEBands[i+c*m->nbEBands],max_decay))); cannam@154: #else cannam@154: f = x-coef*oldE-prev[c]; cannam@154: /* Rounding to nearest integer here is really important! */ cannam@154: qi = (int)floor(.5f+f); cannam@154: decay_bound = MAX16(-QCONST16(28.f,DB_SHIFT), oldEBands[i+c*m->nbEBands]) - max_decay; cannam@154: #endif cannam@154: /* Prevent the energy from going down too quickly (e.g. for bands cannam@154: that have just one bin) */ cannam@154: if (qi < 0 && x < decay_bound) cannam@154: { cannam@154: qi += (int)SHR16(SUB16(decay_bound,x), DB_SHIFT); cannam@154: if (qi > 0) cannam@154: qi = 0; cannam@154: } cannam@154: qi0 = qi; cannam@154: /* If we don't have enough bits to encode all the energy, just assume cannam@154: something safe. */ cannam@154: tell = ec_tell(enc); cannam@154: bits_left = budget-tell-3*C*(end-i); cannam@154: if (i!=start && bits_left < 30) cannam@154: { cannam@154: if (bits_left < 24) cannam@154: qi = IMIN(1, qi); cannam@154: if (bits_left < 16) cannam@154: qi = IMAX(-1, qi); cannam@154: } cannam@154: if (lfe && i>=2) cannam@154: qi = IMIN(qi, 0); cannam@154: if (budget-tell >= 15) cannam@154: { cannam@154: int pi; cannam@154: pi = 2*IMIN(i,20); cannam@154: ec_laplace_encode(enc, &qi, cannam@154: prob_model[pi]<<7, prob_model[pi+1]<<6); cannam@154: } cannam@154: else if(budget-tell >= 2) cannam@154: { cannam@154: qi = IMAX(-1, IMIN(qi, 1)); cannam@154: ec_enc_icdf(enc, 2*qi^-(qi<0), small_energy_icdf, 2); cannam@154: } cannam@154: else if(budget-tell >= 1) cannam@154: { cannam@154: qi = IMIN(0, qi); cannam@154: ec_enc_bit_logp(enc, -qi, 1); cannam@154: } cannam@154: else cannam@154: qi = -1; cannam@154: error[i+c*m->nbEBands] = PSHR32(f,7) - SHL16(qi,DB_SHIFT); cannam@154: badness += abs(qi0-qi); cannam@154: q = (opus_val32)SHL32(EXTEND32(qi),DB_SHIFT); cannam@154: cannam@154: tmp = PSHR32(MULT16_16(coef,oldE),8) + prev[c] + SHL32(q,7); cannam@154: #ifdef FIXED_POINT cannam@154: tmp = MAX32(-QCONST32(28.f, DB_SHIFT+7), tmp); cannam@154: #endif cannam@154: oldEBands[i+c*m->nbEBands] = PSHR32(tmp, 7); cannam@154: prev[c] = prev[c] + SHL32(q,7) - MULT16_16(beta,PSHR32(q,8)); cannam@154: } while (++c < C); cannam@154: } cannam@154: return lfe ? 0 : badness; cannam@154: } cannam@154: cannam@154: void quant_coarse_energy(const CELTMode *m, int start, int end, int effEnd, cannam@154: const opus_val16 *eBands, opus_val16 *oldEBands, opus_uint32 budget, cannam@154: opus_val16 *error, ec_enc *enc, int C, int LM, int nbAvailableBytes, cannam@154: int force_intra, opus_val32 *delayedIntra, int two_pass, int loss_rate, int lfe) cannam@154: { cannam@154: int intra; cannam@154: opus_val16 max_decay; cannam@154: VARDECL(opus_val16, oldEBands_intra); cannam@154: VARDECL(opus_val16, error_intra); cannam@154: ec_enc enc_start_state; cannam@154: opus_uint32 tell; cannam@154: int badness1=0; cannam@154: opus_int32 intra_bias; cannam@154: opus_val32 new_distortion; cannam@154: SAVE_STACK; cannam@154: cannam@154: intra = force_intra || (!two_pass && *delayedIntra>2*C*(end-start) && nbAvailableBytes > (end-start)*C); cannam@154: intra_bias = (opus_int32)((budget**delayedIntra*loss_rate)/(C*512)); cannam@154: new_distortion = loss_distortion(eBands, oldEBands, start, effEnd, m->nbEBands, C); cannam@154: cannam@154: tell = ec_tell(enc); cannam@154: if (tell+3 > budget) cannam@154: two_pass = intra = 0; cannam@154: cannam@154: max_decay = QCONST16(16.f,DB_SHIFT); cannam@154: if (end-start>10) cannam@154: { cannam@154: #ifdef FIXED_POINT cannam@154: max_decay = MIN32(max_decay, SHL32(EXTEND32(nbAvailableBytes),DB_SHIFT-3)); cannam@154: #else cannam@154: max_decay = MIN32(max_decay, .125f*nbAvailableBytes); cannam@154: #endif cannam@154: } cannam@154: if (lfe) cannam@154: max_decay = QCONST16(3.f,DB_SHIFT); cannam@154: enc_start_state = *enc; cannam@154: cannam@154: ALLOC(oldEBands_intra, C*m->nbEBands, opus_val16); cannam@154: ALLOC(error_intra, C*m->nbEBands, opus_val16); cannam@154: OPUS_COPY(oldEBands_intra, oldEBands, C*m->nbEBands); cannam@154: cannam@154: if (two_pass || intra) cannam@154: { cannam@154: badness1 = quant_coarse_energy_impl(m, start, end, eBands, oldEBands_intra, budget, cannam@154: tell, e_prob_model[LM][1], error_intra, enc, C, LM, 1, max_decay, lfe); cannam@154: } cannam@154: cannam@154: if (!intra) cannam@154: { cannam@154: unsigned char *intra_buf; cannam@154: ec_enc enc_intra_state; cannam@154: opus_int32 tell_intra; cannam@154: opus_uint32 nstart_bytes; cannam@154: opus_uint32 nintra_bytes; cannam@154: opus_uint32 save_bytes; cannam@154: int badness2; cannam@154: VARDECL(unsigned char, intra_bits); cannam@154: cannam@154: tell_intra = ec_tell_frac(enc); cannam@154: cannam@154: enc_intra_state = *enc; cannam@154: cannam@154: nstart_bytes = ec_range_bytes(&enc_start_state); cannam@154: nintra_bytes = ec_range_bytes(&enc_intra_state); cannam@154: intra_buf = ec_get_buffer(&enc_intra_state) + nstart_bytes; cannam@154: save_bytes = nintra_bytes-nstart_bytes; cannam@154: if (save_bytes == 0) cannam@154: save_bytes = ALLOC_NONE; cannam@154: ALLOC(intra_bits, save_bytes, unsigned char); cannam@154: /* Copy bits from intra bit-stream */ cannam@154: OPUS_COPY(intra_bits, intra_buf, nintra_bytes - nstart_bytes); cannam@154: cannam@154: *enc = enc_start_state; cannam@154: cannam@154: badness2 = quant_coarse_energy_impl(m, start, end, eBands, oldEBands, budget, cannam@154: tell, e_prob_model[LM][intra], error, enc, C, LM, 0, max_decay, lfe); cannam@154: cannam@154: if (two_pass && (badness1 < badness2 || (badness1 == badness2 && ((opus_int32)ec_tell_frac(enc))+intra_bias > tell_intra))) cannam@154: { cannam@154: *enc = enc_intra_state; cannam@154: /* Copy intra bits to bit-stream */ cannam@154: OPUS_COPY(intra_buf, intra_bits, nintra_bytes - nstart_bytes); cannam@154: OPUS_COPY(oldEBands, oldEBands_intra, C*m->nbEBands); cannam@154: OPUS_COPY(error, error_intra, C*m->nbEBands); cannam@154: intra = 1; cannam@154: } cannam@154: } else { cannam@154: OPUS_COPY(oldEBands, oldEBands_intra, C*m->nbEBands); cannam@154: OPUS_COPY(error, error_intra, C*m->nbEBands); cannam@154: } cannam@154: cannam@154: if (intra) cannam@154: *delayedIntra = new_distortion; cannam@154: else cannam@154: *delayedIntra = ADD32(MULT16_32_Q15(MULT16_16_Q15(pred_coef[LM], pred_coef[LM]),*delayedIntra), cannam@154: new_distortion); cannam@154: cannam@154: RESTORE_STACK; cannam@154: } cannam@154: cannam@154: void quant_fine_energy(const CELTMode *m, int start, int end, opus_val16 *oldEBands, opus_val16 *error, int *fine_quant, ec_enc *enc, int C) cannam@154: { cannam@154: int i, c; cannam@154: cannam@154: /* Encode finer resolution */ cannam@154: for (i=start;inbEBands]+QCONST16(.5f,DB_SHIFT))>>(DB_SHIFT-fine_quant[i]); cannam@154: #else cannam@154: q2 = (int)floor((error[i+c*m->nbEBands]+.5f)*frac); cannam@154: #endif cannam@154: if (q2 > frac-1) cannam@154: q2 = frac-1; cannam@154: if (q2<0) cannam@154: q2 = 0; cannam@154: ec_enc_bits(enc, q2, fine_quant[i]); cannam@154: #ifdef FIXED_POINT cannam@154: offset = SUB16(SHR32(SHL32(EXTEND32(q2),DB_SHIFT)+QCONST16(.5f,DB_SHIFT),fine_quant[i]),QCONST16(.5f,DB_SHIFT)); cannam@154: #else cannam@154: offset = (q2+.5f)*(1<<(14-fine_quant[i]))*(1.f/16384) - .5f; cannam@154: #endif cannam@154: oldEBands[i+c*m->nbEBands] += offset; cannam@154: error[i+c*m->nbEBands] -= offset; cannam@154: /*printf ("%f ", error[i] - offset);*/ cannam@154: } while (++c < C); cannam@154: } cannam@154: } cannam@154: cannam@154: void quant_energy_finalise(const CELTMode *m, int start, int end, opus_val16 *oldEBands, opus_val16 *error, int *fine_quant, int *fine_priority, int bits_left, ec_enc *enc, int C) cannam@154: { cannam@154: int i, prio, c; cannam@154: cannam@154: /* Use up the remaining bits */ cannam@154: for (prio=0;prio<2;prio++) cannam@154: { cannam@154: for (i=start;i=C ;i++) cannam@154: { cannam@154: if (fine_quant[i] >= MAX_FINE_BITS || fine_priority[i]!=prio) cannam@154: continue; cannam@154: c=0; cannam@154: do { cannam@154: int q2; cannam@154: opus_val16 offset; cannam@154: q2 = error[i+c*m->nbEBands]<0 ? 0 : 1; cannam@154: ec_enc_bits(enc, q2, 1); cannam@154: #ifdef FIXED_POINT cannam@154: offset = SHR16(SHL16(q2,DB_SHIFT)-QCONST16(.5f,DB_SHIFT),fine_quant[i]+1); cannam@154: #else cannam@154: offset = (q2-.5f)*(1<<(14-fine_quant[i]-1))*(1.f/16384); cannam@154: #endif cannam@154: oldEBands[i+c*m->nbEBands] += offset; cannam@154: error[i+c*m->nbEBands] -= offset; cannam@154: bits_left--; cannam@154: } while (++c < C); cannam@154: } cannam@154: } cannam@154: } cannam@154: cannam@154: void unquant_coarse_energy(const CELTMode *m, int start, int end, opus_val16 *oldEBands, int intra, ec_dec *dec, int C, int LM) cannam@154: { cannam@154: const unsigned char *prob_model = e_prob_model[LM][intra]; cannam@154: int i, c; cannam@154: opus_val32 prev[2] = {0, 0}; cannam@154: opus_val16 coef; cannam@154: opus_val16 beta; cannam@154: opus_int32 budget; cannam@154: opus_int32 tell; cannam@154: cannam@154: if (intra) cannam@154: { cannam@154: coef = 0; cannam@154: beta = beta_intra; cannam@154: } else { cannam@154: beta = beta_coef[LM]; cannam@154: coef = pred_coef[LM]; cannam@154: } cannam@154: cannam@154: budget = dec->storage*8; cannam@154: cannam@154: /* Decode at a fixed coarse resolution */ cannam@154: for (i=start;i=15) cannam@154: { cannam@154: int pi; cannam@154: pi = 2*IMIN(i,20); cannam@154: qi = ec_laplace_decode(dec, cannam@154: prob_model[pi]<<7, prob_model[pi+1]<<6); cannam@154: } cannam@154: else if(budget-tell>=2) cannam@154: { cannam@154: qi = ec_dec_icdf(dec, small_energy_icdf, 2); cannam@154: qi = (qi>>1)^-(qi&1); cannam@154: } cannam@154: else if(budget-tell>=1) cannam@154: { cannam@154: qi = -ec_dec_bit_logp(dec, 1); cannam@154: } cannam@154: else cannam@154: qi = -1; cannam@154: q = (opus_val32)SHL32(EXTEND32(qi),DB_SHIFT); cannam@154: cannam@154: oldEBands[i+c*m->nbEBands] = MAX16(-QCONST16(9.f,DB_SHIFT), oldEBands[i+c*m->nbEBands]); cannam@154: tmp = PSHR32(MULT16_16(coef,oldEBands[i+c*m->nbEBands]),8) + prev[c] + SHL32(q,7); cannam@154: #ifdef FIXED_POINT cannam@154: tmp = MAX32(-QCONST32(28.f, DB_SHIFT+7), tmp); cannam@154: #endif cannam@154: oldEBands[i+c*m->nbEBands] = PSHR32(tmp, 7); cannam@154: prev[c] = prev[c] + SHL32(q,7) - MULT16_16(beta,PSHR32(q,8)); cannam@154: } while (++c < C); cannam@154: } cannam@154: } cannam@154: cannam@154: void unquant_fine_energy(const CELTMode *m, int start, int end, opus_val16 *oldEBands, int *fine_quant, ec_dec *dec, int C) cannam@154: { cannam@154: int i, c; cannam@154: /* Decode finer resolution */ cannam@154: for (i=start;inbEBands] += offset; cannam@154: } while (++c < C); cannam@154: } cannam@154: } cannam@154: cannam@154: void unquant_energy_finalise(const CELTMode *m, int start, int end, opus_val16 *oldEBands, int *fine_quant, int *fine_priority, int bits_left, ec_dec *dec, int C) cannam@154: { cannam@154: int i, prio, c; cannam@154: cannam@154: /* Use up the remaining bits */ cannam@154: for (prio=0;prio<2;prio++) cannam@154: { cannam@154: for (i=start;i=C ;i++) cannam@154: { cannam@154: if (fine_quant[i] >= MAX_FINE_BITS || fine_priority[i]!=prio) cannam@154: continue; cannam@154: c=0; cannam@154: do { cannam@154: int q2; cannam@154: opus_val16 offset; cannam@154: q2 = ec_dec_bits(dec, 1); cannam@154: #ifdef FIXED_POINT cannam@154: offset = SHR16(SHL16(q2,DB_SHIFT)-QCONST16(.5f,DB_SHIFT),fine_quant[i]+1); cannam@154: #else cannam@154: offset = (q2-.5f)*(1<<(14-fine_quant[i]-1))*(1.f/16384); cannam@154: #endif cannam@154: oldEBands[i+c*m->nbEBands] += offset; cannam@154: bits_left--; cannam@154: } while (++c < C); cannam@154: } cannam@154: } cannam@154: } cannam@154: cannam@154: void amp2Log2(const CELTMode *m, int effEnd, int end, cannam@154: celt_ener *bandE, opus_val16 *bandLogE, int C) cannam@154: { cannam@154: int c, i; cannam@154: c=0; cannam@154: do { cannam@154: for (i=0;inbEBands] = cannam@154: celt_log2(bandE[i+c*m->nbEBands]) cannam@154: - SHL16((opus_val16)eMeans[i],6); cannam@154: #ifdef FIXED_POINT cannam@154: /* Compensate for bandE[] being Q12 but celt_log2() taking a Q14 input. */ cannam@154: bandLogE[i+c*m->nbEBands] += QCONST16(2.f, DB_SHIFT); cannam@154: #endif cannam@154: } cannam@154: for (i=effEnd;inbEBands+i] = -QCONST16(14.f,DB_SHIFT); cannam@154: } while (++c < C); cannam@154: }