Commit 15d146c9 authored by Justin Ruggles's avatar Justin Ruggles

shorten: merge decoding of FN_DIFF* subblocks into decode_subframe_lpc()

parent 034f42df
...@@ -267,37 +267,54 @@ static int16_t * interleave_buffer(int16_t *samples, int nchan, int blocksize, i ...@@ -267,37 +267,54 @@ static int16_t * interleave_buffer(int16_t *samples, int nchan, int blocksize, i
return samples; return samples;
} }
static int decode_subframe_lpc(ShortenContext *s, int channel, static const int fixed_coeffs[3][3] = {
{ 1, 0, 0 },
{ 2, -1, 0 },
{ 3, -3, 1 }
};
static int decode_subframe_lpc(ShortenContext *s, int command, int channel,
int residual_size, int32_t coffset) int residual_size, int32_t coffset)
{ {
int pred_order, sum, i, j; int pred_order, sum, qshift, init_sum, i, j;
int *coeffs = s->coeffs; const int *coeffs;
/* read/validate prediction order */ if (command == FN_QLPC) {
pred_order = get_ur_golomb_shorten(&s->gb, LPCQSIZE); /* read/validate prediction order */
if (pred_order > s->nwrap) { pred_order = get_ur_golomb_shorten(&s->gb, LPCQSIZE);
av_log(s->avctx, AV_LOG_ERROR, "invalid pred_order %d\n", pred_order); if (pred_order > s->nwrap) {
return AVERROR(EINVAL); av_log(s->avctx, AV_LOG_ERROR, "invalid pred_order %d\n", pred_order);
return AVERROR(EINVAL);
}
/* read LPC coefficients */
for (i=0; i<pred_order; i++)
s->coeffs[i] = get_sr_golomb_shorten(&s->gb, LPCQUANT);
coeffs = s->coeffs;
qshift = LPCQUANT;
} else {
/* fixed LPC coeffs */
pred_order = command;
coeffs = fixed_coeffs[pred_order-1];
qshift = 0;
} }
/* read LPC coefficients */
for (i=0; i<pred_order; i++)
coeffs[i] = get_sr_golomb_shorten(&s->gb, LPCQUANT);
/* subtract offset from previous samples to use in prediction */ /* subtract offset from previous samples to use in prediction */
if (coffset) if (command == FN_QLPC && coffset)
for (i = -pred_order; i < 0; i++) for (i = -pred_order; i < 0; i++)
s->decoded[channel][i] -= coffset; s->decoded[channel][i] -= coffset;
/* decode residual and do LPC prediction */ /* decode residual and do LPC prediction */
init_sum = pred_order ? (command == FN_QLPC ? s->lpcqoffset : 0) : coffset;
for (i=0; i < s->blocksize; i++) { for (i=0; i < s->blocksize; i++) {
sum = s->lpcqoffset; sum = init_sum;
for (j=0; j<pred_order; j++) for (j=0; j<pred_order; j++)
sum += coeffs[j] * s->decoded[channel][i-j-1]; sum += coeffs[j] * s->decoded[channel][i-j-1];
s->decoded[channel][i] = get_sr_golomb_shorten(&s->gb, residual_size) + (sum >> LPCQUANT); s->decoded[channel][i] = get_sr_golomb_shorten(&s->gb, residual_size) + (sum >> qshift);
} }
/* add offset to current samples */ /* add offset to current samples */
if (coffset != 0) if (command == FN_QLPC && coffset)
for (i = 0; i < s->blocksize; i++) for (i = 0; i < s->blocksize; i++)
s->decoded[channel][i] += coffset; s->decoded[channel][i] += coffset;
...@@ -497,34 +514,12 @@ static int shorten_decode_frame(AVCodecContext *avctx, ...@@ -497,34 +514,12 @@ static int shorten_decode_frame(AVCodecContext *avctx,
if (s->version >= 2) if (s->version >= 2)
coffset >>= FFMIN(1, s->bitshift); coffset >>= FFMIN(1, s->bitshift);
} }
switch (cmd) { if (cmd == FN_ZERO) {
case FN_ZERO: for (i=0; i<s->blocksize; i++)
for (i=0; i<s->blocksize; i++) s->decoded[channel][i] = 0;
s->decoded[channel][i] = 0; } else {
break; if ((ret = decode_subframe_lpc(s, cmd, channel, residual_size, coffset)) < 0)
case FN_DIFF0: return ret;
for (i=0; i<s->blocksize; i++)
s->decoded[channel][i] = get_sr_golomb_shorten(&s->gb, residual_size) + coffset;
break;
case FN_DIFF1:
for (i=0; i<s->blocksize; i++)
s->decoded[channel][i] = get_sr_golomb_shorten(&s->gb, residual_size) + s->decoded[channel][i - 1];
break;
case FN_DIFF2:
for (i=0; i<s->blocksize; i++)
s->decoded[channel][i] = get_sr_golomb_shorten(&s->gb, residual_size) + 2*s->decoded[channel][i-1]
- s->decoded[channel][i-2];
break;
case FN_DIFF3:
for (i=0; i<s->blocksize; i++)
s->decoded[channel][i] = get_sr_golomb_shorten(&s->gb, residual_size) + 3*s->decoded[channel][i-1]
- 3*s->decoded[channel][i-2]
+ s->decoded[channel][i-3];
break;
case FN_QLPC:
if ((ret = decode_subframe_lpc(s, channel, residual_size, coffset)) < 0)
return ret;
break;
} }
if (s->nmean > 0) { if (s->nmean > 0) {
int32_t sum = (s->version < 2) ? 0 : s->blocksize / 2; int32_t sum = (s->version < 2) ? 0 : s->blocksize / 2;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment