source: 3DVCSoftware/branches/HTM-DEV-0.3-dev2a/source/Lib/TLibEncoder/TEncCu.cpp @ 465

Last change on this file since 465 was 465, checked in by tech, 11 years ago

Fixed some minor memory leaks.

  • Property svn:eol-style set to native
File size: 70.8 KB
Line 
1/* The copyright in this software is being made available under the BSD
2 * License, included below. This software may be subject to other third party
3 * and contributor rights, including patent rights, and no such rights are
4 * granted under this license. 
5 *
6 * Copyright (c) 2010-2013, ITU/ISO/IEC
7 * All rights reserved.
8 *
9 * Redistribution and use in source and binary forms, with or without
10 * modification, are permitted provided that the following conditions are met:
11 *
12 *  * Redistributions of source code must retain the above copyright notice,
13 *    this list of conditions and the following disclaimer.
14 *  * Redistributions in binary form must reproduce the above copyright notice,
15 *    this list of conditions and the following disclaimer in the documentation
16 *    and/or other materials provided with the distribution.
17 *  * Neither the name of the ITU/ISO/IEC nor the names of its contributors may
18 *    be used to endorse or promote products derived from this software without
19 *    specific prior written permission.
20 *
21 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
22 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
23 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
24 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS
25 * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
26 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
27 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
28 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
29 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
30 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
31 * THE POSSIBILITY OF SUCH DAMAGE.
32 */
33
34/** \file     TEncCu.cpp
35    \brief    Coding Unit (CU) encoder class
36*/
37
38#include <stdio.h>
39#include "TEncTop.h"
40#include "TEncCu.h"
41#include "TEncAnalyze.h"
42
43#include <cmath>
44#include <algorithm>
45using namespace std;
46
47//! \ingroup TLibEncoder
48//! \{
49
50// ====================================================================================================================
51// Constructor / destructor / create / destroy
52// ====================================================================================================================
53
54/**
55 \param    uiTotalDepth  total number of allowable depth
56 \param    uiMaxWidth    largest CU width
57 \param    uiMaxHeight   largest CU height
58 */
59Void TEncCu::create(UChar uhTotalDepth, UInt uiMaxWidth, UInt uiMaxHeight)
60{
61  Int i;
62 
63  m_uhTotalDepth   = uhTotalDepth + 1;
64  m_ppcBestCU      = new TComDataCU*[m_uhTotalDepth-1];
65  m_ppcTempCU      = new TComDataCU*[m_uhTotalDepth-1];
66   
67#if H_3D_ARP
68  m_ppcWeightedTempCU = new TComDataCU*[m_uhTotalDepth-1];
69#endif
70
71  m_ppcPredYuvBest = new TComYuv*[m_uhTotalDepth-1];
72  m_ppcResiYuvBest = new TComYuv*[m_uhTotalDepth-1];
73  m_ppcRecoYuvBest = new TComYuv*[m_uhTotalDepth-1];
74  m_ppcPredYuvTemp = new TComYuv*[m_uhTotalDepth-1];
75  m_ppcResiYuvTemp = new TComYuv*[m_uhTotalDepth-1];
76  m_ppcRecoYuvTemp = new TComYuv*[m_uhTotalDepth-1];
77  m_ppcOrigYuv     = new TComYuv*[m_uhTotalDepth-1];
78 
79  UInt uiNumPartitions;
80  for( i=0 ; i<m_uhTotalDepth-1 ; i++)
81  {
82    uiNumPartitions = 1<<( ( m_uhTotalDepth - i - 1 )<<1 );
83    UInt uiWidth  = uiMaxWidth  >> i;
84    UInt uiHeight = uiMaxHeight >> i;
85   
86    m_ppcBestCU[i] = new TComDataCU; m_ppcBestCU[i]->create( uiNumPartitions, uiWidth, uiHeight, false, uiMaxWidth >> (m_uhTotalDepth - 1) );
87    m_ppcTempCU[i] = new TComDataCU; m_ppcTempCU[i]->create( uiNumPartitions, uiWidth, uiHeight, false, uiMaxWidth >> (m_uhTotalDepth - 1) );
88   
89#if H_3D_ARP
90    m_ppcWeightedTempCU[i] = new TComDataCU; m_ppcWeightedTempCU[i]->create( uiNumPartitions, uiWidth, uiHeight, false, uiMaxWidth >> (m_uhTotalDepth - 1) );
91#endif 
92
93    m_ppcPredYuvBest[i] = new TComYuv; m_ppcPredYuvBest[i]->create(uiWidth, uiHeight);
94    m_ppcResiYuvBest[i] = new TComYuv; m_ppcResiYuvBest[i]->create(uiWidth, uiHeight);
95    m_ppcRecoYuvBest[i] = new TComYuv; m_ppcRecoYuvBest[i]->create(uiWidth, uiHeight);
96   
97    m_ppcPredYuvTemp[i] = new TComYuv; m_ppcPredYuvTemp[i]->create(uiWidth, uiHeight);
98    m_ppcResiYuvTemp[i] = new TComYuv; m_ppcResiYuvTemp[i]->create(uiWidth, uiHeight);
99    m_ppcRecoYuvTemp[i] = new TComYuv; m_ppcRecoYuvTemp[i]->create(uiWidth, uiHeight);
100   
101    m_ppcOrigYuv    [i] = new TComYuv; m_ppcOrigYuv    [i]->create(uiWidth, uiHeight);
102  }
103 
104  m_bEncodeDQP = false;
105#if RATE_CONTROL_LAMBDA_DOMAIN
106  m_LCUPredictionSAD = 0;
107  m_addSADDepth      = 0;
108  m_temporalSAD      = 0;
109#endif
110
111  // initialize partition order.
112  UInt* piTmp = &g_auiZscanToRaster[0];
113  initZscanToRaster( m_uhTotalDepth, 1, 0, piTmp);
114  initRasterToZscan( uiMaxWidth, uiMaxHeight, m_uhTotalDepth );
115 
116  // initialize conversion matrix from partition index to pel
117  initRasterToPelXY( uiMaxWidth, uiMaxHeight, m_uhTotalDepth );
118}
119
120Void TEncCu::destroy()
121{
122  Int i;
123 
124  for( i=0 ; i<m_uhTotalDepth-1 ; i++)
125  {
126    if(m_ppcBestCU[i])
127    {
128      m_ppcBestCU[i]->destroy();      delete m_ppcBestCU[i];      m_ppcBestCU[i] = NULL;
129    }
130    if(m_ppcTempCU[i])
131    {
132      m_ppcTempCU[i]->destroy();      delete m_ppcTempCU[i];      m_ppcTempCU[i] = NULL;
133    }
134#if H_3D_ARP
135    if(m_ppcWeightedTempCU[i])
136    {
137      m_ppcWeightedTempCU[i]->destroy(); delete m_ppcWeightedTempCU[i]; m_ppcWeightedTempCU[i] = NULL;
138    }
139#endif
140    if(m_ppcPredYuvBest[i])
141    {
142      m_ppcPredYuvBest[i]->destroy(); delete m_ppcPredYuvBest[i]; m_ppcPredYuvBest[i] = NULL;
143    }
144    if(m_ppcResiYuvBest[i])
145    {
146      m_ppcResiYuvBest[i]->destroy(); delete m_ppcResiYuvBest[i]; m_ppcResiYuvBest[i] = NULL;
147    }
148    if(m_ppcRecoYuvBest[i])
149    {
150      m_ppcRecoYuvBest[i]->destroy(); delete m_ppcRecoYuvBest[i]; m_ppcRecoYuvBest[i] = NULL;
151    }
152    if(m_ppcPredYuvTemp[i])
153    {
154      m_ppcPredYuvTemp[i]->destroy(); delete m_ppcPredYuvTemp[i]; m_ppcPredYuvTemp[i] = NULL;
155    }
156    if(m_ppcResiYuvTemp[i])
157    {
158      m_ppcResiYuvTemp[i]->destroy(); delete m_ppcResiYuvTemp[i]; m_ppcResiYuvTemp[i] = NULL;
159    }
160    if(m_ppcRecoYuvTemp[i])
161    {
162      m_ppcRecoYuvTemp[i]->destroy(); delete m_ppcRecoYuvTemp[i]; m_ppcRecoYuvTemp[i] = NULL;
163    }
164    if(m_ppcOrigYuv[i])
165    {
166      m_ppcOrigYuv[i]->destroy();     delete m_ppcOrigYuv[i];     m_ppcOrigYuv[i] = NULL;
167    }
168  }
169  if(m_ppcBestCU)
170  {
171    delete [] m_ppcBestCU;
172    m_ppcBestCU = NULL;
173  }
174  if(m_ppcTempCU)
175  {
176    delete [] m_ppcTempCU;
177    m_ppcTempCU = NULL;
178  }
179
180#if H_3D_ARP
181  if(m_ppcWeightedTempCU)
182  {
183    delete [] m_ppcWeightedTempCU; 
184    m_ppcWeightedTempCU = NULL; 
185  }
186#endif
187  if(m_ppcPredYuvBest)
188  {
189    delete [] m_ppcPredYuvBest;
190    m_ppcPredYuvBest = NULL;
191  }
192  if(m_ppcResiYuvBest)
193  {
194    delete [] m_ppcResiYuvBest;
195    m_ppcResiYuvBest = NULL;
196  }
197  if(m_ppcRecoYuvBest)
198  {
199    delete [] m_ppcRecoYuvBest;
200    m_ppcRecoYuvBest = NULL;
201  }
202  if(m_ppcPredYuvTemp)
203  {
204    delete [] m_ppcPredYuvTemp;
205    m_ppcPredYuvTemp = NULL;
206  }
207  if(m_ppcResiYuvTemp)
208  {
209    delete [] m_ppcResiYuvTemp;
210    m_ppcResiYuvTemp = NULL;
211  }
212  if(m_ppcRecoYuvTemp)
213  {
214    delete [] m_ppcRecoYuvTemp;
215    m_ppcRecoYuvTemp = NULL;
216  }
217  if(m_ppcOrigYuv)
218  {
219    delete [] m_ppcOrigYuv;
220    m_ppcOrigYuv = NULL;
221  }
222}
223
224/** \param    pcEncTop      pointer of encoder class
225 */
226Void TEncCu::init( TEncTop* pcEncTop )
227{
228  m_pcEncCfg           = pcEncTop;
229  m_pcPredSearch       = pcEncTop->getPredSearch();
230  m_pcTrQuant          = pcEncTop->getTrQuant();
231  m_pcBitCounter       = pcEncTop->getBitCounter();
232  m_pcRdCost           = pcEncTop->getRdCost();
233 
234  m_pcEntropyCoder     = pcEncTop->getEntropyCoder();
235  m_pcCavlcCoder       = pcEncTop->getCavlcCoder();
236  m_pcSbacCoder       = pcEncTop->getSbacCoder();
237  m_pcBinCABAC         = pcEncTop->getBinCABAC();
238 
239  m_pppcRDSbacCoder   = pcEncTop->getRDSbacCoder();
240  m_pcRDGoOnSbacCoder = pcEncTop->getRDGoOnSbacCoder();
241 
242  m_bUseSBACRD        = pcEncTop->getUseSBACRD();
243  m_pcRateCtrl        = pcEncTop->getRateCtrl();
244}
245
246// ====================================================================================================================
247// Public member functions
248// ====================================================================================================================
249
250/** \param  rpcCU pointer of CU data class
251 */
252Void TEncCu::compressCU( TComDataCU*& rpcCU )
253{
254  // initialize CU data
255  m_ppcBestCU[0]->initCU( rpcCU->getPic(), rpcCU->getAddr() );
256  m_ppcTempCU[0]->initCU( rpcCU->getPic(), rpcCU->getAddr() );
257
258#if RATE_CONTROL_LAMBDA_DOMAIN
259  m_addSADDepth      = 0;
260  m_LCUPredictionSAD = 0;
261  m_temporalSAD      = 0;
262#endif
263
264  // analysis of CU
265  xCompressCU( m_ppcBestCU[0], m_ppcTempCU[0], 0 );
266
267#if ADAPTIVE_QP_SELECTION
268  if( m_pcEncCfg->getUseAdaptQpSelect() )
269  {
270    if(rpcCU->getSlice()->getSliceType()!=I_SLICE) //IIII
271    {
272      xLcuCollectARLStats( rpcCU);
273    }
274  }
275#endif
276}
277/** \param  pcCU  pointer of CU data class
278 */
279Void TEncCu::encodeCU ( TComDataCU* pcCU )
280{
281  if ( pcCU->getSlice()->getPPS()->getUseDQP() )
282  {
283    setdQPFlag(true);
284  }
285
286  // Encode CU data
287  xEncodeCU( pcCU, 0, 0 );
288}
289
290// ====================================================================================================================
291// Protected member functions
292// ====================================================================================================================
293/** Derive small set of test modes for AMP encoder speed-up
294 *\param   rpcBestCU
295 *\param   eParentPartSize
296 *\param   bTestAMP_Hor
297 *\param   bTestAMP_Ver
298 *\param   bTestMergeAMP_Hor
299 *\param   bTestMergeAMP_Ver
300 *\returns Void
301*/
302#if AMP_ENC_SPEEDUP
303#if AMP_MRG
304Void TEncCu::deriveTestModeAMP (TComDataCU *&rpcBestCU, PartSize eParentPartSize, Bool &bTestAMP_Hor, Bool &bTestAMP_Ver, Bool &bTestMergeAMP_Hor, Bool &bTestMergeAMP_Ver)
305#else
306Void TEncCu::deriveTestModeAMP (TComDataCU *&rpcBestCU, PartSize eParentPartSize, Bool &bTestAMP_Hor, Bool &bTestAMP_Ver)
307#endif
308{
309  if ( rpcBestCU->getPartitionSize(0) == SIZE_2NxN )
310  {
311    bTestAMP_Hor = true;
312  }
313  else if ( rpcBestCU->getPartitionSize(0) == SIZE_Nx2N )
314  {
315    bTestAMP_Ver = true;
316  }
317  else if ( rpcBestCU->getPartitionSize(0) == SIZE_2Nx2N && rpcBestCU->getMergeFlag(0) == false && rpcBestCU->isSkipped(0) == false )
318  {
319    bTestAMP_Hor = true;         
320    bTestAMP_Ver = true;         
321  }
322
323#if AMP_MRG
324  //! Utilizing the partition size of parent PU   
325  if ( eParentPartSize >= SIZE_2NxnU && eParentPartSize <= SIZE_nRx2N )
326  { 
327    bTestMergeAMP_Hor = true;
328    bTestMergeAMP_Ver = true;
329  }
330
331  if ( eParentPartSize == SIZE_NONE ) //! if parent is intra
332  {
333    if ( rpcBestCU->getPartitionSize(0) == SIZE_2NxN )
334    {
335      bTestMergeAMP_Hor = true;
336    }
337    else if ( rpcBestCU->getPartitionSize(0) == SIZE_Nx2N )
338    {
339      bTestMergeAMP_Ver = true;
340    }
341  }
342
343  if ( rpcBestCU->getPartitionSize(0) == SIZE_2Nx2N && rpcBestCU->isSkipped(0) == false )
344  {
345    bTestMergeAMP_Hor = true;         
346    bTestMergeAMP_Ver = true;         
347  }
348
349  if ( rpcBestCU->getWidth(0) == 64 )
350  { 
351    bTestAMP_Hor = false;
352    bTestAMP_Ver = false;
353  }   
354#else
355  //! Utilizing the partition size of parent PU       
356  if ( eParentPartSize >= SIZE_2NxnU && eParentPartSize <= SIZE_nRx2N )
357  { 
358    bTestAMP_Hor = true;
359    bTestAMP_Ver = true;
360  }
361
362  if ( eParentPartSize == SIZE_2Nx2N )
363  { 
364    bTestAMP_Hor = false;
365    bTestAMP_Ver = false;
366  }     
367#endif
368}
369#endif
370
371// ====================================================================================================================
372// Protected member functions
373// ====================================================================================================================
374/** Compress a CU block recursively with enabling sub-LCU-level delta QP
375 *\param   rpcBestCU
376 *\param   rpcTempCU
377 *\param   uiDepth
378 *\returns Void
379 *
380 *- for loop of QP value to compress the current CU with all possible QP
381*/
382#if AMP_ENC_SPEEDUP
383Void TEncCu::xCompressCU( TComDataCU*& rpcBestCU, TComDataCU*& rpcTempCU, UInt uiDepth, PartSize eParentPartSize )
384#else
385Void TEncCu::xCompressCU( TComDataCU*& rpcBestCU, TComDataCU*& rpcTempCU, UInt uiDepth )
386#endif
387{
388  TComPic* pcPic = rpcBestCU->getPic();
389
390  // get Original YUV data from picture
391  m_ppcOrigYuv[uiDepth]->copyFromPicYuv( pcPic->getPicYuvOrg(), rpcBestCU->getAddr(), rpcBestCU->getZorderIdxInCU() );
392
393  // variables for fast encoder decision
394  Bool    bEarlySkip  = false;
395  Bool    bTrySplit    = true;
396  Double  fRD_Skip    = MAX_DOUBLE;
397
398  // variable for Early CU determination
399  Bool    bSubBranch = true;
400
401  // variable for Cbf fast mode PU decision
402  Bool    doNotBlockPu = true;
403  Bool earlyDetectionSkipMode = false;
404
405  Bool    bTrySplitDQP  = true;
406
407  static  Double  afCost[ MAX_CU_DEPTH ];
408  static  Int      aiNum [ MAX_CU_DEPTH ];
409
410  if ( rpcBestCU->getAddr() == 0 )
411  {
412    ::memset( afCost, 0, sizeof( afCost ) );
413    ::memset( aiNum,  0, sizeof( aiNum  ) );
414  }
415
416  Bool bBoundary = false;
417  UInt uiLPelX   = rpcBestCU->getCUPelX();
418  UInt uiRPelX   = uiLPelX + rpcBestCU->getWidth(0)  - 1;
419  UInt uiTPelY   = rpcBestCU->getCUPelY();
420  UInt uiBPelY   = uiTPelY + rpcBestCU->getHeight(0) - 1;
421
422  Int iBaseQP = xComputeQP( rpcBestCU, uiDepth );
423  Int iMinQP;
424  Int iMaxQP;
425  Bool isAddLowestQP = false;
426  Int lowestQP = -rpcTempCU->getSlice()->getSPS()->getQpBDOffsetY();
427
428  if( (g_uiMaxCUWidth>>uiDepth) >= rpcTempCU->getSlice()->getPPS()->getMinCuDQPSize() )
429  {
430    Int idQP = m_pcEncCfg->getMaxDeltaQP();
431    iMinQP = Clip3( -rpcTempCU->getSlice()->getSPS()->getQpBDOffsetY(), MAX_QP, iBaseQP-idQP );
432    iMaxQP = Clip3( -rpcTempCU->getSlice()->getSPS()->getQpBDOffsetY(), MAX_QP, iBaseQP+idQP );
433    if ( (rpcTempCU->getSlice()->getSPS()->getUseLossless()) && (lowestQP < iMinQP) && rpcTempCU->getSlice()->getPPS()->getUseDQP() )
434    {
435      isAddLowestQP = true; 
436      iMinQP = iMinQP - 1;
437    }
438  }
439  else
440  {
441    iMinQP = rpcTempCU->getQP(0);
442    iMaxQP = rpcTempCU->getQP(0);
443  }
444
445#if RATE_CONTROL_LAMBDA_DOMAIN
446  if ( m_pcEncCfg->getUseRateCtrl() )
447  {
448    iMinQP = m_pcRateCtrl->getRCQP();
449    iMaxQP = m_pcRateCtrl->getRCQP();
450  }
451#else
452  if(m_pcEncCfg->getUseRateCtrl())
453  {
454    Int qp = m_pcRateCtrl->getUnitQP();
455    iMinQP  = Clip3( MIN_QP, MAX_QP, qp);
456    iMaxQP  = Clip3( MIN_QP, MAX_QP, qp);
457  }
458#endif
459
460  // If slice start or slice end is within this cu...
461  TComSlice * pcSlice = rpcTempCU->getPic()->getSlice(rpcTempCU->getPic()->getCurrSliceIdx());
462  Bool bSliceStart = pcSlice->getSliceSegmentCurStartCUAddr()>rpcTempCU->getSCUAddr()&&pcSlice->getSliceSegmentCurStartCUAddr()<rpcTempCU->getSCUAddr()+rpcTempCU->getTotalNumPart();
463  Bool bSliceEnd = (pcSlice->getSliceSegmentCurEndCUAddr()>rpcTempCU->getSCUAddr()&&pcSlice->getSliceSegmentCurEndCUAddr()<rpcTempCU->getSCUAddr()+rpcTempCU->getTotalNumPart());
464  Bool bInsidePicture = ( uiRPelX < rpcBestCU->getSlice()->getSPS()->getPicWidthInLumaSamples() ) && ( uiBPelY < rpcBestCU->getSlice()->getSPS()->getPicHeightInLumaSamples() );
465  // We need to split, so don't try these modes.
466  if(!bSliceEnd && !bSliceStart && bInsidePicture )
467  {
468    for (Int iQP=iMinQP; iQP<=iMaxQP; iQP++)
469    {
470      if (isAddLowestQP && (iQP == iMinQP))
471      {
472        iQP = lowestQP;
473      }
474      // variables for fast encoder decision
475      bEarlySkip  = false;
476      bTrySplit    = true;
477      fRD_Skip    = MAX_DOUBLE;
478
479      rpcTempCU->initEstData( uiDepth, iQP );
480#if H_3D_NBDV
481      DisInfo DvInfo; 
482      DvInfo.bDV = false;
483      DvInfo.m_acNBDV.setZero();
484      DvInfo.m_aVIdxCan = 0;
485#if H_3D_NBDV_REF
486      DvInfo.m_acDoNBDV.setZero();
487#endif
488
489      if( rpcTempCU->getSlice()->getSliceType() != I_SLICE )
490      {
491#if H_3D_ARP
492        if( rpcTempCU->getSlice()->getVPS()->getUseAdvRP(rpcTempCU->getSlice()->getLayerId()) )
493#else
494        if(rpcTempCU->getSlice()->getViewIndex() && !rpcTempCU->getSlice()->getIsDepth()) //Notes from QC: this condition shall be changed once the configuration is completed, e.g. in pcSlice->getSPS()->getMultiviewMvPredMode() || ARP in prev. HTM. Remove this comment once it is done. Remove this comment once it is done.
495#endif
496        {
497          PartSize ePartTemp = rpcTempCU->getPartitionSize(0);
498          rpcTempCU->setPartSizeSubParts( SIZE_2Nx2N, 0, uiDepth );     
499#if H_3D_NBDV_REF
500          if(m_pcSlice->getSPS()->getUseDVPRefine())  //Notes from QC: please check the condition for DoNBDV. Remove this comment once it is done.
501            DvInfo.bDV = rpcTempCU->getDisMvpCandNBDV(&DvInfo, true);
502          else
503#endif
504            DvInfo.bDV = rpcTempCU->getDisMvpCandNBDV(&DvInfo);
505
506          rpcTempCU->setDvInfoSubParts(DvInfo, 0, uiDepth);
507          rpcBestCU->setDvInfoSubParts(DvInfo, 0, uiDepth);
508          rpcTempCU->setPartSizeSubParts( ePartTemp, 0, uiDepth );
509        }
510      }
511#endif
512      // do inter modes, SKIP and 2Nx2N
513      if( rpcBestCU->getSlice()->getSliceType() != I_SLICE )
514      {
515        // 2Nx2N
516        if(m_pcEncCfg->getUseEarlySkipDetection())
517        {
518          xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_2Nx2N );  rpcTempCU->initEstData( uiDepth, iQP );//by Competition for inter_2Nx2N
519        }
520        // SKIP
521        xCheckRDCostMerge2Nx2N( rpcBestCU, rpcTempCU, &earlyDetectionSkipMode );//by Merge for inter_2Nx2N
522        rpcTempCU->initEstData( uiDepth, iQP );
523
524        // fast encoder decision for early skip
525        if ( m_pcEncCfg->getUseFastEnc() )
526        {
527          Int iIdx = g_aucConvertToBit[ rpcBestCU->getWidth(0) ];
528          if ( aiNum [ iIdx ] > 5 && fRD_Skip < EARLY_SKIP_THRES*afCost[ iIdx ]/aiNum[ iIdx ] )
529          {
530            bEarlySkip = true;
531            bTrySplit  = false;
532          }
533        }
534
535        if(!m_pcEncCfg->getUseEarlySkipDetection())
536        {
537          // 2Nx2N, NxN
538          if ( !bEarlySkip )
539          {
540            xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_2Nx2N );  rpcTempCU->initEstData( uiDepth, iQP );
541            if(m_pcEncCfg->getUseCbfFastMode())
542            {
543              doNotBlockPu = rpcBestCU->getQtRootCbf( 0 ) != 0;
544            }
545          }
546        }
547      }
548
549      if( (g_uiMaxCUWidth>>uiDepth) >= rpcTempCU->getSlice()->getPPS()->getMinCuDQPSize() )
550      {
551        if(iQP == iBaseQP)
552        {
553          bTrySplitDQP = bTrySplit;
554        }
555      }
556      else
557      {
558        bTrySplitDQP = bTrySplit;
559      }
560      if (isAddLowestQP && (iQP == lowestQP))
561      {
562        iQP = iMinQP;
563      }
564    }
565
566#if RATE_CONTROL_LAMBDA_DOMAIN
567    if ( uiDepth <= m_addSADDepth )
568    {
569      m_LCUPredictionSAD += m_temporalSAD;
570      m_addSADDepth = uiDepth;
571    }
572#endif
573
574    if(!earlyDetectionSkipMode)
575    {
576      for (Int iQP=iMinQP; iQP<=iMaxQP; iQP++)
577      {
578        if (isAddLowestQP && (iQP == iMinQP))
579        {
580          iQP = lowestQP;
581        }
582        rpcTempCU->initEstData( uiDepth, iQP );
583
584        // do inter modes, NxN, 2NxN, and Nx2N
585        if( rpcBestCU->getSlice()->getSliceType() != I_SLICE )
586        {
587          // 2Nx2N, NxN
588          if ( !bEarlySkip )
589          {
590            if(!( (rpcBestCU->getWidth(0)==8) && (rpcBestCU->getHeight(0)==8) ))
591            {
592              if( uiDepth == g_uiMaxCUDepth - g_uiAddCUDepth && doNotBlockPu)
593              {
594                xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_NxN   );
595                rpcTempCU->initEstData( uiDepth, iQP );
596              }
597            }
598          }
599
600          // 2NxN, Nx2N
601          if(doNotBlockPu)
602          {
603            xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_Nx2N  );
604            rpcTempCU->initEstData( uiDepth, iQP );
605            if(m_pcEncCfg->getUseCbfFastMode() && rpcBestCU->getPartitionSize(0) == SIZE_Nx2N )
606            {
607              doNotBlockPu = rpcBestCU->getQtRootCbf( 0 ) != 0;
608            }
609          }
610          if(doNotBlockPu)
611          {
612            xCheckRDCostInter      ( rpcBestCU, rpcTempCU, SIZE_2NxN  );
613            rpcTempCU->initEstData( uiDepth, iQP );
614            if(m_pcEncCfg->getUseCbfFastMode() && rpcBestCU->getPartitionSize(0) == SIZE_2NxN)
615            {
616              doNotBlockPu = rpcBestCU->getQtRootCbf( 0 ) != 0;
617            }
618          }
619
620#if 1
621          //! Try AMP (SIZE_2NxnU, SIZE_2NxnD, SIZE_nLx2N, SIZE_nRx2N)
622          if( pcPic->getSlice(0)->getSPS()->getAMPAcc(uiDepth) )
623          {
624#if AMP_ENC_SPEEDUP       
625            Bool bTestAMP_Hor = false, bTestAMP_Ver = false;
626
627#if AMP_MRG
628            Bool bTestMergeAMP_Hor = false, bTestMergeAMP_Ver = false;
629
630            deriveTestModeAMP (rpcBestCU, eParentPartSize, bTestAMP_Hor, bTestAMP_Ver, bTestMergeAMP_Hor, bTestMergeAMP_Ver);
631#else
632            deriveTestModeAMP (rpcBestCU, eParentPartSize, bTestAMP_Hor, bTestAMP_Ver);
633#endif
634
635            //! Do horizontal AMP
636            if ( bTestAMP_Hor )
637            {
638              if(doNotBlockPu)
639              {
640                xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_2NxnU );
641                rpcTempCU->initEstData( uiDepth, iQP );
642                if(m_pcEncCfg->getUseCbfFastMode() && rpcBestCU->getPartitionSize(0) == SIZE_2NxnU )
643                {
644                  doNotBlockPu = rpcBestCU->getQtRootCbf( 0 ) != 0;
645                }
646              }
647              if(doNotBlockPu)
648              {
649                xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_2NxnD );
650                rpcTempCU->initEstData( uiDepth, iQP );
651                if(m_pcEncCfg->getUseCbfFastMode() && rpcBestCU->getPartitionSize(0) == SIZE_2NxnD )
652                {
653                  doNotBlockPu = rpcBestCU->getQtRootCbf( 0 ) != 0;
654                }
655              }
656            }
657#if AMP_MRG
658            else if ( bTestMergeAMP_Hor ) 
659            {
660              if(doNotBlockPu)
661              {
662                xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_2NxnU, true );
663                rpcTempCU->initEstData( uiDepth, iQP );
664                if(m_pcEncCfg->getUseCbfFastMode() && rpcBestCU->getPartitionSize(0) == SIZE_2NxnU )
665                {
666                  doNotBlockPu = rpcBestCU->getQtRootCbf( 0 ) != 0;
667                }
668              }
669              if(doNotBlockPu)
670              {
671                xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_2NxnD, true );
672                rpcTempCU->initEstData( uiDepth, iQP );
673                if(m_pcEncCfg->getUseCbfFastMode() && rpcBestCU->getPartitionSize(0) == SIZE_2NxnD )
674                {
675                  doNotBlockPu = rpcBestCU->getQtRootCbf( 0 ) != 0;
676                }
677              }
678            }
679#endif
680
681            //! Do horizontal AMP
682            if ( bTestAMP_Ver )
683            {
684              if(doNotBlockPu)
685              {
686                xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_nLx2N );
687                rpcTempCU->initEstData( uiDepth, iQP );
688                if(m_pcEncCfg->getUseCbfFastMode() && rpcBestCU->getPartitionSize(0) == SIZE_nLx2N )
689                {
690                  doNotBlockPu = rpcBestCU->getQtRootCbf( 0 ) != 0;
691                }
692              }
693              if(doNotBlockPu)
694              {
695                xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_nRx2N );
696                rpcTempCU->initEstData( uiDepth, iQP );
697              }
698            }
699#if AMP_MRG
700            else if ( bTestMergeAMP_Ver )
701            {
702              if(doNotBlockPu)
703              {
704                xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_nLx2N, true );
705                rpcTempCU->initEstData( uiDepth, iQP );
706                if(m_pcEncCfg->getUseCbfFastMode() && rpcBestCU->getPartitionSize(0) == SIZE_nLx2N )
707                {
708                  doNotBlockPu = rpcBestCU->getQtRootCbf( 0 ) != 0;
709                }
710              }
711              if(doNotBlockPu)
712              {
713                xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_nRx2N, true );
714                rpcTempCU->initEstData( uiDepth, iQP );
715              }
716            }
717#endif
718
719#else
720            xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_2NxnU );
721            rpcTempCU->initEstData( uiDepth, iQP );
722            xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_2NxnD );
723            rpcTempCU->initEstData( uiDepth, iQP );
724            xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_nLx2N );
725            rpcTempCU->initEstData( uiDepth, iQP );
726
727            xCheckRDCostInter( rpcBestCU, rpcTempCU, SIZE_nRx2N );
728            rpcTempCU->initEstData( uiDepth, iQP );
729
730#endif
731          }   
732#endif
733        }
734
735        // do normal intra modes
736        if ( !bEarlySkip )
737        {
738          // speedup for inter frames
739          if( rpcBestCU->getSlice()->getSliceType() == I_SLICE || 
740            rpcBestCU->getCbf( 0, TEXT_LUMA     ) != 0   ||
741            rpcBestCU->getCbf( 0, TEXT_CHROMA_U ) != 0   ||
742            rpcBestCU->getCbf( 0, TEXT_CHROMA_V ) != 0     ) // avoid very complex intra if it is unlikely
743          {
744            xCheckRDCostIntra( rpcBestCU, rpcTempCU, SIZE_2Nx2N );
745            rpcTempCU->initEstData( uiDepth, iQP );
746            if( uiDepth == g_uiMaxCUDepth - g_uiAddCUDepth )
747            {
748              if( rpcTempCU->getWidth(0) > ( 1 << rpcTempCU->getSlice()->getSPS()->getQuadtreeTULog2MinSize() ) )
749              {
750                xCheckRDCostIntra( rpcBestCU, rpcTempCU, SIZE_NxN   );
751                rpcTempCU->initEstData( uiDepth, iQP );
752              }
753            }
754          }
755        }
756
757        // test PCM
758        if(pcPic->getSlice(0)->getSPS()->getUsePCM()
759          && rpcTempCU->getWidth(0) <= (1<<pcPic->getSlice(0)->getSPS()->getPCMLog2MaxSize())
760          && rpcTempCU->getWidth(0) >= (1<<pcPic->getSlice(0)->getSPS()->getPCMLog2MinSize()) )
761        {
762          UInt uiRawBits = (2 * g_bitDepthY + g_bitDepthC) * rpcBestCU->getWidth(0) * rpcBestCU->getHeight(0) / 2;
763          UInt uiBestBits = rpcBestCU->getTotalBits();
764#if H_3D_VSO // M7
765          Double dRDCostTemp = m_pcRdCost->getUseVSO() ? m_pcRdCost->calcRdCostVSO(uiRawBits, 0) : m_pcRdCost->calcRdCost(uiRawBits, 0);
766          if((uiBestBits > uiRawBits) || (rpcBestCU->getTotalCost() > dRDCostTemp ))
767#else
768          if((uiBestBits > uiRawBits) || (rpcBestCU->getTotalCost() > m_pcRdCost->calcRdCost(uiRawBits, 0)))
769#endif
770          {
771            xCheckIntraPCM (rpcBestCU, rpcTempCU);
772            rpcTempCU->initEstData( uiDepth, iQP );
773          }
774        }
775        if (isAddLowestQP && (iQP == lowestQP))
776        {
777          iQP = iMinQP;
778        }
779      }
780    }
781
782    m_pcEntropyCoder->resetBits();
783    m_pcEntropyCoder->encodeSplitFlag( rpcBestCU, 0, uiDepth, true );
784    rpcBestCU->getTotalBits() += m_pcEntropyCoder->getNumberOfWrittenBits(); // split bits
785    if(m_pcEncCfg->getUseSBACRD())
786    {
787      rpcBestCU->getTotalBins() += ((TEncBinCABAC *)((TEncSbac*)m_pcEntropyCoder->m_pcEntropyCoderIf)->getEncBinIf())->getBinsCoded();
788    }
789
790#if H_3D_VSO // M8
791    if ( m_pcRdCost->getUseVSO() )   
792      rpcBestCU->getTotalCost()  = m_pcRdCost->calcRdCostVSO( rpcBestCU->getTotalBits(), rpcBestCU->getTotalDistortion() );   
793    else
794#endif
795    rpcBestCU->getTotalCost()  = m_pcRdCost->calcRdCost( rpcBestCU->getTotalBits(), rpcBestCU->getTotalDistortion() );
796
797    // accumulate statistics for early skip
798    if ( m_pcEncCfg->getUseFastEnc() )
799    {
800      if ( rpcBestCU->isSkipped(0) )
801      {
802        Int iIdx = g_aucConvertToBit[ rpcBestCU->getWidth(0) ];
803        afCost[ iIdx ] += rpcBestCU->getTotalCost();
804        aiNum [ iIdx ] ++;
805      }
806    }
807
808    // Early CU determination
809    if( m_pcEncCfg->getUseEarlyCU() && rpcBestCU->isSkipped(0) )
810    {
811      bSubBranch = false;
812    }
813    else
814    {
815      bSubBranch = true;
816    }
817  }
818  else if(!(bSliceEnd && bInsidePicture))
819  {
820    bBoundary = true;
821#if RATE_CONTROL_LAMBDA_DOMAIN
822    m_addSADDepth++;
823#endif
824  }
825
826  // copy orginal YUV samples to PCM buffer
827  if( rpcBestCU->isLosslessCoded(0) && (rpcBestCU->getIPCMFlag(0) == false))
828  {
829    xFillPCMBuffer(rpcBestCU, m_ppcOrigYuv[uiDepth]);
830  }
831  if( (g_uiMaxCUWidth>>uiDepth) == rpcTempCU->getSlice()->getPPS()->getMinCuDQPSize() )
832  {
833    Int idQP = m_pcEncCfg->getMaxDeltaQP();
834    iMinQP = Clip3( -rpcTempCU->getSlice()->getSPS()->getQpBDOffsetY(), MAX_QP, iBaseQP-idQP );
835    iMaxQP = Clip3( -rpcTempCU->getSlice()->getSPS()->getQpBDOffsetY(), MAX_QP, iBaseQP+idQP );
836    if ( (rpcTempCU->getSlice()->getSPS()->getUseLossless()) && (lowestQP < iMinQP) && rpcTempCU->getSlice()->getPPS()->getUseDQP() )
837    {
838      isAddLowestQP = true;
839      iMinQP = iMinQP - 1;     
840    }
841  }
842  else if( (g_uiMaxCUWidth>>uiDepth) > rpcTempCU->getSlice()->getPPS()->getMinCuDQPSize() )
843  {
844    iMinQP = iBaseQP;
845    iMaxQP = iBaseQP;
846  }
847  else
848  {
849    Int iStartQP;
850    if( pcPic->getCU( rpcTempCU->getAddr() )->getSliceSegmentStartCU(rpcTempCU->getZorderIdxInCU()) == pcSlice->getSliceSegmentCurStartCUAddr())
851    {
852      iStartQP = rpcTempCU->getQP(0);
853    }
854    else
855    {
856      UInt uiCurSliceStartPartIdx = pcSlice->getSliceSegmentCurStartCUAddr() % pcPic->getNumPartInCU() - rpcTempCU->getZorderIdxInCU();
857      iStartQP = rpcTempCU->getQP(uiCurSliceStartPartIdx);
858    }
859    iMinQP = iStartQP;
860    iMaxQP = iStartQP;
861  }
862#if RATE_CONTROL_LAMBDA_DOMAIN
863  if ( m_pcEncCfg->getUseRateCtrl() )
864  {
865    iMinQP = m_pcRateCtrl->getRCQP();
866    iMaxQP = m_pcRateCtrl->getRCQP();
867  }
868#else
869  if(m_pcEncCfg->getUseRateCtrl())
870  {
871    Int qp = m_pcRateCtrl->getUnitQP();
872    iMinQP  = Clip3( MIN_QP, MAX_QP, qp);
873    iMaxQP  = Clip3( MIN_QP, MAX_QP, qp);
874  }
875#endif
876  for (Int iQP=iMinQP; iQP<=iMaxQP; iQP++)
877  {
878    if (isAddLowestQP && (iQP == iMinQP))
879    {
880      iQP = lowestQP;
881    }
882    rpcTempCU->initEstData( uiDepth, iQP );
883
884    // further split
885    if( bSubBranch && bTrySplitDQP && uiDepth < g_uiMaxCUDepth - g_uiAddCUDepth )
886    {
887#if H_3D_VSO // M9
888      // reset Model
889      if( m_pcRdCost->getUseRenModel() )
890      {
891        UInt  uiWidth     = m_ppcOrigYuv[uiDepth]->getWidth ( );
892        UInt  uiHeight    = m_ppcOrigYuv[uiDepth]->getHeight( );
893        Pel*  piSrc       = m_ppcOrigYuv[uiDepth]->getLumaAddr( 0 );
894        UInt  uiSrcStride = m_ppcOrigYuv[uiDepth]->getStride();
895        m_pcRdCost->setRenModelData( m_ppcBestCU[uiDepth], 0, piSrc, uiSrcStride, uiWidth, uiHeight );
896      }
897#endif
898
899      UChar       uhNextDepth         = uiDepth+1;
900      TComDataCU* pcSubBestPartCU     = m_ppcBestCU[uhNextDepth];
901      TComDataCU* pcSubTempPartCU     = m_ppcTempCU[uhNextDepth];
902
903      for ( UInt uiPartUnitIdx = 0; uiPartUnitIdx < 4; uiPartUnitIdx++ )
904      {
905        pcSubBestPartCU->initSubCU( rpcTempCU, uiPartUnitIdx, uhNextDepth, iQP );           // clear sub partition datas or init.
906        pcSubTempPartCU->initSubCU( rpcTempCU, uiPartUnitIdx, uhNextDepth, iQP );           // clear sub partition datas or init.
907
908        Bool bInSlice = pcSubBestPartCU->getSCUAddr()+pcSubBestPartCU->getTotalNumPart()>pcSlice->getSliceSegmentCurStartCUAddr()&&pcSubBestPartCU->getSCUAddr()<pcSlice->getSliceSegmentCurEndCUAddr();
909        if(bInSlice && ( pcSubBestPartCU->getCUPelX() < pcSlice->getSPS()->getPicWidthInLumaSamples() ) && ( pcSubBestPartCU->getCUPelY() < pcSlice->getSPS()->getPicHeightInLumaSamples() ) )
910        {
911          if( m_bUseSBACRD )
912          {
913            if ( 0 == uiPartUnitIdx) //initialize RD with previous depth buffer
914            {
915              m_pppcRDSbacCoder[uhNextDepth][CI_CURR_BEST]->load(m_pppcRDSbacCoder[uiDepth][CI_CURR_BEST]);
916            }
917            else
918            {
919              m_pppcRDSbacCoder[uhNextDepth][CI_CURR_BEST]->load(m_pppcRDSbacCoder[uhNextDepth][CI_NEXT_BEST]);
920            }
921          }
922
923#if AMP_ENC_SPEEDUP
924          if ( rpcBestCU->isIntra(0) )
925          {
926            xCompressCU( pcSubBestPartCU, pcSubTempPartCU, uhNextDepth, SIZE_NONE );
927          }
928          else
929          {
930            xCompressCU( pcSubBestPartCU, pcSubTempPartCU, uhNextDepth, rpcBestCU->getPartitionSize(0) );
931          }
932#else
933          xCompressCU( pcSubBestPartCU, pcSubTempPartCU, uhNextDepth );
934#endif
935
936          rpcTempCU->copyPartFrom( pcSubBestPartCU, uiPartUnitIdx, uhNextDepth );         // Keep best part data to current temporary data.
937          xCopyYuv2Tmp( pcSubBestPartCU->getTotalNumPart()*uiPartUnitIdx, uhNextDepth );
938        }
939        else if (bInSlice)
940        {
941          pcSubBestPartCU->copyToPic( uhNextDepth );
942          rpcTempCU->copyPartFrom( pcSubBestPartCU, uiPartUnitIdx, uhNextDepth );
943        }
944      }
945
946      if( !bBoundary )
947      {
948        m_pcEntropyCoder->resetBits();
949        m_pcEntropyCoder->encodeSplitFlag( rpcTempCU, 0, uiDepth, true );
950
951        rpcTempCU->getTotalBits() += m_pcEntropyCoder->getNumberOfWrittenBits(); // split bits
952        if(m_pcEncCfg->getUseSBACRD())
953        {
954          rpcTempCU->getTotalBins() += ((TEncBinCABAC *)((TEncSbac*)m_pcEntropyCoder->m_pcEntropyCoderIf)->getEncBinIf())->getBinsCoded();
955        }
956      }
957
958#if H_3D_VSO // M10
959      if ( m_pcRdCost->getUseVSO() )
960        rpcTempCU->getTotalCost()  = m_pcRdCost->calcRdCostVSO( rpcTempCU->getTotalBits(), rpcTempCU->getTotalDistortion() );
961      else
962#endif
963      rpcTempCU->getTotalCost()  = m_pcRdCost->calcRdCost( rpcTempCU->getTotalBits(), rpcTempCU->getTotalDistortion() );
964
965      if( (g_uiMaxCUWidth>>uiDepth) == rpcTempCU->getSlice()->getPPS()->getMinCuDQPSize() && rpcTempCU->getSlice()->getPPS()->getUseDQP())
966      {
967        Bool hasResidual = false;
968        for( UInt uiBlkIdx = 0; uiBlkIdx < rpcTempCU->getTotalNumPart(); uiBlkIdx ++)
969        {
970          if( ( pcPic->getCU( rpcTempCU->getAddr() )->getSliceSegmentStartCU(uiBlkIdx+rpcTempCU->getZorderIdxInCU()) == rpcTempCU->getSlice()->getSliceSegmentCurStartCUAddr() ) && 
971              ( rpcTempCU->getCbf( uiBlkIdx, TEXT_LUMA ) || rpcTempCU->getCbf( uiBlkIdx, TEXT_CHROMA_U ) || rpcTempCU->getCbf( uiBlkIdx, TEXT_CHROMA_V ) ) )
972          {
973            hasResidual = true;
974            break;
975          }
976        }
977
978        UInt uiTargetPartIdx;
979        if ( pcPic->getCU( rpcTempCU->getAddr() )->getSliceSegmentStartCU(rpcTempCU->getZorderIdxInCU()) != pcSlice->getSliceSegmentCurStartCUAddr() )
980        {
981          uiTargetPartIdx = pcSlice->getSliceSegmentCurStartCUAddr() % pcPic->getNumPartInCU() - rpcTempCU->getZorderIdxInCU();
982        }
983        else
984        {
985          uiTargetPartIdx = 0;
986        }
987        if ( hasResidual )
988        {
989#if !RDO_WITHOUT_DQP_BITS
990          m_pcEntropyCoder->resetBits();
991          m_pcEntropyCoder->encodeQP( rpcTempCU, uiTargetPartIdx, false );
992          rpcTempCU->getTotalBits() += m_pcEntropyCoder->getNumberOfWrittenBits(); // dQP bits
993          if(m_pcEncCfg->getUseSBACRD())
994          {
995            rpcTempCU->getTotalBins() += ((TEncBinCABAC *)((TEncSbac*)m_pcEntropyCoder->m_pcEntropyCoderIf)->getEncBinIf())->getBinsCoded();
996          }
997#if H_3D_VSO // M11
998          if ( m_pcRdCost->getUseLambdaScaleVSO())         
999            rpcTempCU->getTotalCost()  = m_pcRdCost->calcRdCostVSO( rpcTempCU->getTotalBits(), rpcTempCU->getTotalDistortion() );         
1000          else
1001#endif
1002          rpcTempCU->getTotalCost()  = m_pcRdCost->calcRdCost( rpcTempCU->getTotalBits(), rpcTempCU->getTotalDistortion() );
1003#endif
1004
1005          Bool foundNonZeroCbf = false;
1006          rpcTempCU->setQPSubCUs( rpcTempCU->getRefQP( uiTargetPartIdx ), rpcTempCU, 0, uiDepth, foundNonZeroCbf );
1007          assert( foundNonZeroCbf );
1008        }
1009        else
1010        {
1011          rpcTempCU->setQPSubParts( rpcTempCU->getRefQP( uiTargetPartIdx ), 0, uiDepth ); // set QP to default QP
1012        }
1013      }
1014
1015      if( m_bUseSBACRD )
1016      {
1017        m_pppcRDSbacCoder[uhNextDepth][CI_NEXT_BEST]->store(m_pppcRDSbacCoder[uiDepth][CI_TEMP_BEST]);
1018      }
1019      Bool isEndOfSlice        = rpcBestCU->getSlice()->getSliceMode()==FIXED_NUMBER_OF_BYTES
1020                                 && (rpcBestCU->getTotalBits()>rpcBestCU->getSlice()->getSliceArgument()<<3);
1021      Bool isEndOfSliceSegment = rpcBestCU->getSlice()->getSliceSegmentMode()==FIXED_NUMBER_OF_BYTES
1022                                 && (rpcBestCU->getTotalBits()>rpcBestCU->getSlice()->getSliceSegmentArgument()<<3);
1023      if(isEndOfSlice||isEndOfSliceSegment)
1024      {
1025        rpcBestCU->getTotalCost()=rpcTempCU->getTotalCost()+1;
1026      }
1027      xCheckBestMode( rpcBestCU, rpcTempCU, uiDepth);                                  // RD compare current larger prediction
1028    }                                                                                  // with sub partitioned prediction.
1029    if (isAddLowestQP && (iQP == lowestQP))
1030    {
1031      iQP = iMinQP;
1032    }
1033  }
1034
1035
1036#if H_3D_VSO // M12
1037  if( m_pcRdCost->getUseRenModel() )
1038  {
1039    UInt  uiWidth     = m_ppcRecoYuvBest[uiDepth]->getWidth   ( );
1040    UInt  uiHeight    = m_ppcRecoYuvBest[uiDepth]->getHeight  ( );
1041    Pel*  piSrc       = m_ppcRecoYuvBest[uiDepth]->getLumaAddr( 0 );
1042    UInt  uiSrcStride = m_ppcRecoYuvBest[uiDepth]->getStride  ( );
1043    m_pcRdCost->setRenModelData( rpcBestCU, 0, piSrc, uiSrcStride, uiWidth, uiHeight );
1044  }
1045#endif
1046
1047  rpcBestCU->copyToPic(uiDepth);                                                     // Copy Best data to Picture for next partition prediction.
1048
1049  xCopyYuv2Pic( rpcBestCU->getPic(), rpcBestCU->getAddr(), rpcBestCU->getZorderIdxInCU(), uiDepth, uiDepth, rpcBestCU, uiLPelX, uiTPelY );   // Copy Yuv data to picture Yuv
1050  if( bBoundary ||(bSliceEnd && bInsidePicture))
1051  {
1052    return;
1053  }
1054
1055  // Assert if Best prediction mode is NONE
1056  // Selected mode's RD-cost must be not MAX_DOUBLE.
1057  assert( rpcBestCU->getPartitionSize ( 0 ) != SIZE_NONE  );
1058  assert( rpcBestCU->getPredictionMode( 0 ) != MODE_NONE  );
1059  assert( rpcBestCU->getTotalCost     (   ) != MAX_DOUBLE );
1060}
1061
1062/** finish encoding a cu and handle end-of-slice conditions
1063 * \param pcCU
1064 * \param uiAbsPartIdx
1065 * \param uiDepth
1066 * \returns Void
1067 */
1068Void TEncCu::finishCU( TComDataCU* pcCU, UInt uiAbsPartIdx, UInt uiDepth )
1069{
1070  TComPic* pcPic = pcCU->getPic();
1071  TComSlice * pcSlice = pcCU->getPic()->getSlice(pcCU->getPic()->getCurrSliceIdx());
1072
1073  //Calculate end address
1074  UInt uiCUAddr = pcCU->getSCUAddr()+uiAbsPartIdx;
1075
1076  UInt uiInternalAddress = pcPic->getPicSym()->getPicSCUAddr(pcSlice->getSliceSegmentCurEndCUAddr()-1) % pcPic->getNumPartInCU();
1077  UInt uiExternalAddress = pcPic->getPicSym()->getPicSCUAddr(pcSlice->getSliceSegmentCurEndCUAddr()-1) / pcPic->getNumPartInCU();
1078  UInt uiPosX = ( uiExternalAddress % pcPic->getFrameWidthInCU() ) * g_uiMaxCUWidth+ g_auiRasterToPelX[ g_auiZscanToRaster[uiInternalAddress] ];
1079  UInt uiPosY = ( uiExternalAddress / pcPic->getFrameWidthInCU() ) * g_uiMaxCUHeight+ g_auiRasterToPelY[ g_auiZscanToRaster[uiInternalAddress] ];
1080  UInt uiWidth = pcSlice->getSPS()->getPicWidthInLumaSamples();
1081  UInt uiHeight = pcSlice->getSPS()->getPicHeightInLumaSamples();
1082  while(uiPosX>=uiWidth||uiPosY>=uiHeight)
1083  {
1084    uiInternalAddress--;
1085    uiPosX = ( uiExternalAddress % pcPic->getFrameWidthInCU() ) * g_uiMaxCUWidth+ g_auiRasterToPelX[ g_auiZscanToRaster[uiInternalAddress] ];
1086    uiPosY = ( uiExternalAddress / pcPic->getFrameWidthInCU() ) * g_uiMaxCUHeight+ g_auiRasterToPelY[ g_auiZscanToRaster[uiInternalAddress] ];
1087  }
1088  uiInternalAddress++;
1089  if(uiInternalAddress==pcCU->getPic()->getNumPartInCU())
1090  {
1091    uiInternalAddress = 0;
1092    uiExternalAddress = pcPic->getPicSym()->getCUOrderMap(pcPic->getPicSym()->getInverseCUOrderMap(uiExternalAddress)+1);
1093  }
1094  UInt uiRealEndAddress = pcPic->getPicSym()->getPicSCUEncOrder(uiExternalAddress*pcPic->getNumPartInCU()+uiInternalAddress);
1095
1096  // Encode slice finish
1097  Bool bTerminateSlice = false;
1098  if (uiCUAddr+(pcCU->getPic()->getNumPartInCU()>>(uiDepth<<1)) == uiRealEndAddress)
1099  {
1100    bTerminateSlice = true;
1101  }
1102  UInt uiGranularityWidth = g_uiMaxCUWidth;
1103  uiPosX = pcCU->getCUPelX() + g_auiRasterToPelX[ g_auiZscanToRaster[uiAbsPartIdx] ];
1104  uiPosY = pcCU->getCUPelY() + g_auiRasterToPelY[ g_auiZscanToRaster[uiAbsPartIdx] ];
1105  Bool granularityBoundary=((uiPosX+pcCU->getWidth(uiAbsPartIdx))%uiGranularityWidth==0||(uiPosX+pcCU->getWidth(uiAbsPartIdx)==uiWidth))
1106    &&((uiPosY+pcCU->getHeight(uiAbsPartIdx))%uiGranularityWidth==0||(uiPosY+pcCU->getHeight(uiAbsPartIdx)==uiHeight));
1107 
1108  if(granularityBoundary)
1109  {
1110    // The 1-terminating bit is added to all streams, so don't add it here when it's 1.
1111    if (!bTerminateSlice)
1112      m_pcEntropyCoder->encodeTerminatingBit( bTerminateSlice ? 1 : 0 );
1113  }
1114 
1115  Int numberOfWrittenBits = 0;
1116  if (m_pcBitCounter)
1117  {
1118    numberOfWrittenBits = m_pcEntropyCoder->getNumberOfWrittenBits();
1119  }
1120 
1121  // Calculate slice end IF this CU puts us over slice bit size.
1122  UInt iGranularitySize = pcCU->getPic()->getNumPartInCU();
1123  Int iGranularityEnd = ((pcCU->getSCUAddr()+uiAbsPartIdx)/iGranularitySize)*iGranularitySize;
1124  if(iGranularityEnd<=pcSlice->getSliceSegmentCurStartCUAddr()) 
1125  {
1126    iGranularityEnd+=max(iGranularitySize,(pcCU->getPic()->getNumPartInCU()>>(uiDepth<<1)));
1127  }
1128  // Set slice end parameter
1129  if(pcSlice->getSliceMode()==FIXED_NUMBER_OF_BYTES&&!pcSlice->getFinalized()&&pcSlice->getSliceBits()+numberOfWrittenBits>pcSlice->getSliceArgument()<<3) 
1130  {
1131    pcSlice->setSliceSegmentCurEndCUAddr(iGranularityEnd);
1132    pcSlice->setSliceCurEndCUAddr(iGranularityEnd);
1133    return;
1134  }
1135  // Set dependent slice end parameter
1136  if(pcSlice->getSliceSegmentMode()==FIXED_NUMBER_OF_BYTES&&!pcSlice->getFinalized()&&pcSlice->getSliceSegmentBits()+numberOfWrittenBits > pcSlice->getSliceSegmentArgument()<<3) 
1137  {
1138    pcSlice->setSliceSegmentCurEndCUAddr(iGranularityEnd);
1139    return;
1140  }
1141  if(granularityBoundary)
1142  {
1143    pcSlice->setSliceBits( (UInt)(pcSlice->getSliceBits() + numberOfWrittenBits) );
1144    pcSlice->setSliceSegmentBits(pcSlice->getSliceSegmentBits()+numberOfWrittenBits);
1145    if (m_pcBitCounter)
1146    {
1147      m_pcEntropyCoder->resetBits();     
1148    }
1149  }
1150}
1151
1152/** Compute QP for each CU
1153 * \param pcCU Target CU
1154 * \param uiDepth CU depth
1155 * \returns quantization parameter
1156 */
1157Int TEncCu::xComputeQP( TComDataCU* pcCU, UInt uiDepth )
1158{
1159  Int iBaseQp = pcCU->getSlice()->getSliceQp();
1160  Int iQpOffset = 0;
1161  if ( m_pcEncCfg->getUseAdaptiveQP() )
1162  {
1163    TEncPic* pcEPic = dynamic_cast<TEncPic*>( pcCU->getPic() );
1164    UInt uiAQDepth = min( uiDepth, pcEPic->getMaxAQDepth()-1 );
1165    TEncPicQPAdaptationLayer* pcAQLayer = pcEPic->getAQLayer( uiAQDepth );
1166    UInt uiAQUPosX = pcCU->getCUPelX() / pcAQLayer->getAQPartWidth();
1167    UInt uiAQUPosY = pcCU->getCUPelY() / pcAQLayer->getAQPartHeight();
1168    UInt uiAQUStride = pcAQLayer->getAQPartStride();
1169    TEncQPAdaptationUnit* acAQU = pcAQLayer->getQPAdaptationUnit();
1170
1171    Double dMaxQScale = pow(2.0, m_pcEncCfg->getQPAdaptationRange()/6.0);
1172    Double dAvgAct = pcAQLayer->getAvgActivity();
1173    Double dCUAct = acAQU[uiAQUPosY * uiAQUStride + uiAQUPosX].getActivity();
1174    Double dNormAct = (dMaxQScale*dCUAct + dAvgAct) / (dCUAct + dMaxQScale*dAvgAct);
1175    Double dQpOffset = log(dNormAct) / log(2.0) * 6.0;
1176    iQpOffset = Int(floor( dQpOffset + 0.49999 ));
1177  }
1178  return Clip3(-pcCU->getSlice()->getSPS()->getQpBDOffsetY(), MAX_QP, iBaseQp+iQpOffset );
1179}
1180
1181/** encode a CU block recursively
1182 * \param pcCU
1183 * \param uiAbsPartIdx
1184 * \param uiDepth
1185 * \returns Void
1186 */
1187Void TEncCu::xEncodeCU( TComDataCU* pcCU, UInt uiAbsPartIdx, UInt uiDepth )
1188{
1189  TComPic* pcPic = pcCU->getPic();
1190 
1191  Bool bBoundary = false;
1192  UInt uiLPelX   = pcCU->getCUPelX() + g_auiRasterToPelX[ g_auiZscanToRaster[uiAbsPartIdx] ];
1193  UInt uiRPelX   = uiLPelX + (g_uiMaxCUWidth>>uiDepth)  - 1;
1194  UInt uiTPelY   = pcCU->getCUPelY() + g_auiRasterToPelY[ g_auiZscanToRaster[uiAbsPartIdx] ];
1195  UInt uiBPelY   = uiTPelY + (g_uiMaxCUHeight>>uiDepth) - 1;
1196 
1197  TComSlice * pcSlice = pcCU->getPic()->getSlice(pcCU->getPic()->getCurrSliceIdx());
1198  // If slice start is within this cu...
1199  Bool bSliceStart = pcSlice->getSliceSegmentCurStartCUAddr() > pcPic->getPicSym()->getInverseCUOrderMap(pcCU->getAddr())*pcCU->getPic()->getNumPartInCU()+uiAbsPartIdx && 
1200    pcSlice->getSliceSegmentCurStartCUAddr() < pcPic->getPicSym()->getInverseCUOrderMap(pcCU->getAddr())*pcCU->getPic()->getNumPartInCU()+uiAbsPartIdx+( pcPic->getNumPartInCU() >> (uiDepth<<1) );
1201  // We need to split, so don't try these modes.
1202  if(!bSliceStart&&( uiRPelX < pcSlice->getSPS()->getPicWidthInLumaSamples() ) && ( uiBPelY < pcSlice->getSPS()->getPicHeightInLumaSamples() ) )
1203  {
1204    m_pcEntropyCoder->encodeSplitFlag( pcCU, uiAbsPartIdx, uiDepth );
1205  }
1206  else
1207  {
1208    bBoundary = true;
1209  }
1210 
1211  if( ( ( uiDepth < pcCU->getDepth( uiAbsPartIdx ) ) && ( uiDepth < (g_uiMaxCUDepth-g_uiAddCUDepth) ) ) || bBoundary )
1212  {
1213    UInt uiQNumParts = ( pcPic->getNumPartInCU() >> (uiDepth<<1) )>>2;
1214    if( (g_uiMaxCUWidth>>uiDepth) == pcCU->getSlice()->getPPS()->getMinCuDQPSize() && pcCU->getSlice()->getPPS()->getUseDQP())
1215    {
1216      setdQPFlag(true);
1217    }
1218    for ( UInt uiPartUnitIdx = 0; uiPartUnitIdx < 4; uiPartUnitIdx++, uiAbsPartIdx+=uiQNumParts )
1219    {
1220      uiLPelX   = pcCU->getCUPelX() + g_auiRasterToPelX[ g_auiZscanToRaster[uiAbsPartIdx] ];
1221      uiTPelY   = pcCU->getCUPelY() + g_auiRasterToPelY[ g_auiZscanToRaster[uiAbsPartIdx] ];
1222      Bool bInSlice = pcCU->getSCUAddr()+uiAbsPartIdx+uiQNumParts>pcSlice->getSliceSegmentCurStartCUAddr()&&pcCU->getSCUAddr()+uiAbsPartIdx<pcSlice->getSliceSegmentCurEndCUAddr();
1223      if(bInSlice&&( uiLPelX < pcSlice->getSPS()->getPicWidthInLumaSamples() ) && ( uiTPelY < pcSlice->getSPS()->getPicHeightInLumaSamples() ) )
1224      {
1225        xEncodeCU( pcCU, uiAbsPartIdx, uiDepth+1 );
1226      }
1227    }
1228    return;
1229  }
1230 
1231  if( (g_uiMaxCUWidth>>uiDepth) >= pcCU->getSlice()->getPPS()->getMinCuDQPSize() && pcCU->getSlice()->getPPS()->getUseDQP())
1232  {
1233    setdQPFlag(true);
1234  }
1235  if (pcCU->getSlice()->getPPS()->getTransquantBypassEnableFlag())
1236  {
1237    m_pcEntropyCoder->encodeCUTransquantBypassFlag( pcCU, uiAbsPartIdx );
1238  }
1239  if( !pcCU->getSlice()->isIntra() )
1240  {
1241    m_pcEntropyCoder->encodeSkipFlag( pcCU, uiAbsPartIdx );
1242  }
1243 
1244  if( pcCU->isSkipped( uiAbsPartIdx ) )
1245  {
1246    m_pcEntropyCoder->encodeMergeIndex( pcCU, uiAbsPartIdx );
1247#if H_3D_ARP
1248    m_pcEntropyCoder->encodeARPW( pcCU , uiAbsPartIdx );
1249#endif
1250    finishCU(pcCU,uiAbsPartIdx,uiDepth);
1251    return;
1252  }
1253  m_pcEntropyCoder->encodePredMode( pcCU, uiAbsPartIdx );
1254 
1255  m_pcEntropyCoder->encodePartSize( pcCU, uiAbsPartIdx, uiDepth );
1256 
1257  if (pcCU->isIntra( uiAbsPartIdx ) && pcCU->getPartitionSize( uiAbsPartIdx ) == SIZE_2Nx2N )
1258  {
1259    m_pcEntropyCoder->encodeIPCMInfo( pcCU, uiAbsPartIdx );
1260
1261    if(pcCU->getIPCMFlag(uiAbsPartIdx))
1262    {
1263      // Encode slice finish
1264      finishCU(pcCU,uiAbsPartIdx,uiDepth);
1265      return;
1266    }
1267  }
1268
1269  // prediction Info ( Intra : direction mode, Inter : Mv, reference idx )
1270  m_pcEntropyCoder->encodePredInfo( pcCU, uiAbsPartIdx );
1271#if H_3D_ARP
1272  m_pcEntropyCoder->encodeARPW( pcCU , uiAbsPartIdx );
1273#endif
1274
1275  // Encode Coefficients
1276  Bool bCodeDQP = getdQPFlag();
1277  m_pcEntropyCoder->encodeCoeff( pcCU, uiAbsPartIdx, uiDepth, pcCU->getWidth (uiAbsPartIdx), pcCU->getHeight(uiAbsPartIdx), bCodeDQP );
1278  setdQPFlag( bCodeDQP );
1279
1280  // --- write terminating bit ---
1281  finishCU(pcCU,uiAbsPartIdx,uiDepth);
1282}
1283
1284/** check RD costs for a CU block encoded with merge
1285 * \param rpcBestCU
1286 * \param rpcTempCU
1287 * \returns Void
1288 */
1289Void TEncCu::xCheckRDCostMerge2Nx2N( TComDataCU*& rpcBestCU, TComDataCU*& rpcTempCU, Bool *earlyDetectionSkipMode )
1290{
1291  assert( rpcTempCU->getSlice()->getSliceType() != I_SLICE );
1292  TComMvField  cMvFieldNeighbours[MRG_MAX_NUM_CANDS << 1]; // double length for mv of both lists
1293  UChar uhInterDirNeighbours[MRG_MAX_NUM_CANDS];
1294  Int numValidMergeCand = 0;
1295
1296  for( UInt ui = 0; ui < rpcTempCU->getSlice()->getMaxNumMergeCand(); ++ui )
1297  {
1298    uhInterDirNeighbours[ui] = 0;
1299  }
1300  UChar uhDepth = rpcTempCU->getDepth( 0 );
1301
1302#if H_3D_VSO // M1  //nececcary here?
1303  if( m_pcRdCost->getUseRenModel() )
1304  {
1305    UInt  uiWidth     = m_ppcOrigYuv[uhDepth]->getWidth ( );
1306    UInt  uiHeight    = m_ppcOrigYuv[uhDepth]->getHeight( );
1307    Pel*  piSrc       = m_ppcOrigYuv[uhDepth]->getLumaAddr( );
1308    UInt  uiSrcStride = m_ppcOrigYuv[uhDepth]->getStride();
1309    m_pcRdCost->setRenModelData( rpcTempCU, 0, piSrc, uiSrcStride, uiWidth, uiHeight );
1310  }
1311#endif
1312
1313  rpcTempCU->setPartSizeSubParts( SIZE_2Nx2N, 0, uhDepth ); // interprets depth relative to LCU level
1314  rpcTempCU->setCUTransquantBypassSubParts( m_pcEncCfg->getCUTransquantBypassFlagValue(), 0, uhDepth );
1315  rpcTempCU->getInterMergeCandidates( 0, 0, cMvFieldNeighbours,uhInterDirNeighbours, numValidMergeCand );
1316
1317  Int mergeCandBuffer[MRG_MAX_NUM_CANDS];
1318  for( UInt ui = 0; ui < rpcTempCU->getSlice()->getMaxNumMergeCand(); ++ui )
1319  {
1320    mergeCandBuffer[ui] = 0;
1321  }
1322
1323  Bool bestIsSkip = false;
1324
1325  UInt iteration;
1326  if ( rpcTempCU->isLosslessCoded(0))
1327  {
1328    iteration = 1;
1329  }
1330  else 
1331  {
1332    iteration = 2;
1333  }
1334
1335#if H_3D_ARP
1336  Int nARPWMax = rpcTempCU->getSlice()->getARPStepNum() - 1;
1337  if( nARPWMax < 0 || !rpcTempCU->getDvInfo(0).bDV )
1338  {
1339    nARPWMax = 0;
1340  }
1341  for( Int nARPW=nARPWMax; nARPW >= 0 ; nARPW-- )
1342  {
1343    memset( mergeCandBuffer, 0, MRG_MAX_NUM_CANDS*sizeof(Int) );
1344#endif
1345  for( UInt uiNoResidual = 0; uiNoResidual < iteration; ++uiNoResidual )
1346  {
1347    for( UInt uiMergeCand = 0; uiMergeCand < numValidMergeCand; ++uiMergeCand )
1348    {
1349      {
1350        if(!(uiNoResidual==1 && mergeCandBuffer[uiMergeCand]==1))
1351        {
1352
1353        if( !(bestIsSkip && uiNoResidual == 0) )
1354        {
1355          // set MC parameters
1356          rpcTempCU->setPredModeSubParts( MODE_INTER, 0, uhDepth ); // interprets depth relative to LCU level
1357          rpcTempCU->setCUTransquantBypassSubParts( m_pcEncCfg->getCUTransquantBypassFlagValue(),     0, uhDepth );
1358          rpcTempCU->setPartSizeSubParts( SIZE_2Nx2N, 0, uhDepth ); // interprets depth relative to LCU level
1359#if H_3D_ARP
1360          rpcTempCU->setARPWSubParts( (UChar)nARPW , 0 , uhDepth );
1361#endif
1362          rpcTempCU->setMergeFlagSubParts( true, 0, 0, uhDepth ); // interprets depth relative to LCU level
1363          rpcTempCU->setMergeIndexSubParts( uiMergeCand, 0, 0, uhDepth ); // interprets depth relative to LCU level
1364          rpcTempCU->setInterDirSubParts( uhInterDirNeighbours[uiMergeCand], 0, 0, uhDepth ); // interprets depth relative to LCU level
1365          rpcTempCU->getCUMvField( REF_PIC_LIST_0 )->setAllMvField( cMvFieldNeighbours[0 + 2*uiMergeCand], SIZE_2Nx2N, 0, 0 ); // interprets depth relative to rpcTempCU level
1366          rpcTempCU->getCUMvField( REF_PIC_LIST_1 )->setAllMvField( cMvFieldNeighbours[1 + 2*uiMergeCand], SIZE_2Nx2N, 0, 0 ); // interprets depth relative to rpcTempCU level
1367
1368#if H_3D_ARP
1369          if( nARPW )
1370          {
1371            Bool bSignalflag[2] = { true, true };
1372            for( UInt uiRefListIdx = 0; uiRefListIdx < 2; uiRefListIdx ++ )
1373            {
1374              Int iRefIdx = cMvFieldNeighbours[uiRefListIdx + 2*uiMergeCand].getRefIdx();
1375              RefPicList eRefList = uiRefListIdx ? REF_PIC_LIST_1 : REF_PIC_LIST_0;
1376              if( iRefIdx < 0 || rpcTempCU->getSlice()->getPOC() == rpcTempCU->getSlice()->getRefPOC(eRefList, iRefIdx) )
1377              {
1378                bSignalflag[uiRefListIdx] = false;
1379              }
1380            }
1381            if( !bSignalflag[0] && !bSignalflag[1] )
1382            {
1383              rpcTempCU->setARPWSubParts( 0 , 0 , uhDepth );
1384            }
1385          }
1386#endif
1387       // do MC
1388       m_pcPredSearch->motionCompensation ( rpcTempCU, m_ppcPredYuvTemp[uhDepth] );
1389       // estimate residual and encode everything
1390
1391#if H_3D_VSO //M2
1392       if( m_pcRdCost->getUseRenModel() )
1393       { //Reset
1394         UInt  uiWidth     = m_ppcOrigYuv[uhDepth]->getWidth    ();
1395         UInt  uiHeight    = m_ppcOrigYuv[uhDepth]->getHeight   ();
1396         Pel*  piSrc       = m_ppcOrigYuv[uhDepth]->getLumaAddr ();
1397         UInt  uiSrcStride = m_ppcOrigYuv[uhDepth]->getStride   ();
1398         m_pcRdCost->setRenModelData( rpcTempCU, 0, piSrc, uiSrcStride, uiWidth, uiHeight );
1399       }
1400#endif
1401
1402       m_pcPredSearch->encodeResAndCalcRdInterCU( rpcTempCU,
1403         m_ppcOrigYuv    [uhDepth],
1404         m_ppcPredYuvTemp[uhDepth],
1405         m_ppcResiYuvTemp[uhDepth],
1406         m_ppcResiYuvBest[uhDepth],
1407         m_ppcRecoYuvTemp[uhDepth],
1408         (uiNoResidual? true:false));
1409
1410
1411       if(uiNoResidual==0)
1412       {
1413         if(rpcTempCU->getQtRootCbf(0) == 0)
1414         {
1415           mergeCandBuffer[uiMergeCand] = 1;
1416         }
1417       }
1418
1419       rpcTempCU->setSkipFlagSubParts( rpcTempCU->getQtRootCbf(0) == 0, 0, uhDepth );
1420          Int orgQP = rpcTempCU->getQP( 0 );
1421          xCheckDQP( rpcTempCU );
1422          xCheckBestMode(rpcBestCU, rpcTempCU, uhDepth);
1423          rpcTempCU->initEstData( uhDepth, orgQP );
1424
1425
1426      if( m_pcEncCfg->getUseFastDecisionForMerge() && !bestIsSkip )
1427      {
1428        bestIsSkip = rpcBestCU->getQtRootCbf(0) == 0;
1429      }
1430
1431    }
1432    }
1433   }
1434  }
1435
1436  if(uiNoResidual == 0 && m_pcEncCfg->getUseEarlySkipDetection())
1437  {
1438    if(rpcBestCU->getQtRootCbf( 0 ) == 0)
1439    {
1440      if( rpcBestCU->getMergeFlag( 0 ))
1441      {
1442        *earlyDetectionSkipMode = true;
1443      }
1444      else
1445      {
1446        Int absoulte_MV=0;
1447        for ( UInt uiRefListIdx = 0; uiRefListIdx < 2; uiRefListIdx++ )
1448        {
1449          if ( rpcBestCU->getSlice()->getNumRefIdx( RefPicList( uiRefListIdx ) ) > 0 )
1450          {
1451            TComCUMvField* pcCUMvField = rpcBestCU->getCUMvField(RefPicList( uiRefListIdx ));
1452            Int iHor = pcCUMvField->getMvd( 0 ).getAbsHor();
1453            Int iVer = pcCUMvField->getMvd( 0 ).getAbsVer();
1454            absoulte_MV+=iHor+iVer;
1455          }
1456        }
1457
1458        if(absoulte_MV == 0)
1459        {
1460          *earlyDetectionSkipMode = true;
1461        }
1462      }
1463    }
1464  }
1465 }
1466#if H_3D_ARP
1467 }
1468#endif
1469}
1470
1471
1472#if AMP_MRG
1473Void TEncCu::xCheckRDCostInter( TComDataCU*& rpcBestCU, TComDataCU*& rpcTempCU, PartSize ePartSize, Bool bUseMRG)
1474#else
1475Void TEncCu::xCheckRDCostInter( TComDataCU*& rpcBestCU, TComDataCU*& rpcTempCU, PartSize ePartSize )
1476#endif
1477{
1478  UChar uhDepth = rpcTempCU->getDepth( 0 );
1479#if H_3D_ARP
1480  Int iLayerId    = rpcTempCU->getSlice()->getLayerId();
1481  Bool bFirstTime = true;
1482  Int nARPWMax    = rpcTempCU->getSlice()->getARPStepNum() - 1;
1483
1484  if( nARPWMax < 0 || ePartSize != SIZE_2Nx2N || !rpcTempCU->getDvInfo(0).bDV  )
1485  {
1486    nARPWMax = 0;
1487  }
1488
1489  for( Int nARPW = 0; nARPW <= nARPWMax; nARPW++ )
1490  {
1491    if( bFirstTime == false && rpcTempCU->getSlice()->getVPS()->getUseAdvRP( iLayerId ) )
1492    {
1493      rpcTempCU->initEstData( rpcTempCU->getDepth(0), rpcTempCU->getQP(0) );
1494    }
1495#endif
1496#if H_3D_VSO // M3
1497  if( m_pcRdCost->getUseRenModel() )
1498  {
1499    UInt  uiWidth     = m_ppcOrigYuv[uhDepth]->getWidth ( );
1500    UInt  uiHeight    = m_ppcOrigYuv[uhDepth]->getHeight( );
1501    Pel*  piSrc       = m_ppcOrigYuv[uhDepth]->getLumaAddr( );
1502    UInt  uiSrcStride = m_ppcOrigYuv[uhDepth]->getStride();
1503    m_pcRdCost->setRenModelData( rpcTempCU, 0, piSrc, uiSrcStride, uiWidth, uiHeight );
1504  }
1505#endif
1506
1507  rpcTempCU->setDepthSubParts( uhDepth, 0 );
1508 
1509  rpcTempCU->setSkipFlagSubParts( false, 0, uhDepth );
1510
1511  rpcTempCU->setPartSizeSubParts  ( ePartSize,  0, uhDepth );
1512  rpcTempCU->setPredModeSubParts  ( MODE_INTER, 0, uhDepth );
1513  rpcTempCU->setCUTransquantBypassSubParts  ( m_pcEncCfg->getCUTransquantBypassFlagValue(),      0, uhDepth );
1514 
1515#if H_3D_ARP
1516  rpcTempCU->setARPWSubParts( (UChar)nARPW , 0 , uhDepth );
1517#endif
1518
1519#if H_3D_ARP
1520  if( bFirstTime == false && nARPWMax )
1521  {
1522    rpcTempCU->copyPartFrom( m_ppcWeightedTempCU[uhDepth] , 0 , uhDepth );
1523    rpcTempCU->setARPWSubParts( (UChar)nARPW , 0 , uhDepth );
1524
1525    m_pcPredSearch->motionCompensation( rpcTempCU , m_ppcPredYuvTemp[uhDepth] );
1526
1527    if(rpcTempCU->getPartitionSize(0)==SIZE_2Nx2N)
1528    {
1529      Bool bSignalflag[2] = { true, true };
1530      for(UInt uiRefListIdx = 0; uiRefListIdx < 2; uiRefListIdx ++ )
1531      {
1532        RefPicList eRefList = uiRefListIdx ? REF_PIC_LIST_1 : REF_PIC_LIST_0;
1533        Int iRefIdx = rpcTempCU->getCUMvField(eRefList)->getRefIdx(0);
1534        if( iRefIdx < 0 || rpcTempCU->getSlice()->getPOC() == rpcTempCU->getSlice()->getRefPOC(eRefList, iRefIdx) )
1535        {
1536          bSignalflag[uiRefListIdx] = false;
1537        }
1538      }
1539      if( !bSignalflag[0] && !bSignalflag[1] )
1540      {
1541        rpcTempCU->setARPWSubParts( 0 , 0 , uhDepth );
1542      }
1543    }
1544  }
1545  else
1546  {
1547    bFirstTime = false;
1548#endif
1549#if AMP_MRG
1550  rpcTempCU->setMergeAMP (true);
1551  m_pcPredSearch->predInterSearch ( rpcTempCU, m_ppcOrigYuv[uhDepth], m_ppcPredYuvTemp[uhDepth], m_ppcResiYuvTemp[uhDepth], m_ppcRecoYuvTemp[uhDepth], false, bUseMRG );
1552#else 
1553  m_pcPredSearch->predInterSearch ( rpcTempCU, m_ppcOrigYuv[uhDepth], m_ppcPredYuvTemp[uhDepth], m_ppcResiYuvTemp[uhDepth], m_ppcRecoYuvTemp[uhDepth] );
1554#endif
1555#if H_3D_ARP
1556   if( nARPWMax )
1557   {
1558     m_ppcWeightedTempCU[uhDepth]->copyPartFrom( rpcTempCU , 0 , uhDepth );
1559
1560     Bool bSignalflag[2] = { true, true };
1561     for(UInt uiRefListIdx = 0; uiRefListIdx < 2; uiRefListIdx ++ )
1562     {
1563       RefPicList eRefList = uiRefListIdx ? REF_PIC_LIST_1 : REF_PIC_LIST_0;
1564       Int iRefIdx = rpcTempCU->getCUMvField(eRefList)->getRefIdx(0);
1565       if( iRefIdx < 0 || rpcTempCU->getSlice()->getPOC() == rpcTempCU->getSlice()->getRefPOC(eRefList, iRefIdx) )
1566       {
1567         bSignalflag[uiRefListIdx] = false;
1568       }
1569     }
1570     if( !bSignalflag[0] && !bSignalflag[1])
1571     {
1572       rpcTempCU->setARPWSubParts( 0 , 0 , uhDepth );
1573     }
1574   }
1575  }
1576#endif
1577
1578#if AMP_MRG
1579  if ( !rpcTempCU->getMergeAMP() )
1580  {
1581#if H_3D_ARP
1582    if( nARPWMax )
1583    {
1584      continue;
1585    }
1586    else
1587#endif
1588    return;
1589  }
1590#endif
1591
1592#if RATE_CONTROL_LAMBDA_DOMAIN
1593  if ( m_pcEncCfg->getUseRateCtrl() && m_pcEncCfg->getLCULevelRC() && ePartSize == SIZE_2Nx2N && uhDepth <= m_addSADDepth )
1594  {
1595    UInt SAD = m_pcRdCost->getSADPart( g_bitDepthY, m_ppcPredYuvTemp[uhDepth]->getLumaAddr(), m_ppcPredYuvTemp[uhDepth]->getStride(),
1596      m_ppcOrigYuv[uhDepth]->getLumaAddr(), m_ppcOrigYuv[uhDepth]->getStride(),
1597      rpcTempCU->getWidth(0), rpcTempCU->getHeight(0) );
1598    m_temporalSAD = (Int)SAD;
1599  }
1600#endif
1601
1602  m_pcPredSearch->encodeResAndCalcRdInterCU( rpcTempCU, m_ppcOrigYuv[uhDepth], m_ppcPredYuvTemp[uhDepth], m_ppcResiYuvTemp[uhDepth], m_ppcResiYuvBest[uhDepth], m_ppcRecoYuvTemp[uhDepth], false );
1603
1604
1605#if H_3D_VSO // M4
1606  if( m_pcRdCost->getUseLambdaScaleVSO() )
1607    rpcTempCU->getTotalCost()  = m_pcRdCost->calcRdCostVSO( rpcTempCU->getTotalBits(), rpcTempCU->getTotalDistortion() );
1608  else
1609#endif
1610  rpcTempCU->getTotalCost()  = m_pcRdCost->calcRdCost( rpcTempCU->getTotalBits(), rpcTempCU->getTotalDistortion() );
1611
1612  xCheckDQP( rpcTempCU );
1613  xCheckBestMode(rpcBestCU, rpcTempCU, uhDepth);
1614#if H_3D_ARP
1615  }
1616#endif
1617}
1618
1619Void TEncCu::xCheckRDCostIntra( TComDataCU*& rpcBestCU, TComDataCU*& rpcTempCU, PartSize eSize )
1620{
1621  UInt uiDepth = rpcTempCU->getDepth( 0 );
1622 
1623#if H_3D_VSO // M5
1624  if( m_pcRdCost->getUseRenModel() )
1625  {
1626    UInt  uiWidth     = m_ppcOrigYuv[uiDepth]->getWidth   ();
1627    UInt  uiHeight    = m_ppcOrigYuv[uiDepth]->getHeight  ();
1628    Pel*  piSrc       = m_ppcOrigYuv[uiDepth]->getLumaAddr();
1629    UInt  uiSrcStride = m_ppcOrigYuv[uiDepth]->getStride  ();
1630    m_pcRdCost->setRenModelData( rpcTempCU, 0, piSrc, uiSrcStride, uiWidth, uiHeight );
1631  }
1632#endif
1633
1634  rpcTempCU->setSkipFlagSubParts( false, 0, uiDepth );
1635
1636  rpcTempCU->setPartSizeSubParts( eSize, 0, uiDepth );
1637  rpcTempCU->setPredModeSubParts( MODE_INTRA, 0, uiDepth );
1638  rpcTempCU->setCUTransquantBypassSubParts( m_pcEncCfg->getCUTransquantBypassFlagValue(), 0, uiDepth );
1639 
1640  Bool bSeparateLumaChroma = true; // choose estimation mode
1641  UInt uiPreCalcDistC      = 0;
1642  if( !bSeparateLumaChroma )
1643  {
1644    m_pcPredSearch->preestChromaPredMode( rpcTempCU, m_ppcOrigYuv[uiDepth], m_ppcPredYuvTemp[uiDepth] );
1645  }
1646  m_pcPredSearch  ->estIntraPredQT      ( rpcTempCU, m_ppcOrigYuv[uiDepth], m_ppcPredYuvTemp[uiDepth], m_ppcResiYuvTemp[uiDepth], m_ppcRecoYuvTemp[uiDepth], uiPreCalcDistC, bSeparateLumaChroma );
1647
1648  m_ppcRecoYuvTemp[uiDepth]->copyToPicLuma(rpcTempCU->getPic()->getPicYuvRec(), rpcTempCU->getAddr(), rpcTempCU->getZorderIdxInCU() );
1649 
1650  m_pcPredSearch  ->estIntraPredChromaQT( rpcTempCU, m_ppcOrigYuv[uiDepth], m_ppcPredYuvTemp[uiDepth], m_ppcResiYuvTemp[uiDepth], m_ppcRecoYuvTemp[uiDepth], uiPreCalcDistC );
1651 
1652  m_pcEntropyCoder->resetBits();
1653  if ( rpcTempCU->getSlice()->getPPS()->getTransquantBypassEnableFlag())
1654  {
1655    m_pcEntropyCoder->encodeCUTransquantBypassFlag( rpcTempCU, 0,          true );
1656  }
1657  m_pcEntropyCoder->encodeSkipFlag ( rpcTempCU, 0,          true );
1658  m_pcEntropyCoder->encodePredMode( rpcTempCU, 0,          true );
1659  m_pcEntropyCoder->encodePartSize( rpcTempCU, 0, uiDepth, true );
1660  m_pcEntropyCoder->encodePredInfo( rpcTempCU, 0,          true );
1661  m_pcEntropyCoder->encodeIPCMInfo(rpcTempCU, 0, true );
1662
1663  // Encode Coefficients
1664  Bool bCodeDQP = getdQPFlag();
1665  m_pcEntropyCoder->encodeCoeff( rpcTempCU, 0, uiDepth, rpcTempCU->getWidth (0), rpcTempCU->getHeight(0), bCodeDQP );
1666  setdQPFlag( bCodeDQP );
1667 
1668  if( m_bUseSBACRD ) m_pcRDGoOnSbacCoder->store(m_pppcRDSbacCoder[uiDepth][CI_TEMP_BEST]);
1669 
1670  rpcTempCU->getTotalBits() = m_pcEntropyCoder->getNumberOfWrittenBits();
1671  if(m_pcEncCfg->getUseSBACRD())
1672  {
1673    rpcTempCU->getTotalBins() = ((TEncBinCABAC *)((TEncSbac*)m_pcEntropyCoder->m_pcEntropyCoderIf)->getEncBinIf())->getBinsCoded();
1674  }
1675
1676#if H_3D_VSO // M6
1677  if( m_pcRdCost->getUseLambdaScaleVSO()) 
1678    rpcTempCU->getTotalCost() = m_pcRdCost->calcRdCostVSO( rpcTempCU->getTotalBits(), rpcTempCU->getTotalDistortion() ); 
1679  else
1680#endif
1681  rpcTempCU->getTotalCost() = m_pcRdCost->calcRdCost( rpcTempCU->getTotalBits(), rpcTempCU->getTotalDistortion() );
1682 
1683  xCheckDQP( rpcTempCU );
1684  xCheckBestMode(rpcBestCU, rpcTempCU, uiDepth);
1685}
1686
1687/** Check R-D costs for a CU with PCM mode.
1688 * \param rpcBestCU pointer to best mode CU data structure
1689 * \param rpcTempCU pointer to testing mode CU data structure
1690 * \returns Void
1691 *
1692 * \note Current PCM implementation encodes sample values in a lossless way. The distortion of PCM mode CUs are zero. PCM mode is selected if the best mode yields bits greater than that of PCM mode.
1693 */
1694Void TEncCu::xCheckIntraPCM( TComDataCU*& rpcBestCU, TComDataCU*& rpcTempCU )
1695{
1696  UInt uiDepth = rpcTempCU->getDepth( 0 );
1697
1698  rpcTempCU->setSkipFlagSubParts( false, 0, uiDepth );
1699
1700  rpcTempCU->setIPCMFlag(0, true);
1701  rpcTempCU->setIPCMFlagSubParts (true, 0, rpcTempCU->getDepth(0));
1702  rpcTempCU->setPartSizeSubParts( SIZE_2Nx2N, 0, uiDepth );
1703  rpcTempCU->setPredModeSubParts( MODE_INTRA, 0, uiDepth );
1704  rpcTempCU->setTrIdxSubParts ( 0, 0, uiDepth );
1705  rpcTempCU->setCUTransquantBypassSubParts( m_pcEncCfg->getCUTransquantBypassFlagValue(), 0, uiDepth );
1706
1707  m_pcPredSearch->IPCMSearch( rpcTempCU, m_ppcOrigYuv[uiDepth], m_ppcPredYuvTemp[uiDepth], m_ppcResiYuvTemp[uiDepth], m_ppcRecoYuvTemp[uiDepth]);
1708
1709  if( m_bUseSBACRD ) m_pcRDGoOnSbacCoder->load(m_pppcRDSbacCoder[uiDepth][CI_CURR_BEST]);
1710
1711  m_pcEntropyCoder->resetBits();
1712  if ( rpcTempCU->getSlice()->getPPS()->getTransquantBypassEnableFlag())
1713  {
1714    m_pcEntropyCoder->encodeCUTransquantBypassFlag( rpcTempCU, 0,          true );
1715  }
1716  m_pcEntropyCoder->encodeSkipFlag ( rpcTempCU, 0,          true );
1717  m_pcEntropyCoder->encodePredMode ( rpcTempCU, 0,          true );
1718  m_pcEntropyCoder->encodePartSize ( rpcTempCU, 0, uiDepth, true );
1719  m_pcEntropyCoder->encodeIPCMInfo ( rpcTempCU, 0, true );
1720
1721  if( m_bUseSBACRD ) m_pcRDGoOnSbacCoder->store(m_pppcRDSbacCoder[uiDepth][CI_TEMP_BEST]);
1722
1723  rpcTempCU->getTotalBits() = m_pcEntropyCoder->getNumberOfWrittenBits();
1724  if(m_pcEncCfg->getUseSBACRD())
1725  {
1726    rpcTempCU->getTotalBins() = ((TEncBinCABAC *)((TEncSbac*)m_pcEntropyCoder->m_pcEntropyCoderIf)->getEncBinIf())->getBinsCoded();
1727  }
1728#if H_3D_VSO // M44
1729  if ( m_pcRdCost->getUseVSO() )
1730    rpcTempCU->getTotalCost() = m_pcRdCost->calcRdCostVSO( rpcTempCU->getTotalBits(), rpcTempCU->getTotalDistortion() );
1731  else
1732#endif
1733  rpcTempCU->getTotalCost() = m_pcRdCost->calcRdCost( rpcTempCU->getTotalBits(), rpcTempCU->getTotalDistortion() );
1734
1735  xCheckDQP( rpcTempCU );
1736  xCheckBestMode( rpcBestCU, rpcTempCU, uiDepth );
1737}
1738
1739/** check whether current try is the best with identifying the depth of current try
1740 * \param rpcBestCU
1741 * \param rpcTempCU
1742 * \returns Void
1743 */
1744Void TEncCu::xCheckBestMode( TComDataCU*& rpcBestCU, TComDataCU*& rpcTempCU, UInt uiDepth )
1745{
1746  if( rpcTempCU->getTotalCost() < rpcBestCU->getTotalCost() )
1747  {
1748    TComYuv* pcYuv;
1749    // Change Information data
1750    TComDataCU* pcCU = rpcBestCU;
1751    rpcBestCU = rpcTempCU;
1752    rpcTempCU = pcCU;
1753
1754    // Change Prediction data
1755    pcYuv = m_ppcPredYuvBest[uiDepth];
1756    m_ppcPredYuvBest[uiDepth] = m_ppcPredYuvTemp[uiDepth];
1757    m_ppcPredYuvTemp[uiDepth] = pcYuv;
1758
1759    // Change Reconstruction data
1760    pcYuv = m_ppcRecoYuvBest[uiDepth];
1761    m_ppcRecoYuvBest[uiDepth] = m_ppcRecoYuvTemp[uiDepth];
1762    m_ppcRecoYuvTemp[uiDepth] = pcYuv;
1763
1764    pcYuv = NULL;
1765    pcCU  = NULL;
1766
1767    if( m_bUseSBACRD )  // store temp best CI for next CU coding
1768      m_pppcRDSbacCoder[uiDepth][CI_TEMP_BEST]->store(m_pppcRDSbacCoder[uiDepth][CI_NEXT_BEST]);
1769  }
1770}
1771
1772Void TEncCu::xCheckDQP( TComDataCU* pcCU )
1773{
1774  UInt uiDepth = pcCU->getDepth( 0 );
1775
1776  if( pcCU->getSlice()->getPPS()->getUseDQP() && (g_uiMaxCUWidth>>uiDepth) >= pcCU->getSlice()->getPPS()->getMinCuDQPSize() )
1777  {
1778    if ( pcCU->getCbf( 0, TEXT_LUMA, 0 ) || pcCU->getCbf( 0, TEXT_CHROMA_U, 0 ) || pcCU->getCbf( 0, TEXT_CHROMA_V, 0 ) )
1779    {
1780#if !RDO_WITHOUT_DQP_BITS
1781      m_pcEntropyCoder->resetBits();
1782      m_pcEntropyCoder->encodeQP( pcCU, 0, false );
1783      pcCU->getTotalBits() += m_pcEntropyCoder->getNumberOfWrittenBits(); // dQP bits
1784      if(m_pcEncCfg->getUseSBACRD())
1785      {
1786        pcCU->getTotalBins() += ((TEncBinCABAC *)((TEncSbac*)m_pcEntropyCoder->m_pcEntropyCoderIf)->getEncBinIf())->getBinsCoded();
1787      }
1788#if H_3D_VSO // M45
1789      if ( m_pcRdCost->getUseVSO() )     
1790        pcCU->getTotalCost() = m_pcRdCost->calcRdCostVSO( pcCU->getTotalBits(), pcCU->getTotalDistortion() );     
1791      else
1792#endif
1793      pcCU->getTotalCost() = m_pcRdCost->calcRdCost( pcCU->getTotalBits(), pcCU->getTotalDistortion() );
1794#endif
1795    }
1796    else
1797    {
1798      pcCU->setQPSubParts( pcCU->getRefQP( 0 ), 0, uiDepth ); // set QP to default QP
1799    }
1800  }
1801}
1802
1803Void TEncCu::xCopyAMVPInfo (AMVPInfo* pSrc, AMVPInfo* pDst)
1804{
1805  pDst->iN = pSrc->iN;
1806  for (Int i = 0; i < pSrc->iN; i++)
1807  {
1808    pDst->m_acMvCand[i] = pSrc->m_acMvCand[i];
1809  }
1810}
1811Void TEncCu::xCopyYuv2Pic(TComPic* rpcPic, UInt uiCUAddr, UInt uiAbsPartIdx, UInt uiDepth, UInt uiSrcDepth, TComDataCU* pcCU, UInt uiLPelX, UInt uiTPelY )
1812{
1813  UInt uiRPelX   = uiLPelX + (g_uiMaxCUWidth>>uiDepth)  - 1;
1814  UInt uiBPelY   = uiTPelY + (g_uiMaxCUHeight>>uiDepth) - 1;
1815  TComSlice * pcSlice = pcCU->getPic()->getSlice(pcCU->getPic()->getCurrSliceIdx());
1816  Bool bSliceStart = pcSlice->getSliceSegmentCurStartCUAddr() > rpcPic->getPicSym()->getInverseCUOrderMap(pcCU->getAddr())*pcCU->getPic()->getNumPartInCU()+uiAbsPartIdx && 
1817    pcSlice->getSliceSegmentCurStartCUAddr() < rpcPic->getPicSym()->getInverseCUOrderMap(pcCU->getAddr())*pcCU->getPic()->getNumPartInCU()+uiAbsPartIdx+( pcCU->getPic()->getNumPartInCU() >> (uiDepth<<1) );
1818  Bool bSliceEnd   = pcSlice->getSliceSegmentCurEndCUAddr() > rpcPic->getPicSym()->getInverseCUOrderMap(pcCU->getAddr())*pcCU->getPic()->getNumPartInCU()+uiAbsPartIdx && 
1819    pcSlice->getSliceSegmentCurEndCUAddr() < rpcPic->getPicSym()->getInverseCUOrderMap(pcCU->getAddr())*pcCU->getPic()->getNumPartInCU()+uiAbsPartIdx+( pcCU->getPic()->getNumPartInCU() >> (uiDepth<<1) );
1820  if(!bSliceEnd && !bSliceStart && ( uiRPelX < pcSlice->getSPS()->getPicWidthInLumaSamples() ) && ( uiBPelY < pcSlice->getSPS()->getPicHeightInLumaSamples() ) )
1821  {
1822    UInt uiAbsPartIdxInRaster = g_auiZscanToRaster[uiAbsPartIdx];
1823    UInt uiSrcBlkWidth = rpcPic->getNumPartInWidth() >> (uiSrcDepth);
1824    UInt uiBlkWidth    = rpcPic->getNumPartInWidth() >> (uiDepth);
1825    UInt uiPartIdxX = ( ( uiAbsPartIdxInRaster % rpcPic->getNumPartInWidth() ) % uiSrcBlkWidth) / uiBlkWidth;
1826    UInt uiPartIdxY = ( ( uiAbsPartIdxInRaster / rpcPic->getNumPartInWidth() ) % uiSrcBlkWidth) / uiBlkWidth;
1827    UInt uiPartIdx = uiPartIdxY * ( uiSrcBlkWidth / uiBlkWidth ) + uiPartIdxX;
1828    m_ppcRecoYuvBest[uiSrcDepth]->copyToPicYuv( rpcPic->getPicYuvRec (), uiCUAddr, uiAbsPartIdx, uiDepth - uiSrcDepth, uiPartIdx);
1829  }
1830  else
1831  {
1832    UInt uiQNumParts = ( pcCU->getPic()->getNumPartInCU() >> (uiDepth<<1) )>>2;
1833
1834    for ( UInt uiPartUnitIdx = 0; uiPartUnitIdx < 4; uiPartUnitIdx++, uiAbsPartIdx+=uiQNumParts )
1835    {
1836      UInt uiSubCULPelX   = uiLPelX + ( g_uiMaxCUWidth >>(uiDepth+1) )*( uiPartUnitIdx &  1 );
1837      UInt uiSubCUTPelY   = uiTPelY + ( g_uiMaxCUHeight>>(uiDepth+1) )*( uiPartUnitIdx >> 1 );
1838
1839      Bool bInSlice = rpcPic->getPicSym()->getInverseCUOrderMap(pcCU->getAddr())*pcCU->getPic()->getNumPartInCU()+uiAbsPartIdx+uiQNumParts > pcSlice->getSliceSegmentCurStartCUAddr() && 
1840        rpcPic->getPicSym()->getInverseCUOrderMap(pcCU->getAddr())*pcCU->getPic()->getNumPartInCU()+uiAbsPartIdx < pcSlice->getSliceSegmentCurEndCUAddr();
1841      if(bInSlice&&( uiSubCULPelX < pcSlice->getSPS()->getPicWidthInLumaSamples() ) && ( uiSubCUTPelY < pcSlice->getSPS()->getPicHeightInLumaSamples() ) )
1842      {
1843        xCopyYuv2Pic( rpcPic, uiCUAddr, uiAbsPartIdx, uiDepth+1, uiSrcDepth, pcCU, uiSubCULPelX, uiSubCUTPelY );   // Copy Yuv data to picture Yuv
1844      }
1845    }
1846  }
1847}
1848
1849Void TEncCu::xCopyYuv2Tmp( UInt uiPartUnitIdx, UInt uiNextDepth )
1850{
1851  UInt uiCurrDepth = uiNextDepth - 1;
1852  m_ppcRecoYuvBest[uiNextDepth]->copyToPartYuv( m_ppcRecoYuvTemp[uiCurrDepth], uiPartUnitIdx );
1853}
1854
1855/** Function for filling the PCM buffer of a CU using its original sample array
1856 * \param pcCU pointer to current CU
1857 * \param pcOrgYuv pointer to original sample array
1858 * \returns Void
1859 */
1860Void TEncCu::xFillPCMBuffer     ( TComDataCU*& pCU, TComYuv* pOrgYuv )
1861{
1862
1863  UInt   width        = pCU->getWidth(0);
1864  UInt   height       = pCU->getHeight(0);
1865
1866  Pel*   pSrcY = pOrgYuv->getLumaAddr(0, width); 
1867  Pel*   pDstY = pCU->getPCMSampleY();
1868  UInt   srcStride = pOrgYuv->getStride();
1869
1870  for(Int y = 0; y < height; y++ )
1871  {
1872    for(Int x = 0; x < width; x++ )
1873    {
1874      pDstY[x] = pSrcY[x];
1875    }
1876    pDstY += width;
1877    pSrcY += srcStride;
1878  }
1879
1880  Pel* pSrcCb       = pOrgYuv->getCbAddr();
1881  Pel* pSrcCr       = pOrgYuv->getCrAddr();;
1882
1883  Pel* pDstCb       = pCU->getPCMSampleCb();
1884  Pel* pDstCr       = pCU->getPCMSampleCr();;
1885
1886  UInt srcStrideC = pOrgYuv->getCStride();
1887  UInt heightC   = height >> 1;
1888  UInt widthC    = width  >> 1;
1889
1890  for(Int y = 0; y < heightC; y++ )
1891  {
1892    for(Int x = 0; x < widthC; x++ )
1893    {
1894      pDstCb[x] = pSrcCb[x];
1895      pDstCr[x] = pSrcCr[x];
1896    }
1897    pDstCb += widthC;
1898    pDstCr += widthC;
1899    pSrcCb += srcStrideC;
1900    pSrcCr += srcStrideC;
1901  }
1902}
1903
1904#if ADAPTIVE_QP_SELECTION
1905/** Collect ARL statistics from one block
1906  */
1907Int TEncCu::xTuCollectARLStats(TCoeff* rpcCoeff, Int* rpcArlCoeff, Int NumCoeffInCU, Double* cSum, UInt* numSamples )
1908{
1909  for( Int n = 0; n < NumCoeffInCU; n++ )
1910  {
1911    Int u = abs( rpcCoeff[ n ] );
1912    Int absc = rpcArlCoeff[ n ];
1913
1914    if( u != 0 )
1915    {
1916      if( u < LEVEL_RANGE )
1917      {
1918        cSum[ u ] += ( Double )absc;
1919        numSamples[ u ]++;
1920      }
1921      else 
1922      {
1923        cSum[ LEVEL_RANGE ] += ( Double )absc - ( Double )( u << ARL_C_PRECISION );
1924        numSamples[ LEVEL_RANGE ]++;
1925      }
1926    }
1927  }
1928
1929  return 0;
1930}
1931
1932/** Collect ARL statistics from one LCU
1933 * \param pcCU
1934 */
1935Void TEncCu::xLcuCollectARLStats(TComDataCU* rpcCU )
1936{
1937  Double cSum[ LEVEL_RANGE + 1 ];     //: the sum of DCT coefficients corresponding to datatype and quantization output
1938  UInt numSamples[ LEVEL_RANGE + 1 ]; //: the number of coefficients corresponding to datatype and quantization output
1939
1940  TCoeff* pCoeffY = rpcCU->getCoeffY();
1941  Int* pArlCoeffY = rpcCU->getArlCoeffY();
1942
1943  UInt uiMinCUWidth = g_uiMaxCUWidth >> g_uiMaxCUDepth;
1944  UInt uiMinNumCoeffInCU = 1 << uiMinCUWidth;
1945
1946  memset( cSum, 0, sizeof( Double )*(LEVEL_RANGE+1) );
1947  memset( numSamples, 0, sizeof( UInt )*(LEVEL_RANGE+1) );
1948
1949  // Collect stats to cSum[][] and numSamples[][]
1950  for(Int i = 0; i < rpcCU->getTotalNumPart(); i ++ )
1951  {
1952    UInt uiTrIdx = rpcCU->getTransformIdx(i);
1953
1954    if(rpcCU->getPredictionMode(i) == MODE_INTER)
1955    if( rpcCU->getCbf( i, TEXT_LUMA, uiTrIdx ) )
1956    {
1957      xTuCollectARLStats(pCoeffY, pArlCoeffY, uiMinNumCoeffInCU, cSum, numSamples);
1958    }//Note that only InterY is processed. QP rounding is based on InterY data only.
1959   
1960    pCoeffY  += uiMinNumCoeffInCU;
1961    pArlCoeffY  += uiMinNumCoeffInCU;
1962  }
1963
1964  for(Int u=1; u<LEVEL_RANGE;u++)
1965  {
1966    m_pcTrQuant->getSliceSumC()[u] += cSum[ u ] ;
1967    m_pcTrQuant->getSliceNSamples()[u] += numSamples[ u ] ;
1968  }
1969  m_pcTrQuant->getSliceSumC()[LEVEL_RANGE] += cSum[ LEVEL_RANGE ] ;
1970  m_pcTrQuant->getSliceNSamples()[LEVEL_RANGE] += numSamples[ LEVEL_RANGE ] ;
1971}
1972#endif
1973//! \}
Note: See TracBrowser for help on using the repository browser.