49 #define OFFSET(x) offsetof(FSPPContext, x)
50 #define FLAGS AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_VIDEO_PARAM
62 { 0, 48, 12, 60, 3, 51, 15, 63, },
63 { 32, 16, 44, 28, 35, 19, 47, 31, },
64 { 8, 56, 4, 52, 11, 59, 7, 55, },
65 { 40, 24, 36, 20, 43, 27, 39, 23, },
66 { 2, 50, 14, 62, 1, 49, 13, 61, },
67 { 34, 18, 46, 30, 33, 17, 45, 29, },
68 { 10, 58, 6, 54, 9, 57, 5, 53, },
69 { 42, 26, 38, 22, 41, 25, 37, 21, },
76 71, 296, 295, 237, 71, 40, 38, 19,
77 245, 193, 185, 121, 102, 73, 53, 27,
78 158, 129, 141, 107, 97, 73, 50, 26,
79 102, 116, 109, 98, 82, 66, 45, 23,
80 71, 94, 95, 81, 70, 56, 38, 20,
81 56, 77, 74, 66, 56, 44, 30, 15,
82 38, 53, 50, 45, 38, 30, 21, 11,
83 20, 27, 26, 23, 20, 15, 11, 5
88 ptrdiff_t dst_stride, ptrdiff_t src_stride,
89 ptrdiff_t
width, ptrdiff_t
height, ptrdiff_t log2_scale)
93 temp = (src[x + pos] + (d[pos] >> log2_scale)) >> (6 - log2_scale); \
94 src[x + pos] = src[x + pos - 8 * src_stride] = 0; \
95 if (temp & 0x100) temp = ~(temp >> 31); \
98 for (y = 0; y <
height; y++) {
100 for (x = 0; x <
width; x += 8) {
118 ptrdiff_t dst_stride, ptrdiff_t src_stride,
119 ptrdiff_t
width, ptrdiff_t
height, ptrdiff_t log2_scale)
122 #define STORE2(pos) \
123 temp = (src[x + pos] + src[x + pos + 16 * src_stride] + (d[pos] >> log2_scale)) >> (6 - log2_scale); \
124 src[x + pos + 16 * src_stride] = 0; \
125 if (temp & 0x100) temp = ~(temp >> 31); \
128 for (y = 0; y <
height; y++) {
130 for (x = 0; x <
width; x += 8) {
146 static void mul_thrmat_c(int16_t *thr_adr_noq, int16_t *thr_adr,
int q)
149 for (
a = 0;
a < 64;
a++)
150 thr_adr[
a] = q * thr_adr_noq[
a];
154 int dst_stride,
int src_stride,
156 uint8_t *qp_store,
int qp_stride,
int is_luma)
158 int x, x0, y, es, qy, t;
162 const int qpsh = 4 - p->
hsub * !is_luma;
163 const int qpsv = 4 - p->
vsub * !is_luma;
166 int16_t *
block = (int16_t *)block_align;
167 int16_t *block3 = (int16_t *)(block_align + 4 * 8 *
BLOCKSZ);
169 memset(block3, 0, 4 * 8 *
BLOCKSZ);
171 if (!
src || !dst)
return;
173 for (y = 0; y <
height; y++) {
176 for (x = 0; x < 8; x++) {
182 for (y = 0; y < 8; y++) {
188 for (y = 8; y < 24; y++)
192 const int y1 = y - 8 +
step;
198 qy = (qy >> qpsv) * qp_stride;
207 for (x = 0; x < 8 * (
BLOCKSZ - 1); x += 8) {
212 t = qp_store[qy + (t >> qpsh)];
220 memmove(block3, block3 + (
BLOCKSZ - 1) * 64, 6 * 8 *
sizeof(int16_t));
231 if (!(y1 & 7) && y1) {
253 int_simd16_t tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
266 for (; cnt > 0; cnt -= 2) {
267 threshold = (int16_t *)thr_adr;
268 for (ctr =
DCTSIZE; ctr > 0; ctr--) {
303 tmp10 = (tmp0 + tmp2) >> 2;
304 tmp11 = (tmp0 - tmp2) >> 2;
306 tmp13 = (tmp1 + tmp3) >>2;
309 tmp0 = tmp10 + tmp13;
310 tmp3 = tmp10 - tmp13;
311 tmp1 = tmp11 + tmp12;
312 tmp2 = tmp11 - tmp12;
342 z10 = (tmp6 - tmp5) << 1;
344 z12 = (tmp4 - tmp7) << 1;
346 tmp7 = (z11 + z13) >> 2;
356 wsptr[
DCTSIZE * 0] += (tmp0 + tmp7);
357 wsptr[
DCTSIZE * 1] += (tmp1 + tmp6);
358 wsptr[
DCTSIZE * 2] += (tmp2 + tmp5);
359 wsptr[
DCTSIZE * 3] += (tmp3 - tmp4);
360 wsptr[
DCTSIZE * 4] += (tmp3 + tmp4);
361 wsptr[
DCTSIZE * 5] += (tmp2 - tmp5);
362 wsptr[
DCTSIZE * 6] = (tmp1 - tmp6);
363 wsptr[
DCTSIZE * 7] = (tmp0 - tmp7);
374 static void row_idct_c(int16_t *workspace, int16_t *output_adr, ptrdiff_t output_stride,
int cnt)
376 int_simd16_t tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
385 for (; cnt > 0; cnt--) {
388 tmp10 = wsptr[2] + wsptr[3];
389 tmp11 = wsptr[2] - wsptr[3];
391 tmp13 = wsptr[0] + wsptr[1];
394 tmp0 = tmp10 + tmp13;
395 tmp3 = tmp10 - tmp13;
396 tmp1 = tmp11 + tmp12;
397 tmp2 = tmp11 - tmp12;
405 z13 = wsptr[4] + wsptr[5];
406 z10 = wsptr[4] - wsptr[5];
407 z11 = wsptr[6] + wsptr[7];
408 z12 = wsptr[6] - wsptr[7];
417 tmp6 = (tmp12 << 3) - tmp7;
418 tmp5 = (tmp11 << 3) - tmp6;
419 tmp4 = (tmp10 << 3) + tmp5;
422 outptr[0 * output_stride] +=
DESCALE(tmp0 + tmp7, 3);
423 outptr[1 * output_stride] +=
DESCALE(tmp1 + tmp6, 3);
424 outptr[2 * output_stride] +=
DESCALE(tmp2 + tmp5, 3);
425 outptr[3 * output_stride] +=
DESCALE(tmp3 - tmp4, 3);
426 outptr[4 * output_stride] +=
DESCALE(tmp3 + tmp4, 3);
427 outptr[5 * output_stride] +=
DESCALE(tmp2 - tmp5, 3);
428 outptr[6 * output_stride] +=
DESCALE(tmp1 - tmp6, 3);
429 outptr[7 * output_stride] +=
DESCALE(tmp0 - tmp7, 3);
436 static void row_fdct_c(int16_t *
data,
const uint8_t *pixels, ptrdiff_t line_size,
int cnt)
438 int_simd16_t tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
447 for (; cnt > 0; cnt--) {
448 tmp0 = pixels[line_size * 0] + pixels[line_size * 7];
449 tmp7 = pixels[line_size * 0] - pixels[line_size * 7];
450 tmp1 = pixels[line_size * 1] + pixels[line_size * 6];
451 tmp6 = pixels[line_size * 1] - pixels[line_size * 6];
452 tmp2 = pixels[line_size * 2] + pixels[line_size * 5];
453 tmp5 = pixels[line_size * 2] - pixels[line_size * 5];
454 tmp3 = pixels[line_size * 3] + pixels[line_size * 4];
455 tmp4 = pixels[line_size * 3] - pixels[line_size * 4];
466 dataptr[2] = tmp10 + tmp11;
467 dataptr[3] = tmp10 - tmp11;
470 dataptr[0] = tmp13 + z1;
471 dataptr[1] = tmp13 - z1;
475 tmp10 = (tmp4 + tmp5) << 2;
476 tmp11 = (tmp5 + tmp6) << 2;
477 tmp12 = (tmp6 + tmp7) << 2;
487 dataptr[4] = z13 + z2;
488 dataptr[5] = z13 - z2;
489 dataptr[6] = z11 + z4;
490 dataptr[7] = z11 - z4;
546 int8_t *qp_table =
NULL;
549 int custom_threshold_m[64];
556 for (
i = 0;
i < 8;
i++) {
558 |(((uint64_t)custom_threshold_m[
i * 8 + 6]) << 16)
559 |(((uint64_t)custom_threshold_m[
i * 8 + 0]) << 32)
560 |(((uint64_t)custom_threshold_m[
i * 8 + 4]) << 48);
563 |(((uint64_t)custom_threshold_m[
i * 8 + 3]) << 16)
564 |(((uint64_t)custom_threshold_m[
i * 8 + 1]) << 32)
565 |(((uint64_t)custom_threshold_m[
i * 8 + 7]) << 48);
595 if (qp_table || fspp->
qp) {
619 cw, ch, qp_table, qp_stride, 0);
621 cw, ch, qp_table, qp_stride, 0);
665 .priv_class = &fspp_class,