FFmpeg  4.0
mlpdsp.c
Go to the documentation of this file.
1 /*
2  * Copyright (c) 2007-2008 Ian Caulfield
3  * 2009 Ramiro Polla
4  *
5  * This file is part of FFmpeg.
6  *
7  * FFmpeg is free software; you can redistribute it and/or
8  * modify it under the terms of the GNU Lesser General Public
9  * License as published by the Free Software Foundation; either
10  * version 2.1 of the License, or (at your option) any later version.
11  *
12  * FFmpeg is distributed in the hope that it will be useful,
13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15  * Lesser General Public License for more details.
16  *
17  * You should have received a copy of the GNU Lesser General Public
18  * License along with FFmpeg; if not, write to the Free Software
19  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20  */
21 
22 #include "config.h"
23 #include "libavutil/attributes.h"
24 #include "mlpdsp.h"
25 #include "mlp.h"
26 
28  int firorder, int iirorder,
29  unsigned int filter_shift, int32_t mask,
30  int blocksize, int32_t *sample_buffer)
31 {
32  int32_t *firbuf = state;
33  int32_t *iirbuf = state + MAX_BLOCKSIZE + MAX_FIR_ORDER;
34  const int32_t *fircoeff = coeff;
35  const int32_t *iircoeff = coeff + MAX_FIR_ORDER;
36  int i;
37 
38  for (i = 0; i < blocksize; i++) {
39  int32_t residual = *sample_buffer;
40  unsigned int order;
41  int64_t accum = 0;
42  int32_t result;
43 
44  for (order = 0; order < firorder; order++)
45  accum += (int64_t) firbuf[order] * fircoeff[order];
46  for (order = 0; order < iirorder; order++)
47  accum += (int64_t) iirbuf[order] * iircoeff[order];
48 
49  accum = accum >> filter_shift;
50  result = (accum + residual) & mask;
51 
52  *--firbuf = result;
53  *--iirbuf = result - accum;
54 
55  *sample_buffer = result;
56  sample_buffer += MAX_CHANNELS;
57  }
58 }
59 
61  const int32_t *coeffs,
62  const uint8_t *bypassed_lsbs,
63  const int8_t *noise_buffer,
64  int index,
65  unsigned int dest_ch,
66  uint16_t blockpos,
67  unsigned int maxchan,
68  int matrix_noise_shift,
69  int access_unit_size_pow2,
70  int32_t mask)
71 {
72  unsigned int src_ch, i;
73  int index2 = 2 * index + 1;
74  for (i = 0; i < blockpos; i++) {
75  int64_t accum = 0;
76 
77  for (src_ch = 0; src_ch <= maxchan; src_ch++)
78  accum += (int64_t) samples[src_ch] * coeffs[src_ch];
79 
80  if (matrix_noise_shift) {
81  index &= access_unit_size_pow2 - 1;
82  accum += noise_buffer[index] << (matrix_noise_shift + 7);
83  index += index2;
84  }
85 
86  samples[dest_ch] = ((accum >> 14) & mask) + *bypassed_lsbs;
87  bypassed_lsbs += MAX_CHANNELS;
88  samples += MAX_CHANNELS;
89  }
90 }
91 
93  int8_t *output_shift,
94  uint8_t max_matrix_channel,
95  int is32))(int32_t, uint16_t, int32_t (*)[], void *, uint8_t*, int8_t *, uint8_t, int)
96 {
97  return ff_mlp_pack_output;
98 }
99 
100 int32_t ff_mlp_pack_output(int32_t lossless_check_data,
101  uint16_t blockpos,
102  int32_t (*sample_buffer)[MAX_CHANNELS],
103  void *data,
104  uint8_t *ch_assign,
105  int8_t *output_shift,
106  uint8_t max_matrix_channel,
107  int is32)
108 {
109  unsigned int i, out_ch = 0;
110  int32_t *data_32 = data;
111  int16_t *data_16 = data;
112 
113  for (i = 0; i < blockpos; i++) {
114  for (out_ch = 0; out_ch <= max_matrix_channel; out_ch++) {
115  int mat_ch = ch_assign[out_ch];
116  int32_t sample = sample_buffer[i][mat_ch] *
117  (1U << output_shift[mat_ch]);
118  lossless_check_data ^= (sample & 0xffffff) << mat_ch;
119  if (is32)
120  *data_32++ = sample * 256U;
121  else
122  *data_16++ = sample >> 8;
123  }
124  }
125  return lossless_check_data;
126 }
127 
129 {
134  if (ARCH_ARM)
136  if (ARCH_X86)
138 }
void(* mlp_rematrix_channel)(int32_t *samples, const int32_t *coeffs, const uint8_t *bypassed_lsbs, const int8_t *noise_buffer, int index, unsigned int dest_ch, uint16_t blockpos, unsigned int maxchan, int matrix_noise_shift, int access_unit_size_pow2, int32_t mask)
Definition: mlpdsp.h:54
static void mlp_filter_channel(int32_t *state, const int32_t *coeff, int firorder, int iirorder, unsigned int filter_shift, int32_t mask, int blocksize, int32_t *sample_buffer)
Definition: mlpdsp.c:27
int32_t(*(* mlp_select_pack_output)(uint8_t *ch_assign, int8_t *output_shift, uint8_t max_matrix_channel, int is32))(int32_t
Definition: mlpdsp.h:65
void ff_mlp_rematrix_channel(int32_t *samples, const int32_t *coeffs, const uint8_t *bypassed_lsbs, const int8_t *noise_buffer, int index, unsigned int dest_ch, uint16_t blockpos, unsigned int maxchan, int matrix_noise_shift, int access_unit_size_pow2, int32_t mask)
Definition: mlpdsp.c:60
static int32_t(*)(int32_t, uint16_t, int32_t(*)[], void *, uint8_t *, int8_t *, uint8_t, int) mlp_select_pack_output(uint8_t *ch_assign, int8_t *output_shift, uint8_t max_matrix_channel, int is32)
Definition: mlpdsp.c:92
#define sample
Macro definitions for various function/variable attributes.
uint8_t
#define av_cold
Definition: attributes.h:82
#define MAX_FIR_ORDER
The maximum number of taps in IIR and FIR filters.
Definition: mlp.h:64
const char data[16]
Definition: mxf.c:90
int32_t ff_mlp_pack_output(int32_t lossless_check_data, uint16_t blockpos, int32_t(*sample_buffer)[MAX_CHANNELS], void *data, uint8_t *ch_assign, int8_t *output_shift, uint8_t max_matrix_channel, int is32)
Definition: mlpdsp.c:100
#define U(x)
Definition: vp56_arith.h:37
static const uint16_t mask[17]
Definition: lzw.c:38
#define ARCH_X86
Definition: config.h:38
int32_t(* mlp_pack_output)(int32_t lossless_check_data, uint16_t blockpos, int32_t(*sample_buffer)[MAX_CHANNELS], void *data, uint8_t *ch_assign, int8_t *output_shift, uint8_t max_matrix_channel, int is32)
Definition: mlpdsp.h:69
#define MAX_BLOCKSIZE
Definition: diracdec.c:54
uint64_t residual
Definition: dirac_vlc.h:29
void(* mlp_filter_channel)(int32_t *state, const int32_t *coeff, int firorder, int iirorder, unsigned int filter_shift, int32_t mask, int blocksize, int32_t *sample_buffer)
Definition: mlpdsp.h:50
#define ARCH_ARM
Definition: config.h:19
int32_t
static struct @271 state
av_cold void ff_mlpdsp_init_arm(MLPDSPContext *c)
av_cold void ff_mlpdsp_init(MLPDSPContext *c)
Definition: mlpdsp.c:128
int index
Definition: gxfenc.c:89
void ff_mlpdsp_init_x86(MLPDSPContext *c)
Definition: mlpdsp_init.c:193
#define MAX_CHANNELS
Definition: aac.h:47
int
static double c[64]
static const int16_t coeffs[]
static const double coeff[2][5]
Definition: vf_owdenoise.c:72