Making the MDCT produce interleaved data
[opus.git] / libcelt / rate.h
index 40f8dbf..e1346ee 100644 (file)
@@ -1,21 +1,18 @@
-/* (C) 2007-2008 Jean-Marc Valin, CSIRO
-*/
+/* Copyright (c) 2007-2008 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
    are met:
-   
+
    - Redistributions of source code must retain the above copyright
    notice, this list of conditions and the following disclaimer.
-   
+
    - Redistributions in binary form must reproduce the above copyright
    notice, this list of conditions and the following disclaimer in the
    documentation and/or other materials provided with the distribution.
-   
-   - Neither the name of the Xiph.org Foundation nor the names of its
-   contributors may be used to endorse or promote products derived from
-   this software without specific prior written permission.
-   
+
    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
    LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
 #ifndef RATE_H
 #define RATE_H
 
+#define MAX_PSEUDO 40
+#define LOG_MAX_PSEUDO 6
+
 #define MAX_PULSES 128
-#define LOG_MAX_PULSES 7
 
-#define BITRES 4
-#define BITROUND 8
-#define BITOVERFLOW 30000
+#define MAX_FINE_BITS 8
+
+#define FINE_OFFSET 21
+#define QTHETA_OFFSET 4
+#define QTHETA_OFFSET_TWOPHASE 16
 
 #include "cwrs.h"
+#include "modes.h"
 
-static inline int bits2pulses(const CELTMode *m, const celt_int16_t *cache, int N, int bits)
+void compute_pulse_cache(CELTMode *m, int LM);
+
+static inline int get_pulses(int i)
+{
+   return i<8 ? i : (8 + (i&7)) << ((i>>3)-1);
+}
+
+static inline int bits2pulses(const CELTMode *m, int band, int LM, int bits)
 {
    int i;
    int lo, hi;
+   const unsigned char *cache;
+
+   LM++;
+   cache = m->cache.bits + m->cache.index[LM*m->nbEBands+band];
+
    lo = 0;
-   hi = MAX_PULSES-1;
-   
-#if 0 /* Disabled until we can make that useful */
-   /* Use of more than MAX_PULSES is disabled until we are able to cwrs that decently */
-   if (bits > cache[MAX_PULSES-1] && N<=4)
-   {
-      /*int pulses;
-      pulses = 127;
-      while (16 + log2_frac(2*(pulses+1)*(pulses+1) + 1, 4) <= bits && pulses < 32767)
-         pulses++;*/
-      lo = 127;
-      switch (N)
-      {
-         case 3:
-            hi = 1024;
-            for (i=0;i<10;i++)
-            {
-               int pulses = (lo+hi)>>1;
-               if (log2_frac(((UMUL16_16(pulses,pulses)>>1)+1)>>1, 4) > bits)
-                  hi = pulses;
-               else
-                  lo = pulses;
-            }
-            break;
-         case 4:
-            hi = 1024;
-            for (i=0;i<10;i++)
-            {
-               int pulses = (lo+hi)>>1;
-               if (log2_frac((UMUL32(UMUL16_16(pulses,pulses)+2,pulses))/3<<3, 4) > bits)
-                  hi = pulses;
-               else
-                  lo = pulses;
-            }
-            break;
-      }
-      return lo;
-   }
-#endif
-   /* 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++)
+   hi = cache[0];
+   bits--;
+   for (i=0;i<LOG_MAX_PSEUDO;i++)
    {
-      int mid = (lo+hi)>>1;
+      int mid = (lo+hi+1)>>1;
       /* OPT: Make sure this is implemented with a conditional move */
       if (cache[mid] >= bits)
          hi = mid;
       else
          lo = mid;
    }
-   if (bits-cache[lo] <= cache[hi]-bits)
+   if (bits- (lo == 0 ? -1 : cache[lo]) <= cache[hi]-bits)
       return lo;
    else
       return hi;
 }
 
-
-static inline int pulses2bits(const celt_int16_t *cache, int N, int pulses)
+static inline int pulses2bits(const CELTMode *m, int band, int LM, int pulses)
 {
-#if 0 /* Use of more than MAX_PULSES is disabled until we are able to cwrs that decently */
-   if (pulses > 127)
-   {
-      int bits;
-      switch (N)
-      {
-         case 3:
-            bits = log2_frac(((UMUL16_16(pulses,pulses)>>1)+1)>>1, 4);
-            break;
-         case 4:
-            bits = log2_frac((UMUL32(UMUL16_16(pulses,pulses)+2,pulses))/3<<3, 4);
-            break;
-      }
-      /*printf ("%d <- %d\n", bits, pulses);*/
-      return bits;
-   }
-#endif
-   return cache[pulses];
+   const unsigned char *cache;
+
+   LM++;
+   cache = m->cache.bits + m->cache.index[LM*m->nbEBands+band];
+   return pulses == 0 ? 0 : cache[pulses]+1;
 }
 
 /** Computes a cache of the pulses->bits mapping in each band */
-celt_int16_t **compute_alloc_cache(CELTMode *m, int C);
+opus_int16 **compute_alloc_cache(CELTMode *m, int M);
 
 /** Compute the pulse allocation, i.e. how many pulses will go in each
-  * band. 
+  * band.
  @param m mode
  @param offsets Requested increase or decrease in the number of bits for
                 each band
@@ -138,7 +98,7 @@ celt_int16_t **compute_alloc_cache(CELTMode *m, int C);
  @param pulses Number of pulses per band (returned)
  @return Total number of bits allocated
 */
-void compute_allocation(const CELTMode *m, int *offsets, const int *stereo_mode, int total, int *pulses, int *ebits);
-
+int compute_allocation(const CELTMode *m, int start, int end, const int *offsets, const int *cap, int alloc_trim, int *intensity, int *dual_stero,
+      opus_int32 total, opus_int32 *balance, int *pulses, int *ebits, int *fine_priority, int _C, int LM, ec_ctx *ec, int encode, int prev);
 
 #endif