FFmpeg coverage

Directory: ../../../ffmpeg/
File: src/libavcodec/h264dsp_template.c
Date: 2022-11-26 13:19:19
Exec Total Coverage
Lines: 142 142 100.0%
Branches: 70 70 100.0%
Line Branch Exec Source
1 /*
2 * H.26L/H.264/AVC/JVT/14496-10/... encoder/decoder
3 * Copyright (c) 2003-2011 Michael Niedermayer <michaelni@gmx.at>
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 /**
23 * @file
24 * H.264 / AVC / MPEG-4 part10 DSP functions.
25 * @author Michael Niedermayer <michaelni@gmx.at>
26 */
27
28 #include "bit_depth_template.c"
29
30 #define op_scale1(x) block[x] = av_clip_pixel( (block[x]*weight + offset) >> log2_denom )
31 #define op_scale2(x) dst[x] = av_clip_pixel( (src[x]*weights + dst[x]*weightd + offset) >> (log2_denom+1))
32 #define H264_WEIGHT(W) \
33 static void FUNCC(weight_h264_pixels ## W)(uint8_t *_block, ptrdiff_t stride, int height, \
34 int log2_denom, int weight, int offset) \
35 { \
36 int y; \
37 pixel *block = (pixel*)_block; \
38 stride >>= sizeof(pixel)-1; \
39 offset = (unsigned)offset << (log2_denom + (BIT_DEPTH-8)); \
40 if(log2_denom) offset += 1<<(log2_denom-1); \
41 for (y = 0; y < height; y++, block += stride) { \
42 op_scale1(0); \
43 op_scale1(1); \
44 if(W==2) continue; \
45 op_scale1(2); \
46 op_scale1(3); \
47 if(W==4) continue; \
48 op_scale1(4); \
49 op_scale1(5); \
50 op_scale1(6); \
51 op_scale1(7); \
52 if(W==8) continue; \
53 op_scale1(8); \
54 op_scale1(9); \
55 op_scale1(10); \
56 op_scale1(11); \
57 op_scale1(12); \
58 op_scale1(13); \
59 op_scale1(14); \
60 op_scale1(15); \
61 } \
62 } \
63 static void FUNCC(biweight_h264_pixels ## W)(uint8_t *_dst, uint8_t *_src, ptrdiff_t stride, int height, \
64 int log2_denom, int weightd, int weights, int offset) \
65 { \
66 int y; \
67 pixel *dst = (pixel*)_dst; \
68 pixel *src = (pixel*)_src; \
69 stride >>= sizeof(pixel)-1; \
70 offset = (unsigned)offset << (BIT_DEPTH-8); \
71 offset = (unsigned)((offset + 1) | 1) << log2_denom; \
72 for (y = 0; y < height; y++, dst += stride, src += stride) { \
73 op_scale2(0); \
74 op_scale2(1); \
75 if(W==2) continue; \
76 op_scale2(2); \
77 op_scale2(3); \
78 if(W==4) continue; \
79 op_scale2(4); \
80 op_scale2(5); \
81 op_scale2(6); \
82 op_scale2(7); \
83 if(W==8) continue; \
84 op_scale2(8); \
85 op_scale2(9); \
86 op_scale2(10); \
87 op_scale2(11); \
88 op_scale2(12); \
89 op_scale2(13); \
90 op_scale2(14); \
91 op_scale2(15); \
92 } \
93 }
94
95
4/4
✓ Branch 0 taken 3013961 times.
✓ Branch 1 taken 636371 times.
✓ Branch 2 taken 9534024 times.
✓ Branch 3 taken 628642 times.
 26368712 H264_WEIGHT(16)
96
4/4
✓ Branch 0 taken 3639700 times.
✓ Branch 1 taken 458331 times.
✓ Branch 2 taken 3453424 times.
✓ Branch 3 taken 462250 times.
 7551455 H264_WEIGHT(8)
97
4/4
✓ Branch 0 taken 1253658 times.
✓ Branch 1 taken 149092 times.
✓ Branch 2 taken 3073392 times.
✓ Branch 3 taken 690606 times.
 4476142 H264_WEIGHT(4)
