index 57c7709..5d85f55 100644 (file)
@@ -1,5 +1,6 @@
-/* (C) 2007-2008 Jean-Marc Valin, CSIRO
-*/
+   Copyright (c) 2007-2009 Xiph.Org Foundation
+   Written by Jean-Marc Valin */
/*
Redistribution and use in source and binary forms, with or without
modification, are permitted provided that the following conditions
/*
Redistribution and use in source and binary forms, with or without
modification, are permitted provided that the following conditions
#include "entcode.h"
#include "rate.h"

#include "entcode.h"
#include "rate.h"

-#define BITRES 4
-#define BITROUND 8
-#define BITOVERFLOW 10000
+
+static const unsigned char LOG2_FRAC_TABLE={
+   0,
+   8,13,
+  16,19,21,23,
+  24,26,27,28,29,30,31,32,
+  32,33,34,34,35,36,36,37,37
+};

#ifndef STATIC_MODES

#ifndef STATIC_MODES
-static int log2_frac(ec_uint32 val, int frac)
+
+/*Determines if V(N,K) fits in a 32-bit unsigned integer.
+  N and K are themselves limited to 15 bits.*/
+static int fits_in32(int _n, int _k)
{
{
-   int i;
-   /* EC_ILOG() actually returns log2()+1, go figure */
-   int L = EC_ILOG(val)-1;
-   /*printf ("in: %d %d ", val, L);*/
-   if (L>14)
-      val >>= L-14;
-   else if (L<14)
-      val <<= 14-L;
-   L <<= frac;
-   /*printf ("%d\n", val);*/
-   for (i=0;i<frac;i++)
+   static const celt_int16 maxN = {
+      32767, 32767, 32767, 1476, 283, 109,  60,  40,
+       29,  24,  20,  18,  16,  14,  13};
+   static const celt_int16 maxK = {
+      32767, 32767, 32767, 32767, 1172, 238,  95,  53,
+       36,  27,  22,  18,  16,  15,  13};
+   if (_n>=14)
{
{
-      val = (val*val) >> 15;
-      /*printf ("%d\n", val);*/
-      if (val > 16384)
-         L |= (1<<(frac-i-1));
-      else
-         val <<= 1;
+      if (_k>=14)
+         return 0;
+      else
+         return _n <= maxN[_k];
+   } else {
+      return _k <= maxK[_n];
}
}
-   return L;
}

}

