git.sesse.net Git - ffmpeg/blob - libavcodec/bfin/mpegvideo_bfin.c

   1 /*
   2  * BlackFin MPEGVIDEO OPTIMIZATIONS
   3  *
   4  * Copyright (C) 2007 Marc Hoffman <mmh@pleasantst.com>
   5  *
   6  * This file is part of Libav.
   7  *
   8  * Libav is free software; you can redistribute it and/or
   9  * modify it under the terms of the GNU Lesser General Public
  10  * License as published by the Free Software Foundation; either
  11  * version 2.1 of the License, or (at your option) any later version.
  12  *
  13  * Libav is distributed in the hope that it will be useful,
  14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  16  * Lesser General Public License for more details.
  17  *
  18  * You should have received a copy of the GNU Lesser General Public
  19  * License along with Libav; if not, write to the Free Software
  20  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  21  */
  22
  23 #include "libavutil/attributes.h"
  24 #include "libavcodec/avcodec.h"
  25 #include "libavcodec/dsputil.h"
  26 #include "libavcodec/mpegvideo.h"
  27 #include "dsputil_bfin.h"
  28
  29 static int dct_quantize_bfin (MpegEncContext *s,
  30                               int16_t *block, int n,
  31                               int qscale, int *overflow)
  32 {
  33     int last_non_zero, q, start_i;
  34     const short *qmat;
  35     short *bias;
  36     const uint8_t *scantable= s->intra_scantable.scantable;
  37     short dc;
  38     int   max=0;
  39
  40     PROF("fdct",0);
  41     s->dsp.fdct(block);
  42     EPROF();
  43
  44     PROF("denoise",1);
  45     if(s->dct_error_sum)
  46         s->denoise_dct(s, block);
  47     EPROF();
  48
  49     PROF("quant-init",2);
  50     if (s->mb_intra) {
  51         if (!s->h263_aic) {
  52             if (n < 4)
  53                 q = s->y_dc_scale;
  54             else
  55                 q = s->c_dc_scale;
  56             q = q << 3;
  57         } else
  58             /* For AIC we skip quant/dequant of INTRADC */
  59             q = 1 << 3;
  60
  61         /* note: block[0] is assumed to be positive */
  62         dc = block[0] = (block[0] + (q >> 1)) / q;
  63         start_i = 1;
  64         last_non_zero = 0;
  65         bias = s->q_intra_matrix16[qscale][1];
  66         qmat = s->q_intra_matrix16[qscale][0];
  67
  68     } else {
  69         start_i = 0;
  70         last_non_zero = -1;
  71         bias = s->q_inter_matrix16[qscale][1];
  72         qmat = s->q_inter_matrix16[qscale][0];
  73
  74     }
  75     EPROF();
  76
  77     PROF("quantize",4);
  78
  79     /*  for(i=start_i; i<64; i++) {                           */
  80     /*      sign     = (block[i]>>15)|1;                      */
  81     /*      level    = ((abs(block[i])+bias[0])*qmat[i])>>16; */
  82     /*      if (level < 0) level = 0;                         */
  83     /*      max     |= level;                                 */
  84     /*      level    = level * sign;                          */
  85     /*      block[i] = level;                                 */
  86     /*  } */
  87
  88     __asm__ volatile
  89         ("i2=%1;\n\t"
  90          "r1=[%1++];                                                         \n\t"
  91          "r0=r1>>>15 (v);                                                    \n\t"
  92          "lsetup (0f,1f) lc0=%3;                                             \n\t"
  93          "0:   r0=r0|%4;                                                     \n\t"
  94          "     r1=abs r1 (v)                                    || r2=[%2++];\n\t"
  95          "     r1=r1+|+%5;                                                   \n\t"
  96          "     r1=max(r1,%6) (v);                                            \n\t"
  97          "     r1.h=(a1 =r1.h*r2.h), r1.l=(a0 =r1.l*r2.l) (tfu);             \n\t"
  98          "     %0=%0|r1;                                                     \n\t"
  99          "     r0.h=(a1 =r1.h*r0.h), r0.l=(a0 =r1.l*r0.l) (is)  || r1=[%1++];\n\t"
 100          "1:   r0=r1>>>15 (v)                                   || [i2++]=r0;\n\t"
 101          "r1=%0>>16;                                                         \n\t"
 102          "%0=%0|r1;                                                          \n\t"
 103          "%0.h=0;                                                            \n\t"
 104          : "=&d" (max)
 105          : "b" (block), "b" (qmat), "a" (32), "d" (0x00010001), "d" (bias[0]*0x10001), "d" (0)
 106          : "R0","R1","R2", "I2");
 107     if (start_i == 1) block[0] = dc;
 108
 109     EPROF();
 110
 111
 112     PROF("zzscan",5);
 113
 114     __asm__ volatile
 115         ("r0=b[%1--] (x);         \n\t"
 116          "lsetup (0f,1f) lc0=%3;  \n\t"     /*    for(i=63; i>=start_i; i--) { */
 117          "0: p0=r0;               \n\t"     /*        j = scantable[i];        */
 118          "   p0=%2+(p0<<1);       \n\t"     /*        if (block[j]) {          */
 119          "   r0=w[p0];            \n\t"     /*           last_non_zero = i;    */
 120          "   cc=r0==0;            \n\t"     /*           break;                */
 121          "   if !cc jump 2f;      \n\t"     /*        }                        */
 122          "1: r0=b[%1--] (x);      \n\t"     /*    }                            */
 123          "   %0=%4;               \n\t"
 124          "   jump 3f;             \n\t"
 125          "2: %0=lc0;              \n\t"
 126          "3:\n\t"
 127
 128          : "=d" (last_non_zero)
 129          : "a" (scantable+63), "a" (block), "a" (63), "d" (last_non_zero)
 130          : "P0","R0");
 131
 132     EPROF();
 133
 134     *overflow= s->max_qcoeff < max; //overflow might have happened
 135
 136     bfprof();
 137
 138     /* we need this permutation so that we correct the IDCT, we only permute the !=0 elements */
 139     if (s->dsp.idct_permutation_type != FF_NO_IDCT_PERM)
 140         ff_block_permute(block, s->dsp.idct_permutation, scantable, last_non_zero);
 141
 142     return last_non_zero;
 143 }
 144
 145 av_cold void ff_MPV_common_init_bfin (MpegEncContext *s)
 146 {
 147 /*     s->dct_quantize= dct_quantize_bfin; */
 148 }