98
4/4
✓ Branch 0 taken 642584 times.
✓ Branch 1 taken 98228 times.
✓ Branch 2 taken 1261776 times.
✓ Branch 3 taken 447120 times.
 2002588 H264_WEIGHT(2)
99
100 #undef op_scale1
101 #undef op_scale2
102 #undef H264_WEIGHT
103
104 73681320 static av_always_inline av_flatten void FUNCC(h264_loop_filter_luma)(uint8_t *p_pix, ptrdiff_t xstride, ptrdiff_t ystride, int inner_iters, int alpha, int beta, int8_t *tc0)
105 {
106 73681320 pixel *pix = (pixel*)p_pix;
107 int i, d;
108 73681320 xstride >>= sizeof(pixel)-1;
109 73681320 ystride >>= sizeof(pixel)-1;
110 73681320 alpha <<= BIT_DEPTH - 8;
111 73681320 beta <<= BIT_DEPTH - 8;
112
2/2
✓ Branch 0 taken 147362640 times.
✓ Branch 1 taken 36840660 times.
 368406600 for( i = 0; i < 4; i++ ) {
113 294725280 const int tc_orig = tc0[i] * (1 << (BIT_DEPTH - 8));
114
2/2
✓ Branch 0 taken 24447326 times.
✓ Branch 1 taken 122915314 times.
 294725280 if( tc_orig < 0 ) {
115 48894652 pix += inner_iters*ystride;
116 48894652 continue;
117 }
118
2/2
✓ Branch 0 taken 489259944 times.
✓ Branch 1 taken 122915314 times.
 1224350516 for( d = 0; d < inner_iters; d++ ) {
119 978519888 const int p0 = pix[-1*xstride];
120 978519888 const int p1 = pix[-2*xstride];
121 978519888 const int p2 = pix[-3*xstride];
122 978519888 const int q0 = pix[0];
123 978519888 const int q1 = pix[1*xstride];
124 978519888 const int q2 = pix[2*xstride];
125
126
2/2
✓ Branch 0 taken 419142803 times.
✓ Branch 1 taken 70117141 times.
 978519888 if( FFABS( p0 - q0 ) < alpha &&
127
2/2
✓ Branch 0 taken 315636926 times.
✓ Branch 1 taken 103505877 times.
 838285606 FFABS( p1 - p0 ) < beta &&
128
2/2
✓ Branch 0 taken 264637583 times.
✓ Branch 1 taken 50999343 times.
 631273852 FFABS( q1 - q0 ) < beta ) {
129
130 529275166 int tc = tc_orig;
131 int i_delta;
132
133
2/2
✓ Branch 0 taken 224593857 times.
✓ Branch 1 taken 40043726 times.
 529275166 if( FFABS( p2 - p0 ) < beta ) {
134
2/2
✓ Branch 0 taken 210096452 times.
✓ Branch 1 taken 14497405 times.
 449187714 if(tc_orig)
135 420192904 pix[-2*xstride] = p1 + av_clip( (( p2 + ( ( p0 + q0 + 1 ) >> 1 ) ) >> 1) - p1, -tc_orig, tc_orig );
136 449187714 tc++;
137 }
138
2/2
✓ Branch 0 taken 223441817 times.
✓ Branch 1 taken 41195766 times.
 529275166 if( FFABS( q2 - q0 ) < beta ) {
139
2/2
✓ Branch 0 taken 208969783 times.
✓ Branch 1 taken 14472034 times.
 446883634 if(tc_orig)
140 417939566 pix[ xstride] = q1 + av_clip( (( q2 + ( ( p0 + q0 + 1 ) >> 1 ) ) >> 1) - q1, -tc_orig, tc_orig );
141 446883634 tc++;
142 }
143
144 529275166 i_delta = av_clip( (((q0 - p0 ) * 4) + (p1 - q1) + 4) >> 3, -tc, tc );
145 529275166 pix[-xstride] = av_clip_pixel( p0 + i_delta ); /* p0' */
146 529275166 pix[0] = av_clip_pixel( q0 - i_delta ); /* q0' */
147 }
148 978519888 pix += ystride;
149 }
150 }
151 }
152 18691820 static void FUNCC(h264_v_loop_filter_luma)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta, int8_t *tc0)
153 {
154 37383640 FUNCC(h264_loop_filter_luma)(pix, stride, sizeof(pixel), 4, alpha, beta, tc0);
155 }
156 17848676 static void FUNCC(h264_h_loop_filter_luma)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta, int8_t *tc0)
157 {
158 35697352 FUNCC(h264_loop_filter_luma)(pix, sizeof(pixel), stride, 4, alpha, beta, tc0);
159 }
160 300164 static void FUNCC(h264_h_loop_filter_luma_mbaff)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta, int8_t *tc0)
161 {
162 600328 FUNCC(h264_loop_filter_luma)(pix, sizeof(pixel), stride, 2, alpha, beta, tc0);
163 }
164
165 11165114 static av_always_inline av_flatten void FUNCC(h264_loop_filter_luma_intra)(uint8_t *p_pix, ptrdiff_t xstride, ptrdiff_t ystride, int inner_iters, int alpha, int beta)
166 {
167 11165114 pixel *pix = (pixel*)p_pix;
168 int d;
169 11165114 xstride >>= sizeof(pixel)-1;
170 11165114 ystride >>= sizeof(pixel)-1;
171 11165114 alpha <<= BIT_DEPTH - 8;
172 11165114 beta <<= BIT_DEPTH - 8;
173
2/2
✓ Branch 0 taken 86322416 times.
✓ Branch 1 taken 5582557 times.
 183809946 for( d = 0; d < 4 * inner_iters; d++ ) {
174 172644832 const int p2 = pix[-3*xstride];
175 172644832 const int p1 = pix[-2*xstride];
176 172644832 const int p0 = pix[-1*xstride];
177
178 172644832 const int q0 = pix[ 0*xstride];
179 172644832 const int q1 = pix[ 1*xstride];
180 172644832 const int q2 = pix[ 2*xstride];
181
182
2/2
✓ Branch 0 taken 75926580 times.
✓ Branch 1 taken 10395836 times.
 172644832 if( FFABS( p0 - q0 ) < alpha &&
183
2/2
✓ Branch 0 taken 60891609 times.
✓ Branch 1 taken 15034971 times.
 151853160 FFABS( p1 - p0 ) < beta &&
184
2/2
✓ Branch 0 taken 52455112 times.
✓ Branch 1 taken 8436497 times.
 121783218 FFABS( q1 - q0 ) < beta ) {
185
186
2/2
✓ Branch 0 taken 43573317 times.
✓ Branch 1 taken 8881795 times.
 104910224 if(FFABS( p0 - q0 ) < (( alpha >> 2 ) + 2 )){
187
2/2
✓ Branch 0 taken 39240448 times.
✓ Branch 1 taken 4332869 times.
 87146634 if( FFABS( p2 - p0 ) < beta)
188 {
189 78480896 const int p3 = pix[-4*xstride];
190 /* p0', p1', p2' */
191 78480896 pix[-1*xstride] = ( p2 + 2*p1 + 2*p0 + 2*q0 + q1 + 4 ) >> 3;
192 78480896 pix[-2*xstride] = ( p2 + p1 + p0 + q0 + 2 ) >> 2;
193 78480896 pix[-3*xstride] = ( 2*p3 + 3*p2 + p1 + p0 + q0 + 4 ) >> 3;
194 } else {
195 /* p0' */
196 8665738 pix[-1*xstride] = ( 2*p1 + p0 + q1 + 2 ) >> 2;
197 }
198
2/2
✓ Branch 0 taken 39012074 times.
✓ Branch 1 taken 4561243 times.
 87146634 if( FFABS( q2 - q0 ) < beta)
199 {
200 78024148 const int q3 = pix[3*xstride];
201 /* q0', q1', q2' */
202 78024148 pix[0*xstride] = ( p1 + 2*p0 + 2*q0 + 2*q1 + q2 + 4 ) >> 3;
203 78024148 pix[1*xstride] = ( p0 + q0 + q1 + q2 + 2 ) >> 2;
204 78024148 pix[2*xstride] = ( 2*q3 + 3*q2 + q1 + q0 + p0 + 4 ) >> 3;
205 } else {
206 /* q0' */
207 9122486 pix[0*xstride] = ( 2*q1 + q0 + p1 + 2 ) >> 2;
208 }
209 }else{
210 /* p0', q0' */
211 17763590 pix[-1*xstride] = ( 2*p1 + p0 + q1 + 2 ) >> 2;
212 17763590 pix[ 0*xstride] = ( 2*q1 + q0 + p1 + 2 ) >> 2;
213 }
214 }
215 172644832 pix += ystride;
216 }
217 }
218 2306607 static void FUNCC(h264_v_loop_filter_luma_intra)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta)
219 {
220 4613214 FUNCC(h264_loop_filter_luma_intra)(pix, stride, sizeof(pixel), 4, alpha, beta);
221 }
222 2901138 static void FUNCC(h264_h_loop_filter_luma_intra)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta)
223 {
224 5802276 FUNCC(h264_loop_filter_luma_intra)(pix, sizeof(pixel), stride, 4, alpha, beta);
225 }
226 374812 static void FUNCC(h264_h_loop_filter_luma_mbaff_intra)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta)
227 {
228 749624 FUNCC(h264_loop_filter_luma_intra)(pix, sizeof(pixel), stride, 2, alpha, beta);
229 }
230
231 87311614 static av_always_inline av_flatten void FUNCC(h264_loop_filter_chroma)(uint8_t *p_pix, ptrdiff_t xstride, ptrdiff_t ystride, int inner_iters, int alpha, int beta, int8_t *tc0)
232 {
233 87311614 pixel *pix = (pixel*)p_pix;
234 int i, d;
235 87311614 alpha <<= BIT_DEPTH - 8;
236 87311614 beta <<= BIT_DEPTH - 8;
237 87311614 xstride >>= sizeof(pixel)-1;
238 87311614 ystride >>= sizeof(pixel)-1;
239
2/2
✓ Branch 0 taken 174623228 times.
✓ Branch 1 taken 43655807 times.
 436558070 for( i = 0; i < 4; i++ ) {
240 349246456 const int tc = ((tc0[i] - 1U) << (BIT_DEPTH - 8)) + 1;
241
2/2
✓ Branch 0 taken 22849220 times.
✓ Branch 1 taken 151774008 times.
 349246456 if( tc <= 0 ) {
242 45698440 pix += inner_iters*ystride;
243 45698440 continue;
244 }
245
2/2
✓ Branch 0 taken 312017810 times.
✓ Branch 1 taken 151774008 times.
 927583636 for( d = 0; d < inner_iters; d++ ) {
246 624035620 const int p0 = pix[-1*xstride];
247 624035620 const int p1 = pix[-2*xstride];
248 624035620 const int q0 = pix[0];
249 624035620 const int q1 = pix[1*xstride];
250
251
2/2
✓ Branch 0 taken 305536068 times.
✓ Branch 1 taken 6481742 times.
 624035620 if( FFABS( p0 - q0 ) < alpha &&
252
2/2
✓ Branch 0 taken 287278987 times.
✓ Branch 1 taken 18257081 times.
 611072136 FFABS( p1 - p0 ) < beta &&
253
2/2
✓ Branch 0 taken 275096273 times.
✓ Branch 1 taken 12182714 times.
 574557974 FFABS( q1 - q0 ) < beta ) {
254
255 550192546 int delta = av_clip( ((q0 - p0) * 4 + (p1 - q1) + 4) >> 3, -tc, tc );
256
257 550192546 pix[-xstride] = av_clip_pixel( p0 + delta ); /* p0' */
258 550192546 pix[0] = av_clip_pixel( q0 - delta ); /* q0' */
259 }
260 624035620 pix += ystride;
261 }
262 }
263 }
264 23724842 static void FUNCC(h264_v_loop_filter_chroma)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta, int8_t *tc0)
265 {
266 47449684 FUNCC(h264_loop_filter_chroma)(pix, stride, sizeof(pixel), 2, alpha, beta, tc0);
267 }
268 17894535 static void FUNCC(h264_h_loop_filter_chroma)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta, int8_t *tc0)
269 {
270 35789070 FUNCC(h264_loop_filter_chroma)(pix, sizeof(pixel), stride, 2, alpha, beta, tc0);
271 }
272 649578 static void FUNCC(h264_h_loop_filter_chroma_mbaff)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta, int8_t *tc0)
273 {
274 1299156 FUNCC(h264_loop_filter_chroma)(pix, sizeof(pixel), stride, 1, alpha, beta, tc0);
275 }
276 1386636 static void FUNCC(h264_h_loop_filter_chroma422)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta, int8_t *tc0)
277 {
278 2773272 FUNCC(h264_loop_filter_chroma)(pix, sizeof(pixel), stride, 4, alpha, beta, tc0);
279 }
280 216 static void FUNCC(h264_h_loop_filter_chroma422_mbaff)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta, int8_t *tc0)
281 {
282 432 FUNCC(h264_loop_filter_chroma)(pix, sizeof(pixel), stride, 2, alpha, beta, tc0);
283 }
284
285 21766146 static av_always_inline av_flatten void FUNCC(h264_loop_filter_chroma_intra)(uint8_t *p_pix, ptrdiff_t xstride, ptrdiff_t ystride, int inner_iters, int alpha, int beta)
286 {
287 21766146 pixel *pix = (pixel*)p_pix;
288 int d;
289 21766146 xstride >>= sizeof(pixel)-1;
290 21766146 ystride >>= sizeof(pixel)-1;
291 21766146 alpha <<= BIT_DEPTH - 8;
292 21766146 beta <<= BIT_DEPTH - 8;
293
2/2
✓ Branch 0 taken 94942480 times.
✓ Branch 1 taken 10883073 times.
 211651106 for( d = 0; d < 4 * inner_iters; d++ ) {
294 189884960 const int p0 = pix[-1*xstride];
295 189884960 const int p1 = pix[-2*xstride];
296 189884960 const int q0 = pix[0];
297 189884960 const int q1 = pix[1*xstride];
298
299
2/2
✓ Branch 0 taken 92087612 times.
✓ Branch 1 taken 2854868 times.
 189884960 if( FFABS( p0 - q0 ) < alpha &&
300
2/2
✓ Branch 0 taken 87086386 times.
✓ Branch 1 taken 5001226 times.
 184175224 FFABS( p1 - p0 ) < beta &&
301
2/2
✓ Branch 0 taken 83570896 times.
✓ Branch 1 taken 3515490 times.
 174172772 FFABS( q1 - q0 ) < beta ) {
302
303 167141792 pix[-xstride] = ( 2*p1 + p0 + q1 + 2 ) >> 2; /* p0' */
304 167141792 pix[0] = ( 2*q1 + q0 + p1 + 2 ) >> 2; /* q0' */
305 }
306 189884960 pix += ystride;
307 }
308 }
309 4462478 static void FUNCC(h264_v_loop_filter_chroma_intra)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta)
310 {
311 8924956 FUNCC(h264_loop_filter_chroma_intra)(pix, stride, sizeof(pixel), 2, alpha, beta);
312 }
313 4421523 static void FUNCC(h264_h_loop_filter_chroma_intra)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta)
314 {
315 8843046 FUNCC(h264_loop_filter_chroma_intra)(pix, sizeof(pixel), stride, 2, alpha, beta);
316 }
317 529690 static void FUNCC(h264_h_loop_filter_chroma_mbaff_intra)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta)
318 {
319 1059380 FUNCC(h264_loop_filter_chroma_intra)(pix, sizeof(pixel), stride, 1, alpha, beta);
320 }
321 1249582 static void FUNCC(h264_h_loop_filter_chroma422_intra)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta)
322 {
323 2499164 FUNCC(h264_loop_filter_chroma_intra)(pix, sizeof(pixel), stride, 4, alpha, beta);
324 }
325 219800 static void FUNCC(h264_h_loop_filter_chroma422_mbaff_intra)(uint8_t *pix, ptrdiff_t stride, int alpha, int beta)
326 {
327 439600 FUNCC(h264_loop_filter_chroma_intra)(pix, sizeof(pixel), stride, 2, alpha, beta);
328 }
329