-static int log2_frac64(ec_uint64 val, int frac)
+void compute_pulse_cache(CELTMode *m, int LM)
{
int i;
{
int i;
-   /* EC_ILOG64() actually returns log2()+1, go figure */
-   int L = EC_ILOG64(val)-1;
-   /*printf ("in: %d %d ", val, L);*/
-   if (L>14)
-      val >>= L-14;
-   else if (L<14)
-      val <<= 14-L;
-   L <<= frac;
-   /*printf ("%d\n", val);*/
-   for (i=0;i<frac;i++)
-   {
-      val = (val*val) >> 15;
-      /*printf ("%d\n", val);*/
-      if (val > 16384)
-         L |= (1<<(frac-i-1));
-      else
-         val <<= 1;
-   }
-   return L;
-}
+   int curr=0;
+   int nbEntries=0;
+   int entryN, entryK, entryI;
+   const celt_int16 *eBands = m->eBands;
+   PulseCache *cache = &m->cache;
+   celt_int16 *cindex;
+   unsigned char *bits;

-void compute_alloc_cache(CELTMode *m)
-{
-   int i, prevN, BC;
-   celt_int16_t **bits;
-   const celt_int16_t *eBands = m->eBands;
+   cindex = celt_alloc(sizeof(cache->index)*m->nbEBands*(LM+2));
+   cache->index = cindex;

-   bits = celt_alloc(m->nbEBands*sizeof(celt_int16_t*));
-
-   BC = m->nbChannels;
-   prevN = -1;
-   for (i=0;i<m->nbEBands;i++)
+   /* Scan for all unique band sizes */
+   for (i=0;i<=LM+1;i++)
{
{
-      int N = BC*(eBands[i+1]-eBands[i]);
-      if (N == prevN && eBands[i] < m->pitchEnd)
+      int j;
+      for (j=0;j<m->nbEBands;j++)
{
{
-         bits[i] = bits[i-1];
-      } else {
-         int j;
-         VARDECL(celt_uint64_t, u);
-         SAVE_STACK;
-         ALLOC(u, N, celt_uint64_t);
-         /* FIXME: We could save memory here */
-         bits[i] = celt_alloc(MAX_PULSES*sizeof(celt_int16_t));
-         for (j=0;j<MAX_PULSES;j++)
+         int k;
+         int N = (eBands[j+1]-eBands[j])<<i>>1;
+         cindex[i*m->nbEBands+j] = -1;
+         /* Find other bands that have the same size */
+         for (k=0;k<=i;k++)
{
{
-            int done = 0;
-            int pulses = j;
-            /* For bands where there's no pitch, id 1 corresponds to intra prediction
-            with no pulse. id 2 means intra prediction with one pulse, and so on.*/
-            if (eBands[i] >= m->pitchEnd)
-               pulses -= 1;
-            if (pulses < 0)
-               bits[i][j] = 0;
-            else {
-               celt_uint64_t nc;
-               nc=pulses?ncwrs_unext64(N, u):ncwrs_u64(N, 0, u);
-               bits[i][j] = log2_frac64(nc,BITRES);
-               /* FIXME: Could there be a better test for the max number of pulses that fit in 64 bits? */
-               if (bits[i][j] > (60<<BITRES))
-                  done = 1;
-               /* Add the intra-frame prediction bits */
-               if (eBands[i] >= m->pitchEnd)
+            int n;
+            for (n=0;n<m->nbEBands && (k!=i || n<j);n++)
+            {
+               if (N == (eBands[n+1]-eBands[n])<<k>>1)
{
{
-                  int max_pos = 2*eBands[i]-eBands[i+1];
-                  if (max_pos > 32)
-                     max_pos = 32;
-                  bits[i][j] += (1<<BITRES) + log2_frac(max_pos,BITRES);
+                  cindex[i*m->nbEBands+j] = cindex[k*m->nbEBands+n];
+                  break;
}
}
}
}
-            if (done)
-               break;
}
}
-         for (;j<MAX_PULSES;j++)
-            bits[i][j] = BITOVERFLOW;
-         prevN = N;
-         RESTORE_STACK;
+         if (cache->index[i*m->nbEBands+j] == -1 && N!=0)
+         {
+            int K;
+            entryN[nbEntries] = N;
+            K = 0;
+            while (fits_in32(N,get_pulses(K+1)) && K<MAX_PSEUDO)
+               K++;
+            entryK[nbEntries] = K;
+            cindex[i*m->nbEBands+j] = curr;
+            entryI[nbEntries] = curr;
+
+            curr += K+1;
+            nbEntries++;
+         }
}
}
}
}
-   m->bits = (const celt_int16_t * const *)bits;
-}
-
-#endif /* !STATIC_MODES */
-
-static inline int bits2pulses(const CELTMode *m, int band, int bits)
-{
-   int i;
-   int lo, hi;
-   lo = 0;
-   hi = MAX_PULSES-1;
-
-   /* Instead of using the "bisection condition" we use a fixed number of
-      iterations because it should be faster */
-   /*while (hi-lo != 1)*/
-   for (i=0;i<LOG_MAX_PULSES;i++)
+   bits = celt_alloc(sizeof(unsigned char)*curr);
+   cache->bits = bits;
+   cache->size = curr;
+   /* Compute the cache for all unique sizes */
+   for (i=0;i<nbEntries;i++)
{
{
-      int mid = (lo+hi)>>1;
-      /* OPT: Make sure this is implemented with a conditional move */
-      if (m->bits[band][mid] >= bits)
-         hi = mid;
-      else
-         lo = mid;
+      int j;
+      unsigned char *ptr = bits+entryI[i];
+      celt_int16 tmp[MAX_PULSES+1];
+      get_required_bits(tmp, entryN[i], get_pulses(entryK[i]), BITRES);
+      for (j=1;j<=entryK[i];j++)
+         ptr[j] = tmp[get_pulses(j)]-1;
+      ptr = entryK[i];
}
}
-   if (bits-m->bits[band][lo] <= m->bits[band][hi]-bits)
-      return lo;
-   else
-      return hi;
}

}

