← Back to branch summary

~oah-dev/oah/gst-plugins-bad

« back to all changes in this revision

Viewing changes to gst/mve/mveaudioenc.c

  • Committer: Haakon Sporsheim
  • Date: 2009-03-12 13:52:03 UTC
  • Revision ID: haakon.sporsheim@tandberg.com-20090312135203-i5k294hgkushb0mt
Initial import of git repository: git://anongit.freedesktop.org/gstreamer/gst-plugins-bad (tag: RELEASE-0_10_10)

Show diffs side-by-side

added added

removed removed

Lines of Context:
gst/mve/mveaudioenc.c
 
1
/* 
 
2
 * Interplay MVE audio compressor
 
3
 * Copyright (C) 2003, 2004 Alexander Belyakov <abel@krasu.ru>
 
4
 * Copyright (C) 2006 Jens Granseuer <jensgr@gmx.net>
 
5
 *
 
6
 * This library is free software; you can redistribute it and/or
 
7
 * modify it under the terms of the GNU Library General Public
 
8
 * License as published by the Free Software Foundation; either
 
9
 * version 2 of the License, or (at your option) any later version.
 
10
 *
 
11
 * This library is distributed in the hope that it will be useful,
 
12
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 
13
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 
14
 * Library General Public License for more details.
 
15
 *
 
16
 * You should have received a copy of the GNU Library General Public
 
17
 * License along with this library; if not, write to the
 
18
 * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
 
19
 * Boston, MA 02111-1307, USA.
 
20
 */
 
21
 
 
22
#include <math.h>
 
23
#include <stdlib.h>
 
24
#include <gst/gst.h>
 
25
 
 
26
static const gint32 dec_table[256] = {
 
27
  0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,
 
28
  16, 17, 18, 19,
 
29
  20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
 
30
  32, 33, 34, 35, 36, 37,
 
31
  38, 39, 40, 41, 42, 43, 47, 51, 56, 61,
 
32
  66, 72, 79, 86, 94, 102, 112,
 
33
  122, 133, 145, 158, 173, 189, 206, 225, 245,
 
34
  267, 292, 318, 348, 379,
 
35
  414, 452, 493, 538, 587, 640, 699, 763, 832, 908, 991,
 
36
  1081, 1180, 1288,
 
37
  1405, 1534, 1673, 1826, 1993, 2175, 2373, 2590, 2826, 3084, 3365, 3672,
 
38
  4008,
 
39
  4373, 4772, 5208, 5683, 6202, 6767, 7385, 8059, 8794, 9597, 10472,
 
40
  11428, 12471, 13609, 14851, 16206,
 
41
  17685, 19298, 21060, 22981, 25078,
 
42
  27367, 29864, 32589, 35563, 38808, 42350, 46214, 50431, 55033, 60055,
 
43
  65535,
 
44
  1, -65535, -60055, -55033, -50431, -46214, -42350, -38808, -35563,
 
45
  -32589, -29864, -27367, -25078, -22981, -21060, -19298,
 
46
  -17685, -16206,
 
47
  -14851, -13609, -12471, -11428, -10472, -9597, -8794, -8059, -7385, -6767,
 
48
  -6202, -5683, -5208, -4772,
 
49
  -4373, -4008, -3672, -3365, -3084, -2826,
 
50
  -2590, -2373, -2175, -1993, -1826, -1673, -1534, -1405, -1288, -1180,
 
51
 
 
52
  -1081, -991, -908, -832, -763, -699, -640, -587, -538, -493, -452, -414,
 
53
  -379, -348, -318, -292,
 
54
  -267, -245, -225, -206, -189, -173, -158, -145,
 
55
  -133, -122, -112, -102, -94, -86, -79, -72,
 
56
  -66, -61, -56, -51, -47, -43,
 
57
  -42, -41, -40, -39, -38, -37, -36, -35, -34, -33,
 
58
  -32, -31, -30, -29,
 
59
  -28, -27, -26, -25, -24, -23, -22, -21, -20, -19, -18, -17,
 
60
  -16, -15,
 
61
  -14, -13, -12, -11, -10, -9, -8, -7, -6, -5, -4, -3, -2, -1
 
62
};
 
63
 
 
64
 
 
65
 
 
66
/* This value could be non-optimal. Without knowledge of the value
 
67
   distribution in the real signal, the actual optimum cannot be evaluated.
 
68
   Should be somewhere between 11.458 and 11.542. */
 
69
static const gdouble DPCM_SCALE = 11.5131;
 
70
 
 
71
static gint8
 
72
mve_enc_delta (guint n)
 
73
{
 
74
  if (n < 44)
 
75
    return n;
 
76
  return floor (DPCM_SCALE * log (n));
 
77
}
 
78
 
 
79
gint
 
80
mve_compress_audio (guint8 * dest, const guint8 * src, guint16 len,
 
81
    guint8 channels)
 
82
{
 
83
  gint16 prev[2], s;
 
84
  gint delta, real_res;
 
85
  gint cur_chan;
 
86
  guint8 v;
 
87
 
 
88
  for (cur_chan = 0; cur_chan < channels; ++cur_chan) {
 
89
    prev[cur_chan] = GST_READ_UINT16_LE (src);
 
90
    GST_WRITE_UINT16_LE (dest, prev[cur_chan]);
 
91
    src += 2;
 
92
    dest += 2;
 
93
    len -= 2;
 
94
  }
 
95
 
 
96
  cur_chan = 0;
 
97
  while (len > 0) {
 
98
    s = GST_READ_UINT16_LE (src);
 
99
    src += 2;
 
100
 
 
101
    delta = s - prev[cur_chan];
 
102
 
 
103
    if (delta >= 0)
 
104
 
 
105
      v = mve_enc_delta (delta);
 
106
 
 
107
    else
 
108
 
 
109
      v = 256 - mve_enc_delta (-delta);
 
110
 
 
111
 
 
112
    real_res = dec_table[v] + prev[cur_chan];
 
113
 
 
114
    if (real_res < -32768 || real_res > 32767) {
 
115
 
 
116
      /* correct overflow */
 
117
      /* GST_DEBUG ("co:%d + %d = %d -> new v:%d, dec_table:%d will be %d",
 
118
         prev[cur_chan], dec_table[v], real_res,
 
119
         v, dec_table[v], prev[cur_chan]+dec_table[v]); */
 
120
      if (s > 0) {
 
121
 
 
122
        if (real_res > 32767)
 
123
          --v;
 
124
 
 
125
      } else {
 
126
 
 
127
        if (real_res < -32768)
 
128
          ++v;
 
129
 
 
130
      }
 
131
 
 
132
      real_res = dec_table[v] + prev[cur_chan];
 
133
 
 
134
    }
 
135
 
 
136
    if (G_UNLIKELY (abs (real_res - s) > 32767)) {
 
137
      GST_ERROR ("sign loss left unfixed in audio stream, deviation:%d",
 
138
          real_res - s);
 
139
      return -1;
 
140
    }
 
141
 
 
142
 
 
143
    *dest++ = v;
 
144
 
 
145
    --len;
 
146
    /* use previous output instead of input. That way output will not go too far from input. */
 
147
    prev[cur_chan] += dec_table[v];
 
148
    cur_chan = channels - 1 - cur_chan;
 
149
 
 
150
  }
 
151
 
 
152
  return 0;
 
153
}