开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > DVD > 查看源码
getblk.cc
资源名称：dvd-munitions.tar.gz　[点击查看]
上传用户：aoeyumen
上传日期：2007-01-06
资源大小：3329k
文件大小：9k
源码类别：

DVD

开发平台：

Unix_Linux

getblk.cc：源码内容
  1. /*
  2.   File: getblk.cc
  5. */
  8. #include "all.hh"
  11. /* defined in getvlc.hh */
  12. typedef struct {
  13.   char run, level, len;
  14. } DCTtab;
  16. extern DCTtab DCTtabfirst[],DCTtabnext[],DCTtab0[],DCTtab1[];
  17. extern DCTtab DCTtab2[],DCTtab3[],DCTtab4[],DCTtab5[],DCTtab6[];
  18. extern DCTtab DCTtab0a[],DCTtab1a[];
  21. // decode one intra coded MPEG-1 block
  23. void LayerData::getintrablock(int comp, int dc_dct_pred[]){
  24.   int val, j, sign;
  25.   unsigned int code;
  26.   DCTtab *tab;
  28.   short *bp = block[comp];
  30.   /* decode DC coefficients */
  31.   if (comp<4)         bp[0] = (dc_dct_pred[0]+=getDClum()) << 3;
  32.   else if (comp==4)   bp[0] = (dc_dct_pred[1]+=getDCchrom()) << 3;
  33.   else                bp[0] = (dc_dct_pred[2]+=getDCchrom()) << 3;
  34. #ifdef HAVE_MMX
  35.   bp[0] <<= 4;
  36. #endif
  38.   if (fault) return;
  40.   /* decode AC coefficients */
  41.   for (int i=1; ; i++){
  42.     code = (vs->showbits(32) >> 16);
  43.     if (code>=16384) tab = &DCTtabnext[(code>>12)-4];
  44.     else if (code>=1024) tab = &DCTtab0[(code>>8)-4];
  45.     else if (code>=512) tab = &DCTtab1[(code>>6)-8];
  46.     else if (code>=256) tab = &DCTtab2[(code>>4)-16];
  47.     else if (code>=128) tab = &DCTtab3[(code>>3)-16];
  48.     else if (code>=64) tab = &DCTtab4[(code>>2)-16];
  49.     else if (code>=32) tab = &DCTtab5[(code>>1)-16];
  50.     else if (code>=16) tab = &DCTtab6[code-16];
  51.     else {
  52.       if (!quiet) fprintf(stderr,"invalid Huffman code in getintrablock()n");
  53.       fault = 1;
  54.       return;
  55.     }
  57.     vs->flushbits(tab->len);
  59.     if (tab->run==64) break;  // end_of_block
  61.     if (tab->run==65){         // escape
  62.       i+= vs->getbits(6);
  63.       if ((val=vs->getbits(8))==0) val=vs->getbits(8);
  64.       else if (val==128)         val=vs->getbits(8) - 256;
  65.       else if (val>128)          val-= 256;
  66.       if ((sign=(val<0))!=0)        val= -val;
  67.     }
  68.     else {
  69.       i+= tab->run;
  70.       val = tab->level;
  71.       sign = vs->getbits(1);
  72.     }
  74.     j = zig_zag_scan[i];
  75. #ifdef HAVE_MMX
  76.     val = (val*quant_scale*intra_quantizer_matrix[j]) << 1;
  77.     /* if (val!=0) should always be true */
  78.     val = (val-16) | 16;
  79. #else
  80.     val = (val*quant_scale*intra_quantizer_matrix[j]) >> 3;
  81.     /* if (val!=0) should always be true */
  82.     val = (val-1) | 1;
  83. #endif
  84.     bp[j] = sign ? -val : val;
  85.   }
  86.   if (j != 0) {
  87.      /* not a sparse matrice ! */
  88.      sparse[comp] = 0;
  89.   }
  91. }
  94. // decode one non-intra coded MPEG-1 block
  96. void LayerData::getinterblock(int comp){
  97.   int val, j, sign;
  98.   unsigned int code;
  99.   DCTtab *tab;
  101.   short *bp = block[comp];
  103.   /* decode AC coefficients */
  104.   for (int i=0; ; i++){
  105.     code = (vs->showbits(32) >> 16);
  106.     if (code>=16384){
  107.       if (i==0) tab = &DCTtabfirst[(code>>12)-4];
  108.       else      tab = &DCTtabnext[(code>>12)-4];
  109.     }
  110.     else if (code>=1024) tab = &DCTtab0[(code>>8)-4];
  111.     else if (code>=512)  tab = &DCTtab1[(code>>6)-8];
  112.     else if (code>=256)  tab = &DCTtab2[(code>>4)-16];
  113.     else if (code>=128)  tab = &DCTtab3[(code>>3)-16];
  114.     else if (code>=64)   tab = &DCTtab4[(code>>2)-16];
  115.     else if (code>=32)   tab = &DCTtab5[(code>>1)-16];
  116.     else if (code>=16)   tab = &DCTtab6[code-16];
  117.     else {
  118.       if (!quiet)
  119.         fprintf(stderr,"invalid Huffman code in getinterblock()n");
  120.       fault = 1;
  121.       return;
  122.     }
  124.     vs->flushbits(tab->len);
  126.     if (tab->run==64)
  127.        break;   // end of block
  129.     if (tab->run==65){          // escape 
  130.       i+= vs->getbits(6);
  131.       if ((val=vs->getbits(8))==0) val=vs->getbits(8);
  132.       else if (val==128)       val=vs->getbits(8) - 256;
  133.       else if (val>128)        val-= 256;
  134.       if ((sign=(val<0))!=0)      val = -val;
  135.     }
  136.     else {
  137.       i+= tab->run;
  138.       val = tab->level;
  139.       sign = vs->getbits(1);
  140.     }
  142.     j = zig_zag_scan[i];
  143. #ifdef HAVE_MMX
  144.     val = (((val<<1)+1)*quant_scale*non_intra_quantizer_matrix[j]);
  145.     /* if (val!=0) should always be true */
  146.     val = (val-16) | 16;
  147. #else
  148.     val = (((val<<1)+1)*quant_scale*non_intra_quantizer_matrix[j]) >> 4;
  149.     /* if (val!=0) should always be true */
  150.     val = (val-1) | 1;
  151. #endif
  152.     bp[j] = sign ? -val : val;
  153.   }
  154.   if (j != 0) {
  155.      /* not a sparse matrice ! */
  156.      sparse[comp] = 0;
  157.   }
  159. }
  162. // decode one intra coded MPEG-2 block
  164. void LayerData::getmpg2intrablock(int comp, int dc_dct_pred[]){
  165.   int val, j, sign, nc;
  166.   unsigned int code;
  167.   DCTtab *tab;
  168.   short *bp;
  169.   int *qmat;
  171.   /* with data partitioning, data always goes to base layer */
  172. //  ld1 = (ld->scalable_mode==SC_DP) ? &base : ld;
  173. //  ld1=ld;
  174.   bp = block[comp];
  176. #ifdef ENHANCEMENT
  177.   if (base.scalable_mode==SC_DP)
  178.     if (base.pri_brk<64) ld = &enhan;
  179.     else ld = &base;
  180. #endif
  182.   qmat = (comp<4 || chroma_format==CHROMA420)
  183.          ? intra_quantizer_matrix
  184.          : chroma_intra_quantizer_matrix;
  186.   /* decode DC coefficients */
  187.   if (comp<4)           val = (dc_dct_pred[0]+= getDClum());
  188.   else if ((comp&1)==0) val = (dc_dct_pred[1]+= getDCchrom());
  189.   else                  val = (dc_dct_pred[2]+= getDCchrom());
  191.   if (fault) return;
  192. #ifdef HAVE_MMX
  193.   bp[0] = val << (7-dc_prec);
  194. #else
  195.   bp[0] = val << (3-dc_prec);
  196. #endif
  197.   nc=0;
  199.   /* decode AC coefficients */
  200.   for (int i=1; ; i++){
  201. //    code = vs->showbits(32) >> 16;
  202.     code = vs->showbits(16);
  203.     if (code>=16384 && !intravlc) tab = &DCTtabnext[(code>>12)-4];
  204.     else if (code>=1024){
  205.       if (intravlc)     tab = &DCTtab0a[(code>>8)-4];
  206.       else              tab = &DCTtab0[(code>>8)-4];
  207.     }
  208.     else if (code>=512){
  209.       if (intravlc)     tab = &DCTtab1a[(code>>6)-8];
  210.       else              tab = &DCTtab1[(code>>6)-8];
  211.     }
  212.     else if (code>=256) tab = &DCTtab2[(code>>4)-16];
  213.     else if (code>=128) tab = &DCTtab3[(code>>3)-16];
  214.     else if (code>=64)  tab = &DCTtab4[(code>>2)-16];
  215.     else if (code>=32)  tab = &DCTtab5[(code>>1)-16];
  216.     else if (code>=16)  tab = &DCTtab6[code-16];
  217.     else {
  218.       if (!quiet) fprintf(stderr,"invalid Huffman code in getmpg2intrablock()n");
  219.       fault = 1;
  220.       return;
  221.     }
  223.     vs->flushbits(tab->len);
  225.     if (tab->run==64)
  226.        break; // end_of_block
  228.     if (tab->run==65){        // escape
  229.       i+= vs->getbits(6);
  231.       val = vs->getbits(12);
  232.       if ((val&2047)==0){
  233.         if (!quiet) fprintf(stderr,"invalid signed_level (escape) in getmpg2intrablock()n");
  234.         fault = 1;
  235.         return;
  236.       }
  237.       if ((sign=(val>=2048))!=0) val = 4096 - val;
  238.     }
  239.     else {
  240.       i+= tab->run;
  241.       val = tab->level;
  242.       sign = vs->getbits(1);
  243.     }
  245.     j = (altscan ? alternate_scan : zig_zag_scan)[i];
  247. #ifdef HAVE_MMX
  248.     val = (val * quant_scale * qmat[j]);
  249. #else
  250.     val = (val * quant_scale * qmat[j]) >> 4;
  251. #endif
  253.     bp[j] = sign ? -val : val;
  254.     nc++;
  256. #ifdef ENHANCEMENT
  257.     if (base.scalable_mode==SC_DP && nc==base.pri_brk-63){
  258.       ld = &enhan;
  259.       incnt = ld->incnt;
  260.       bfr = vs->showbits(32);
  261.     }
  262. #endif
  263.   }
  264.   if (j != 0) {
  265.      /* not a sparse matrice ! */
  266.      sparse[comp] = 0;
  267.   }
  269. }
  272. // decode one non-intra coded MPEG-2 block
  274. void LayerData::getmpg2interblock(int comp){
  275.   int val, j, sign, nc;
  276.   unsigned int code;
  277.   DCTtab *tab;
  278.   short *bp;
  279.   int *qmat;
  281.   /* with data partitioning, data always goes to base layer */
  282. //  ld1 = (ld->scalable_mode==SC_DP) ? &base : ld;
  283. //  ld1=ld;
  284.   bp = block[comp];
  286. #ifdef ENHANCEMENT
  287.   if (base.scalable_mode==SC_DP)
  288.     if (base.pri_brk<64) ld = &enhan;
  289.     else                 ld = &base;
  290. #endif
  292.   qmat = (comp<4 || chroma_format==CHROMA420)
  293.          ? non_intra_quantizer_matrix
  294.          : chroma_non_intra_quantizer_matrix;
  296.   nc = 0;
  298.   /* decode AC coefficients */
  299.   for (int i=0; ; i++){
  300. //    code = vs->showbits(32) >> 16;
  301.     code = vs->showbits(16);
  302.     if (code>=16384){
  303.       if (i==0) tab = &DCTtabfirst[(code>>12)-4];
  304.       else      tab = &DCTtabnext[(code>>12)-4];
  305.     }
  306.     else if (code>=1024) tab = &DCTtab0[(code>>8)-4];
  307.     else if (code>=512)  tab = &DCTtab1[(code>>6)-8];
  308.     else if (code>=256)  tab = &DCTtab2[(code>>4)-16];
  309.     else if (code>=128)  tab = &DCTtab3[(code>>3)-16];
  310.     else if (code>=64)   tab = &DCTtab4[(code>>2)-16];
  311.     else if (code>=32)   tab = &DCTtab5[(code>>1)-16];
  312.     else if (code>=16)   tab = &DCTtab6[code-16];
  313.     else {
  314.       if (!quiet) fprintf(stderr,"invalid Huffman code in getmpg2interblock()n");
  315.       fault = 1;
  316.       return;
  317.     }
  319.     vs->flushbits(tab->len);
  321.     if (tab->run==64)
  322.        break;          // end_of_block
  324.     if (tab->run==65){                 // escape
  325.       i+= vs->getbits(6);
  326.       val = vs->getbits(12);
  327.       if ((val&2047)==0){
  328.         if (!quiet) fprintf(stderr,"invalid signed_level (escape) in getmpg2intrablock()n");
  329.         fault = 1;
  330.         return;
  331.       }
  332.       if ((sign=(val>=2048))!=0) val = 4096 - val;
  333.     }
  334.     else {
  335.       i+= tab->run;
  336.       val = tab->level;
  337.       sign = vs->getbits(1);
  338.     }
  340.     j = (altscan ? alternate_scan : zig_zag_scan)[i];
  341. #ifdef HAVE_MMX
  342.      val = (((val<<1)+1) * quant_scale * qmat[j]) >> 1;
  343. #else
  344.      val = (((val<<1)+1) * quant_scale * qmat[j]) >> 5;
  345. #endif
  347.     bp[j] = sign ? (-val) : val ;
  348.     nc++;
  350. /*
  351.     if (base.scalable_mode==SC_DP && nc==base.pri_brk-63){
  352.       ld = &enhan;
  353.       incnt = ld->incnt;
  354.       bfr = vs->showbits(32);
  355.     }
  356. */
  357.   }
  358.   if (j != 0) {
  359.     /* not a sparse matrice ! */
  360.     sparse[comp] = 0;
  361.   }
  362. }