40 "pshuflw $0, "a", "a" \n\t"\
41 "punpcklwd "a", "a" \n\t"
43 "movhlps "a", "b" \n\t"\
44 "pmaxsw "b", "a" \n\t"\
45 "pshuflw $0x0E, "a", "b" \n\t"\
46 "pmaxsw "b", "a" \n\t"\
47 "pshuflw $0x01, "a", "b" \n\t"\
48 "pmaxsw "b", "a" \n\t"
50 #if COMPILE_TEMPLATE_SSSE3
51 #define SAVE_SIGN(a,b) \
52 "movdqa "b", "a" \n\t"\
54 #define RESTORE_SIGN(a,b) \
55 "psignw "a", "b" \n\t"
57 #define SAVE_SIGN(a,b) \
59 "pcmpgtw "b", "a" \n\t" \
62 #define RESTORE_SIGN(a,b) \
64 "psubw "a", "b" \n\t" // out=((ABS(block[i])*qmat[0] - bias[0]*qmat[0])>>16)*sign(block[i])
68 int16_t *
block,
int n,
79 if (
s->dct_error_sum) {
80 const int intra =
s->c.mb_intra;
81 s->dct_count[intra]++;
82 #if HAVE_SSE2_EXTERNAL
85 s->mpvencdsp.denoise_dct(
block,
s->dct_error_sum[intra],
s->dct_offset[intra]);
93 qmat =
s->q_intra_matrix16[qscale][0];
96 qmat =
s->q_chroma_intra_matrix16[qscale][0];
111 last_non_zero_p1 = 1;
113 last_non_zero_p1 = 0;
114 qmat =
s->q_inter_matrix16[qscale][0];
119 "movd %%"FF_REG_a
", %%xmm3 \n\t"
121 "pxor %%xmm2, %%xmm2 \n\t"
122 "pxor %%xmm4, %%xmm4 \n\t"
123 "movdqa (%2), %%xmm5 \n\t"
124 "pxor %%xmm6, %%xmm6 \n\t"
125 "psubw 128(%2), %%xmm6 \n\t"
126 "mov $-128, %%"FF_REG_a
" \n\t"
129 "movdqa (%1, %%"FF_REG_a
"), %%xmm0 \n\t"
131 "psubusw %%xmm6, %%xmm0 \n\t"
132 "pmulhw %%xmm5, %%xmm0 \n\t"
133 "por %%xmm0, %%xmm4 \n\t"
135 "movdqa %%xmm0, (%4, %0) \n\t"
136 "pcmpeqw %%xmm2, %%xmm0 \n\t"
137 "movdqa (%3, %0), %%xmm1 \n\t"
138 "movdqa %%xmm2, (%1, %%"FF_REG_a
") \n\t"
139 "pandn %%xmm1, %%xmm0 \n\t"
140 "pmaxsw %%xmm0, %%xmm3 \n\t"
141 "add $16, %%"FF_REG_a
" \n\t"
143 PMAX(
"%%xmm3",
"%%xmm0")
144 "movd %%xmm3, %%"FF_REG_a
" \n\t"
145 "movzbl %%al, %%eax \n\t"
146 :
"+a" (last_non_zero_p1)
147 :
"r" (
block+64),
"r" (qmat),
150 "%xmm4",
"%xmm5",
"%xmm6")
154 "movd %%"FF_REG_a
", %%xmm3 \n\t"
156 "pxor %%xmm2, %%xmm2 \n\t"
157 "pxor %%xmm4, %%xmm4 \n\t"
158 "mov $-128, %%"FF_REG_a
" \n\t"
161 "movdqa (%1, %%"FF_REG_a
"), %%xmm0 \n\t"
163 "movdqa 128(%2, %0), %%xmm6 \n\t"
164 "paddusw %%xmm6, %%xmm0 \n\t"
165 "movdqa (%2, %%"FF_REG_a
"), %%xmm5 \n\t"
166 "pmulhw %%xmm5, %%xmm0 \n\t"
167 "por %%xmm0, %%xmm4 \n\t"
169 "movdqa %%xmm0, (%4, %0) \n\t"
170 "pcmpeqw %%xmm2, %%xmm0 \n\t"
171 "movdqa (%3, %0), %%xmm1 \n\t"
172 "movdqa %%xmm2, (%1, %%"FF_REG_a
") \n\t"
173 "pandn %%xmm1, %%xmm0 \n\t"
174 "pmaxsw %%xmm0, %%xmm3 \n\t"
175 "add $16, %%"FF_REG_a
" \n\t"
177 PMAX(
"%%xmm3",
"%%xmm0")
178 "movd %%xmm3, %%"FF_REG_a
" \n\t"
179 "movzbl %%al, %%eax \n\t"
180 :
"+a" (last_non_zero_p1)
181 :
"r" (
block+64),
"r" (qmat+64),
184 "%xmm4",
"%xmm5",
"%xmm6")
188 "movd %1, %%xmm1 \n\t"
190 "psubusw %%xmm1, %%xmm4 \n\t"
191 "packuswb %%xmm4, %%xmm4 \n\t"
192 "packsswb %%xmm4, %%xmm4 \n\t"
193 "movd %%xmm4, %0 \n\t"
195 :
"g" (
s->max_qcoeff)
202 if(last_non_zero_p1 <= 1)
goto end;
203 block[0x08] = temp_block[0x01];
block[0x10] = temp_block[0x08];
204 block[0x20] = temp_block[0x10];
205 if(last_non_zero_p1 <= 4)
goto end;
206 block[0x18] = temp_block[0x09];
block[0x04] = temp_block[0x02];
207 block[0x09] = temp_block[0x03];
208 if(last_non_zero_p1 <= 7)
goto end;
209 block[0x14] = temp_block[0x0A];
block[0x28] = temp_block[0x11];
210 block[0x12] = temp_block[0x18];
block[0x02] = temp_block[0x20];
211 if(last_non_zero_p1 <= 11)
goto end;
212 block[0x1A] = temp_block[0x19];
block[0x24] = temp_block[0x12];
213 block[0x19] = temp_block[0x0B];
block[0x01] = temp_block[0x04];
214 block[0x0C] = temp_block[0x05];
215 if(last_non_zero_p1 <= 16)
goto end;
216 block[0x11] = temp_block[0x0C];
block[0x29] = temp_block[0x13];
217 block[0x16] = temp_block[0x1A];
block[0x0A] = temp_block[0x21];
218 block[0x30] = temp_block[0x28];
block[0x22] = temp_block[0x30];
219 block[0x38] = temp_block[0x29];
block[0x06] = temp_block[0x22];
220 if(last_non_zero_p1 <= 24)
goto end;
221 block[0x1B] = temp_block[0x1B];
block[0x21] = temp_block[0x14];
222 block[0x1C] = temp_block[0x0D];
block[0x05] = temp_block[0x06];
223 block[0x0D] = temp_block[0x07];
block[0x15] = temp_block[0x0E];
224 block[0x2C] = temp_block[0x15];
block[0x13] = temp_block[0x1C];
225 if(last_non_zero_p1 <= 32)
goto end;
226 block[0x0B] = temp_block[0x23];
block[0x34] = temp_block[0x2A];
227 block[0x2A] = temp_block[0x31];
block[0x32] = temp_block[0x38];
228 block[0x3A] = temp_block[0x39];
block[0x26] = temp_block[0x32];
229 block[0x39] = temp_block[0x2B];
block[0x03] = temp_block[0x24];
230 if(last_non_zero_p1 <= 40)
goto end;
231 block[0x1E] = temp_block[0x1D];
block[0x25] = temp_block[0x16];
232 block[0x1D] = temp_block[0x0F];
block[0x2D] = temp_block[0x17];
233 block[0x17] = temp_block[0x1E];
block[0x0E] = temp_block[0x25];
234 block[0x31] = temp_block[0x2C];
block[0x2B] = temp_block[0x33];
235 if(last_non_zero_p1 <= 48)
goto end;
236 block[0x36] = temp_block[0x3A];
block[0x3B] = temp_block[0x3B];
237 block[0x23] = temp_block[0x34];
block[0x3C] = temp_block[0x2D];
238 block[0x07] = temp_block[0x26];
block[0x1F] = temp_block[0x1F];
239 block[0x0F] = temp_block[0x27];
block[0x35] = temp_block[0x2E];
240 if(last_non_zero_p1 <= 56)
goto end;
241 block[0x2E] = temp_block[0x35];
block[0x33] = temp_block[0x3C];
242 block[0x3E] = temp_block[0x3D];
block[0x27] = temp_block[0x36];
243 block[0x3D] = temp_block[0x2F];
block[0x2F] = temp_block[0x37];
244 block[0x37] = temp_block[0x3E];
block[0x3F] = temp_block[0x3F];
246 if(last_non_zero_p1 <= 1)
goto end;
247 block[0x04] = temp_block[0x01];
248 block[0x08] = temp_block[0x08];
block[0x10] = temp_block[0x10];
249 if(last_non_zero_p1 <= 4)
goto end;
250 block[0x0C] = temp_block[0x09];
block[0x01] = temp_block[0x02];
251 block[0x05] = temp_block[0x03];
252 if(last_non_zero_p1 <= 7)
goto end;
253 block[0x09] = temp_block[0x0A];
block[0x14] = temp_block[0x11];
254 block[0x18] = temp_block[0x18];
block[0x20] = temp_block[0x20];
255 if(last_non_zero_p1 <= 11)
goto end;
256 block[0x1C] = temp_block[0x19];
257 block[0x11] = temp_block[0x12];
block[0x0D] = temp_block[0x0B];
258 block[0x02] = temp_block[0x04];
block[0x06] = temp_block[0x05];
259 if(last_non_zero_p1 <= 16)
goto end;
260 block[0x0A] = temp_block[0x0C];
block[0x15] = temp_block[0x13];
261 block[0x19] = temp_block[0x1A];
block[0x24] = temp_block[0x21];
262 block[0x28] = temp_block[0x28];
block[0x30] = temp_block[0x30];
263 block[0x2C] = temp_block[0x29];
block[0x21] = temp_block[0x22];
264 if(last_non_zero_p1 <= 24)
goto end;
265 block[0x1D] = temp_block[0x1B];
block[0x12] = temp_block[0x14];
266 block[0x0E] = temp_block[0x0D];
block[0x03] = temp_block[0x06];
267 block[0x07] = temp_block[0x07];
block[0x0B] = temp_block[0x0E];
268 block[0x16] = temp_block[0x15];
block[0x1A] = temp_block[0x1C];
269 if(last_non_zero_p1 <= 32)
goto end;
270 block[0x25] = temp_block[0x23];
block[0x29] = temp_block[0x2A];
271 block[0x34] = temp_block[0x31];
block[0x38] = temp_block[0x38];
272 block[0x3C] = temp_block[0x39];
block[0x31] = temp_block[0x32];
273 block[0x2D] = temp_block[0x2B];
block[0x22] = temp_block[0x24];
274 if(last_non_zero_p1 <= 40)
goto end;
275 block[0x1E] = temp_block[0x1D];
block[0x13] = temp_block[0x16];
276 block[0x0F] = temp_block[0x0F];
block[0x17] = temp_block[0x17];
277 block[0x1B] = temp_block[0x1E];
block[0x26] = temp_block[0x25];
278 block[0x2A] = temp_block[0x2C];
block[0x35] = temp_block[0x33];
279 if(last_non_zero_p1 <= 48)
goto end;
280 block[0x39] = temp_block[0x3A];
block[0x3D] = temp_block[0x3B];
281 block[0x32] = temp_block[0x34];
block[0x2E] = temp_block[0x2D];
282 block[0x23] = temp_block[0x26];
block[0x1F] = temp_block[0x1F];
283 block[0x27] = temp_block[0x27];
block[0x2B] = temp_block[0x2E];
284 if(last_non_zero_p1 <= 56)
goto end;
285 block[0x36] = temp_block[0x35];
block[0x3A] = temp_block[0x3C];
286 block[0x3E] = temp_block[0x3D];
block[0x33] = temp_block[0x36];
287 block[0x2F] = temp_block[0x2F];
block[0x37] = temp_block[0x37];
288 block[0x3B] = temp_block[0x3E];
block[0x3F] = temp_block[0x3F];
290 if(last_non_zero_p1 <= 1)
goto end;
291 block[0x01] = temp_block[0x01];
292 block[0x08] = temp_block[0x08];
block[0x10] = temp_block[0x10];
293 if(last_non_zero_p1 <= 4)
goto end;
294 block[0x09] = temp_block[0x09];
block[0x02] = temp_block[0x02];
295 block[0x03] = temp_block[0x03];
296 if(last_non_zero_p1 <= 7)
goto end;
297 block[0x0A] = temp_block[0x0A];
block[0x11] = temp_block[0x11];
298 block[0x18] = temp_block[0x18];
block[0x20] = temp_block[0x20];
299 if(last_non_zero_p1 <= 11)
goto end;
300 block[0x19] = temp_block[0x19];
301 block[0x12] = temp_block[0x12];
block[0x0B] = temp_block[0x0B];
302 block[0x04] = temp_block[0x04];
block[0x05] = temp_block[0x05];
303 if(last_non_zero_p1 <= 16)
goto end;
304 block[0x0C] = temp_block[0x0C];
block[0x13] = temp_block[0x13];
305 block[0x1A] = temp_block[0x1A];
block[0x21] = temp_block[0x21];
306 block[0x28] = temp_block[0x28];
block[0x30] = temp_block[0x30];
307 block[0x29] = temp_block[0x29];
block[0x22] = temp_block[0x22];
308 if(last_non_zero_p1 <= 24)
goto end;
309 block[0x1B] = temp_block[0x1B];
block[0x14] = temp_block[0x14];
310 block[0x0D] = temp_block[0x0D];
block[0x06] = temp_block[0x06];
311 block[0x07] = temp_block[0x07];
block[0x0E] = temp_block[0x0E];
312 block[0x15] = temp_block[0x15];
block[0x1C] = temp_block[0x1C];
313 if(last_non_zero_p1 <= 32)
goto end;
314 block[0x23] = temp_block[0x23];
block[0x2A] = temp_block[0x2A];
315 block[0x31] = temp_block[0x31];
block[0x38] = temp_block[0x38];
316 block[0x39] = temp_block[0x39];
block[0x32] = temp_block[0x32];
317 block[0x2B] = temp_block[0x2B];
block[0x24] = temp_block[0x24];
318 if(last_non_zero_p1 <= 40)
goto end;
319 block[0x1D] = temp_block[0x1D];
block[0x16] = temp_block[0x16];
320 block[0x0F] = temp_block[0x0F];
block[0x17] = temp_block[0x17];
321 block[0x1E] = temp_block[0x1E];
block[0x25] = temp_block[0x25];
322 block[0x2C] = temp_block[0x2C];
block[0x33] = temp_block[0x33];
323 if(last_non_zero_p1 <= 48)
goto end;
324 block[0x3A] = temp_block[0x3A];
block[0x3B] = temp_block[0x3B];
325 block[0x34] = temp_block[0x34];
block[0x2D] = temp_block[0x2D];
326 block[0x26] = temp_block[0x26];
block[0x1F] = temp_block[0x1F];
327 block[0x27] = temp_block[0x27];
block[0x2E] = temp_block[0x2E];
328 if(last_non_zero_p1 <= 56)
goto end;
329 block[0x35] = temp_block[0x35];
block[0x3C] = temp_block[0x3C];
330 block[0x3D] = temp_block[0x3D];
block[0x36] = temp_block[0x36];
331 block[0x2F] = temp_block[0x2F];
block[0x37] = temp_block[0x37];
332 block[0x3E] = temp_block[0x3E];
block[0x3F] = temp_block[0x3F];
334 if(last_non_zero_p1 <= 1)
goto end;
335 block[0x08] = temp_block[0x01];
336 block[0x01] = temp_block[0x08];
block[0x02] = temp_block[0x10];
337 if(last_non_zero_p1 <= 4)
goto end;
338 block[0x09] = temp_block[0x09];
block[0x10] = temp_block[0x02];
339 block[0x18] = temp_block[0x03];
340 if(last_non_zero_p1 <= 7)
goto end;
341 block[0x11] = temp_block[0x0A];
block[0x0A] = temp_block[0x11];
342 block[0x03] = temp_block[0x18];
block[0x04] = temp_block[0x20];
343 if(last_non_zero_p1 <= 11)
goto end;
344 block[0x0B] = temp_block[0x19];
345 block[0x12] = temp_block[0x12];
block[0x19] = temp_block[0x0B];
346 block[0x20] = temp_block[0x04];
block[0x28] = temp_block[0x05];
347 if(last_non_zero_p1 <= 16)
goto end;
348 block[0x21] = temp_block[0x0C];
block[0x1A] = temp_block[0x13];
349 block[0x13] = temp_block[0x1A];
block[0x0C] = temp_block[0x21];
350 block[0x05] = temp_block[0x28];
block[0x06] = temp_block[0x30];
351 block[0x0D] = temp_block[0x29];
block[0x14] = temp_block[0x22];
352 if(last_non_zero_p1 <= 24)
goto end;
353 block[0x1B] = temp_block[0x1B];
block[0x22] = temp_block[0x14];
354 block[0x29] = temp_block[0x0D];
block[0x30] = temp_block[0x06];
355 block[0x38] = temp_block[0x07];
block[0x31] = temp_block[0x0E];
356 block[0x2A] = temp_block[0x15];
block[0x23] = temp_block[0x1C];
357 if(last_non_zero_p1 <= 32)
goto end;
358 block[0x1C] = temp_block[0x23];
block[0x15] = temp_block[0x2A];
359 block[0x0E] = temp_block[0x31];
block[0x07] = temp_block[0x38];
360 block[0x0F] = temp_block[0x39];
block[0x16] = temp_block[0x32];
361 block[0x1D] = temp_block[0x2B];
block[0x24] = temp_block[0x24];
362 if(last_non_zero_p1 <= 40)
goto end;
363 block[0x2B] = temp_block[0x1D];
block[0x32] = temp_block[0x16];
364 block[0x39] = temp_block[0x0F];
block[0x3A] = temp_block[0x17];
365 block[0x33] = temp_block[0x1E];
block[0x2C] = temp_block[0x25];
366 block[0x25] = temp_block[0x2C];
block[0x1E] = temp_block[0x33];
367 if(last_non_zero_p1 <= 48)
goto end;
368 block[0x17] = temp_block[0x3A];
block[0x1F] = temp_block[0x3B];
369 block[0x26] = temp_block[0x34];
block[0x2D] = temp_block[0x2D];
370 block[0x34] = temp_block[0x26];
block[0x3B] = temp_block[0x1F];
371 block[0x3C] = temp_block[0x27];
block[0x35] = temp_block[0x2E];
372 if(last_non_zero_p1 <= 56)
goto end;
373 block[0x2E] = temp_block[0x35];
block[0x27] = temp_block[0x3C];
374 block[0x2F] = temp_block[0x3D];
block[0x36] = temp_block[0x36];
375 block[0x3D] = temp_block[0x2F];
block[0x3E] = temp_block[0x37];
376 block[0x37] = temp_block[0x3E];
block[0x3F] = temp_block[0x3F];
379 (
int)
s->c.idsp.perm_type);
386 return last_non_zero_p1 - 1;