FFmpeg  4.0
vp3dsp.c
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2004 The FFmpeg project
3  *
4  * This file is part of FFmpeg.
5  *
6  * FFmpeg is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU Lesser General Public
8  * License as published by the Free Software Foundation; either
9  * version 2.1 of the License, or (at your option) any later version.
10  *
11  * FFmpeg 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  * Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with FFmpeg; if not, write to the Free Software
18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19  */
20 
21 /**
22  * @file
23  * Standard C DSP-oriented functions cribbed from the original VP3
24  * source code.
25  */
26 
27 #include "libavutil/attributes.h"
28 #include "libavutil/common.h"
29 #include "libavutil/intreadwrite.h"
30 
31 #include "avcodec.h"
32 #include "rnd_avg.h"
33 #include "vp3dsp.h"
34 
35 #define IdctAdjustBeforeShift 8
36 #define xC1S7 64277
37 #define xC2S6 60547
38 #define xC3S5 54491
39 #define xC4S4 46341
40 #define xC5S3 36410
41 #define xC6S2 25080
42 #define xC7S1 12785
43 
44 #define M(a, b) ((int)((SUINT)(a) * (b)) >> 16)
45 
46 static av_always_inline void idct(uint8_t *dst, ptrdiff_t stride,
47  int16_t *input, int type)
48 {
49  int16_t *ip = input;
50 
51  int A, B, C, D, Ad, Bd, Cd, Dd, E, F, G, H;
52  int Ed, Gd, Add, Bdd, Fd, Hd;
53 
54  int i;
55 
56  /* Inverse DCT on the rows now */
57  for (i = 0; i < 8; i++) {
58  /* Check for non-zero values */
59  if (ip[0 * 8] | ip[1 * 8] | ip[2 * 8] | ip[3 * 8] |
60  ip[4 * 8] | ip[5 * 8] | ip[6 * 8] | ip[7 * 8]) {
61  A = M(xC1S7, ip[1 * 8]) + M(xC7S1, ip[7 * 8]);
62  B = M(xC7S1, ip[1 * 8]) - M(xC1S7, ip[7 * 8]);
63  C = M(xC3S5, ip[3 * 8]) + M(xC5S3, ip[5 * 8]);
64  D = M(xC3S5, ip[5 * 8]) - M(xC5S3, ip[3 * 8]);
65 
66  Ad = M(xC4S4, (A - C));
67  Bd = M(xC4S4, (B - D));
68 
69  Cd = A + C;
70  Dd = B + D;
71 
72  E = M(xC4S4, (ip[0 * 8] + ip[4 * 8]));
73  F = M(xC4S4, (ip[0 * 8] - ip[4 * 8]));
74 
75  G = M(xC2S6, ip[2 * 8]) + M(xC6S2, ip[6 * 8]);
76  H = M(xC6S2, ip[2 * 8]) - M(xC2S6, ip[6 * 8]);
77 
78  Ed = E - G;
79  Gd = E + G;
80 
81  Add = F + Ad;
82  Bdd = Bd - H;
83 
84  Fd = F - Ad;
85  Hd = Bd + H;
86 
87  /* Final sequence of operations over-write original inputs. */
88  ip[0 * 8] = Gd + Cd;
89  ip[7 * 8] = Gd - Cd;
90 
91  ip[1 * 8] = Add + Hd;
92  ip[2 * 8] = Add - Hd;
93 
94  ip[3 * 8] = Ed + Dd;
95  ip[4 * 8] = Ed - Dd;
96 
97  ip[5 * 8] = Fd + Bdd;
98  ip[6 * 8] = Fd - Bdd;
99  }
100 
101  ip += 1; /* next row */
102  }
103 
104  ip = input;
105 
106  for (i = 0; i < 8; i++) {
107  /* Check for non-zero values (bitwise or faster than ||) */
108  if (ip[1] | ip[2] | ip[3] |
109  ip[4] | ip[5] | ip[6] | ip[7]) {
110  A = M(xC1S7, ip[1]) + M(xC7S1, ip[7]);
111  B = M(xC7S1, ip[1]) - M(xC1S7, ip[7]);
112  C = M(xC3S5, ip[3]) + M(xC5S3, ip[5]);
113  D = M(xC3S5, ip[5]) - M(xC5S3, ip[3]);
114 
115  Ad = M(xC4S4, (A - C));
116  Bd = M(xC4S4, (B - D));
117 
118  Cd = A + C;
119  Dd = B + D;
120 
121  E = M(xC4S4, (ip[0] + ip[4])) + 8;
122  F = M(xC4S4, (ip[0] - ip[4])) + 8;
123 
124  if (type == 1) { // HACK
125  E += 16 * 128;
126  F += 16 * 128;
127  }
128 
129  G = M(xC2S6, ip[2]) + M(xC6S2, ip[6]);
130  H = M(xC6S2, ip[2]) - M(xC2S6, ip[6]);
131 
132  Ed = E - G;
133  Gd = E + G;
134 
135  Add = F + Ad;
136  Bdd = Bd - H;
137 
138  Fd = F - Ad;
139  Hd = Bd + H;
140 
141  /* Final sequence of operations over-write original inputs. */
142  if (type == 1) {
143  dst[0 * stride] = av_clip_uint8((Gd + Cd) >> 4);
144  dst[7 * stride] = av_clip_uint8((Gd - Cd) >> 4);
145 
146  dst[1 * stride] = av_clip_uint8((Add + Hd) >> 4);
147  dst[2 * stride] = av_clip_uint8((Add - Hd) >> 4);
148 
149  dst[3 * stride] = av_clip_uint8((Ed + Dd) >> 4);
150  dst[4 * stride] = av_clip_uint8((Ed - Dd) >> 4);
151 
152  dst[5 * stride] = av_clip_uint8((Fd + Bdd) >> 4);
153  dst[6 * stride] = av_clip_uint8((Fd - Bdd) >> 4);
154  } else {
155  dst[0 * stride] = av_clip_uint8(dst[0 * stride] + ((Gd + Cd) >> 4));
156  dst[7 * stride] = av_clip_uint8(dst[7 * stride] + ((Gd - Cd) >> 4));
157 
158  dst[1 * stride] = av_clip_uint8(dst[1 * stride] + ((Add + Hd) >> 4));
159  dst[2 * stride] = av_clip_uint8(dst[2 * stride] + ((Add - Hd) >> 4));
160 
161  dst[3 * stride] = av_clip_uint8(dst[3 * stride] + ((Ed + Dd) >> 4));
162  dst[4 * stride] = av_clip_uint8(dst[4 * stride] + ((Ed - Dd) >> 4));
163 
164  dst[5 * stride] = av_clip_uint8(dst[5 * stride] + ((Fd + Bdd) >> 4));
165  dst[6 * stride] = av_clip_uint8(dst[6 * stride] + ((Fd - Bdd) >> 4));
166  }
167  } else {
168  if (type == 1) {
169  dst[0*stride] =
170  dst[1*stride] =
171  dst[2*stride] =
172  dst[3*stride] =
173  dst[4*stride] =
174  dst[5*stride] =
175  dst[6*stride] =
176  dst[7*stride] = av_clip_uint8(128 + ((xC4S4 * ip[0] + (IdctAdjustBeforeShift << 16)) >> 20));
177  } else {
178  if (ip[0]) {
179  int v = (xC4S4 * ip[0] + (IdctAdjustBeforeShift << 16)) >> 20;
180  dst[0 * stride] = av_clip_uint8(dst[0 * stride] + v);
181  dst[1 * stride] = av_clip_uint8(dst[1 * stride] + v);
182  dst[2 * stride] = av_clip_uint8(dst[2 * stride] + v);
183  dst[3 * stride] = av_clip_uint8(dst[3 * stride] + v);
184  dst[4 * stride] = av_clip_uint8(dst[4 * stride] + v);
185  dst[5 * stride] = av_clip_uint8(dst[5 * stride] + v);
186  dst[6 * stride] = av_clip_uint8(dst[6 * stride] + v);
187  dst[7 * stride] = av_clip_uint8(dst[7 * stride] + v);
188  }
189  }
190  }
191 
192  ip += 8; /* next column */
193  dst++;
194  }
195 }
196 
197 static void vp3_idct_put_c(uint8_t *dest /* align 8 */, ptrdiff_t stride,
198  int16_t *block /* align 16 */)
199 {
200  idct(dest, stride, block, 1);
201  memset(block, 0, sizeof(*block) * 64);
202 }
203 
204 static void vp3_idct_add_c(uint8_t *dest /* align 8 */, ptrdiff_t stride,
205  int16_t *block /* align 16 */)
206 {
207  idct(dest, stride, block, 2);
208  memset(block, 0, sizeof(*block) * 64);
209 }
210 
211 static void vp3_idct_dc_add_c(uint8_t *dest /* align 8 */, ptrdiff_t stride,
212  int16_t *block /* align 16 */)
213 {
214  int i, dc = (block[0] + 15) >> 5;
215 
216  for (i = 0; i < 8; i++) {
217  dest[0] = av_clip_uint8(dest[0] + dc);
218  dest[1] = av_clip_uint8(dest[1] + dc);
219  dest[2] = av_clip_uint8(dest[2] + dc);
220  dest[3] = av_clip_uint8(dest[3] + dc);
221  dest[4] = av_clip_uint8(dest[4] + dc);
222  dest[5] = av_clip_uint8(dest[5] + dc);
223  dest[6] = av_clip_uint8(dest[6] + dc);
224  dest[7] = av_clip_uint8(dest[7] + dc);
225  dest += stride;
226  }
227  block[0] = 0;
228 }
229 
230 static void vp3_v_loop_filter_c(uint8_t *first_pixel, ptrdiff_t stride,
231  int *bounding_values)
232 {
233  unsigned char *end;
234  int filter_value;
235  const ptrdiff_t nstride = -stride;
236 
237  for (end = first_pixel + 8; first_pixel < end; first_pixel++) {
238  filter_value = (first_pixel[2 * nstride] - first_pixel[stride]) +
239  (first_pixel[0] - first_pixel[nstride]) * 3;
240  filter_value = bounding_values[(filter_value + 4) >> 3];
241 
242  first_pixel[nstride] = av_clip_uint8(first_pixel[nstride] + filter_value);
243  first_pixel[0] = av_clip_uint8(first_pixel[0] - filter_value);
244  }
245 }
246 
247 static void vp3_h_loop_filter_c(uint8_t *first_pixel, ptrdiff_t stride,
248  int *bounding_values)
249 {
250  unsigned char *end;
251  int filter_value;
252 
253  for (end = first_pixel + 8 * stride; first_pixel != end; first_pixel += stride) {
254  filter_value = (first_pixel[-2] - first_pixel[1]) +
255  (first_pixel[ 0] - first_pixel[-1]) * 3;
256  filter_value = bounding_values[(filter_value + 4) >> 3];
257 
258  first_pixel[-1] = av_clip_uint8(first_pixel[-1] + filter_value);
259  first_pixel[ 0] = av_clip_uint8(first_pixel[ 0] - filter_value);
260  }
261 }
262 
263 static void put_no_rnd_pixels_l2(uint8_t *dst, const uint8_t *src1,
264  const uint8_t *src2, ptrdiff_t stride, int h)
265 {
266  int i;
267 
268  for (i = 0; i < h; i++) {
269  uint32_t a, b;
270 
271  a = AV_RN32(&src1[i * stride]);
272  b = AV_RN32(&src2[i * stride]);
273  AV_WN32A(&dst[i * stride], no_rnd_avg32(a, b));
274  a = AV_RN32(&src1[i * stride + 4]);
275  b = AV_RN32(&src2[i * stride + 4]);
276  AV_WN32A(&dst[i * stride + 4], no_rnd_avg32(a, b));
277  }
278 }
279 
281 {
283 
289 
290  if (ARCH_ARM)
291  ff_vp3dsp_init_arm(c, flags);
292  if (ARCH_PPC)
293  ff_vp3dsp_init_ppc(c, flags);
294  if (ARCH_X86)
295  ff_vp3dsp_init_x86(c, flags);
296 }
static void vp3_idct_put_c(uint8_t *dest, ptrdiff_t stride, int16_t *block)
Definition: vp3dsp.c:197
#define C
#define ARCH_PPC
Definition: config.h:29
void(* put_no_rnd_pixels_l2)(uint8_t *dst, const uint8_t *a, const uint8_t *b, ptrdiff_t stride, int h)
Copy 8xH pixels from source to destination buffer using a bilinear filter with no rounding (i...
Definition: vp3dsp.h:36
const char * b
Definition: vf_curves.c:113
#define xC4S4
Definition: vp3dsp.c:39
void(* v_loop_filter)(uint8_t *src, ptrdiff_t stride, int *bounding_values)
Definition: vp3dsp.h:44
static void put_no_rnd_pixels_l2(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, ptrdiff_t stride, int h)
Definition: vp3dsp.c:263
int stride
Definition: mace.c:144
#define AV_WN32A(p, v)
Definition: intreadwrite.h:538
Macro definitions for various function/variable attributes.
static int16_t block[64]
Definition: dct.c:115
#define xC1S7
Definition: vp3dsp.c:36
#define M(a, b)
Definition: vp3dsp.c:44
uint8_t
#define av_cold
Definition: attributes.h:82
static av_cold int end(AVCodecContext *avctx)
Definition: avrndec.c:90
static int flags
Definition: log.c:55
#define A(x)
Definition: vp56_arith.h:28
static void vp3_idct_add_c(uint8_t *dest, ptrdiff_t stride, int16_t *block)
Definition: vp3dsp.c:204
#define xC6S2
Definition: vp3dsp.c:41
#define ARCH_X86
Definition: config.h:38
static void vp3_v_loop_filter_c(uint8_t *first_pixel, ptrdiff_t stride, int *bounding_values)
Definition: vp3dsp.c:230
#define B
Definition: huffyuvdsp.h:32
void(* idct_add)(uint8_t *dest, ptrdiff_t stride, int16_t *block)
Definition: vp3dsp.h:42
#define xC2S6
Definition: vp3dsp.c:37
void(* h_loop_filter)(uint8_t *src, ptrdiff_t stride, int *bounding_values)
Definition: vp3dsp.h:45
#define xC7S1
Definition: vp3dsp.c:42
static uint32_t no_rnd_avg32(uint32_t a, uint32_t b)
Definition: rnd_avg.h:36
#define E
Definition: avdct.c:32
static int filter_value(int in, int rrp[8], int v[9])
#define ARCH_ARM
Definition: config.h:19
static void vp3_h_loop_filter_c(uint8_t *first_pixel, ptrdiff_t stride, int *bounding_values)
Definition: vp3dsp.c:247
void(* idct_dc_add)(uint8_t *dest, ptrdiff_t stride, int16_t *block)
Definition: vp3dsp.h:43
#define src1
Definition: h264pred.c:139
Libavcodec external API header.
av_cold void ff_vp3dsp_init_arm(VP3DSPContext *c, int flags)
static void vp3_idct_dc_add_c(uint8_t *dest, ptrdiff_t stride, int16_t *block)
Definition: vp3dsp.c:211
cl_device_type type
void ff_vp3dsp_init_x86(VP3DSPContext *c, int flags)
Definition: vp3dsp_init.c:46
#define AV_RN32(p)
Definition: intreadwrite.h:364
uint8_t pi<< 24) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_U8, uint8_t,(*(const uint8_t *) pi - 0x80) *(1.0f/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_U8, uint8_t,(*(const uint8_t *) pi - 0x80) *(1.0/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S16, int16_t,(*(const int16_t *) pi >> 8)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S16, int16_t, *(const int16_t *) pi *(1.0f/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S16, int16_t, *(const int16_t *) pi *(1.0/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S32, int32_t,(*(const int32_t *) pi >> 24)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S32, int32_t, *(const int32_t *) pi *(1.0f/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S32, int32_t, *(const int32_t *) pi *(1.0/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_FLT, float, av_clip_uint8(lrintf(*(const float *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_FLT, float, av_clip_int16(lrintf(*(const float *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_FLT, float, av_clipl_int32(llrintf(*(const float *) pi *(1U<< 31)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_DBL, double, av_clip_uint8(lrint(*(const double *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_DBL, double, av_clip_int16(lrint(*(const double *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_DBL, double, av_clipl_int32(llrint(*(const double *) pi *(1U<< 31)))) #define SET_CONV_FUNC_GROUP(ofmt, ifmt) static void set_generic_function(AudioConvert *ac) { } void ff_audio_convert_free(AudioConvert **ac) { if(! *ac) return;ff_dither_free(&(*ac) ->dc);av_freep(ac);} AudioConvert *ff_audio_convert_alloc(AVAudioResampleContext *avr, enum AVSampleFormat out_fmt, enum AVSampleFormat in_fmt, int channels, int sample_rate, int apply_map) { AudioConvert *ac;int in_planar, out_planar;ac=av_mallocz(sizeof(*ac));if(!ac) return NULL;ac->avr=avr;ac->out_fmt=out_fmt;ac->in_fmt=in_fmt;ac->channels=channels;ac->apply_map=apply_map;if(avr->dither_method !=AV_RESAMPLE_DITHER_NONE &&av_get_packed_sample_fmt(out_fmt)==AV_SAMPLE_FMT_S16 &&av_get_bytes_per_sample(in_fmt) > 2) { ac->dc=ff_dither_alloc(avr, out_fmt, in_fmt, channels, sample_rate, apply_map);if(!ac->dc) { av_free(ac);return NULL;} return ac;} in_planar=ff_sample_fmt_is_planar(in_fmt, channels);out_planar=ff_sample_fmt_is_planar(out_fmt, channels);if(in_planar==out_planar) { ac->func_type=CONV_FUNC_TYPE_FLAT;ac->planes=in_planar ? ac->channels :1;} else if(in_planar) ac->func_type=CONV_FUNC_TYPE_INTERLEAVE;else ac->func_type=CONV_FUNC_TYPE_DEINTERLEAVE;set_generic_function(ac);if(ARCH_AARCH64) ff_audio_convert_init_aarch64(ac);if(ARCH_ARM) ff_audio_convert_init_arm(ac);if(ARCH_X86) ff_audio_convert_init_x86(ac);return ac;} int ff_audio_convert(AudioConvert *ac, AudioData *out, AudioData *in) { int use_generic=1;int len=in->nb_samples;int p;if(ac->dc) { av_log(ac->avr, AV_LOG_TRACE, "%d samples - audio_convert: %s to %s (dithered)\", len, av_get_sample_fmt_name(ac->in_fmt), av_get_sample_fmt_name(ac->out_fmt));return ff_convert_dither(ac-> dc
#define xC3S5
Definition: vp3dsp.c:38
common internal and external API header
#define G
Definition: huffyuvdsp.h:33
D(D(float, sse)
Definition: rematrix_init.c:28
static double c[64]
av_cold void ff_vp3dsp_init_ppc(VP3DSPContext *c, int flags)
static av_always_inline void idct(uint8_t *dst, ptrdiff_t stride, int16_t *input, int type)
Definition: vp3dsp.c:46
#define F(x)
#define H
Definition: pixlet.c:39
#define av_always_inline
Definition: attributes.h:39
#define IdctAdjustBeforeShift
Definition: vp3dsp.c:35
void(* idct_put)(uint8_t *dest, ptrdiff_t stride, int16_t *block)
Definition: vp3dsp.h:41
av_cold void ff_vp3dsp_init(VP3DSPContext *c, int flags)
Definition: vp3dsp.c:280
#define xC5S3
Definition: vp3dsp.c:40