-static int vec_bits2pulses(const CELTMode *m, int *bits, int *pulses, int len)
-{
-   int i;
-   int sum=0;
+#endif /* !STATIC_MODES */

-   for (i=0;i<len;i++)
-   {
-      pulses[i] = bits2pulses(m, i, bits[i]);
-      sum += m->bits[i][pulses[i]];
-   }
-   /*printf ("sum = %d\n", sum);*/
-   return sum;
-}

-static int interp_bits2pulses(const CELTMode *m, int *bits1, int *bits2, int total, int *pulses, int len)
+#define ALLOC_STEPS 6
+
+static inline int interp_bits2pulses(const CELTMode *m, int start, int end, int skip_start,
+      const int *bits1, const int *bits2, const int *thresh, int total, int skip_rsv,
+      int *intensity, int intensity_rsv, int *dual_stereo, int dual_stereo_rsv, int *bits,
+      int *ebits, int *fine_priority, int _C, int LM, void *ec, int encode, int prev)
{
{
-   int lo, hi, out;
-   int j;
-   VARDECL(int, bits);
+   int psum;
+   int lo, hi;
+   int i, j;
+   int logM;
+   const int C = CHANNELS(_C);
+   int stereo;
+   int codedBands=-1;
+   int alloc_floor;
+   int left, percoeff;
+   int done;
+   int balance;
SAVE_STACK;
SAVE_STACK;
-   ALLOC(bits, len, int);
+
+   alloc_floor = C<<BITRES;
+   stereo = C>1;
+
+   logM = LM<<BITRES;
lo = 0;
lo = 0;
-   hi = 1<<BITRES;
-   while (hi-lo != 1)
+   hi = 1<<ALLOC_STEPS;
+   for (i=0;i<ALLOC_STEPS;i++)
{
int mid = (lo+hi)>>1;
{
int mid = (lo+hi)>>1;
-      for (j=0;j<len;j++)
-         bits[j] = ((1<<BITRES)-mid)*bits1[j] + mid*bits2[j];
-      if (vec_bits2pulses(m, bits, pulses, len) > total<<BITRES)
+      psum = 0;
+      done = 0;
+      for (j=end;j-->start;)
+      {
+         int tmp = bits1[j] + (mid*bits2[j]>>ALLOC_STEPS);
+         if (tmp >= thresh[j] || done)
+         {
+            done = 1;
+            /* Don't allocate more than we can actually use */
+            psum += IMIN(tmp, 64*C<<BITRES<<LM);
+         } else {
+            if (tmp >= alloc_floor)
+               psum += alloc_floor;
+         }
+      }
+      if (psum > total)
hi = mid;
else
lo = mid;
}
hi = mid;
else
lo = mid;
}
+   psum = 0;
/*printf ("interp bisection gave %d\n", lo);*/
/*printf ("interp bisection gave %d\n", lo);*/
-   for (j=0;j<len;j++)
-      bits[j] = ((1<<BITRES)-lo)*bits1[j] + lo*bits2[j];
-   out = vec_bits2pulses(m, bits, pulses, len);
-   /* Do some refinement to use up all bits. In the first pass, we can only add pulses to
-      bands that are under their allocated budget. In the second pass, anything goes */
-   for (j=0;j<len;j++)
+   done = 0;
+   for (j=end;j-->start;)
+   {
+      int tmp = bits1[j] + (lo*bits2[j]>>ALLOC_STEPS);
+      if (tmp < thresh[j] && !done)
+      {
+         if (tmp >= alloc_floor)
+            tmp = alloc_floor;
+         else
+            tmp = 0;
+      } else
+         done = 1;
+      /* Don't allocate more than we can actually use */
+      tmp = IMIN(tmp, 64*C<<BITRES<<LM);
+      bits[j] = tmp;
+      psum += tmp;
+   }
+
+   /* Decide which bands to skip, working backwards from the end. */
+   for (codedBands=end;;codedBands--)
{
{
-      if (m->bits[j][pulses[j]] < bits[j] && pulses[j]<MAX_PULSES-1)
+      int band_width;
+      int band_bits;
+      int rem;
+      j = codedBands-1;
+      /* Never skip the first band, nor a band that has been boosted by
+          dynalloc.
+         In the first case, we'd be coding a bit to signal we're going to waste
+          all the other bits.
+         In the second case, we'd be coding a bit to redistribute all the bits
+          we just signaled should be cocentrated in this band. */
+      if (j<=skip_start)
+      {
+         /* Give the bit we reserved to end skipping back. */
+         total += skip_rsv;
+         break;
+      }
+      /*Figure out how many left-over bits we would be adding to this band.
+        This can include bits we've stolen back from higher, skipped bands.*/
+      left = total-psum;
+      percoeff = left/(m->eBands[codedBands]-m->eBands[start]);
+      left -= (m->eBands[codedBands]-m->eBands[start])*percoeff;
+      rem = IMAX(left-(m->eBands[j]-m->eBands[start]),0);
+      band_width = m->eBands[codedBands]-m->eBands[j];
+      band_bits = bits[j] + percoeff*band_width + rem;
+      /*Only code a skip decision if we're above the threshold for this band.
+        Otherwise it is force-skipped.
+        This ensures that we have enough bits to code the skip flag.*/
+      if (band_bits >= IMAX(thresh[j], alloc_floor+(1<<BITRES)))
{
{
-         if (out+m->bits[j][pulses[j]+1]-m->bits[j][pulses[j]] <= total<<BITRES)
+         if (encode)
{
{
-            out = out+m->bits[j][pulses[j]+1]-m->bits[j][pulses[j]];
-            pulses[j] += 1;
+            /*This if() block is the only part of the allocation function that
+               is not a mandatory part of the bitstream: any bands we choose to
+               skip here must be explicitly signaled.*/
+            /*Choose a threshold with some hysteresis to keep bands from
+               fluctuating in and out.*/
+            if (band_bits > ((j<prev?7:9)*band_width<<LM<<BITRES)>>4)
+            {
+               ec_enc_bit_logp((ec_enc *)ec, 1, 1);
+               break;
+            }
+            ec_enc_bit_logp((ec_enc *)ec, 0, 1);
+         } else if (ec_dec_bit_logp((ec_dec *)ec, 1)) {
+            break;
}
}
+         /*We used a bit to skip this band.*/
+         psum += 1<<BITRES;
+         band_bits -= 1<<BITRES;
}
}
+      /*Reclaim the bits originally allocated to this band.*/
+      psum -= bits[j]+intensity_rsv;
+      if (intensity_rsv > 0)
+         intensity_rsv = LOG2_FRAC_TABLE[j-start];
+      psum += intensity_rsv;
+      if (band_bits >= alloc_floor)
+      {
+         /*If we have enough for a fine energy bit per channel, use it.*/
+         psum += alloc_floor;
+         bits[j] = alloc_floor;
+      } else {
+         /*Otherwise this band gets nothing at all.*/
+         bits[j] = 0;
+      }
+   }
+
+   celt_assert(codedBands > start);
+   /* Code the intensity and dual stereo parameters. */
+   if (intensity_rsv > 0)
+   {
+      if (encode)
+      {
+         *intensity = IMIN(*intensity, codedBands);
+         ec_enc_uint((ec_enc *)ec, *intensity-start, codedBands+1-start);
+      }
+      else
+         *intensity = start+ec_dec_uint((ec_dec *)ec, codedBands+1-start);
+   }
+   else
+      *intensity = 0;
+   if (*intensity <= start)
+   {
+      total += dual_stereo_rsv;
+      dual_stereo_rsv = 0;
+   }
+   if (dual_stereo_rsv > 0)
+   {
+      if (encode)
+         ec_enc_bit_logp((ec_enc *)ec, *dual_stereo, 1);
+      else
+         *dual_stereo = ec_dec_bit_logp((ec_dec *)ec, 1);
}
}
-   while(1)
+   else
+      *dual_stereo = 0;
+
+   /* Allocate the remaining bits */
+   left = total-psum;
+   percoeff = left/(m->eBands[codedBands]-m->eBands[start]);
+   left -= (m->eBands[codedBands]-m->eBands[start])*percoeff;
+   for (j=start;j<codedBands;j++)
+      bits[j] += percoeff*(m->eBands[j+1]-m->eBands[j]);
+   for (j=start;j<codedBands;j++)
{
{
-      int incremented = 0;
-      for (j=0;j<len;j++)
+      int tmp = IMIN(left, m->eBands[j+1]-m->eBands[j]);
+      bits[j] += tmp;
+      left -= tmp;
+   }
+   /*for (j=0;j<end;j++)printf("%d ", bits[j]);printf("\n");*/
+
+   balance = 0;
+   for (j=start;j<codedBands;j++)
+   {
+      int N0, N, den;
+      int offset;
+      int NClogN;
+
+      celt_assert(bits[j] >= 0);
+      N0 = m->eBands[j+1]-m->eBands[j];
+      N=N0<<LM;
+
+      if (N>1)
{
{
-         if (pulses[j]<MAX_PULSES-1)
+         NClogN = N*C*(m->logN[j] + logM);
+
+         /* Compensate for the extra DoF in stereo */
+         den=(C*N+ ((C==2 && N>2) ? 1 : 0));
+
+         /* Offset for the number of fine bits by log2(N)/2 + FINE_OFFSET
+            compared to their "fair share" of total/N */
+         offset = (NClogN>>1)-N*C*FINE_OFFSET;
+
+         /* N=2 is the only point that doesn't match the curve */
+         if (N==2)
+            offset += N*C<<BITRES>>2;
+
+         /* Changing the offset for allocating the second and third
+             fine energy bit */
+         if (bits[j] + offset < den*2<<BITRES)
+            offset += NClogN>>2;
+         else if (bits[j] + offset < den*3<<BITRES)
+            offset += NClogN>>3;
+
+         /* Divide with rounding */
+         ebits[j] = IMAX(0, (bits[j] + offset + (den<<(BITRES-1))) / (den<<BITRES));
+
+         /* Make sure not to bust */
+         if (C*ebits[j] > (bits[j]>>BITRES))
+            ebits[j] = bits[j] >> stereo >> BITRES;
+
+         /* More than that is useless because that's about as far as PVQ can go */
+         ebits[j] = IMIN(ebits[j], MAX_FINE_BITS);
+
+         /* If we rounded down or capped this band, make it a candidate for the
+             final fine energy pass */
+         fine_priority[j] = ebits[j]*(den<<BITRES) >= bits[j]+offset;
+
+      } else {
+         /* For N=1, all bits go to fine energy except for a single sign bit */
+         ebits[j] = IMIN(IMAX(0,(bits[j] >> stereo >> BITRES)-1),MAX_FINE_BITS);
+         fine_priority[j] = (ebits[j]+1)*C<<BITRES >= (bits[j]-balance);
+         /* N=1 bands can't take advantage of the re-balancing in
+             quant_all_bands() because they don't have shape, only fine energy.
+            Instead, do the re-balancing here.*/
+         balance = IMAX(0,bits[j] - ((ebits[j]+1)*C<<BITRES));
+         if (j+1<codedBands)
{
{
-            if (out+m->bits[j][pulses[j]+1]-m->bits[j][pulses[j]] <= total<<BITRES)
-            {
-               out = out+m->bits[j][pulses[j]+1]-m->bits[j][pulses[j]];
-               pulses[j] += 1;
-               incremented = 1;
-            }
+            bits[j] -= balance;
+            bits[j+1] += balance;
}
}
}
}
-      if (!incremented)
-            break;
+
+      /* The other bits are assigned to PVQ */
+      bits[j] -= C*ebits[j]<<BITRES;
+      celt_assert(bits[j] >= 0);
+      celt_assert(ebits[j] >= 0);
+   }
+   /* The skipped bands use all their bits for fine energy. */
+   for (;j<end;j++)
+   {
+      ebits[j] = bits[j] >> stereo >> BITRES;
+      celt_assert(C*ebits[j]<<BITRES == bits[j]);
+      bits[j] = 0;
+      fine_priority[j] = ebits[j]<1;
}
RESTORE_STACK;
}
RESTORE_STACK;
-   return (out+BITROUND) >> BITRES;
+   return codedBands;
}

}

