LCOV - code coverage report
Current view: top level - libavcodec - hq_hqadsp.c (source / functions) Hit Total Coverage
Test: coverage.info Lines: 75 75 100.0 %
Date: 2017-12-13 02:34:56 Functions: 4 4 100.0 %

          Line data    Source code
       1             : /*
       2             :  * Canopus HQ/HQA decoder
       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             : #include <stdint.h>
      22             : 
      23             : #include "libavutil/attributes.h"
      24             : #include "libavutil/common.h"
      25             : 
      26             : #include "hq_hqadsp.h"
      27             : 
      28             : #define FIX_1_082 17734
      29             : #define FIX_1_847 30274
      30             : #define FIX_1_414 23170
      31             : #define FIX_2_613 21407 // divided by two to fit the range
      32             : 
      33             : #define IDCTMUL(a, b) ((int)((a) * (unsigned)(b)) >> 16)
      34             : 
      35     3544320 : static inline void idct_row(int16_t *blk)
      36             : {
      37             :     int tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7, tmp8, tmp9, tmpA;
      38             :     int tmpB, tmpC, tmpD, tmpE, tmpF, tmp10, tmp11, tmp12, tmp13, tmp14;
      39             : 
      40     3544320 :     tmp0 = blk[5] - blk[3];
      41     3544320 :     tmp1 = blk[5] + blk[3];
      42     3544320 :     tmp2 = blk[1] - blk[7];
      43     3544320 :     tmp3 = blk[1] + blk[7];
      44     3544320 :     tmp4 = tmp3 - tmp1;
      45     3544320 :     tmp5 = IDCTMUL(tmp0 + tmp2, FIX_1_847);
      46     3544320 :     tmp6 = IDCTMUL(tmp2,        FIX_1_082) - tmp5;
      47     3544320 :     tmp7 = tmp5 - IDCTMUL(tmp0, FIX_2_613) * 2;
      48     3544320 :     tmp8 = tmp3 + tmp1;
      49     3544320 :     tmp9 = tmp7 * 4 - tmp8;
      50     3544320 :     tmpA = IDCTMUL(tmp4, FIX_1_414) * 4 - tmp9;
      51     3544320 :     tmpB = tmp6 * 4 + tmpA;
      52     3544320 :     tmpC = blk[2] + blk[6];
      53     3544320 :     tmpD = blk[2] - blk[6];
      54     3544320 :     tmpE = blk[0] - blk[4];
      55     3544320 :     tmpF = blk[0] + blk[4];
      56             : 
      57     3544320 :     tmp10 = IDCTMUL(tmpD, FIX_1_414) * 4 - tmpC;
      58     3544320 :     tmp11 = tmpE - tmp10;
      59     3544320 :     tmp12 = tmpF - tmpC;
      60     3544320 :     tmp13 = tmpE + tmp10;
      61     3544320 :     tmp14 = tmpF + tmpC;
      62             : 
      63     3544320 :     blk[0] = tmp14 + tmp8;
      64     3544320 :     blk[1] = tmp13 + tmp9;
      65     3544320 :     blk[2] = tmp11 + tmpA;
      66     3544320 :     blk[3] = tmp12 - tmpB;
      67     3544320 :     blk[4] = tmp12 + tmpB;
      68     3544320 :     blk[5] = tmp11 - tmpA;
      69     3544320 :     blk[6] = tmp13 - tmp9;
      70     3544320 :     blk[7] = tmp14 - tmp8;
      71     3544320 : }
      72             : 
      73     3544320 : static inline void idct_col(int16_t *blk)
      74             : {
      75             :     int tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7, tmp8, tmp9, tmpA;
      76             :     int tmpB, tmpC, tmpD, tmpE, tmpF, tmp10, tmp11, tmp12, tmp13, tmp14;
      77             : 
      78     3544320 :     tmp0 = blk[5 * 8] - blk[3 * 8];
      79     3544320 :     tmp1 = blk[5 * 8] + blk[3 * 8];
      80     3544320 :     tmp2 = blk[1 * 8] * 2 - (blk[7 * 8] >> 2);
      81     3544320 :     tmp3 = blk[1 * 8] * 2 + (blk[7 * 8] >> 2);
      82     3544320 :     tmp4 = tmp3 - tmp1;
      83     3544320 :     tmp5 = IDCTMUL(tmp0 + tmp2, FIX_1_847);
      84     3544320 :     tmp6 = IDCTMUL(tmp2,        FIX_1_082) - tmp5;
      85     3544320 :     tmp7 = tmp5 - IDCTMUL(tmp0, FIX_2_613) * 2;
      86     3544320 :     tmp8 = (tmp3 + tmp1) >> 1;
      87     3544320 :     tmp9 = tmp7 * 2 - tmp8;
      88     3544320 :     tmpA = IDCTMUL(tmp4, FIX_1_414) * 2 - tmp9;
      89     3544320 :     tmpB = tmp6 * 2 + tmpA;
      90     3544320 :     tmpC =  blk[2 * 8] + (blk[6 * 8] >> 1) >> 1;
      91     3544320 :     tmpD =  blk[2 * 8] - (blk[6 * 8] >> 1);
      92     3544320 :     tmpE = (blk[0 * 8] >> 1) - (blk[4 * 8] >> 1) + 0x2020;
      93     3544320 :     tmpF = (blk[0 * 8] >> 1) + (blk[4 * 8] >> 1) + 0x2020;
      94             : 
      95     3544320 :     tmp10 = IDCTMUL(tmpD, FIX_1_414) * 2 - tmpC;
      96     3544320 :     tmp11 = tmpE - tmp10;
      97     3544320 :     tmp12 = tmpF - tmpC;
      98     3544320 :     tmp13 = tmpE + tmp10;
      99     3544320 :     tmp14 = tmpF + tmpC;
     100             : 
     101     3544320 :     blk[0 * 8] = (tmp14 + tmp8) >> 6;
     102     3544320 :     blk[1 * 8] = (tmp13 + tmp9) >> 6;
     103     3544320 :     blk[2 * 8] = (tmp11 + tmpA) >> 6;
     104     3544320 :     blk[3 * 8] = (tmp12 - tmpB) >> 6;
     105     3544320 :     blk[4 * 8] = (tmp12 + tmpB) >> 6;
     106     3544320 :     blk[5 * 8] = (tmp11 - tmpA) >> 6;
     107     3544320 :     blk[6 * 8] = (tmp13 - tmp9) >> 6;
     108     3544320 :     blk[7 * 8] = (tmp14 - tmp8) >> 6;
     109     3544320 : }
     110             : 
     111      443040 : static void hq_idct_put(uint8_t *dst, int stride, int16_t *block)
     112             : {
     113             :     int i, j;
     114             : 
     115     3987360 :     for (i = 0; i < 8; i++)
     116     3544320 :         idct_row(block + i * 8);
     117     3987360 :     for (i = 0; i < 8; i++)
     118     3544320 :         idct_col(block + i);
     119             : 
     120     3987360 :     for (i = 0; i < 8; i++) {
     121    31898880 :         for (j = 0; j < 8; j++)
     122    28354560 :             dst[j] = av_clip_uint8(block[j + i * 8]);
     123     3544320 :         dst += stride;
     124             :     }
     125      443040 : }
     126             : 
     127           6 : av_cold void ff_hqdsp_init(HQDSPContext *c)
     128             : {
     129           6 :     c->idct_put = hq_idct_put;
     130           6 : }

Generated by: LCOV version 1.13