1/* 2 * This file is part of FFmpeg. 3 * 4 * FFmpeg is free software; you can redistribute it and/or 5 * modify it under the terms of the GNU Lesser General Public 6 * License as published by the Free Software Foundation; either 7 * version 2.1 of the License, or (at your option) any later version. 8 * 9 * FFmpeg is distributed in the hope that it will be useful, 10 * but WITHOUT ANY WARRANTY; without even the implied warranty of 11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 12 * Lesser General Public License for more details. 13 * 14 * You should have received a copy of the GNU Lesser General Public 15 * License along with FFmpeg; if not, write to the Free Software 16 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA 17 */ 18 19#include "libavutil/attributes.h" 20#include "ttadsp.h" 21#include "config.h" 22 23static void tta_filter_process_c(int32_t *qmi, int32_t *dx, int32_t *dl, 24 int32_t *error, int32_t *in, int32_t shift, 25 int32_t round) { 26 uint32_t *qm = qmi; 27 28 if (*error < 0) { 29 qm[0] -= dx[0]; qm[1] -= dx[1]; qm[2] -= dx[2]; qm[3] -= dx[3]; 30 qm[4] -= dx[4]; qm[5] -= dx[5]; qm[6] -= dx[6]; qm[7] -= dx[7]; 31 } else if (*error > 0) { 32 qm[0] += dx[0]; qm[1] += dx[1]; qm[2] += dx[2]; qm[3] += dx[3]; 33 qm[4] += dx[4]; qm[5] += dx[5]; qm[6] += dx[6]; qm[7] += dx[7]; 34 } 35 36 round += dl[0] * qm[0] + dl[1] * qm[1] + dl[2] * qm[2] + dl[3] * qm[3] + 37 dl[4] * qm[4] + dl[5] * qm[5] + dl[6] * qm[6] + dl[7] * qm[7]; 38 39 dx[0] = dx[1]; dx[1] = dx[2]; dx[2] = dx[3]; dx[3] = dx[4]; 40 dl[0] = dl[1]; dl[1] = dl[2]; dl[2] = dl[3]; dl[3] = dl[4]; 41 42 dx[4] = ((dl[4] >> 30) | 1); 43 dx[5] = ((dl[5] >> 30) | 2) & ~1; 44 dx[6] = ((dl[6] >> 30) | 2) & ~1; 45 dx[7] = ((dl[7] >> 30) | 4) & ~3; 46 47 *error = *in; 48 *in += (round >> shift); 49 50 dl[4] = -(unsigned)dl[5]; dl[5] = -(unsigned)dl[6]; 51 dl[6] = *in -(unsigned)dl[7]; dl[7] = *in; 52 dl[5] += (unsigned)dl[6]; dl[4] += (unsigned)dl[5]; 53} 54 55av_cold void ff_ttadsp_init(TTADSPContext *c) 56{ 57 c->filter_process = tta_filter_process_c; 58 59#if ARCH_X86 60 ff_ttadsp_init_x86(c); 61#endif 62} 63