-int compute_allocation(const CELTMode *m, int *offsets, int total, int *pulses)
+int compute_allocation(const CELTMode *m, int start, int end, const int *offsets, int alloc_trim, int *intensity, int *dual_stereo,
+      int total, int *pulses, int *ebits, int *fine_priority, int _C, int LM, void *ec, int encode, int prev)
{
{
-   int lo, hi, len, ret;
+   int lo, hi, len, j;
+   const int C = CHANNELS(_C);
+   int codedBands;
+   int skip_start;
+   int skip_rsv;
+   int intensity_rsv;
+   int dual_stereo_rsv;
VARDECL(int, bits1);
VARDECL(int, bits2);
VARDECL(int, bits1);
VARDECL(int, bits2);
+   VARDECL(int, thresh);
+   VARDECL(int, trim_offset);
SAVE_STACK;

SAVE_STACK;

+   total = IMAX(total, 0);
len = m->nbEBands;
len = m->nbEBands;
+   skip_start = start;
+   /* Reserve a bit to signal the end of manually skipped bands. */
+   skip_rsv = total >= 1<<BITRES ? 1<<BITRES : 0;
+   total -= skip_rsv;
+   /* Reserve bits for the intensity and dual stereo parameters. */
+   intensity_rsv = dual_stereo_rsv = 0;
+   if (C==2)
+   {
+      intensity_rsv = LOG2_FRAC_TABLE[end-start];
+      if (intensity_rsv>total)
+         intensity_rsv = 0;
+      else
+      {
+         total -= intensity_rsv;
+         dual_stereo_rsv = total>=1<<BITRES ? 1<<BITRES : 0;
+         total -= dual_stereo_rsv;
+      }
+   }
ALLOC(bits1, len, int);
ALLOC(bits2, len, int);
ALLOC(bits1, len, int);
ALLOC(bits2, len, int);
-   lo = 0;
-   hi = m->nbAllocVectors - 1;
-   while (hi-lo != 1)
+   ALLOC(thresh, len, int);
+   ALLOC(trim_offset, len, int);
+
+   for (j=start;j<end;j++)
{
{
-      int j;
+      /* Below this threshold, we're sure not to allocate any PVQ bits */
+      thresh[j] = IMAX((C)<<BITRES, (3*(m->eBands[j+1]-m->eBands[j])<<LM<<BITRES)>>4);
+      /* Tilt of the allocation curve */
+      trim_offset[j] = C*(m->eBands[j+1]-m->eBands[j])*(alloc_trim-5-LM)*(m->nbEBands-j-1)
+            <<(LM+BITRES)>>6;
+      /* Giving less resolution to single-coefficient bands because they get
+         more benefit from having one coarse value per coefficient*/
+      if ((m->eBands[j+1]-m->eBands[j])<<LM==1)
+         trim_offset[j] -= C<<BITRES;
+   }
+   lo = 1;
+   hi = m->nbAllocVectors - 2;
+   do
+   {
+      int done = 0;
+      int psum = 0;
int mid = (lo+hi) >> 1;
int mid = (lo+hi) >> 1;
-      for (j=0;j<len;j++)
+      for (j=end;j-->start;)
{
{
-         bits1[j] = (m->allocVectors[mid*len+j] + offsets[j])<<BITRES;
-         if (bits1[j] < 0)
-            bits1[j] = 0;
-         /*printf ("%d ", bits[j]);*/
+         int N = m->eBands[j+1]-m->eBands[j];
+         bits1[j] = C*N*m->allocVectors[mid*len+j]<<LM>>2;
+         if (bits1[j] > 0)
+            bits1[j] = IMAX(0, bits1[j] + trim_offset[j]);
+         bits1[j] += offsets[j];
+         if (bits1[j] >= thresh[j] || done)
+         {
+            done = 1;
+            /* Don't allocate more than we can actually use */
+            psum += IMIN(bits1[j], 64*C<<BITRES<<LM);
+         } else {
+            if (bits1[j] >= C<<BITRES)
+               psum += C<<BITRES;
+         }
}
}
-      /*printf ("\n");*/
-      if (vec_bits2pulses(m, bits1, pulses, len) > total<<BITRES)
-         hi = mid;
+      if (psum > total)
+         hi = mid - 1;
else
else
-         lo = mid;
+         lo = mid + 1;
/*printf ("lo = %d, hi = %d\n", lo, hi);*/
}
/*printf ("lo = %d, hi = %d\n", lo, hi);*/
}
+   while (lo <= hi);
+   hi = lo--;
+   /*printf ("interp between %d and %d\n", lo, hi);*/
+   for (j=start;j<end;j++)
{
{
-      int j;
-      for (j=0;j<len;j++)
-      {
-         bits1[j] = m->allocVectors[lo*len+j] + offsets[j];
-         bits2[j] = m->allocVectors[hi*len+j] + offsets[j];
-         if (bits1[j] < 0)
-            bits1[j] = 0;
-         if (bits2[j] < 0)
-            bits2[j] = 0;
-      }
-      ret = interp_bits2pulses(m, bits1, bits2, total, pulses, len);
-      RESTORE_STACK;
-      return ret;
+      int N = m->eBands[j+1]-m->eBands[j];
+      bits1[j] = C*N*m->allocVectors[lo*len+j]<<LM>>2;
+      bits2[j] = C*N*m->allocVectors[hi*len+j]<<LM>>2;
+      if (bits1[j] > 0)
+         bits1[j] = IMAX(0, bits1[j] + trim_offset[j]);
+      if (bits2[j] > 0)
+         bits2[j] = IMAX(0, bits2[j] + trim_offset[j]);
+      if (lo > 0)
+         bits1[j] += offsets[j];
+      bits2[j] += offsets[j];
+      if (offsets[j]>0)
+         skip_start = j;
+      bits2[j] -= bits1[j];
}
}
+   codedBands = interp_bits2pulses(m, start, end, skip_start, bits1, bits2, thresh,
+         total, skip_rsv, intensity, intensity_rsv, dual_stereo, dual_stereo_rsv,
+         pulses, ebits, fine_priority, C, LM, ec, encode, prev);
+   RESTORE_STACK;
+   return codedBands;
}

}