FFmpeg
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
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) (((a) * (b)) >> 16)
45 
46 static av_always_inline void idct(uint8_t *dst, int 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 */, int line_size,
198  int16_t *block /* align 16 */)
199 {
200  idct(dest, line_size, block, 1);
201  memset(block, 0, sizeof(*block) * 64);
202 }
203 
204 static void vp3_idct_add_c(uint8_t *dest /* align 8 */, int line_size,
205  int16_t *block /* align 16 */)
206 {
207  idct(dest, line_size, block, 2);
208  memset(block, 0, sizeof(*block) * 64);
209 }
210 
211 static void vp3_idct_dc_add_c(uint8_t *dest /* align 8 */, int line_size,
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 += line_size;
226  }
227  block[0] = 0;
228 }
229 
230 static void vp3_v_loop_filter_c(uint8_t *first_pixel, int stride,
231  int *bounding_values)
232 {
233  unsigned char *end;
234  int filter_value;
235  const int 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, int 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 }
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
float v
void(* idct_put)(uint8_t *dest, int line_size, int16_t *block)
Definition: vp3dsp.h:41
void(* h_loop_filter)(uint8_t *src, int stride, int *bounding_values)
Definition: vp3dsp.h:45
#define C
const char * b
Definition: vf_curves.c:109
#define xC4S4
Definition: vp3dsp.c:39
static void vp3_idct_put_c(uint8_t *dest, int line_size, int16_t *block)
Definition: vp3dsp.c:197
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
#define AV_WN32A(p, v)
Definition: intreadwrite.h:538
Macro definitions for various function/variable attributes.
void(* idct_dc_add)(uint8_t *dest, int line_size, int16_t *block)
Definition: vp3dsp.h:43
#define xC1S7
Definition: vp3dsp.c:36
#define M(a, b)
Definition: vp3dsp.c:44
uint8_t
#define av_cold
Definition: attributes.h:74
#define H
Definition: swscale-test.c:344
static av_cold int end(AVCodecContext *avctx)
Definition: avrndec.c:67
static void vp3_v_loop_filter_c(uint8_t *first_pixel, int stride, int *bounding_values)
Definition: vp3dsp.c:230
#define A(x)
Definition: vp56_arith.h:28
#define xC6S2
Definition: vp3dsp.c:41
#define F
Definition: af_aformat.c:50
#define xC2S6
Definition: vp3dsp.c:37
Libavcodec external API header.
#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 av_always_inline void idct(uint8_t *dst, int stride, int16_t *input, int type)
Definition: vp3dsp.c:46
static int filter_value(int in, int rrp[8], int v[9])
void(* idct_add)(uint8_t *dest, int line_size, int16_t *block)
Definition: vp3dsp.h:42
static void vp3_idct_dc_add_c(uint8_t *dest, int line_size, int16_t *block)
Definition: vp3dsp.c:211
Definition: vf_geq.c:45
av_cold void ff_vp3dsp_init_arm(VP3DSPContext *c, int flags)
GLint GLenum type
Definition: opengl_enc.c:105
static void vp3_h_loop_filter_c(uint8_t *first_pixel, int stride, int *bounding_values)
Definition: vp3dsp.c:247
static void vp3_idct_add_c(uint8_t *dest, int line_size, int16_t *block)
Definition: vp3dsp.c:204
void ff_vp3dsp_init_x86(VP3DSPContext *c, int flags)
Definition: vp3dsp_init.c:47
#define AV_RN32(p)
Definition: intreadwrite.h:364
static int flags
Definition: cpu.c:47
#define xC3S5
Definition: vp3dsp.c:38
GLint GLenum GLboolean GLsizei stride
Definition: opengl_enc.c:105
common internal and external API header
void(* v_loop_filter)(uint8_t *src, int stride, int *bounding_values)
Definition: vp3dsp.h:44
D(D(float, sse)
Definition: rematrix_init.c:28
static double c[64]
av_cold void ff_vp3dsp_init_ppc(VP3DSPContext *c, int flags)
uint8_t pi<< 24) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0f/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S16, int16_t,(*(constint16_t *) pi >>8)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0f/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S32, int32_t,(*(constint32_t *) pi >>24)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0f/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_FLT, float, av_clip_uint8(lrintf(*(constfloat *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_FLT, float, av_clip_int16(lrintf(*(constfloat *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_FLT, float, av_clipl_int32(llrintf(*(constfloat *) pi *(1U<< 31)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_DBL, double, av_clip_uint8(lrint(*(constdouble *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_DBL, double, av_clip_int16(lrint(*(constdouble *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_DBL, double, av_clipl_int32(llrint(*(constdouble *) pi *(1U<< 31))))#defineSET_CONV_FUNC_GROUP(ofmt, ifmt) staticvoidset_generic_function(AudioConvert *ac){}voidff_audio_convert_free(AudioConvert **ac){if(!*ac) return;ff_dither_free(&(*ac) ->dc);av_freep(ac);}AudioConvert *ff_audio_convert_alloc(AVAudioResampleContext *avr, enumAVSampleFormatout_fmt, enumAVSampleFormatin_fmt, intchannels, intsample_rate, intapply_map){AudioConvert *ac;intin_planar, out_planar;ac=av_mallocz(sizeof(*ac));if(!ac) returnNULL;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);returnNULL;}returnac;}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;}elseif(in_planar) ac->func_type=CONV_FUNC_TYPE_INTERLEAVE;elseac->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);returnac;}intff_audio_convert(AudioConvert *ac, AudioData *out, AudioData *in){intuse_generic=1;intlen=in->nb_samples;intp;if(ac->dc){av_log(ac->avr, AV_LOG_TRACE,"%dsamples-audio_convert:%sto%s(dithered)\n", len, av_get_sample_fmt_name(ac->in_fmt), av_get_sample_fmt_name(ac->out_fmt));returnff_convert_dither(ac-> dc
#define av_always_inline
Definition: attributes.h:37
#define IdctAdjustBeforeShift
Definition: vp3dsp.c:35
#define stride
av_cold void ff_vp3dsp_init(VP3DSPContext *c, int flags)
Definition: vp3dsp.c:280
#define xC5S3
Definition: vp3dsp.c:40
Definition: vf_geq.c:45
static int16_t block[64]
Definition: dct-test.c:110