source: SHVCSoftware/branches/SHM-2.1-dev/source/Lib/TLibEncoder/TEncGOP.cpp @ 291

Last change on this file since 291 was 291, checked in by nokia, 11 years ago

JCTVC-M0457: Co-located picture signaling

File size: 117.6 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     TEncGOP.cpp
35    \brief    GOP encoder class
36*/
37
38#include <list>
39#include <algorithm>
40#include <functional>
41
42#include "TEncTop.h"
43#include "TEncGOP.h"
44#include "TEncAnalyze.h"
45#include "libmd5/MD5.h"
46#include "TLibCommon/SEI.h"
47#include "TLibCommon/NAL.h"
48#include "NALwrite.h"
49#include <time.h>
50#include <math.h>
51
52using namespace std;
53//! \ingroup TLibEncoder
54//! \{
55
56// ====================================================================================================================
57// Constructor / destructor / initialization / destroy
58// ====================================================================================================================
59Int getLSB(Int poc, Int maxLSB)
60{
61  if (poc >= 0)
62  {
63    return poc % maxLSB;
64  }
65  else
66  {
67    return (maxLSB - ((-poc) % maxLSB)) % maxLSB;
68  }
69}
70
71TEncGOP::TEncGOP()
72{
73  m_iLastIDR            = 0;
74  m_iGopSize            = 0;
75  m_iNumPicCoded        = 0; //Niko
76  m_bFirst              = true;
77 
78  m_pcCfg               = NULL;
79  m_pcSliceEncoder      = NULL;
80  m_pcListPic           = NULL;
81 
82  m_pcEntropyCoder      = NULL;
83  m_pcCavlcCoder        = NULL;
84  m_pcSbacCoder         = NULL;
85  m_pcBinCABAC          = NULL;
86 
87  m_bSeqFirst           = true;
88 
89  m_bRefreshPending     = 0;
90  m_pocCRA            = 0;
91  m_numLongTermRefPicSPS = 0;
92  ::memset(m_ltRefPicPocLsbSps, 0, sizeof(m_ltRefPicPocLsbSps));
93  ::memset(m_ltRefPicUsedByCurrPicFlag, 0, sizeof(m_ltRefPicUsedByCurrPicFlag));
94  m_cpbRemovalDelay   = 0;
95  m_lastBPSEI         = 0;
96#if L0045_NON_NESTED_SEI_RESTRICTIONS
97  xResetNonNestedSEIPresentFlags();
98#if K0180_SCALABLE_NESTING_SEI
99  xResetNestedSEIPresentFlags();
100#endif
101#endif
102#if SVC_UPSAMPLING
103  m_pcPredSearch        = NULL;
104#endif
105  return;
106}
107
108TEncGOP::~TEncGOP()
109{
110}
111
112/** Create list to contain pointers to LCU start addresses of slice.
113 */
114#if SVC_EXTENSION
115Void  TEncGOP::create( UInt layerId )
116{
117  m_bLongtermTestPictureHasBeenCoded = 0;
118  m_bLongtermTestPictureHasBeenCoded2 = 0;
119  m_layerId = layerId;
120}
121#else
122Void  TEncGOP::create()
123{
124  m_bLongtermTestPictureHasBeenCoded = 0;
125  m_bLongtermTestPictureHasBeenCoded2 = 0;
126}
127#endif
128
129Void  TEncGOP::destroy()
130{
131}
132
133Void TEncGOP::init ( TEncTop* pcTEncTop )
134{
135  m_pcEncTop     = pcTEncTop;
136  m_pcCfg                = pcTEncTop;
137  m_pcSliceEncoder       = pcTEncTop->getSliceEncoder();
138  m_pcListPic            = pcTEncTop->getListPic(); 
139 
140  m_pcEntropyCoder       = pcTEncTop->getEntropyCoder();
141  m_pcCavlcCoder         = pcTEncTop->getCavlcCoder();
142  m_pcSbacCoder          = pcTEncTop->getSbacCoder();
143  m_pcBinCABAC           = pcTEncTop->getBinCABAC();
144  m_pcLoopFilter         = pcTEncTop->getLoopFilter();
145  m_pcBitCounter         = pcTEncTop->getBitCounter();
146 
147  //--Adaptive Loop filter
148  m_pcSAO                = pcTEncTop->getSAO();
149  m_pcRateCtrl           = pcTEncTop->getRateCtrl();
150  m_lastBPSEI          = 0;
151  m_totalCoded         = 0;
152
153#if SVC_EXTENSION
154  m_ppcTEncTop           = pcTEncTop->getLayerEnc();
155#endif
156#if SVC_UPSAMPLING
157  m_pcPredSearch         = pcTEncTop->getPredSearch();                       ///< encoder search class
158#endif
159}
160
161SEIActiveParameterSets* TEncGOP::xCreateSEIActiveParameterSets (TComSPS *sps)
162{
163  SEIActiveParameterSets *seiActiveParameterSets = new SEIActiveParameterSets(); 
164  seiActiveParameterSets->activeVPSId = m_pcCfg->getVPS()->getVPSId(); 
165#if L0047_APS_FLAGS
166  seiActiveParameterSets->m_fullRandomAccessFlag = false;
167  seiActiveParameterSets->m_noParamSetUpdateFlag = false;
168#endif
169  seiActiveParameterSets->numSpsIdsMinus1 = 0;
170  seiActiveParameterSets->activeSeqParamSetId.resize(seiActiveParameterSets->numSpsIdsMinus1 + 1); 
171  seiActiveParameterSets->activeSeqParamSetId[0] = sps->getSPSId();
172  return seiActiveParameterSets;
173}
174
175SEIFramePacking* TEncGOP::xCreateSEIFramePacking()
176{
177  SEIFramePacking *seiFramePacking = new SEIFramePacking();
178  seiFramePacking->m_arrangementId = m_pcCfg->getFramePackingArrangementSEIId();
179  seiFramePacking->m_arrangementCancelFlag = 0;
180  seiFramePacking->m_arrangementType = m_pcCfg->getFramePackingArrangementSEIType();
181#if L0444_FPA_TYPE
182  assert((seiFramePacking->m_arrangementType > 2) && (seiFramePacking->m_arrangementType < 6) );
183#endif
184  seiFramePacking->m_quincunxSamplingFlag = m_pcCfg->getFramePackingArrangementSEIQuincunx();
185  seiFramePacking->m_contentInterpretationType = m_pcCfg->getFramePackingArrangementSEIInterpretation();
186  seiFramePacking->m_spatialFlippingFlag = 0;
187  seiFramePacking->m_frame0FlippedFlag = 0;
188  seiFramePacking->m_fieldViewsFlag = (seiFramePacking->m_arrangementType == 2);
189  seiFramePacking->m_currentFrameIsFrame0Flag = ((seiFramePacking->m_arrangementType == 5) && m_iNumPicCoded&1);
190  seiFramePacking->m_frame0SelfContainedFlag = 0;
191  seiFramePacking->m_frame1SelfContainedFlag = 0;
192  seiFramePacking->m_frame0GridPositionX = 0;
193  seiFramePacking->m_frame0GridPositionY = 0;
194  seiFramePacking->m_frame1GridPositionX = 0;
195  seiFramePacking->m_frame1GridPositionY = 0;
196  seiFramePacking->m_arrangementReservedByte = 0;
197#if L0045_PERSISTENCE_FLAGS
198  seiFramePacking->m_arrangementPersistenceFlag = true;
199#else
200  seiFramePacking->m_arrangementRepetetionPeriod = 1;
201#endif
202  seiFramePacking->m_upsampledAspectRatio = 0;
203  return seiFramePacking;
204}
205
206SEIDisplayOrientation* TEncGOP::xCreateSEIDisplayOrientation()
207{
208  SEIDisplayOrientation *seiDisplayOrientation = new SEIDisplayOrientation();
209  seiDisplayOrientation->cancelFlag = false;
210  seiDisplayOrientation->horFlip = false;
211  seiDisplayOrientation->verFlip = false;
212  seiDisplayOrientation->anticlockwiseRotation = m_pcCfg->getDisplayOrientationSEIAngle();
213  return seiDisplayOrientation;
214}
215
216#if J0149_TONE_MAPPING_SEI
217SEIToneMappingInfo*  TEncGOP::xCreateSEIToneMappingInfo()
218{
219  SEIToneMappingInfo *seiToneMappingInfo = new SEIToneMappingInfo();
220  seiToneMappingInfo->m_toneMapId = m_pcCfg->getTMISEIToneMapId();
221  seiToneMappingInfo->m_toneMapCancelFlag = m_pcCfg->getTMISEIToneMapCancelFlag();
222  seiToneMappingInfo->m_toneMapPersistenceFlag = m_pcCfg->getTMISEIToneMapPersistenceFlag();
223
224  seiToneMappingInfo->m_codedDataBitDepth = m_pcCfg->getTMISEICodedDataBitDepth();
225  assert(seiToneMappingInfo->m_codedDataBitDepth >= 8 && seiToneMappingInfo->m_codedDataBitDepth <= 14);
226  seiToneMappingInfo->m_targetBitDepth = m_pcCfg->getTMISEITargetBitDepth();
227  assert( (seiToneMappingInfo->m_targetBitDepth >= 1 && seiToneMappingInfo->m_targetBitDepth <= 17) || (seiToneMappingInfo->m_targetBitDepth  == 255) );
228  seiToneMappingInfo->m_modelId = m_pcCfg->getTMISEIModelID();
229  assert(seiToneMappingInfo->m_modelId >=0 &&seiToneMappingInfo->m_modelId<=4);
230
231  switch( seiToneMappingInfo->m_modelId)
232  {
233  case 0:
234    {
235      seiToneMappingInfo->m_minValue = m_pcCfg->getTMISEIMinValue();
236      seiToneMappingInfo->m_maxValue = m_pcCfg->getTMISEIMaxValue();
237      break;
238    }
239  case 1:
240    {
241      seiToneMappingInfo->m_sigmoidMidpoint = m_pcCfg->getTMISEISigmoidMidpoint();
242      seiToneMappingInfo->m_sigmoidWidth = m_pcCfg->getTMISEISigmoidWidth();
243      break;
244    }
245  case 2:
246    {
247      UInt num = 1u<<(seiToneMappingInfo->m_targetBitDepth);
248      seiToneMappingInfo->m_startOfCodedInterval.resize(num);
249      Int* ptmp = m_pcCfg->getTMISEIStartOfCodedInterva();
250      if(ptmp)
251      {
252        for(int i=0; i<num;i++)
253        {
254          seiToneMappingInfo->m_startOfCodedInterval[i] = ptmp[i];
255        }
256      }
257      break;
258    }
259  case 3:
260    {
261      seiToneMappingInfo->m_numPivots = m_pcCfg->getTMISEINumPivots();
262      seiToneMappingInfo->m_codedPivotValue.resize(seiToneMappingInfo->m_numPivots);
263      seiToneMappingInfo->m_targetPivotValue.resize(seiToneMappingInfo->m_numPivots);
264      Int* ptmpcoded = m_pcCfg->getTMISEICodedPivotValue();
265      Int* ptmptarget = m_pcCfg->getTMISEITargetPivotValue();
266      if(ptmpcoded&&ptmptarget)
267      {
268        for(int i=0; i<(seiToneMappingInfo->m_numPivots);i++)
269        {
270          seiToneMappingInfo->m_codedPivotValue[i]=ptmpcoded[i];
271          seiToneMappingInfo->m_targetPivotValue[i]=ptmptarget[i];
272         }
273       }
274       break;
275     }
276  case 4:
277     {
278       seiToneMappingInfo->m_cameraIsoSpeedIdc = m_pcCfg->getTMISEICameraIsoSpeedIdc();
279       seiToneMappingInfo->m_cameraIsoSpeedValue = m_pcCfg->getTMISEICameraIsoSpeedValue();
280       assert( seiToneMappingInfo->m_cameraIsoSpeedValue !=0 );
281       seiToneMappingInfo->m_exposureCompensationValueSignFlag = m_pcCfg->getTMISEIExposureCompensationValueSignFlag();
282       seiToneMappingInfo->m_exposureCompensationValueNumerator = m_pcCfg->getTMISEIExposureCompensationValueNumerator();
283       seiToneMappingInfo->m_exposureCompensationValueDenomIdc = m_pcCfg->getTMISEIExposureCompensationValueDenomIdc();
284       seiToneMappingInfo->m_refScreenLuminanceWhite = m_pcCfg->getTMISEIRefScreenLuminanceWhite();
285       seiToneMappingInfo->m_extendedRangeWhiteLevel = m_pcCfg->getTMISEIExtendedRangeWhiteLevel();
286       assert( seiToneMappingInfo->m_extendedRangeWhiteLevel >= 100 );
287       seiToneMappingInfo->m_nominalBlackLevelLumaCodeValue = m_pcCfg->getTMISEINominalBlackLevelLumaCodeValue();
288       seiToneMappingInfo->m_nominalWhiteLevelLumaCodeValue = m_pcCfg->getTMISEINominalWhiteLevelLumaCodeValue();
289       assert( seiToneMappingInfo->m_nominalWhiteLevelLumaCodeValue > seiToneMappingInfo->m_nominalBlackLevelLumaCodeValue );
290       seiToneMappingInfo->m_extendedWhiteLevelLumaCodeValue = m_pcCfg->getTMISEIExtendedWhiteLevelLumaCodeValue();
291       assert( seiToneMappingInfo->m_extendedWhiteLevelLumaCodeValue >= seiToneMappingInfo->m_nominalWhiteLevelLumaCodeValue );
292       break;
293    }
294  default:
295    {
296      assert(!"Undefined SEIToneMapModelId");
297      break;
298    }
299  }
300  return seiToneMappingInfo;
301}
302#endif
303Void TEncGOP::xCreateLeadingSEIMessages (/*SEIMessages seiMessages,*/ AccessUnit &accessUnit, TComSPS *sps)
304{
305  OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI);
306
307  if(m_pcCfg->getActiveParameterSetsSEIEnabled())
308  {
309    SEIActiveParameterSets *sei = xCreateSEIActiveParameterSets (sps);
310
311    //nalu = NALUnit(NAL_UNIT_SEI);
312    m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
313    m_seiWriter.writeSEImessage(nalu.m_Bitstream, *sei, sps); 
314    writeRBSPTrailingBits(nalu.m_Bitstream);
315    accessUnit.push_back(new NALUnitEBSP(nalu));
316    delete sei;
317#if L0045_NON_NESTED_SEI_RESTRICTIONS
318    m_activeParameterSetSEIPresentInAU = true;
319#endif
320  }
321
322  if(m_pcCfg->getFramePackingArrangementSEIEnabled())
323  {
324    SEIFramePacking *sei = xCreateSEIFramePacking ();
325
326    nalu = NALUnit(NAL_UNIT_PREFIX_SEI);
327    m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
328    m_seiWriter.writeSEImessage(nalu.m_Bitstream, *sei, sps);
329    writeRBSPTrailingBits(nalu.m_Bitstream);
330    accessUnit.push_back(new NALUnitEBSP(nalu));
331    delete sei;
332  }
333  if (m_pcCfg->getDisplayOrientationSEIAngle())
334  {
335    SEIDisplayOrientation *sei = xCreateSEIDisplayOrientation();
336
337    nalu = NALUnit(NAL_UNIT_PREFIX_SEI); 
338    m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
339    m_seiWriter.writeSEImessage(nalu.m_Bitstream, *sei, sps); 
340    writeRBSPTrailingBits(nalu.m_Bitstream);
341    accessUnit.push_back(new NALUnitEBSP(nalu));
342    delete sei;
343  }
344#if J0149_TONE_MAPPING_SEI
345  if(m_pcCfg->getToneMappingInfoSEIEnabled())
346  {
347    SEIToneMappingInfo *sei = xCreateSEIToneMappingInfo ();
348     
349    nalu = NALUnit(NAL_UNIT_PREFIX_SEI); 
350    m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
351    m_seiWriter.writeSEImessage(nalu.m_Bitstream, *sei, sps); 
352    writeRBSPTrailingBits(nalu.m_Bitstream);
353    accessUnit.push_back(new NALUnitEBSP(nalu));
354    delete sei;
355  }
356#endif
357}
358
359// ====================================================================================================================
360// Public member functions
361// ====================================================================================================================
362#if SVC_EXTENSION
363Void TEncGOP::compressGOP( Int iPicIdInGOP, Int iPOCLast, Int iNumPicRcvd, TComList<TComPic*>& rcListPic, TComList<TComPicYuv*>& rcListPicYuvRecOut, std::list<AccessUnit>& accessUnitsInGOP)
364#else
365Void TEncGOP::compressGOP( Int iPOCLast, Int iNumPicRcvd, TComList<TComPic*>& rcListPic, TComList<TComPicYuv*>& rcListPicYuvRecOut, std::list<AccessUnit>& accessUnitsInGOP)
366#endif
367{
368  TComPic*        pcPic;
369  TComPicYuv*     pcPicYuvRecOut;
370  TComSlice*      pcSlice;
371  TComOutputBitstream  *pcBitstreamRedirect;
372  pcBitstreamRedirect = new TComOutputBitstream;
373  AccessUnit::iterator  itLocationToPushSliceHeaderNALU; // used to store location where NALU containing slice header is to be inserted
374  UInt                  uiOneBitstreamPerSliceLength = 0;
375  TEncSbac* pcSbacCoders = NULL;
376  TComOutputBitstream* pcSubstreamsOut = NULL;
377
378  xInitGOP( iPOCLast, iNumPicRcvd, rcListPic, rcListPicYuvRecOut );
379
380  m_iNumPicCoded = 0;
381  SEIPictureTiming pictureTimingSEI;
382#if L0208_SOP_DESCRIPTION_SEI
383  Bool writeSOP = m_pcCfg->getSOPDescriptionSEIEnabled();
384#endif
385#if K0180_SCALABLE_NESTING_SEI
386  // Initialize Scalable Nesting SEI with single layer values
387  SEIScalableNesting scalableNestingSEI;
388  scalableNestingSEI.m_bitStreamSubsetFlag           = 1;      // If the nested SEI messages are picture buffereing SEI mesages, picure timing SEI messages or sub-picture timing SEI messages, bitstream_subset_flag shall be equal to 1
389  scalableNestingSEI.m_nestingOpFlag                 = 0;
390  scalableNestingSEI.m_nestingNumOpsMinus1           = 0;      //nesting_num_ops_minus1
391  scalableNestingSEI.m_allLayersFlag                 = 0;
392  scalableNestingSEI.m_nestingNoOpMaxTemporalIdPlus1 = 6 + 1;  //nesting_no_op_max_temporal_id_plus1
393  scalableNestingSEI.m_nestingNumLayersMinus1        = 1 - 1;  //nesting_num_layers_minus1
394  scalableNestingSEI.m_nestingLayerId[0]             = 0;
395  scalableNestingSEI.m_callerOwnsSEIs                = true;
396#endif
397#if L0044_DU_DPB_OUTPUT_DELAY_HRD
398  Int picSptDpbOutputDuDelay = 0;
399#endif
400  UInt *accumBitsDU = NULL;
401  UInt *accumNalsDU = NULL;
402  SEIDecodingUnitInfo decodingUnitInfoSEI;
403#if SVC_EXTENSION
404  for ( Int iGOPid=iPicIdInGOP; iGOPid < iPicIdInGOP+1; iGOPid++ )
405#else
406  for ( Int iGOPid=0; iGOPid < m_iGopSize; iGOPid++ )
407#endif
408  {
409    UInt uiColDir = 1;
410    //-- For time output for each slice
411    long iBeforeTime = clock();
412
413    //select uiColDir
414    Int iCloseLeft=1, iCloseRight=-1;
415    for(Int i = 0; i<m_pcCfg->getGOPEntry(iGOPid).m_numRefPics; i++) 
416    {
417      Int iRef = m_pcCfg->getGOPEntry(iGOPid).m_referencePics[i];
418      if(iRef>0&&(iRef<iCloseRight||iCloseRight==-1))
419      {
420        iCloseRight=iRef;
421      }
422      else if(iRef<0&&(iRef>iCloseLeft||iCloseLeft==1))
423      {
424        iCloseLeft=iRef;
425      }
426    }
427    if(iCloseRight>-1)
428    {
429      iCloseRight=iCloseRight+m_pcCfg->getGOPEntry(iGOPid).m_POC-1;
430    }
431    if(iCloseLeft<1) 
432    {
433      iCloseLeft=iCloseLeft+m_pcCfg->getGOPEntry(iGOPid).m_POC-1;
434      while(iCloseLeft<0)
435      {
436        iCloseLeft+=m_iGopSize;
437      }
438    }
439    Int iLeftQP=0, iRightQP=0;
440    for(Int i=0; i<m_iGopSize; i++)
441    {
442      if(m_pcCfg->getGOPEntry(i).m_POC==(iCloseLeft%m_iGopSize)+1)
443      {
444        iLeftQP= m_pcCfg->getGOPEntry(i).m_QPOffset;
445      }
446      if (m_pcCfg->getGOPEntry(i).m_POC==(iCloseRight%m_iGopSize)+1)
447      {
448        iRightQP=m_pcCfg->getGOPEntry(i).m_QPOffset;
449      }
450    }
451    if(iCloseRight>-1&&iRightQP<iLeftQP)
452    {
453      uiColDir=0;
454    }
455
456    /////////////////////////////////////////////////////////////////////////////////////////////////// Initial to start encoding
457    Int pocCurr = iPOCLast -iNumPicRcvd+ m_pcCfg->getGOPEntry(iGOPid).m_POC;
458    Int iTimeOffset = m_pcCfg->getGOPEntry(iGOPid).m_POC;
459    if(iPOCLast == 0)
460    {
461      pocCurr=0;
462      iTimeOffset = 1;
463    }
464    if(pocCurr>=m_pcCfg->getFramesToBeEncoded())
465    {
466      continue;
467    }
468
469    if( getNalUnitType(pocCurr, m_iLastIDR) == NAL_UNIT_CODED_SLICE_IDR_W_RADL || getNalUnitType(pocCurr, m_iLastIDR) == NAL_UNIT_CODED_SLICE_IDR_N_LP )
470    {
471      m_iLastIDR = pocCurr;
472    }       
473    // start a new access unit: create an entry in the list of output access units
474    accessUnitsInGOP.push_back(AccessUnit());
475    AccessUnit& accessUnit = accessUnitsInGOP.back();
476    xGetBuffer( rcListPic, rcListPicYuvRecOut, iNumPicRcvd, iTimeOffset, pcPic, pcPicYuvRecOut, pocCurr );
477
478    //  Slice data initialization
479    pcPic->clearSliceBuffer();
480    assert(pcPic->getNumAllocatedSlice() == 1);
481    m_pcSliceEncoder->setSliceIdx(0);
482    pcPic->setCurrSliceIdx(0);
483#if SVC_EXTENSION
484    pcPic->setLayerId( m_layerId );
485    m_pcSliceEncoder->initEncSlice ( pcPic, iPOCLast, pocCurr, iNumPicRcvd, iGOPid, pcSlice, m_pcEncTop->getSPS(), m_pcEncTop->getPPS(), m_pcEncTop->getVPS() );
486#else
487    m_pcSliceEncoder->initEncSlice ( pcPic, iPOCLast, pocCurr, iNumPicRcvd, iGOPid, pcSlice, m_pcEncTop->getSPS(), m_pcEncTop->getPPS() );
488#endif
489#if M0457_COL_PICTURE_SIGNALING
490    pcSlice->setMFMEnabledFlag(m_pcEncTop->getMFMEnabledFlag());
491    pcSlice->setAltColIndicationFlag(m_pcEncTop->getNumMotionPredRefLayers() > 0 ? true : false);
492    pcSlice->setColRefLayerIdx(0);
493#endif
494    pcSlice->setLastIDR(m_iLastIDR);
495    pcSlice->setSliceIdx(0);
496    //set default slice level flag to the same as SPS level flag
497    pcSlice->setLFCrossSliceBoundaryFlag(  pcSlice->getPPS()->getLoopFilterAcrossSlicesEnabledFlag()  );
498    pcSlice->setScalingList ( m_pcEncTop->getScalingList()  );
499    pcSlice->getScalingList()->setUseTransformSkip(m_pcEncTop->getPPS()->getUseTransformSkip());
500    if(m_pcEncTop->getUseScalingListId() == SCALING_LIST_OFF)
501    {
502      m_pcEncTop->getTrQuant()->setFlatScalingList();
503      m_pcEncTop->getTrQuant()->setUseScalingList(false);
504      m_pcEncTop->getSPS()->setScalingListPresentFlag(false);
505      m_pcEncTop->getPPS()->setScalingListPresentFlag(false);
506    }
507    else if(m_pcEncTop->getUseScalingListId() == SCALING_LIST_DEFAULT)
508    {
509      pcSlice->setDefaultScalingList ();
510      m_pcEncTop->getSPS()->setScalingListPresentFlag(false);
511      m_pcEncTop->getPPS()->setScalingListPresentFlag(false);
512      m_pcEncTop->getTrQuant()->setScalingList(pcSlice->getScalingList());
513      m_pcEncTop->getTrQuant()->setUseScalingList(true);
514    }
515    else if(m_pcEncTop->getUseScalingListId() == SCALING_LIST_FILE_READ)
516    {
517      if(pcSlice->getScalingList()->xParseScalingList(m_pcCfg->getScalingListFile()))
518      {
519        pcSlice->setDefaultScalingList ();
520      }
521      pcSlice->getScalingList()->checkDcOfMatrix();
522      m_pcEncTop->getSPS()->setScalingListPresentFlag(pcSlice->checkDefaultScalingList());
523      m_pcEncTop->getPPS()->setScalingListPresentFlag(false);
524      m_pcEncTop->getTrQuant()->setScalingList(pcSlice->getScalingList());
525      m_pcEncTop->getTrQuant()->setUseScalingList(true);
526    }
527    else
528    {
529      printf("error : ScalingList == %d no support\n",m_pcEncTop->getUseScalingListId());
530      assert(0);
531    }
532
533    if(pcSlice->getSliceType()==B_SLICE&&m_pcCfg->getGOPEntry(iGOPid).m_sliceType=='P')
534    {
535      pcSlice->setSliceType(P_SLICE);
536    }
537    // Set the nal unit type
538    pcSlice->setNalUnitType(getNalUnitType(pocCurr, m_iLastIDR));
539#if SVC_EXTENSION
540    if (m_layerId > 0)
541    {
542      for( Int i = 0; i < pcSlice->getActiveNumILRRefIdx(); i++ )
543      {
544        UInt refLayerIdc = pcSlice->getInterLayerPredLayerIdc(i);
545#if VPS_EXTN_DIRECT_REF_LAYERS
546        TComList<TComPic*> *cListPic = m_ppcTEncTop[m_layerId]->getRefLayerEnc(refLayerIdc)->getListPic();
547#else
548        TComList<TComPic*> *cListPic = m_ppcTEncTop[m_layerId-1]->getListPic();
549#endif
550        pcSlice->setBaseColPic( *cListPic, refLayerIdc );
551
552#if SIMPLIFIED_MV_POS_SCALING
553#if SCALED_REF_LAYER_OFFSETS
554        const Window &scalEL = m_pcEncTop->getScaledRefLayerWindow(refLayerIdc);
555
556        Int widthBL   = pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec()->getWidth();
557        Int heightBL  = pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec()->getHeight();
558
559        Int widthEL   = pcPic->getPicYuvRec()->getWidth()  - scalEL.getWindowLeftOffset() - scalEL.getWindowRightOffset();
560        Int heightEL  = pcPic->getPicYuvRec()->getHeight() - scalEL.getWindowTopOffset()  - scalEL.getWindowBottomOffset();
561#else
562        const Window &confBL = pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec()->getConformanceWindow();
563        const Window &confEL = pcPic->getPicYuvRec()->getConformanceWindow();
564
565        Int widthBL   = pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec()->getWidth () - confBL.getWindowLeftOffset() - confBL.getWindowRightOffset();
566        Int heightBL  = pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec()->getHeight() - confBL.getWindowTopOffset() - confBL.getWindowBottomOffset();
567
568        Int widthEL   = pcPic->getPicYuvRec()->getWidth() - confEL.getWindowLeftOffset() - confEL.getWindowRightOffset();
569        Int heightEL  = pcPic->getPicYuvRec()->getHeight() - confEL.getWindowTopOffset() - confEL.getWindowBottomOffset();
570#endif
571        g_mvScalingFactor[refLayerIdc][0] = widthEL  == widthBL  ? 4096 : Clip3(-4096, 4095, ((widthEL  << 8) + (widthBL  >> 1)) / widthBL);
572        g_mvScalingFactor[refLayerIdc][1] = heightEL == heightBL ? 4096 : Clip3(-4096, 4095, ((heightEL << 8) + (heightBL >> 1)) / heightBL);
573
574        g_posScalingFactor[refLayerIdc][0] = ((widthBL  << 16) + (widthEL  >> 1)) / widthEL;
575        g_posScalingFactor[refLayerIdc][1] = ((heightBL << 16) + (heightEL >> 1)) / heightEL;
576#endif
577
578#if SVC_UPSAMPLING
579        if( pcPic->isSpatialEnhLayer(refLayerIdc))
580        { 
581#if SCALED_REF_LAYER_OFFSETS
582          m_pcPredSearch->upsampleBasePic( refLayerIdc, pcPic->getFullPelBaseRec(refLayerIdc), pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec(), pcPic->getPicYuvRec(), pcSlice->getSPS()->getScaledRefLayerWindow(refLayerIdc) );
583#else
584          m_pcPredSearch->upsampleBasePic( refLayerIdc, pcPic->getFullPelBaseRec(refLayerIdc), pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec(), pcPic->getPicYuvRec() );
585#endif
586        }
587        else
588        {
589          pcPic->setFullPelBaseRec( refLayerIdc, pcSlice->getBaseColPic(refLayerIdc)->getPicYuvRec() );
590        }
591        pcSlice->setFullPelBaseRec ( refLayerIdc, pcPic->getFullPelBaseRec(refLayerIdc) );
592#endif
593      }
594    }
595#endif
596
597#if REF_IDX_FRAMEWORK
598    if( m_layerId > 0 && (pocCurr % m_pcCfg->getIntraPeriod() == 0) )
599    {
600#if IDR_ALIGNMENT
601      TComList<TComPic*> *cListPic = m_ppcTEncTop[m_layerId]->getRefLayerEnc(0)->getListPic();
602      TComPic* picLayer0 = pcSlice->getRefPic(*cListPic, pcSlice->getPOC() );
603      if( picLayer0->getSlice(0)->getNalUnitType() == NAL_UNIT_CODED_SLICE_IDR_W_RADL || picLayer0->getSlice(0)->getNalUnitType() == NAL_UNIT_CODED_SLICE_IDR_N_LP )
604      {
605        pcSlice->setNalUnitType(picLayer0->getSlice(0)->getNalUnitType());
606      }
607      else
608#endif
609      pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_CRA);
610    }
611#if ZERO_NUM_DIRECT_LAYERS
612    if( m_layerId > 0 && pcSlice->getActiveNumILRRefIdx() == 0 && pcSlice->getNalUnitType() >= NAL_UNIT_CODED_SLICE_BLA_W_LP && pcSlice->getNalUnitType() <= NAL_UNIT_CODED_SLICE_CRA )
613    {
614      pcSlice->setSliceType(I_SLICE);
615    }
616    else
617#endif
618    if( m_layerId > 0 && !m_pcEncTop->getElRapSliceTypeB() )
619    {
620      if( (pcSlice->getNalUnitType() >= NAL_UNIT_CODED_SLICE_BLA_W_LP) &&
621          (pcSlice->getNalUnitType() <= NAL_UNIT_CODED_SLICE_CRA) &&
622           pcSlice->getSliceType() == B_SLICE )
623      {
624        pcSlice->setSliceType(P_SLICE);
625      }
626    }
627#endif
628    if(pcSlice->getTemporalLayerNonReferenceFlag())
629    {
630      if(pcSlice->getNalUnitType()==NAL_UNIT_CODED_SLICE_TRAIL_R)
631      {
632        pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_TRAIL_N);
633    }
634      if(pcSlice->getNalUnitType()==NAL_UNIT_CODED_SLICE_RADL_R)
635      {
636        pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_RADL_N);
637      }
638      if(pcSlice->getNalUnitType()==NAL_UNIT_CODED_SLICE_RASL_R)
639      {
640        pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_RASL_N);
641      }
642    }
643
644    // Do decoding refresh marking if any
645    pcSlice->decodingRefreshMarking(m_pocCRA, m_bRefreshPending, rcListPic);
646    m_pcEncTop->selectReferencePictureSet(pcSlice, pocCurr, iGOPid);
647    pcSlice->getRPS()->setNumberOfLongtermPictures(0);
648
649    if(pcSlice->checkThatAllRefPicsAreAvailable(rcListPic, pcSlice->getRPS(), false) != 0)
650    {
651      pcSlice->createExplicitReferencePictureSetFromReference(rcListPic, pcSlice->getRPS());
652    }
653    pcSlice->applyReferencePictureSet(rcListPic, pcSlice->getRPS());
654
655    if(pcSlice->getTLayer() > 0)
656    {
657      if(pcSlice->isTemporalLayerSwitchingPoint(rcListPic) || pcSlice->getSPS()->getTemporalIdNestingFlag())
658      {
659        if(pcSlice->getTemporalLayerNonReferenceFlag())
660        {
661          pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_TSA_N);
662        }
663        else
664        {
665          pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_TLA_R);
666        }
667      }
668      else if(pcSlice->isStepwiseTemporalLayerSwitchingPointCandidate(rcListPic))
669      {
670        Bool isSTSA=true;
671        for(Int ii=iGOPid+1;(ii<m_pcCfg->getGOPSize() && isSTSA==true);ii++)
672        {
673          Int lTid= m_pcCfg->getGOPEntry(ii).m_temporalId;
674          if(lTid==pcSlice->getTLayer()) 
675          {
676            TComReferencePictureSet* nRPS = pcSlice->getSPS()->getRPSList()->getReferencePictureSet(ii);
677            for(Int jj=0;jj<nRPS->getNumberOfPictures();jj++)
678            {
679              if(nRPS->getUsed(jj)) 
680              {
681                Int tPoc=m_pcCfg->getGOPEntry(ii).m_POC+nRPS->getDeltaPOC(jj);
682                Int kk=0;
683                for(kk=0;kk<m_pcCfg->getGOPSize();kk++)
684                {
685                  if(m_pcCfg->getGOPEntry(kk).m_POC==tPoc)
686                    break;
687                }
688                Int tTid=m_pcCfg->getGOPEntry(kk).m_temporalId;
689                if(tTid >= pcSlice->getTLayer())
690                {
691                  isSTSA=false;
692                  break;
693                }
694              }
695            }
696          }
697        }
698        if(isSTSA==true)
699        {   
700          if(pcSlice->getTemporalLayerNonReferenceFlag())
701          {
702            pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_STSA_N);
703          }
704          else
705          {
706            pcSlice->setNalUnitType(NAL_UNIT_CODED_SLICE_STSA_R);
707          }
708        }
709      }
710    }
711    arrangeLongtermPicturesInRPS(pcSlice, rcListPic);
712    TComRefPicListModification* refPicListModification = pcSlice->getRefPicListModification();
713    refPicListModification->setRefPicListModificationFlagL0(0);
714    refPicListModification->setRefPicListModificationFlagL1(0);
715    pcSlice->setNumRefIdx(REF_PIC_LIST_0,min(m_pcCfg->getGOPEntry(iGOPid).m_numRefPicsActive,pcSlice->getRPS()->getNumberOfPictures()));
716    pcSlice->setNumRefIdx(REF_PIC_LIST_1,min(m_pcCfg->getGOPEntry(iGOPid).m_numRefPicsActive,pcSlice->getRPS()->getNumberOfPictures()));
717
718#if REF_IDX_FRAMEWORK
719#if ZERO_NUM_DIRECT_LAYERS
720    if( m_layerId > 0 && pcSlice->getActiveNumILRRefIdx() )
721#else
722    if(m_layerId > 0)
723#endif
724    {
725#if RESTR_CHK
726      if (pcSlice->getPOC()>0  && pcSlice->isRADL()&& pcPic->getSlice(0)->isRASL())
727      {
728#if JCTVC_M0458_INTERLAYER_RPS_SIG
729        pcSlice->setActiveNumILRRefIdx(0);
730        pcSlice->setInterLayerPredEnabledFlag(0);
731#else
732        pcSlice->setNumILRRefIdx(0);
733#endif
734      }
735#endif
736#if JCTVC_M0458_INTERLAYER_RPS_SIG
737      if( pcSlice->getNalUnitType() >= NAL_UNIT_CODED_SLICE_BLA_W_LP && pcSlice->getNalUnitType() <= NAL_UNIT_CODED_SLICE_CRA )
738      {
739        pcSlice->setNumRefIdx(REF_PIC_LIST_0, pcSlice->getActiveNumILRRefIdx());
740        pcSlice->setNumRefIdx(REF_PIC_LIST_1, pcSlice->getActiveNumILRRefIdx());
741      }
742      else
743      {
744        pcSlice->setNumRefIdx(REF_PIC_LIST_0, pcSlice->getNumRefIdx(REF_PIC_LIST_0)+pcSlice->getActiveNumILRRefIdx());
745        pcSlice->setNumRefIdx(REF_PIC_LIST_1, pcSlice->getNumRefIdx(REF_PIC_LIST_1)+pcSlice->getActiveNumILRRefIdx());
746      }
747#else
748      if( pcSlice->getNalUnitType() >= NAL_UNIT_CODED_SLICE_BLA_W_LP && pcSlice->getNalUnitType() <= NAL_UNIT_CODED_SLICE_CRA )
749      {
750        pcSlice->setNumRefIdx(REF_PIC_LIST_0, pcSlice->getNumILRRefIdx());
751        pcSlice->setNumRefIdx(REF_PIC_LIST_1, pcSlice->getNumILRRefIdx());
752      }
753      else
754      {
755        pcSlice->setNumRefIdx(REF_PIC_LIST_0, pcSlice->getNumRefIdx(REF_PIC_LIST_0)+pcSlice->getNumILRRefIdx());
756        pcSlice->setNumRefIdx(REF_PIC_LIST_1, pcSlice->getNumRefIdx(REF_PIC_LIST_1)+pcSlice->getNumILRRefIdx());
757      }
758#endif
759    }
760#endif
761
762#if ADAPTIVE_QP_SELECTION
763    pcSlice->setTrQuant( m_pcEncTop->getTrQuant() );
764#endif     
765
766#if REF_IDX_FRAMEWORK
767#if M0457_COL_PICTURE_SIGNALING
768    if ( pcSlice->getSliceType() == B_SLICE && !pcSlice->getAltColIndicationFlag() )
769#else
770    if( pcSlice->getSliceType() == B_SLICE )
771#endif
772    {
773      pcSlice->setColFromL0Flag(1-uiColDir);
774    }
775#endif
776
777    //  Set reference list
778#if REF_IDX_FRAMEWORK
779#if ZERO_NUM_DIRECT_LAYERS
780    if(m_layerId ==  0 || ( m_layerId > 0 && pcSlice->getActiveNumILRRefIdx() == 0 ) )
781#else
782    if(m_layerId ==  0)
783#endif
784    {
785      pcSlice->setRefPicList( rcListPic);
786    }
787#else
788    pcSlice->setRefPicList ( rcListPic );
789#endif
790#if REF_IDX_FRAMEWORK
791#if ZERO_NUM_DIRECT_LAYERS
792    if( m_layerId > 0 && pcSlice->getActiveNumILRRefIdx() )
793#else
794    if(m_layerId > 0)
795#endif
796    {
797      m_pcEncTop->setILRPic(pcPic);
798
799#if REF_IDX_MFM
800#if M0457_COL_PICTURE_SIGNALING
801      if( pcSlice->getMFMEnabledFlag() )
802#else
803      if( pcSlice->getSPS()->getMFMEnabledFlag() )
804#endif
805      {
806        pcSlice->setRefPOCListILP(m_pcEncTop->getIlpList(), pcSlice->getBaseColPic());
807#if M0457_COL_PICTURE_SIGNALING
808        pcSlice->setMotionPredIlp(getMotionPredIlp(pcSlice));
809#endif
810      }
811#endif
812      pcSlice->setRefPicListModificationSvc();
813      pcSlice->setRefPicList( rcListPic, false, m_pcEncTop->getIlpList());
814
815#if REF_IDX_MFM
816#if M0457_COL_PICTURE_SIGNALING
817      if( pcSlice->getMFMEnabledFlag() && !pcSlice->getAltColIndicationFlag() )
818#else
819      if( pcSlice->getSPS()->getMFMEnabledFlag() )
820#endif
821      {
822        Bool found         = false;
823        UInt ColFromL0Flag = pcSlice->getColFromL0Flag();
824        UInt ColRefIdx     = pcSlice->getColRefIdx();
825        for(Int colIdx = 0; colIdx < pcSlice->getNumRefIdx( RefPicList(1 - ColFromL0Flag) ); colIdx++) 
826        { 
827          if( pcSlice->getRefPic( RefPicList(1 - ColFromL0Flag), colIdx)->isILR(m_layerId) ) 
828          { 
829            ColRefIdx = colIdx; 
830            found = true;
831            break; 
832          }
833        }
834
835        if( found == false )
836        {
837          ColFromL0Flag = 1 - ColFromL0Flag;
838          for(Int colIdx = 0; colIdx < pcSlice->getNumRefIdx( RefPicList(1 - ColFromL0Flag) ); colIdx++) 
839          { 
840            if( pcSlice->getRefPic( RefPicList(1 - ColFromL0Flag), colIdx)->isILR(m_layerId) ) 
841            { 
842              ColRefIdx = colIdx; 
843              found = true; 
844              break; 
845            } 
846          }
847        }
848
849        if(found == true)
850        {
851          pcSlice->setColFromL0Flag(ColFromL0Flag);
852          pcSlice->setColRefIdx(ColRefIdx);
853        }
854      }
855#endif
856    }
857#endif
858
859    //  Slice info. refinement
860    if ( (pcSlice->getSliceType() == B_SLICE) && (pcSlice->getNumRefIdx(REF_PIC_LIST_1) == 0) )
861    {
862      pcSlice->setSliceType ( P_SLICE );
863    }
864#if !L0034_COMBINED_LIST_CLEANUP
865    if (pcSlice->getSliceType() != B_SLICE || !pcSlice->getSPS()->getUseLComb())
866    {
867      pcSlice->setNumRefIdx(REF_PIC_LIST_C, 0);
868      pcSlice->setRefPicListCombinationFlag(false);
869      pcSlice->setRefPicListModificationFlagLC(false);
870    }
871    else
872    {
873      pcSlice->setRefPicListCombinationFlag(pcSlice->getSPS()->getUseLComb());
874      pcSlice->setNumRefIdx(REF_PIC_LIST_C, pcSlice->getNumRefIdx(REF_PIC_LIST_0));
875    }
876#endif
877
878    if (pcSlice->getSliceType() == B_SLICE)
879    {
880#if !REF_IDX_FRAMEWORK
881      pcSlice->setColFromL0Flag(1-uiColDir);
882#endif
883      Bool bLowDelay = true;
884      Int  iCurrPOC  = pcSlice->getPOC();
885      Int iRefIdx = 0;
886
887      for (iRefIdx = 0; iRefIdx < pcSlice->getNumRefIdx(REF_PIC_LIST_0) && bLowDelay; iRefIdx++)
888      {
889        if ( pcSlice->getRefPic(REF_PIC_LIST_0, iRefIdx)->getPOC() > iCurrPOC )
890        {
891          bLowDelay = false;
892        }
893      }
894      for (iRefIdx = 0; iRefIdx < pcSlice->getNumRefIdx(REF_PIC_LIST_1) && bLowDelay; iRefIdx++)
895      {
896        if ( pcSlice->getRefPic(REF_PIC_LIST_1, iRefIdx)->getPOC() > iCurrPOC )
897        {
898          bLowDelay = false;
899        }
900      }
901
902      pcSlice->setCheckLDC(bLowDelay); 
903    }
904    else
905    {
906      pcSlice->setCheckLDC(true); 
907    }
908
909    uiColDir = 1-uiColDir;
910
911    //-------------------------------------------------------------
912    pcSlice->setRefPOCList();
913
914#if L0034_COMBINED_LIST_CLEANUP
915    pcSlice->setList1IdxToList0Idx();
916#else
917    pcSlice->setNoBackPredFlag( false );
918    if ( pcSlice->getSliceType() == B_SLICE && !pcSlice->getRefPicListCombinationFlag())
919    {
920      if ( pcSlice->getNumRefIdx(RefPicList( 0 ) ) == pcSlice->getNumRefIdx(RefPicList( 1 ) ) )
921      {
922        pcSlice->setNoBackPredFlag( true );
923        Int i;
924        for ( i=0; i < pcSlice->getNumRefIdx(RefPicList( 1 ) ); i++ )
925        {
926          if ( pcSlice->getRefPOC(RefPicList(1), i) != pcSlice->getRefPOC(RefPicList(0), i) ) 
927          {
928            pcSlice->setNoBackPredFlag( false );
929            break;
930          }
931        }
932      }
933    }
934
935    if(pcSlice->getNoBackPredFlag())
936    {
937      pcSlice->setNumRefIdx(REF_PIC_LIST_C, 0);
938    }
939    pcSlice->generateCombinedList();
940#endif
941
942    if (m_pcEncTop->getTMVPModeId() == 2)
943    {
944      if (iGOPid == 0) // first picture in SOP (i.e. forward B)
945      {
946        pcSlice->setEnableTMVPFlag(0);
947      }
948      else
949      {
950        // Note: pcSlice->getColFromL0Flag() is assumed to be always 0 and getcolRefIdx() is always 0.
951        pcSlice->setEnableTMVPFlag(1);
952      }
953      pcSlice->getSPS()->setTMVPFlagsPresent(1);
954    }
955    else if (m_pcEncTop->getTMVPModeId() == 1)
956    {
957      pcSlice->getSPS()->setTMVPFlagsPresent(1);
958#if SVC_EXTENSION
959      if( pcSlice->getIdrPicFlag() )
960      {
961        pcSlice->setEnableTMVPFlag(0);
962      }
963      else
964#endif
965      pcSlice->setEnableTMVPFlag(1);
966    }
967    else
968    {
969      pcSlice->getSPS()->setTMVPFlagsPresent(0);
970      pcSlice->setEnableTMVPFlag(0);
971    }
972    /////////////////////////////////////////////////////////////////////////////////////////////////// Compress a slice
973    //  Slice compression
974    if (m_pcCfg->getUseASR())
975    {
976      m_pcSliceEncoder->setSearchRange(pcSlice);
977    }
978
979    Bool bGPBcheck=false;
980    if ( pcSlice->getSliceType() == B_SLICE)
981    {
982      if ( pcSlice->getNumRefIdx(RefPicList( 0 ) ) == pcSlice->getNumRefIdx(RefPicList( 1 ) ) )
983      {
984        bGPBcheck=true;
985        Int i;
986        for ( i=0; i < pcSlice->getNumRefIdx(RefPicList( 1 ) ); i++ )
987        {
988          if ( pcSlice->getRefPOC(RefPicList(1), i) != pcSlice->getRefPOC(RefPicList(0), i) ) 
989          {
990            bGPBcheck=false;
991            break;
992          }
993        }
994      }
995    }
996    if(bGPBcheck)
997    {
998      pcSlice->setMvdL1ZeroFlag(true);
999    }
1000    else
1001    {
1002      pcSlice->setMvdL1ZeroFlag(false);
1003    }
1004    pcPic->getSlice(pcSlice->getSliceIdx())->setMvdL1ZeroFlag(pcSlice->getMvdL1ZeroFlag());
1005
1006#if RATE_CONTROL_LAMBDA_DOMAIN
1007    Int sliceQP              = pcSlice->getSliceQp();
1008    Double lambda            = 0.0;
1009    Int actualHeadBits       = 0;
1010    Int actualTotalBits      = 0;
1011    Int estimatedBits        = 0;
1012    Int tmpBitsBeforeWriting = 0;
1013    if ( m_pcCfg->getUseRateCtrl() )
1014    {
1015      Int frameLevel = m_pcRateCtrl->getRCSeq()->getGOPID2Level( iGOPid );
1016      if ( pcPic->getSlice(0)->getSliceType() == I_SLICE )
1017      {
1018        frameLevel = 0;
1019      }
1020      m_pcRateCtrl->initRCPic( frameLevel );
1021      estimatedBits = m_pcRateCtrl->getRCPic()->getTargetBits();
1022
1023      if ( ( pcSlice->getPOC() == 0 && m_pcCfg->getInitialQP() > 0 ) || ( frameLevel == 0 && m_pcCfg->getForceIntraQP() ) ) // QP is specified
1024      {
1025        sliceQP              = m_pcCfg->getInitialQP();
1026        Int    NumberBFrames = ( m_pcCfg->getGOPSize() - 1 );
1027        Double dLambda_scale = 1.0 - Clip3( 0.0, 0.5, 0.05*(Double)NumberBFrames );
1028        Double dQPFactor     = 0.57*dLambda_scale;
1029        Int    SHIFT_QP      = 12;
1030        Int    bitdepth_luma_qp_scale = 0;
1031        Double qp_temp = (Double) sliceQP + bitdepth_luma_qp_scale - SHIFT_QP;
1032        lambda = dQPFactor*pow( 2.0, qp_temp/3.0 );
1033      }
1034      else if ( frameLevel == 0 )   // intra case, but use the model
1035      {
1036        if ( m_pcCfg->getIntraPeriod() != 1 )   // do not refine allocated bits for all intra case
1037        {
1038          Int bits = m_pcRateCtrl->getRCSeq()->getLeftAverageBits();
1039          bits = m_pcRateCtrl->getRCSeq()->getRefineBitsForIntra( bits );
1040          if ( bits < 200 )
1041          {
1042            bits = 200;
1043          }
1044          m_pcRateCtrl->getRCPic()->setTargetBits( bits );
1045        }
1046
1047        list<TEncRCPic*> listPreviousPicture = m_pcRateCtrl->getPicList();
1048        lambda  = m_pcRateCtrl->getRCPic()->estimatePicLambda( listPreviousPicture );
1049        sliceQP = m_pcRateCtrl->getRCPic()->estimatePicQP( lambda, listPreviousPicture );
1050      }
1051      else    // normal case
1052      {
1053        list<TEncRCPic*> listPreviousPicture = m_pcRateCtrl->getPicList();
1054        lambda  = m_pcRateCtrl->getRCPic()->estimatePicLambda( listPreviousPicture );
1055        sliceQP = m_pcRateCtrl->getRCPic()->estimatePicQP( lambda, listPreviousPicture );
1056      }
1057
1058      sliceQP = Clip3( -pcSlice->getSPS()->getQpBDOffsetY(), MAX_QP, sliceQP );
1059      m_pcRateCtrl->getRCPic()->setPicEstQP( sliceQP );
1060
1061      m_pcSliceEncoder->resetQP( pcPic, sliceQP, lambda );
1062    }
1063#endif
1064
1065    UInt uiNumSlices = 1;
1066
1067    UInt uiInternalAddress = pcPic->getNumPartInCU()-4;
1068    UInt uiExternalAddress = pcPic->getPicSym()->getNumberOfCUsInFrame()-1;
1069    UInt uiPosX = ( uiExternalAddress % pcPic->getFrameWidthInCU() ) * g_uiMaxCUWidth+ g_auiRasterToPelX[ g_auiZscanToRaster[uiInternalAddress] ];
1070    UInt uiPosY = ( uiExternalAddress / pcPic->getFrameWidthInCU() ) * g_uiMaxCUHeight+ g_auiRasterToPelY[ g_auiZscanToRaster[uiInternalAddress] ];
1071    UInt uiWidth = pcSlice->getSPS()->getPicWidthInLumaSamples();
1072    UInt uiHeight = pcSlice->getSPS()->getPicHeightInLumaSamples();
1073    while(uiPosX>=uiWidth||uiPosY>=uiHeight) 
1074    {
1075      uiInternalAddress--;
1076      uiPosX = ( uiExternalAddress % pcPic->getFrameWidthInCU() ) * g_uiMaxCUWidth+ g_auiRasterToPelX[ g_auiZscanToRaster[uiInternalAddress] ];
1077      uiPosY = ( uiExternalAddress / pcPic->getFrameWidthInCU() ) * g_uiMaxCUHeight+ g_auiRasterToPelY[ g_auiZscanToRaster[uiInternalAddress] ];
1078    }
1079    uiInternalAddress++;
1080    if(uiInternalAddress==pcPic->getNumPartInCU()) 
1081    {
1082      uiInternalAddress = 0;
1083      uiExternalAddress++;
1084    }
1085    UInt uiRealEndAddress = uiExternalAddress*pcPic->getNumPartInCU()+uiInternalAddress;
1086
1087    UInt uiCummulativeTileWidth;
1088    UInt uiCummulativeTileHeight;
1089    Int  p, j;
1090    UInt uiEncCUAddr;
1091
1092    //set NumColumnsMinus1 and NumRowsMinus1
1093    pcPic->getPicSym()->setNumColumnsMinus1( pcSlice->getPPS()->getNumColumnsMinus1() );
1094    pcPic->getPicSym()->setNumRowsMinus1( pcSlice->getPPS()->getNumRowsMinus1() );
1095
1096    //create the TComTileArray
1097    pcPic->getPicSym()->xCreateTComTileArray();
1098
1099    if( pcSlice->getPPS()->getUniformSpacingFlag() == 1 )
1100    {
1101      //set the width for each tile
1102      for(j=0; j < pcPic->getPicSym()->getNumRowsMinus1()+1; j++)
1103      {
1104        for(p=0; p < pcPic->getPicSym()->getNumColumnsMinus1()+1; p++)
1105        {
1106          pcPic->getPicSym()->getTComTile( j * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + p )->
1107            setTileWidth( (p+1)*pcPic->getPicSym()->getFrameWidthInCU()/(pcPic->getPicSym()->getNumColumnsMinus1()+1) 
1108            - (p*pcPic->getPicSym()->getFrameWidthInCU())/(pcPic->getPicSym()->getNumColumnsMinus1()+1) );
1109        }
1110      }
1111
1112      //set the height for each tile
1113      for(j=0; j < pcPic->getPicSym()->getNumColumnsMinus1()+1; j++)
1114      {
1115        for(p=0; p < pcPic->getPicSym()->getNumRowsMinus1()+1; p++)
1116        {
1117          pcPic->getPicSym()->getTComTile( p * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + j )->
1118            setTileHeight( (p+1)*pcPic->getPicSym()->getFrameHeightInCU()/(pcPic->getPicSym()->getNumRowsMinus1()+1) 
1119            - (p*pcPic->getPicSym()->getFrameHeightInCU())/(pcPic->getPicSym()->getNumRowsMinus1()+1) );   
1120        }
1121      }
1122    }
1123    else
1124    {
1125      //set the width for each tile
1126      for(j=0; j < pcPic->getPicSym()->getNumRowsMinus1()+1; j++)
1127      {
1128        uiCummulativeTileWidth = 0;
1129        for(p=0; p < pcPic->getPicSym()->getNumColumnsMinus1(); p++)
1130        {
1131          pcPic->getPicSym()->getTComTile( j * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + p )->setTileWidth( pcSlice->getPPS()->getColumnWidth(p) );
1132          uiCummulativeTileWidth += pcSlice->getPPS()->getColumnWidth(p);
1133        }
1134        pcPic->getPicSym()->getTComTile(j * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + p)->setTileWidth( pcPic->getPicSym()->getFrameWidthInCU()-uiCummulativeTileWidth );
1135      }
1136
1137      //set the height for each tile
1138      for(j=0; j < pcPic->getPicSym()->getNumColumnsMinus1()+1; j++)
1139      {
1140        uiCummulativeTileHeight = 0;
1141        for(p=0; p < pcPic->getPicSym()->getNumRowsMinus1(); p++)
1142        {
1143          pcPic->getPicSym()->getTComTile( p * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + j )->setTileHeight( pcSlice->getPPS()->getRowHeight(p) );
1144          uiCummulativeTileHeight += pcSlice->getPPS()->getRowHeight(p);
1145        }
1146        pcPic->getPicSym()->getTComTile(p * (pcPic->getPicSym()->getNumColumnsMinus1()+1) + j)->setTileHeight( pcPic->getPicSym()->getFrameHeightInCU()-uiCummulativeTileHeight );
1147      }
1148    }
1149    //intialize each tile of the current picture
1150    pcPic->getPicSym()->xInitTiles();
1151
1152    // Allocate some coders, now we know how many tiles there are.
1153    Int iNumSubstreams = pcSlice->getPPS()->getNumSubstreams();
1154
1155    //generate the Coding Order Map and Inverse Coding Order Map
1156    for(p=0, uiEncCUAddr=0; p<pcPic->getPicSym()->getNumberOfCUsInFrame(); p++, uiEncCUAddr = pcPic->getPicSym()->xCalculateNxtCUAddr(uiEncCUAddr))
1157    {
1158      pcPic->getPicSym()->setCUOrderMap(p, uiEncCUAddr);
1159      pcPic->getPicSym()->setInverseCUOrderMap(uiEncCUAddr, p);
1160    }
1161    pcPic->getPicSym()->setCUOrderMap(pcPic->getPicSym()->getNumberOfCUsInFrame(), pcPic->getPicSym()->getNumberOfCUsInFrame());   
1162    pcPic->getPicSym()->setInverseCUOrderMap(pcPic->getPicSym()->getNumberOfCUsInFrame(), pcPic->getPicSym()->getNumberOfCUsInFrame());
1163
1164    // Allocate some coders, now we know how many tiles there are.
1165    m_pcEncTop->createWPPCoders(iNumSubstreams);
1166    pcSbacCoders = m_pcEncTop->getSbacCoders();
1167    pcSubstreamsOut = new TComOutputBitstream[iNumSubstreams];
1168
1169    UInt startCUAddrSliceIdx = 0; // used to index "m_uiStoredStartCUAddrForEncodingSlice" containing locations of slice boundaries
1170    UInt startCUAddrSlice    = 0; // used to keep track of current slice's starting CU addr.
1171    pcSlice->setSliceCurStartCUAddr( startCUAddrSlice ); // Setting "start CU addr" for current slice
1172    m_storedStartCUAddrForEncodingSlice.clear();
1173
1174    UInt startCUAddrSliceSegmentIdx = 0; // used to index "m_uiStoredStartCUAddrForEntropyEncodingSlice" containing locations of slice boundaries
1175    UInt startCUAddrSliceSegment    = 0; // used to keep track of current Dependent slice's starting CU addr.
1176    pcSlice->setSliceSegmentCurStartCUAddr( startCUAddrSliceSegment ); // Setting "start CU addr" for current Dependent slice
1177
1178    m_storedStartCUAddrForEncodingSliceSegment.clear();
1179    UInt nextCUAddr = 0;
1180    m_storedStartCUAddrForEncodingSlice.push_back (nextCUAddr);
1181    startCUAddrSliceIdx++;
1182    m_storedStartCUAddrForEncodingSliceSegment.push_back(nextCUAddr);
1183    startCUAddrSliceSegmentIdx++;
1184#if AVC_BASE
1185    if( m_layerId == 0 && m_pcEncTop->getVPS()->getAvcBaseLayerFlag() )
1186    {
1187      pcPic->getPicYuvOrg()->copyToPic( pcPic->getPicYuvRec() );
1188#if AVC_SYNTAX
1189      pcPic->readBLSyntax( m_ppcTEncTop[0]->getBLSyntaxFile(), SYNTAX_BYTES );
1190#endif
1191      return;
1192    }
1193#endif
1194
1195    while(nextCUAddr<uiRealEndAddress) // determine slice boundaries
1196    {
1197      pcSlice->setNextSlice       ( false );
1198      pcSlice->setNextSliceSegment( false );
1199      assert(pcPic->getNumAllocatedSlice() == startCUAddrSliceIdx);
1200      m_pcSliceEncoder->precompressSlice( pcPic );
1201      m_pcSliceEncoder->compressSlice   ( pcPic );
1202
1203      Bool bNoBinBitConstraintViolated = (!pcSlice->isNextSlice() && !pcSlice->isNextSliceSegment());
1204      if (pcSlice->isNextSlice() || (bNoBinBitConstraintViolated && m_pcCfg->getSliceMode()==FIXED_NUMBER_OF_LCU))
1205      {
1206        startCUAddrSlice = pcSlice->getSliceCurEndCUAddr();
1207        // Reconstruction slice
1208        m_storedStartCUAddrForEncodingSlice.push_back(startCUAddrSlice);
1209        startCUAddrSliceIdx++;
1210        // Dependent slice
1211        if (startCUAddrSliceSegmentIdx>0 && m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx-1] != startCUAddrSlice)
1212        {
1213          m_storedStartCUAddrForEncodingSliceSegment.push_back(startCUAddrSlice);
1214          startCUAddrSliceSegmentIdx++;
1215        }
1216
1217        if (startCUAddrSlice < uiRealEndAddress)
1218        {
1219          pcPic->allocateNewSlice();         
1220          pcPic->setCurrSliceIdx                  ( startCUAddrSliceIdx-1 );
1221          m_pcSliceEncoder->setSliceIdx           ( startCUAddrSliceIdx-1 );
1222          pcSlice = pcPic->getSlice               ( startCUAddrSliceIdx-1 );
1223          pcSlice->copySliceInfo                  ( pcPic->getSlice(0)      );
1224          pcSlice->setSliceIdx                    ( startCUAddrSliceIdx-1 );
1225          pcSlice->setSliceCurStartCUAddr         ( startCUAddrSlice      );
1226          pcSlice->setSliceSegmentCurStartCUAddr  ( startCUAddrSlice      );
1227          pcSlice->setSliceBits(0);
1228          uiNumSlices ++;
1229        }
1230      }
1231      else if (pcSlice->isNextSliceSegment() || (bNoBinBitConstraintViolated && m_pcCfg->getSliceSegmentMode()==FIXED_NUMBER_OF_LCU))
1232      {
1233        startCUAddrSliceSegment                                                     = pcSlice->getSliceSegmentCurEndCUAddr();
1234        m_storedStartCUAddrForEncodingSliceSegment.push_back(startCUAddrSliceSegment);
1235        startCUAddrSliceSegmentIdx++;
1236        pcSlice->setSliceSegmentCurStartCUAddr( startCUAddrSliceSegment );
1237      }
1238      else
1239      {
1240        startCUAddrSlice                                                            = pcSlice->getSliceCurEndCUAddr();
1241        startCUAddrSliceSegment                                                     = pcSlice->getSliceSegmentCurEndCUAddr();
1242      }       
1243
1244      nextCUAddr = (startCUAddrSlice > startCUAddrSliceSegment) ? startCUAddrSlice : startCUAddrSliceSegment;
1245    }
1246    m_storedStartCUAddrForEncodingSlice.push_back( pcSlice->getSliceCurEndCUAddr());
1247    startCUAddrSliceIdx++;
1248    m_storedStartCUAddrForEncodingSliceSegment.push_back(pcSlice->getSliceCurEndCUAddr());
1249    startCUAddrSliceSegmentIdx++;
1250
1251    pcSlice = pcPic->getSlice(0);
1252
1253    // SAO parameter estimation using non-deblocked pixels for LCU bottom and right boundary areas
1254    if( m_pcCfg->getSaoLcuBasedOptimization() && m_pcCfg->getSaoLcuBoundary() )
1255    {
1256      m_pcSAO->resetStats();
1257      m_pcSAO->calcSaoStatsCu_BeforeDblk( pcPic );
1258    }
1259
1260    //-- Loop filter
1261    Bool bLFCrossTileBoundary = pcSlice->getPPS()->getLoopFilterAcrossTilesEnabledFlag();
1262    m_pcLoopFilter->setCfg(bLFCrossTileBoundary);
1263#if L0386_DB_METRIC
1264    if ( m_pcCfg->getDeblockingFilterMetric() )
1265    {
1266      dblMetric(pcPic, uiNumSlices);
1267    }
1268#endif
1269    m_pcLoopFilter->loopFilterPic( pcPic );
1270
1271    pcSlice = pcPic->getSlice(0);
1272    if(pcSlice->getSPS()->getUseSAO())
1273    {
1274      std::vector<Bool> LFCrossSliceBoundaryFlag;
1275      for(Int s=0; s< uiNumSlices; s++)
1276      {
1277        LFCrossSliceBoundaryFlag.push_back(  ((uiNumSlices==1)?true:pcPic->getSlice(s)->getLFCrossSliceBoundaryFlag()) );
1278      }
1279      m_storedStartCUAddrForEncodingSlice.resize(uiNumSlices+1);
1280      pcPic->createNonDBFilterInfo(m_storedStartCUAddrForEncodingSlice, 0, &LFCrossSliceBoundaryFlag ,pcPic->getPicSym()->getNumTiles() ,bLFCrossTileBoundary);
1281    }
1282
1283
1284    pcSlice = pcPic->getSlice(0);
1285
1286    if(pcSlice->getSPS()->getUseSAO())
1287    {
1288      m_pcSAO->createPicSaoInfo(pcPic);
1289    }
1290
1291    /////////////////////////////////////////////////////////////////////////////////////////////////// File writing
1292    // Set entropy coder
1293    m_pcEntropyCoder->setEntropyCoder   ( m_pcCavlcCoder, pcSlice );
1294
1295    /* write various header sets. */
1296    if ( m_bSeqFirst )
1297    {
1298#if SVC_EXTENSION
1299      OutputNALUnit nalu(NAL_UNIT_VPS, 0, m_layerId);
1300#if AVC_BASE
1301      if( ( m_layerId == 1 && m_pcEncTop->getVPS()->getAvcBaseLayerFlag() ) || ( m_layerId == 0 && !m_pcEncTop->getVPS()->getAvcBaseLayerFlag() ) )
1302#else
1303      if( m_layerId == 0 )
1304#endif
1305      {
1306#else
1307      OutputNALUnit nalu(NAL_UNIT_VPS);
1308#endif
1309      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1310      m_pcEntropyCoder->encodeVPS(m_pcEncTop->getVPS());
1311      writeRBSPTrailingBits(nalu.m_Bitstream);
1312      accessUnit.push_back(new NALUnitEBSP(nalu));
1313#if RATE_CONTROL_LAMBDA_DOMAIN
1314      actualTotalBits += UInt(accessUnit.back()->m_nalUnitData.str().size()) * 8;
1315#endif
1316#if SVC_EXTENSION
1317      }
1318#endif
1319
1320#if SVC_EXTENSION
1321      nalu = NALUnit(NAL_UNIT_SPS, 0, m_layerId);
1322#else
1323      nalu = NALUnit(NAL_UNIT_SPS);
1324#endif
1325      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1326      if (m_bSeqFirst)
1327      {
1328        pcSlice->getSPS()->setNumLongTermRefPicSPS(m_numLongTermRefPicSPS);
1329        for (Int k = 0; k < m_numLongTermRefPicSPS; k++)
1330        {
1331          pcSlice->getSPS()->setLtRefPicPocLsbSps(k, m_ltRefPicPocLsbSps[k]);
1332          pcSlice->getSPS()->setUsedByCurrPicLtSPSFlag(k, m_ltRefPicUsedByCurrPicFlag[k]);
1333        }
1334      }
1335      if( m_pcCfg->getPictureTimingSEIEnabled() || m_pcCfg->getDecodingUnitInfoSEIEnabled() )
1336      {
1337        UInt maxCU = m_pcCfg->getSliceArgument() >> ( pcSlice->getSPS()->getMaxCUDepth() << 1);
1338        UInt numDU = ( m_pcCfg->getSliceMode() == 1 ) ? ( pcPic->getNumCUsInFrame() / maxCU ) : ( 0 );
1339        if( pcPic->getNumCUsInFrame() % maxCU != 0 )
1340        {
1341          numDU ++;
1342        }
1343        pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->setNumDU( numDU );
1344        pcSlice->getSPS()->setHrdParameters( m_pcCfg->getFrameRate(), numDU, m_pcCfg->getTargetBitrate(), ( m_pcCfg->getIntraPeriod() > 0 ) );
1345      }
1346      if( m_pcCfg->getBufferingPeriodSEIEnabled() || m_pcCfg->getPictureTimingSEIEnabled() || m_pcCfg->getDecodingUnitInfoSEIEnabled() )
1347      {
1348        pcSlice->getSPS()->getVuiParameters()->setHrdParametersPresentFlag( true );
1349      }
1350      m_pcEntropyCoder->encodeSPS(pcSlice->getSPS());
1351      writeRBSPTrailingBits(nalu.m_Bitstream);
1352      accessUnit.push_back(new NALUnitEBSP(nalu));
1353#if RATE_CONTROL_LAMBDA_DOMAIN
1354      actualTotalBits += UInt(accessUnit.back()->m_nalUnitData.str().size()) * 8;
1355#endif
1356
1357#if SVC_EXTENSION
1358      nalu = NALUnit(NAL_UNIT_PPS, 0, m_layerId);
1359#else
1360      nalu = NALUnit(NAL_UNIT_PPS);
1361#endif
1362      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1363      m_pcEntropyCoder->encodePPS(pcSlice->getPPS());
1364      writeRBSPTrailingBits(nalu.m_Bitstream);
1365      accessUnit.push_back(new NALUnitEBSP(nalu));
1366#if RATE_CONTROL_LAMBDA_DOMAIN
1367      actualTotalBits += UInt(accessUnit.back()->m_nalUnitData.str().size()) * 8;
1368#endif
1369
1370      xCreateLeadingSEIMessages(accessUnit, pcSlice->getSPS());
1371
1372      m_bSeqFirst = false;
1373    }
1374
1375#if L0208_SOP_DESCRIPTION_SEI
1376    if (writeSOP) // write SOP description SEI (if enabled) at the beginning of GOP
1377    {
1378      Int SOPcurrPOC = pocCurr;
1379
1380      OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI);
1381      m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1382      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1383
1384      SEISOPDescription SOPDescriptionSEI;
1385      SOPDescriptionSEI.m_sopSeqParameterSetId = pcSlice->getSPS()->getSPSId();
1386
1387      UInt i = 0;
1388      UInt prevEntryId = iGOPid;
1389      for (j = iGOPid; j < m_iGopSize; j++)
1390      {
1391        Int deltaPOC = m_pcCfg->getGOPEntry(j).m_POC - m_pcCfg->getGOPEntry(prevEntryId).m_POC;
1392        if ((SOPcurrPOC + deltaPOC) < m_pcCfg->getFramesToBeEncoded())
1393        {
1394          SOPcurrPOC += deltaPOC;
1395          SOPDescriptionSEI.m_sopDescVclNaluType[i] = getNalUnitType(SOPcurrPOC, m_iLastIDR);
1396          SOPDescriptionSEI.m_sopDescTemporalId[i] = m_pcCfg->getGOPEntry(j).m_temporalId;
1397          SOPDescriptionSEI.m_sopDescStRpsIdx[i] = m_pcEncTop->getReferencePictureSetIdxForSOP(pcSlice, SOPcurrPOC, j);
1398          SOPDescriptionSEI.m_sopDescPocDelta[i] = deltaPOC;
1399
1400          prevEntryId = j;
1401          i++;
1402        }
1403      }
1404
1405      SOPDescriptionSEI.m_numPicsInSopMinus1 = i - 1;
1406
1407      m_seiWriter.writeSEImessage( nalu.m_Bitstream, SOPDescriptionSEI, pcSlice->getSPS());
1408      writeRBSPTrailingBits(nalu.m_Bitstream);
1409      accessUnit.push_back(new NALUnitEBSP(nalu));
1410
1411      writeSOP = false;
1412    }
1413#endif
1414
1415    if( ( m_pcCfg->getPictureTimingSEIEnabled() || m_pcCfg->getDecodingUnitInfoSEIEnabled() ) &&
1416        ( pcSlice->getSPS()->getVuiParametersPresentFlag() ) &&
1417        ( ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getNalHrdParametersPresentFlag() ) 
1418       || ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getVclHrdParametersPresentFlag() ) ) )
1419    {
1420      if( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getSubPicCpbParamsPresentFlag() )
1421      {
1422        UInt numDU = pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getNumDU();
1423        pictureTimingSEI.m_numDecodingUnitsMinus1     = ( numDU - 1 );
1424        pictureTimingSEI.m_duCommonCpbRemovalDelayFlag = false;
1425
1426        if( pictureTimingSEI.m_numNalusInDuMinus1 == NULL )
1427        {
1428          pictureTimingSEI.m_numNalusInDuMinus1       = new UInt[ numDU ];
1429        }
1430        if( pictureTimingSEI.m_duCpbRemovalDelayMinus1  == NULL )
1431        {
1432          pictureTimingSEI.m_duCpbRemovalDelayMinus1  = new UInt[ numDU ];
1433        }
1434        if( accumBitsDU == NULL )
1435        {
1436          accumBitsDU                                  = new UInt[ numDU ];
1437        }
1438        if( accumNalsDU == NULL )
1439        {
1440          accumNalsDU                                  = new UInt[ numDU ];
1441        }
1442      }
1443      pictureTimingSEI.m_auCpbRemovalDelay = std::max<Int>(1, m_totalCoded - m_lastBPSEI); // Syntax element signalled as minus, hence the .
1444      pictureTimingSEI.m_picDpbOutputDelay = pcSlice->getSPS()->getNumReorderPics(0) + pcSlice->getPOC() - m_totalCoded;
1445#if L0044_DU_DPB_OUTPUT_DELAY_HRD
1446      Int factor = pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getTickDivisorMinus2() + 2;
1447      pictureTimingSEI.m_picDpbOutputDuDelay = factor * pictureTimingSEI.m_picDpbOutputDelay;
1448      if( m_pcCfg->getDecodingUnitInfoSEIEnabled() )
1449      {
1450        picSptDpbOutputDuDelay = factor * pictureTimingSEI.m_picDpbOutputDelay;
1451      }
1452#endif
1453    }
1454
1455    if( ( m_pcCfg->getBufferingPeriodSEIEnabled() ) && ( pcSlice->getSliceType() == I_SLICE ) &&
1456        ( pcSlice->getSPS()->getVuiParametersPresentFlag() ) && 
1457        ( ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getNalHrdParametersPresentFlag() ) 
1458       || ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getVclHrdParametersPresentFlag() ) ) )
1459    {
1460      OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI);
1461      m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1462      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1463
1464      SEIBufferingPeriod sei_buffering_period;
1465     
1466      UInt uiInitialCpbRemovalDelay = (90000/2);                      // 0.5 sec
1467      sei_buffering_period.m_initialCpbRemovalDelay      [0][0]     = uiInitialCpbRemovalDelay;
1468      sei_buffering_period.m_initialCpbRemovalDelayOffset[0][0]     = uiInitialCpbRemovalDelay;
1469      sei_buffering_period.m_initialCpbRemovalDelay      [0][1]     = uiInitialCpbRemovalDelay;
1470      sei_buffering_period.m_initialCpbRemovalDelayOffset[0][1]     = uiInitialCpbRemovalDelay;
1471
1472#if L0043_TIMING_INFO
1473      Double dTmp = (Double)pcSlice->getSPS()->getVuiParameters()->getTimingInfo()->getNumUnitsInTick() / (Double)pcSlice->getSPS()->getVuiParameters()->getTimingInfo()->getTimeScale();
1474#else
1475      Double dTmp = (Double)pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getNumUnitsInTick() / (Double)pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getTimeScale();
1476#endif
1477
1478      UInt uiTmp = (UInt)( dTmp * 90000.0 ); 
1479      uiInitialCpbRemovalDelay -= uiTmp;
1480      uiInitialCpbRemovalDelay -= uiTmp / ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getTickDivisorMinus2() + 2 );
1481      sei_buffering_period.m_initialAltCpbRemovalDelay      [0][0]  = uiInitialCpbRemovalDelay;
1482      sei_buffering_period.m_initialAltCpbRemovalDelayOffset[0][0]  = uiInitialCpbRemovalDelay;
1483      sei_buffering_period.m_initialAltCpbRemovalDelay      [0][1]  = uiInitialCpbRemovalDelay;
1484      sei_buffering_period.m_initialAltCpbRemovalDelayOffset[0][1]  = uiInitialCpbRemovalDelay;
1485
1486      sei_buffering_period.m_rapCpbParamsPresentFlag              = 0;
1487#if L0328_SPLICING
1488      //for the concatenation, it can be set to one during splicing.
1489      sei_buffering_period.m_concatenationFlag = 0;
1490      //since the temporal layer HRD is not ready, we assumed it is fixed
1491      sei_buffering_period.m_auCpbRemovalDelayDelta = 1;
1492#endif
1493#if L0044_CPB_DPB_DELAY_OFFSET
1494      sei_buffering_period.m_cpbDelayOffset = 0;
1495      sei_buffering_period.m_dpbDelayOffset = 0;
1496#endif
1497
1498      m_seiWriter.writeSEImessage( nalu.m_Bitstream, sei_buffering_period, pcSlice->getSPS());
1499      writeRBSPTrailingBits(nalu.m_Bitstream);
1500#if L0045_NON_NESTED_SEI_RESTRICTIONS
1501      {
1502      UInt seiPositionInAu = xGetFirstSeiLocation(accessUnit);
1503      UInt offsetPosition = m_activeParameterSetSEIPresentInAU;   // Insert BP SEI after APS SEI
1504      AccessUnit::iterator it;
1505      for(j = 0, it = accessUnit.begin(); j < seiPositionInAu + offsetPosition; j++)
1506      {
1507        it++;
1508      }
1509      accessUnit.insert(it, new NALUnitEBSP(nalu));
1510      m_bufferingPeriodSEIPresentInAU = true;
1511      }
1512#else
1513      accessUnit.push_back(new NALUnitEBSP(nalu));
1514#endif
1515
1516#if K0180_SCALABLE_NESTING_SEI
1517      if (m_pcCfg->getScalableNestingSEIEnabled())
1518      {
1519        OutputNALUnit naluTmp(NAL_UNIT_PREFIX_SEI);
1520        m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1521        m_pcEntropyCoder->setBitstream(&naluTmp.m_Bitstream);
1522        scalableNestingSEI.m_nestedSEIs.clear();
1523        scalableNestingSEI.m_nestedSEIs.push_back(&sei_buffering_period);
1524        m_seiWriter.writeSEImessage( naluTmp.m_Bitstream, scalableNestingSEI, pcSlice->getSPS());
1525        writeRBSPTrailingBits(naluTmp.m_Bitstream);
1526#if L0045_NON_NESTED_SEI_RESTRICTIONS
1527        UInt seiPositionInAu = xGetFirstSeiLocation(accessUnit);
1528        UInt offsetPosition = m_activeParameterSetSEIPresentInAU + m_bufferingPeriodSEIPresentInAU + m_pictureTimingSEIPresentInAU;   // Insert BP SEI after non-nested APS, BP and PT SEIs
1529        AccessUnit::iterator it;
1530        for(j = 0, it = accessUnit.begin(); j < seiPositionInAu + offsetPosition; j++)
1531        {
1532          it++;
1533        }
1534        accessUnit.insert(it, new NALUnitEBSP(naluTmp));
1535        m_nestedBufferingPeriodSEIPresentInAU = true;
1536#else
1537        accessUnit.push_back(new NALUnitEBSP(naluTmp));
1538#endif
1539      }
1540#endif
1541
1542      m_lastBPSEI = m_totalCoded;
1543      m_cpbRemovalDelay = 0;
1544    }
1545    m_cpbRemovalDelay ++;
1546    if( ( m_pcEncTop->getRecoveryPointSEIEnabled() ) && ( pcSlice->getSliceType() == I_SLICE ) )
1547    {
1548      if( m_pcEncTop->getGradualDecodingRefreshInfoEnabled() && !pcSlice->getRapPicFlag() )
1549      {
1550        // Gradual decoding refresh SEI
1551        OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI);
1552        m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1553        m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1554
1555        SEIGradualDecodingRefreshInfo seiGradualDecodingRefreshInfo;
1556        seiGradualDecodingRefreshInfo.m_gdrForegroundFlag = true; // Indicating all "foreground"
1557
1558        m_seiWriter.writeSEImessage( nalu.m_Bitstream, seiGradualDecodingRefreshInfo, pcSlice->getSPS() );
1559        writeRBSPTrailingBits(nalu.m_Bitstream);
1560        accessUnit.push_back(new NALUnitEBSP(nalu));
1561      }
1562    // Recovery point SEI
1563      OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI);
1564      m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1565      m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1566
1567      SEIRecoveryPoint sei_recovery_point;
1568      sei_recovery_point.m_recoveryPocCnt    = 0;
1569      sei_recovery_point.m_exactMatchingFlag = ( pcSlice->getPOC() == 0 ) ? (true) : (false);
1570      sei_recovery_point.m_brokenLinkFlag    = false;
1571
1572      m_seiWriter.writeSEImessage( nalu.m_Bitstream, sei_recovery_point, pcSlice->getSPS() );
1573      writeRBSPTrailingBits(nalu.m_Bitstream);
1574      accessUnit.push_back(new NALUnitEBSP(nalu));
1575    }
1576
1577    /* use the main bitstream buffer for storing the marshalled picture */
1578    m_pcEntropyCoder->setBitstream(NULL);
1579
1580    startCUAddrSliceIdx = 0;
1581    startCUAddrSlice    = 0; 
1582
1583    startCUAddrSliceSegmentIdx = 0;
1584    startCUAddrSliceSegment    = 0; 
1585    nextCUAddr                 = 0;
1586    pcSlice = pcPic->getSlice(startCUAddrSliceIdx);
1587
1588    Int processingState = (pcSlice->getSPS()->getUseSAO())?(EXECUTE_INLOOPFILTER):(ENCODE_SLICE);
1589    Bool skippedSlice=false;
1590    while (nextCUAddr < uiRealEndAddress) // Iterate over all slices
1591    {
1592      switch(processingState)
1593      {
1594      case ENCODE_SLICE:
1595        {
1596          pcSlice->setNextSlice       ( false );
1597          pcSlice->setNextSliceSegment( false );
1598          if (nextCUAddr == m_storedStartCUAddrForEncodingSlice[startCUAddrSliceIdx])
1599          {
1600            pcSlice = pcPic->getSlice(startCUAddrSliceIdx);
1601            if(startCUAddrSliceIdx > 0 && pcSlice->getSliceType()!= I_SLICE)
1602            {
1603              pcSlice->checkColRefIdx(startCUAddrSliceIdx, pcPic);
1604            }
1605            pcPic->setCurrSliceIdx(startCUAddrSliceIdx);
1606            m_pcSliceEncoder->setSliceIdx(startCUAddrSliceIdx);
1607            assert(startCUAddrSliceIdx == pcSlice->getSliceIdx());
1608            // Reconstruction slice
1609            pcSlice->setSliceCurStartCUAddr( nextCUAddr );  // to be used in encodeSlice() + context restriction
1610            pcSlice->setSliceCurEndCUAddr  ( m_storedStartCUAddrForEncodingSlice[startCUAddrSliceIdx+1 ] );
1611            // Dependent slice
1612            pcSlice->setSliceSegmentCurStartCUAddr( nextCUAddr );  // to be used in encodeSlice() + context restriction
1613            pcSlice->setSliceSegmentCurEndCUAddr  ( m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx+1 ] );
1614
1615            pcSlice->setNextSlice       ( true );
1616
1617            startCUAddrSliceIdx++;
1618            startCUAddrSliceSegmentIdx++;
1619          } 
1620          else if (nextCUAddr == m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx])
1621          {
1622            // Dependent slice
1623            pcSlice->setSliceSegmentCurStartCUAddr( nextCUAddr );  // to be used in encodeSlice() + context restriction
1624            pcSlice->setSliceSegmentCurEndCUAddr  ( m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx+1 ] );
1625
1626            pcSlice->setNextSliceSegment( true );
1627
1628            startCUAddrSliceSegmentIdx++;
1629          }
1630#if REF_IDX_FRAMEWORK && M0457_COL_PICTURE_SIGNALING
1631          pcSlice->setNumMotionPredRefLayers(m_pcEncTop->getNumMotionPredRefLayers());
1632#endif
1633          pcSlice->setRPS(pcPic->getSlice(0)->getRPS());
1634          pcSlice->setRPSidx(pcPic->getSlice(0)->getRPSidx());
1635          UInt uiDummyStartCUAddr;
1636          UInt uiDummyBoundingCUAddr;
1637          m_pcSliceEncoder->xDetermineStartAndBoundingCUAddr(uiDummyStartCUAddr,uiDummyBoundingCUAddr,pcPic,true);
1638
1639          uiInternalAddress = pcPic->getPicSym()->getPicSCUAddr(pcSlice->getSliceSegmentCurEndCUAddr()-1) % pcPic->getNumPartInCU();
1640          uiExternalAddress = pcPic->getPicSym()->getPicSCUAddr(pcSlice->getSliceSegmentCurEndCUAddr()-1) / pcPic->getNumPartInCU();
1641          uiPosX = ( uiExternalAddress % pcPic->getFrameWidthInCU() ) * g_uiMaxCUWidth+ g_auiRasterToPelX[ g_auiZscanToRaster[uiInternalAddress] ];
1642          uiPosY = ( uiExternalAddress / pcPic->getFrameWidthInCU() ) * g_uiMaxCUHeight+ g_auiRasterToPelY[ g_auiZscanToRaster[uiInternalAddress] ];
1643          uiWidth = pcSlice->getSPS()->getPicWidthInLumaSamples();
1644          uiHeight = pcSlice->getSPS()->getPicHeightInLumaSamples();
1645          while(uiPosX>=uiWidth||uiPosY>=uiHeight)
1646          {
1647            uiInternalAddress--;
1648            uiPosX = ( uiExternalAddress % pcPic->getFrameWidthInCU() ) * g_uiMaxCUWidth+ g_auiRasterToPelX[ g_auiZscanToRaster[uiInternalAddress] ];
1649            uiPosY = ( uiExternalAddress / pcPic->getFrameWidthInCU() ) * g_uiMaxCUHeight+ g_auiRasterToPelY[ g_auiZscanToRaster[uiInternalAddress] ];
1650          }
1651          uiInternalAddress++;
1652          if(uiInternalAddress==pcPic->getNumPartInCU())
1653          {
1654            uiInternalAddress = 0;
1655            uiExternalAddress = pcPic->getPicSym()->getCUOrderMap(pcPic->getPicSym()->getInverseCUOrderMap(uiExternalAddress)+1);
1656          }
1657          UInt endAddress = pcPic->getPicSym()->getPicSCUEncOrder(uiExternalAddress*pcPic->getNumPartInCU()+uiInternalAddress);
1658          if(endAddress<=pcSlice->getSliceSegmentCurStartCUAddr()) 
1659          {
1660            UInt boundingAddrSlice, boundingAddrSliceSegment;
1661            boundingAddrSlice          = m_storedStartCUAddrForEncodingSlice[startCUAddrSliceIdx];         
1662            boundingAddrSliceSegment = m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx];         
1663            nextCUAddr               = min(boundingAddrSlice, boundingAddrSliceSegment);
1664            if(pcSlice->isNextSlice())
1665            {
1666              skippedSlice=true;
1667            }
1668            continue;
1669          }
1670          if(skippedSlice) 
1671          {
1672            pcSlice->setNextSlice       ( true );
1673            pcSlice->setNextSliceSegment( false );
1674          }
1675          skippedSlice=false;
1676          pcSlice->allocSubstreamSizes( iNumSubstreams );
1677          for ( UInt ui = 0 ; ui < iNumSubstreams; ui++ )
1678          {
1679            pcSubstreamsOut[ui].clear();
1680          }
1681
1682          m_pcEntropyCoder->setEntropyCoder   ( m_pcCavlcCoder, pcSlice );
1683          m_pcEntropyCoder->resetEntropy      ();
1684          /* start slice NALunit */
1685#if SVC_EXTENSION
1686          OutputNALUnit nalu( pcSlice->getNalUnitType(), pcSlice->getTLayer(), m_layerId );
1687#else
1688          OutputNALUnit nalu( pcSlice->getNalUnitType(), pcSlice->getTLayer() );
1689#endif
1690          Bool sliceSegment = (!pcSlice->isNextSlice());
1691          if (!sliceSegment)
1692          {
1693            uiOneBitstreamPerSliceLength = 0; // start of a new slice
1694          }
1695          m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1696#if RATE_CONTROL_LAMBDA_DOMAIN
1697          tmpBitsBeforeWriting = m_pcEntropyCoder->getNumberOfWrittenBits();
1698#endif
1699          m_pcEntropyCoder->encodeSliceHeader(pcSlice);
1700#if RATE_CONTROL_LAMBDA_DOMAIN
1701          actualHeadBits += ( m_pcEntropyCoder->getNumberOfWrittenBits() - tmpBitsBeforeWriting );
1702#endif
1703
1704          // is it needed?
1705          {
1706            if (!sliceSegment)
1707            {
1708              pcBitstreamRedirect->writeAlignOne();
1709            }
1710            else
1711            {
1712              // We've not completed our slice header info yet, do the alignment later.
1713            }
1714            m_pcSbacCoder->init( (TEncBinIf*)m_pcBinCABAC );
1715            m_pcEntropyCoder->setEntropyCoder ( m_pcSbacCoder, pcSlice );
1716            m_pcEntropyCoder->resetEntropy    ();
1717            for ( UInt ui = 0 ; ui < pcSlice->getPPS()->getNumSubstreams() ; ui++ )
1718            {
1719              m_pcEntropyCoder->setEntropyCoder ( &pcSbacCoders[ui], pcSlice );
1720              m_pcEntropyCoder->resetEntropy    ();
1721            }
1722          }
1723
1724          if(pcSlice->isNextSlice())
1725          {
1726            // set entropy coder for writing
1727            m_pcSbacCoder->init( (TEncBinIf*)m_pcBinCABAC );
1728            {
1729              for ( UInt ui = 0 ; ui < pcSlice->getPPS()->getNumSubstreams() ; ui++ )
1730              {
1731                m_pcEntropyCoder->setEntropyCoder ( &pcSbacCoders[ui], pcSlice );
1732                m_pcEntropyCoder->resetEntropy    ();
1733              }
1734              pcSbacCoders[0].load(m_pcSbacCoder);
1735              m_pcEntropyCoder->setEntropyCoder ( &pcSbacCoders[0], pcSlice );  //ALF is written in substream #0 with CABAC coder #0 (see ALF param encoding below)
1736            }
1737            m_pcEntropyCoder->resetEntropy    ();
1738            // File writing
1739            if (!sliceSegment)
1740            {
1741              m_pcEntropyCoder->setBitstream(pcBitstreamRedirect);
1742            }
1743            else
1744            {
1745              m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1746            }
1747            // for now, override the TILES_DECODER setting in order to write substreams.
1748            m_pcEntropyCoder->setBitstream    ( &pcSubstreamsOut[0] );
1749
1750          }
1751          pcSlice->setFinalized(true);
1752
1753          m_pcSbacCoder->load( &pcSbacCoders[0] );
1754
1755          pcSlice->setTileOffstForMultES( uiOneBitstreamPerSliceLength );
1756            pcSlice->setTileLocationCount ( 0 );
1757          m_pcSliceEncoder->encodeSlice(pcPic, pcSubstreamsOut);
1758
1759          {
1760            // Construct the final bitstream by flushing and concatenating substreams.
1761            // The final bitstream is either nalu.m_Bitstream or pcBitstreamRedirect;
1762            UInt* puiSubstreamSizes = pcSlice->getSubstreamSizes();
1763            UInt uiTotalCodedSize = 0; // for padding calcs.
1764            UInt uiNumSubstreamsPerTile = iNumSubstreams;
1765            if (iNumSubstreams > 1)
1766            {
1767              uiNumSubstreamsPerTile /= pcPic->getPicSym()->getNumTiles();
1768            }
1769            for ( UInt ui = 0 ; ui < iNumSubstreams; ui++ )
1770            {
1771              // Flush all substreams -- this includes empty ones.
1772              // Terminating bit and flush.
1773              m_pcEntropyCoder->setEntropyCoder   ( &pcSbacCoders[ui], pcSlice );
1774              m_pcEntropyCoder->setBitstream      (  &pcSubstreamsOut[ui] );
1775              m_pcEntropyCoder->encodeTerminatingBit( 1 );
1776              m_pcEntropyCoder->encodeSliceFinish();
1777
1778              pcSubstreamsOut[ui].writeByteAlignment();   // Byte-alignment in slice_data() at end of sub-stream
1779              // Byte alignment is necessary between tiles when tiles are independent.
1780              uiTotalCodedSize += pcSubstreamsOut[ui].getNumberOfWrittenBits();
1781
1782              Bool bNextSubstreamInNewTile = ((ui+1) < iNumSubstreams)&& ((ui+1)%uiNumSubstreamsPerTile == 0);
1783              if (bNextSubstreamInNewTile)
1784              {
1785                pcSlice->setTileLocation(ui/uiNumSubstreamsPerTile, pcSlice->getTileOffstForMultES()+(uiTotalCodedSize>>3));
1786              }
1787              if (ui+1 < pcSlice->getPPS()->getNumSubstreams())
1788              {
1789                puiSubstreamSizes[ui] = pcSubstreamsOut[ui].getNumberOfWrittenBits() + (pcSubstreamsOut[ui].countStartCodeEmulations()<<3);
1790              }
1791            }
1792
1793            // Complete the slice header info.
1794            m_pcEntropyCoder->setEntropyCoder   ( m_pcCavlcCoder, pcSlice );
1795            m_pcEntropyCoder->setBitstream(&nalu.m_Bitstream);
1796            m_pcEntropyCoder->encodeTilesWPPEntryPoint( pcSlice );
1797
1798            // Substreams...
1799            TComOutputBitstream *pcOut = pcBitstreamRedirect;
1800          Int offs = 0;
1801          Int nss = pcSlice->getPPS()->getNumSubstreams();
1802          if (pcSlice->getPPS()->getEntropyCodingSyncEnabledFlag())
1803          {
1804            // 1st line present for WPP.
1805            offs = pcSlice->getSliceSegmentCurStartCUAddr()/pcSlice->getPic()->getNumPartInCU()/pcSlice->getPic()->getFrameWidthInCU();
1806            nss  = pcSlice->getNumEntryPointOffsets()+1;
1807          }
1808          for ( UInt ui = 0 ; ui < nss; ui++ )
1809          {
1810            pcOut->addSubstream(&pcSubstreamsOut[ui+offs]);
1811            }
1812          }
1813
1814          UInt boundingAddrSlice, boundingAddrSliceSegment;
1815          boundingAddrSlice        = m_storedStartCUAddrForEncodingSlice[startCUAddrSliceIdx];         
1816          boundingAddrSliceSegment = m_storedStartCUAddrForEncodingSliceSegment[startCUAddrSliceSegmentIdx];         
1817          nextCUAddr               = min(boundingAddrSlice, boundingAddrSliceSegment);
1818          // If current NALU is the first NALU of slice (containing slice header) and more NALUs exist (due to multiple dependent slices) then buffer it.
1819          // If current NALU is the last NALU of slice and a NALU was buffered, then (a) Write current NALU (b) Update an write buffered NALU at approproate location in NALU list.
1820          Bool bNALUAlignedWrittenToList    = false; // used to ensure current NALU is not written more than once to the NALU list.
1821          xAttachSliceDataToNalUnit(nalu, pcBitstreamRedirect);
1822          accessUnit.push_back(new NALUnitEBSP(nalu));
1823#if RATE_CONTROL_LAMBDA_DOMAIN
1824          actualTotalBits += UInt(accessUnit.back()->m_nalUnitData.str().size()) * 8;
1825#endif
1826          bNALUAlignedWrittenToList = true; 
1827          uiOneBitstreamPerSliceLength += nalu.m_Bitstream.getNumberOfWrittenBits(); // length of bitstream after byte-alignment
1828
1829          if (!bNALUAlignedWrittenToList)
1830          {
1831            {
1832              nalu.m_Bitstream.writeAlignZero();
1833            }
1834            accessUnit.push_back(new NALUnitEBSP(nalu));
1835            uiOneBitstreamPerSliceLength += nalu.m_Bitstream.getNumberOfWrittenBits() + 24; // length of bitstream after byte-alignment + 3 byte startcode 0x000001
1836          }
1837
1838          if( ( m_pcCfg->getPictureTimingSEIEnabled() || m_pcCfg->getDecodingUnitInfoSEIEnabled() ) &&
1839              ( pcSlice->getSPS()->getVuiParametersPresentFlag() ) &&
1840              ( ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getNalHrdParametersPresentFlag() ) 
1841             || ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getVclHrdParametersPresentFlag() ) ) &&
1842              ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getSubPicCpbParamsPresentFlag() ) )
1843          {
1844              UInt numNalus = 0;
1845            UInt numRBSPBytes = 0;
1846            for (AccessUnit::const_iterator it = accessUnit.begin(); it != accessUnit.end(); it++)
1847            {
1848              UInt numRBSPBytes_nal = UInt((*it)->m_nalUnitData.str().size());
1849              if ((*it)->m_nalUnitType != NAL_UNIT_PREFIX_SEI && (*it)->m_nalUnitType != NAL_UNIT_SUFFIX_SEI)
1850              {
1851                numRBSPBytes += numRBSPBytes_nal;
1852                numNalus ++;
1853              }
1854            }
1855            accumBitsDU[ pcSlice->getSliceIdx() ] = ( numRBSPBytes << 3 );
1856            accumNalsDU[ pcSlice->getSliceIdx() ] = numNalus;   // SEI not counted for bit count; hence shouldn't be counted for # of NALUs - only for consistency
1857          }
1858          processingState = ENCODE_SLICE;
1859          }
1860          break;
1861        case EXECUTE_INLOOPFILTER:
1862          {
1863            // set entropy coder for RD
1864            m_pcEntropyCoder->setEntropyCoder ( m_pcSbacCoder, pcSlice );
1865            if ( pcSlice->getSPS()->getUseSAO() )
1866            {
1867              m_pcEntropyCoder->resetEntropy();
1868              m_pcEntropyCoder->setBitstream( m_pcBitCounter );
1869              m_pcSAO->startSaoEnc(pcPic, m_pcEntropyCoder, m_pcEncTop->getRDSbacCoder(), m_pcEncTop->getRDGoOnSbacCoder());
1870              SAOParam& cSaoParam = *pcSlice->getPic()->getPicSym()->getSaoParam();
1871
1872#if SAO_CHROMA_LAMBDA
1873#if SAO_ENCODING_CHOICE
1874              m_pcSAO->SAOProcess(&cSaoParam, pcPic->getSlice(0)->getLambdaLuma(), pcPic->getSlice(0)->getLambdaChroma(), pcPic->getSlice(0)->getDepth());
1875#else
1876              m_pcSAO->SAOProcess(&cSaoParam, pcPic->getSlice(0)->getLambdaLuma(), pcPic->getSlice(0)->getLambdaChroma());
1877#endif
1878#else
1879              m_pcSAO->SAOProcess(&cSaoParam, pcPic->getSlice(0)->getLambda());
1880#endif
1881              m_pcSAO->endSaoEnc();
1882              m_pcSAO->PCMLFDisableProcess(pcPic);
1883            }
1884#if SAO_RDO
1885            m_pcEntropyCoder->setEntropyCoder ( m_pcCavlcCoder, pcSlice );
1886#endif
1887            processingState = ENCODE_SLICE;
1888
1889            for(Int s=0; s< uiNumSlices; s++)
1890            {
1891              if (pcSlice->getSPS()->getUseSAO())
1892              {
1893                pcPic->getSlice(s)->setSaoEnabledFlag((pcSlice->getPic()->getPicSym()->getSaoParam()->bSaoFlag[0]==1)?true:false);
1894              }
1895            }
1896          }
1897          break;
1898        default:
1899          {
1900            printf("Not a supported encoding state\n");
1901            assert(0);
1902            exit(-1);
1903          }
1904        }
1905      } // end iteration over slices
1906
1907      if(pcSlice->getSPS()->getUseSAO())
1908      {
1909        if(pcSlice->getSPS()->getUseSAO())
1910        {
1911          m_pcSAO->destroyPicSaoInfo();
1912        }
1913        pcPic->destroyNonDBFilterInfo();
1914      }
1915
1916      pcPic->compressMotion(); 
1917     
1918      //-- For time output for each slice
1919      Double dEncTime = (Double)(clock()-iBeforeTime) / CLOCKS_PER_SEC;
1920
1921      const Char* digestStr = NULL;
1922      if (m_pcCfg->getDecodedPictureHashSEIEnabled())
1923      {
1924        /* calculate MD5sum for entire reconstructed picture */
1925        SEIDecodedPictureHash sei_recon_picture_digest;
1926        if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 1)
1927        {
1928          sei_recon_picture_digest.method = SEIDecodedPictureHash::MD5;
1929          calcMD5(*pcPic->getPicYuvRec(), sei_recon_picture_digest.digest);
1930          digestStr = digestToString(sei_recon_picture_digest.digest, 16);
1931        }
1932        else if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 2)
1933        {
1934          sei_recon_picture_digest.method = SEIDecodedPictureHash::CRC;
1935          calcCRC(*pcPic->getPicYuvRec(), sei_recon_picture_digest.digest);
1936          digestStr = digestToString(sei_recon_picture_digest.digest, 2);
1937        }
1938        else if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 3)
1939        {
1940          sei_recon_picture_digest.method = SEIDecodedPictureHash::CHECKSUM;
1941          calcChecksum(*pcPic->getPicYuvRec(), sei_recon_picture_digest.digest);
1942          digestStr = digestToString(sei_recon_picture_digest.digest, 4);
1943        }
1944#if SVC_EXTENSION
1945        OutputNALUnit nalu(NAL_UNIT_SUFFIX_SEI, pcSlice->getTLayer(), m_layerId);
1946#else
1947        OutputNALUnit nalu(NAL_UNIT_SUFFIX_SEI, pcSlice->getTLayer());
1948#endif
1949
1950        /* write the SEI messages */
1951        m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1952        m_seiWriter.writeSEImessage(nalu.m_Bitstream, sei_recon_picture_digest, pcSlice->getSPS());
1953        writeRBSPTrailingBits(nalu.m_Bitstream);
1954
1955        accessUnit.insert(accessUnit.end(), new NALUnitEBSP(nalu));
1956      }
1957      if (m_pcCfg->getTemporalLevel0IndexSEIEnabled())
1958      {
1959        SEITemporalLevel0Index sei_temporal_level0_index;
1960        if (pcSlice->getRapPicFlag())
1961        {
1962          m_tl0Idx = 0;
1963          m_rapIdx = (m_rapIdx + 1) & 0xFF;
1964        }
1965        else
1966        {
1967          m_tl0Idx = (m_tl0Idx + (pcSlice->getTLayer() ? 0 : 1)) & 0xFF;
1968        }
1969        sei_temporal_level0_index.tl0Idx = m_tl0Idx;
1970        sei_temporal_level0_index.rapIdx = m_rapIdx;
1971
1972        OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI); 
1973
1974        /* write the SEI messages */
1975        m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
1976        m_seiWriter.writeSEImessage(nalu.m_Bitstream, sei_temporal_level0_index, pcSlice->getSPS());
1977        writeRBSPTrailingBits(nalu.m_Bitstream);
1978
1979        /* insert the SEI message NALUnit before any Slice NALUnits */
1980        AccessUnit::iterator it = find_if(accessUnit.begin(), accessUnit.end(), mem_fun(&NALUnit::isSlice));
1981        accessUnit.insert(it, new NALUnitEBSP(nalu));
1982      }
1983
1984      xCalculateAddPSNR( pcPic, pcPic->getPicYuvRec(), accessUnit, dEncTime );
1985
1986      if (digestStr)
1987      {
1988        if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 1)
1989        {
1990          printf(" [MD5:%s]", digestStr);
1991        }
1992        else if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 2)
1993        {
1994          printf(" [CRC:%s]", digestStr);
1995        }
1996        else if(m_pcCfg->getDecodedPictureHashSEIEnabled() == 3)
1997        {
1998          printf(" [Checksum:%s]", digestStr);
1999        }
2000      }
2001#if RATE_CONTROL_LAMBDA_DOMAIN
2002      if ( m_pcCfg->getUseRateCtrl() )
2003      {
2004        Double effectivePercentage = m_pcRateCtrl->getRCPic()->getEffectivePercentage();
2005        Double avgQP     = m_pcRateCtrl->getRCPic()->calAverageQP();
2006        Double avgLambda = m_pcRateCtrl->getRCPic()->calAverageLambda();
2007        if ( avgLambda < 0.0 )
2008        {
2009          avgLambda = lambda;
2010        }
2011        m_pcRateCtrl->getRCPic()->updateAfterPicture( actualHeadBits, actualTotalBits, avgQP, avgLambda, effectivePercentage );
2012        m_pcRateCtrl->getRCPic()->addToPictureLsit( m_pcRateCtrl->getPicList() );
2013
2014        m_pcRateCtrl->getRCSeq()->updateAfterPic( actualTotalBits );
2015        if ( pcSlice->getSliceType() != I_SLICE )
2016        {
2017          m_pcRateCtrl->getRCGOP()->updateAfterPicture( actualTotalBits );
2018        }
2019        else    // for intra picture, the estimated bits are used to update the current status in the GOP
2020        {
2021          m_pcRateCtrl->getRCGOP()->updateAfterPicture( estimatedBits );
2022        }
2023      }
2024#else
2025      if(m_pcCfg->getUseRateCtrl())
2026      {
2027        UInt  frameBits = m_vRVM_RP[m_vRVM_RP.size()-1];
2028        m_pcRateCtrl->updataRCFrameStatus((Int)frameBits, pcSlice->getSliceType());
2029      }
2030#endif
2031      if( ( m_pcCfg->getPictureTimingSEIEnabled() || m_pcCfg->getDecodingUnitInfoSEIEnabled() ) &&
2032          ( pcSlice->getSPS()->getVuiParametersPresentFlag() ) &&
2033          ( ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getNalHrdParametersPresentFlag() ) 
2034         || ( pcSlice->getSPS()->getVuiParameters()->getHrdParameters()->getVclHrdParametersPresentFlag() ) ) )
2035      {
2036        TComVUI *vui = pcSlice->getSPS()->getVuiParameters();
2037        TComHRD *hrd = vui->getHrdParameters();
2038
2039        if( hrd->getSubPicCpbParamsPresentFlag() )
2040        {
2041          Int i;
2042          UInt64 ui64Tmp;
2043          UInt uiPrev = 0;
2044          UInt numDU = ( pictureTimingSEI.m_numDecodingUnitsMinus1 + 1 );
2045          UInt *pCRD = &pictureTimingSEI.m_duCpbRemovalDelayMinus1[0];
2046          UInt maxDiff = ( hrd->getTickDivisorMinus2() + 2 ) - 1;
2047
2048          for( i = 0; i < numDU; i ++ )
2049          {
2050            pictureTimingSEI.m_numNalusInDuMinus1[ i ]       = ( i == 0 ) ? ( accumNalsDU[ i ] - 1 ) : ( accumNalsDU[ i ] - accumNalsDU[ i - 1] - 1 );
2051          }
2052
2053          if( numDU == 1 )
2054          {
2055            pCRD[ 0 ] = 0; /* don't care */
2056          }
2057          else
2058          {
2059            pCRD[ numDU - 1 ] = 0;/* by definition */
2060            UInt tmp = 0;
2061            UInt accum = 0;
2062
2063            for( i = ( numDU - 2 ); i >= 0; i -- )
2064            {
2065#if L0043_TIMING_INFO
2066              ui64Tmp = ( ( ( accumBitsDU[ numDU - 1 ]  - accumBitsDU[ i ] ) * ( vui->getTimingInfo()->getTimeScale() / vui->getTimingInfo()->getNumUnitsInTick() ) * ( hrd->getTickDivisorMinus2() + 2 ) ) / ( m_pcCfg->getTargetBitrate() ) );
2067#else
2068              ui64Tmp = ( ( ( accumBitsDU[ numDU - 1 ]  - accumBitsDU[ i ] ) * ( hrd->getTimeScale() / hrd->getNumUnitsInTick() ) * ( hrd->getTickDivisorMinus2() + 2 ) ) / ( m_pcCfg->getTargetBitrate() ) );
2069#endif
2070              if( (UInt)ui64Tmp > maxDiff )
2071              {
2072                tmp ++;
2073              }
2074            }
2075            uiPrev = 0;
2076
2077            UInt flag = 0;
2078            for( i = ( numDU - 2 ); i >= 0; i -- )
2079            {
2080              flag = 0;
2081#if L0043_TIMING_INFO
2082              ui64Tmp = ( ( ( accumBitsDU[ numDU - 1 ]  - accumBitsDU[ i ] ) * ( vui->getTimingInfo()->getTimeScale() / vui->getTimingInfo()->getNumUnitsInTick() ) * ( hrd->getTickDivisorMinus2() + 2 ) ) / ( m_pcCfg->getTargetBitrate() ) );
2083#else
2084              ui64Tmp = ( ( ( accumBitsDU[ numDU - 1 ]  - accumBitsDU[ i ] ) * ( hrd->getTimeScale() / hrd->getNumUnitsInTick() ) * ( hrd->getTickDivisorMinus2() + 2 ) ) / ( m_pcCfg->getTargetBitrate() ) );
2085#endif
2086
2087              if( (UInt)ui64Tmp > maxDiff )
2088              {
2089                if(uiPrev >= maxDiff - tmp)
2090                {
2091                  ui64Tmp = uiPrev + 1;
2092                  flag = 1;
2093                }
2094                else                            ui64Tmp = maxDiff - tmp + 1;
2095              }
2096              pCRD[ i ] = (UInt)ui64Tmp - uiPrev - 1;
2097              if( (Int)pCRD[ i ] < 0 )
2098              {
2099                pCRD[ i ] = 0;
2100              }
2101              else if (tmp > 0 && flag == 1) 
2102              {
2103                tmp --;
2104              }
2105              accum += pCRD[ i ] + 1;
2106              uiPrev = accum;
2107            }
2108          }
2109        }
2110        if( m_pcCfg->getPictureTimingSEIEnabled() )
2111        {
2112          {
2113            OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI, pcSlice->getTLayer());
2114          m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
2115          m_seiWriter.writeSEImessage(nalu.m_Bitstream, pictureTimingSEI, pcSlice->getSPS());
2116          writeRBSPTrailingBits(nalu.m_Bitstream);
2117#if L0045_NON_NESTED_SEI_RESTRICTIONS
2118          UInt seiPositionInAu = xGetFirstSeiLocation(accessUnit);
2119          UInt offsetPosition = m_activeParameterSetSEIPresentInAU
2120                                    + m_bufferingPeriodSEIPresentInAU;    // Insert PT SEI after APS and BP SEI
2121          AccessUnit::iterator it;
2122          for(j = 0, it = accessUnit.begin(); j < seiPositionInAu + offsetPosition; j++)
2123          {
2124            it++;
2125          }
2126          accessUnit.insert(it, new NALUnitEBSP(nalu));
2127          m_pictureTimingSEIPresentInAU = true;
2128#else
2129          AccessUnit::iterator it = find_if(accessUnit.begin(), accessUnit.end(), mem_fun(&NALUnit::isSlice));
2130          accessUnit.insert(it, new NALUnitEBSP(nalu));
2131#endif
2132        }
2133#if K0180_SCALABLE_NESTING_SEI
2134          if ( m_pcCfg->getScalableNestingSEIEnabled() ) // put picture timing SEI into scalable nesting SEI
2135          {
2136            OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI, pcSlice->getTLayer());
2137            m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
2138            scalableNestingSEI.m_nestedSEIs.clear();
2139            scalableNestingSEI.m_nestedSEIs.push_back(&pictureTimingSEI);
2140            m_seiWriter.writeSEImessage(nalu.m_Bitstream, scalableNestingSEI, pcSlice->getSPS());
2141            writeRBSPTrailingBits(nalu.m_Bitstream);
2142#if L0045_NON_NESTED_SEI_RESTRICTIONS
2143            UInt seiPositionInAu = xGetFirstSeiLocation(accessUnit);
2144            UInt offsetPosition = m_activeParameterSetSEIPresentInAU
2145              + m_bufferingPeriodSEIPresentInAU + m_pictureTimingSEIPresentInAU + m_nestedBufferingPeriodSEIPresentInAU;    // Insert PT SEI after APS and BP SEI
2146            AccessUnit::iterator it;
2147            for(j = 0, it = accessUnit.begin(); j < seiPositionInAu + offsetPosition; j++)
2148            {
2149              it++;
2150            }
2151            accessUnit.insert(it, new NALUnitEBSP(nalu));
2152            m_nestedPictureTimingSEIPresentInAU = true;
2153#else
2154            AccessUnit::iterator it = find_if(accessUnit.begin(), accessUnit.end(), mem_fun(&NALUnit::isSlice));
2155            accessUnit.insert(it, new NALUnitEBSP(nalu));
2156#endif
2157          }
2158#endif
2159
2160        }
2161        if( m_pcCfg->getDecodingUnitInfoSEIEnabled() && hrd->getSubPicCpbParamsPresentFlag() )
2162        {             
2163          m_pcEntropyCoder->setEntropyCoder(m_pcCavlcCoder, pcSlice);
2164          for( Int i = 0; i < ( pictureTimingSEI.m_numDecodingUnitsMinus1 + 1 ); i ++ )
2165          {
2166            OutputNALUnit nalu(NAL_UNIT_PREFIX_SEI, pcSlice->getTLayer());
2167
2168            SEIDecodingUnitInfo tempSEI;
2169            tempSEI.m_decodingUnitIdx = i;
2170            tempSEI.m_duSptCpbRemovalDelay = pictureTimingSEI.m_duCpbRemovalDelayMinus1[i] + 1;
2171#if L0044_DU_DPB_OUTPUT_DELAY_HRD
2172            tempSEI.m_dpbOutputDuDelayPresentFlag = false;
2173            tempSEI.m_picSptDpbOutputDuDelay = picSptDpbOutputDuDelay;
2174#endif
2175
2176            AccessUnit::iterator it;
2177            // Insert the first one in the right location, before the first slice
2178            if(i == 0)
2179            {
2180              // Insert before the first slice.
2181              m_seiWriter.writeSEImessage(nalu.m_Bitstream, tempSEI, pcSlice->getSPS());
2182              writeRBSPTrailingBits(nalu.m_Bitstream);
2183
2184#if L0045_NON_NESTED_SEI_RESTRICTIONS
2185              UInt seiPositionInAu = xGetFirstSeiLocation(accessUnit);
2186              UInt offsetPosition = m_activeParameterSetSEIPresentInAU
2187                                    + m_bufferingPeriodSEIPresentInAU
2188                                    + m_pictureTimingSEIPresentInAU;  // Insert DU info SEI after APS, BP and PT SEI
2189              for(j = 0, it = accessUnit.begin(); j < seiPositionInAu + offsetPosition; j++)
2190              {
2191                it++;
2192              }
2193              accessUnit.insert(it, new NALUnitEBSP(nalu));
2194#else
2195              it = find_if(accessUnit.begin(), accessUnit.end(), mem_fun(&NALUnit::isSlice));
2196              accessUnit.insert(it, new NALUnitEBSP(nalu)); 
2197#endif
2198            }
2199            else
2200            {
2201              Int ctr;
2202              // For the second decoding unit onwards we know how many NALUs are present
2203              for (ctr = 0, it = accessUnit.begin(); it != accessUnit.end(); it++)
2204              {           
2205                if(ctr == accumNalsDU[ i - 1 ])
2206                {
2207                  // Insert before the first slice.
2208                  m_seiWriter.writeSEImessage(nalu.m_Bitstream, tempSEI, pcSlice->getSPS());
2209                  writeRBSPTrailingBits(nalu.m_Bitstream);
2210
2211                  accessUnit.insert(it, new NALUnitEBSP(nalu));
2212                  break;
2213                }
2214                if ((*it)->m_nalUnitType != NAL_UNIT_PREFIX_SEI && (*it)->m_nalUnitType != NAL_UNIT_SUFFIX_SEI)
2215                {
2216                  ctr++;
2217                }
2218              }
2219            }           
2220          }
2221        }
2222      }
2223#if L0045_NON_NESTED_SEI_RESTRICTIONS
2224      xResetNonNestedSEIPresentFlags();
2225#if K0180_SCALABLE_NESTING_SEI
2226      xResetNestedSEIPresentFlags();
2227#endif
2228#endif
2229      pcPic->getPicYuvRec()->copyToPic(pcPicYuvRecOut);
2230
2231      pcPic->setReconMark   ( true );
2232      m_bFirst = false;
2233      m_iNumPicCoded++;
2234      m_totalCoded ++;
2235      /* logging: insert a newline at end of picture period */
2236      printf("\n");
2237      fflush(stdout);
2238
2239      delete[] pcSubstreamsOut;
2240  }
2241#if !RATE_CONTROL_LAMBDA_DOMAIN
2242  if(m_pcCfg->getUseRateCtrl())
2243  {
2244    m_pcRateCtrl->updateRCGOPStatus();
2245  }
2246#endif
2247  delete pcBitstreamRedirect;
2248
2249  if( accumBitsDU != NULL) delete accumBitsDU;
2250  if( accumNalsDU != NULL) delete accumNalsDU;
2251
2252#if SVC_EXTENSION
2253  assert ( m_iNumPicCoded <= 1 );
2254#else
2255  assert ( m_iNumPicCoded == iNumPicRcvd );
2256#endif
2257}
2258
2259#if !SVC_EXTENSION
2260Void TEncGOP::printOutSummary(UInt uiNumAllPicCoded)
2261{
2262  assert (uiNumAllPicCoded == m_gcAnalyzeAll.getNumPic());
2263 
2264   
2265  //--CFG_KDY
2266  m_gcAnalyzeAll.setFrmRate( m_pcCfg->getFrameRate() );
2267  m_gcAnalyzeI.setFrmRate( m_pcCfg->getFrameRate() );
2268  m_gcAnalyzeP.setFrmRate( m_pcCfg->getFrameRate() );
2269  m_gcAnalyzeB.setFrmRate( m_pcCfg->getFrameRate() );
2270 
2271  //-- all
2272  printf( "\n\nSUMMARY --------------------------------------------------------\n" );
2273  m_gcAnalyzeAll.printOut('a');
2274 
2275  printf( "\n\nI Slices--------------------------------------------------------\n" );
2276  m_gcAnalyzeI.printOut('i');
2277 
2278  printf( "\n\nP Slices--------------------------------------------------------\n" );
2279  m_gcAnalyzeP.printOut('p');
2280 
2281  printf( "\n\nB Slices--------------------------------------------------------\n" );
2282  m_gcAnalyzeB.printOut('b');
2283 
2284#if _SUMMARY_OUT_
2285  m_gcAnalyzeAll.printSummaryOut();
2286#endif
2287#if _SUMMARY_PIC_
2288  m_gcAnalyzeI.printSummary('I');
2289  m_gcAnalyzeP.printSummary('P');
2290  m_gcAnalyzeB.printSummary('B');
2291#endif
2292
2293  printf("\nRVM: %.3lf\n" , xCalculateRVM());
2294}
2295#endif
2296
2297Void TEncGOP::preLoopFilterPicAll( TComPic* pcPic, UInt64& ruiDist, UInt64& ruiBits )
2298{
2299  TComSlice* pcSlice = pcPic->getSlice(pcPic->getCurrSliceIdx());
2300  Bool bCalcDist = false;
2301  m_pcLoopFilter->setCfg(m_pcCfg->getLFCrossTileBoundaryFlag());
2302  m_pcLoopFilter->loopFilterPic( pcPic );
2303 
2304  m_pcEntropyCoder->setEntropyCoder ( m_pcEncTop->getRDGoOnSbacCoder(), pcSlice );
2305  m_pcEntropyCoder->resetEntropy    ();
2306  m_pcEntropyCoder->setBitstream    ( m_pcBitCounter );
2307  pcSlice = pcPic->getSlice(0);
2308  if(pcSlice->getSPS()->getUseSAO())
2309  {
2310    std::vector<Bool> LFCrossSliceBoundaryFlag(1, true);
2311    std::vector<Int>  sliceStartAddress;
2312    sliceStartAddress.push_back(0);
2313    sliceStartAddress.push_back(pcPic->getNumCUsInFrame()* pcPic->getNumPartInCU());
2314    pcPic->createNonDBFilterInfo(sliceStartAddress, 0, &LFCrossSliceBoundaryFlag);
2315  }
2316 
2317  if( pcSlice->getSPS()->getUseSAO())
2318  {
2319    pcPic->destroyNonDBFilterInfo();
2320  }
2321 
2322  m_pcEntropyCoder->resetEntropy    ();
2323  ruiBits += m_pcEntropyCoder->getNumberOfWrittenBits();
2324 
2325  if (!bCalcDist)
2326    ruiDist = xFindDistortionFrame(pcPic->getPicYuvOrg(), pcPic->getPicYuvRec());
2327}
2328
2329// ====================================================================================================================
2330// Protected member functions
2331// ====================================================================================================================
2332
2333Void TEncGOP::xInitGOP( Int iPOCLast, Int iNumPicRcvd, TComList<TComPic*>& rcListPic, TComList<TComPicYuv*>& rcListPicYuvRecOut )
2334{
2335  assert( iNumPicRcvd > 0 );
2336  //  Exception for the first frame
2337  if ( iPOCLast == 0 )
2338  {
2339    m_iGopSize    = 1;
2340  }
2341  else
2342    m_iGopSize    = m_pcCfg->getGOPSize();
2343 
2344  assert (m_iGopSize > 0); 
2345
2346  return;
2347}
2348
2349Void TEncGOP::xGetBuffer( TComList<TComPic*>&      rcListPic,
2350                         TComList<TComPicYuv*>&    rcListPicYuvRecOut,
2351                         Int                       iNumPicRcvd,
2352                         Int                       iTimeOffset,
2353                         TComPic*&                 rpcPic,
2354                         TComPicYuv*&              rpcPicYuvRecOut,
2355                         Int                       pocCurr )
2356{
2357  Int i;
2358  //  Rec. output
2359  TComList<TComPicYuv*>::iterator     iterPicYuvRec = rcListPicYuvRecOut.end();
2360  for ( i = 0; i < iNumPicRcvd - iTimeOffset + 1; i++ )
2361  {
2362    iterPicYuvRec--;
2363  }
2364 
2365  rpcPicYuvRecOut = *(iterPicYuvRec);
2366 
2367  //  Current pic.
2368  TComList<TComPic*>::iterator        iterPic       = rcListPic.begin();
2369  while (iterPic != rcListPic.end())
2370  {
2371    rpcPic = *(iterPic);
2372    rpcPic->setCurrSliceIdx(0);
2373    if (rpcPic->getPOC() == pocCurr)
2374    {
2375      break;
2376    }
2377    iterPic++;
2378  }
2379 
2380  assert (rpcPic->getPOC() == pocCurr);
2381 
2382  return;
2383}
2384
2385UInt64 TEncGOP::xFindDistortionFrame (TComPicYuv* pcPic0, TComPicYuv* pcPic1)
2386{
2387  Int     x, y;
2388  Pel*  pSrc0   = pcPic0 ->getLumaAddr();
2389  Pel*  pSrc1   = pcPic1 ->getLumaAddr();
2390  UInt  uiShift = 2 * DISTORTION_PRECISION_ADJUSTMENT(g_bitDepthY-8);
2391  Int   iTemp;
2392 
2393  Int   iStride = pcPic0->getStride();
2394  Int   iWidth  = pcPic0->getWidth();
2395  Int   iHeight = pcPic0->getHeight();
2396 
2397  UInt64  uiTotalDiff = 0;
2398 
2399  for( y = 0; y < iHeight; y++ )
2400  {
2401    for( x = 0; x < iWidth; x++ )
2402    {
2403      iTemp = pSrc0[x] - pSrc1[x]; uiTotalDiff += (iTemp*iTemp) >> uiShift;
2404    }
2405    pSrc0 += iStride;
2406    pSrc1 += iStride;
2407  }
2408 
2409  uiShift = 2 * DISTORTION_PRECISION_ADJUSTMENT(g_bitDepthC-8);
2410  iHeight >>= 1;
2411  iWidth  >>= 1;
2412  iStride >>= 1;
2413 
2414  pSrc0  = pcPic0->getCbAddr();
2415  pSrc1  = pcPic1->getCbAddr();
2416 
2417  for( y = 0; y < iHeight; y++ )
2418  {
2419    for( x = 0; x < iWidth; x++ )
2420    {
2421      iTemp = pSrc0[x] - pSrc1[x]; uiTotalDiff += (iTemp*iTemp) >> uiShift;
2422    }
2423    pSrc0 += iStride;
2424    pSrc1 += iStride;
2425  }
2426 
2427  pSrc0  = pcPic0->getCrAddr();
2428  pSrc1  = pcPic1->getCrAddr();
2429 
2430  for( y = 0; y < iHeight; y++ )
2431  {
2432    for( x = 0; x < iWidth; x++ )
2433    {
2434      iTemp = pSrc0[x] - pSrc1[x]; uiTotalDiff += (iTemp*iTemp) >> uiShift;
2435    }
2436    pSrc0 += iStride;
2437    pSrc1 += iStride;
2438  }
2439 
2440  return uiTotalDiff;
2441}
2442
2443#if VERBOSE_RATE
2444static const Char* nalUnitTypeToString(NalUnitType type)
2445{
2446  switch (type)
2447  {
2448    case NAL_UNIT_CODED_SLICE_TRAIL_R: return "TRAIL_R";
2449    case NAL_UNIT_CODED_SLICE_TRAIL_N: return "TRAIL_N";
2450    case NAL_UNIT_CODED_SLICE_TLA_R:      return "TLA_R";
2451    case NAL_UNIT_CODED_SLICE_TSA_N: return "TSA_N";
2452    case NAL_UNIT_CODED_SLICE_STSA_R: return "STSA_R";
2453    case NAL_UNIT_CODED_SLICE_STSA_N: return "STSA_N";
2454    case NAL_UNIT_CODED_SLICE_BLA_W_LP:   return "BLA_W_LP";
2455    case NAL_UNIT_CODED_SLICE_BLA_W_RADL: return "BLA_W_RADL";
2456    case NAL_UNIT_CODED_SLICE_BLA_N_LP: return "BLA_N_LP";
2457    case NAL_UNIT_CODED_SLICE_IDR_W_RADL: return "IDR_W_RADL";
2458    case NAL_UNIT_CODED_SLICE_IDR_N_LP: return "IDR_N_LP";
2459    case NAL_UNIT_CODED_SLICE_CRA: return "CRA";
2460    case NAL_UNIT_CODED_SLICE_RADL_R:     return "RADL_R";
2461    case NAL_UNIT_CODED_SLICE_RASL_R:     return "RASL_R";
2462    case NAL_UNIT_VPS: return "VPS";
2463    case NAL_UNIT_SPS: return "SPS";
2464    case NAL_UNIT_PPS: return "PPS";
2465    case NAL_UNIT_ACCESS_UNIT_DELIMITER: return "AUD";
2466    case NAL_UNIT_EOS: return "EOS";
2467    case NAL_UNIT_EOB: return "EOB";
2468    case NAL_UNIT_FILLER_DATA: return "FILLER";
2469    case NAL_UNIT_PREFIX_SEI:             return "SEI";
2470    case NAL_UNIT_SUFFIX_SEI:             return "SEI";
2471    default: return "UNK";
2472  }
2473}
2474#endif
2475
2476Void TEncGOP::xCalculateAddPSNR( TComPic* pcPic, TComPicYuv* pcPicD, const AccessUnit& accessUnit, Double dEncTime )
2477{
2478  Int     x, y;
2479  UInt64 uiSSDY  = 0;
2480  UInt64 uiSSDU  = 0;
2481  UInt64 uiSSDV  = 0;
2482 
2483  Double  dYPSNR  = 0.0;
2484  Double  dUPSNR  = 0.0;
2485  Double  dVPSNR  = 0.0;
2486 
2487  //===== calculate PSNR =====
2488  Pel*  pOrg    = pcPic ->getPicYuvOrg()->getLumaAddr();
2489  Pel*  pRec    = pcPicD->getLumaAddr();
2490  Int   iStride = pcPicD->getStride();
2491 
2492  Int   iWidth;
2493  Int   iHeight;
2494 
2495  iWidth  = pcPicD->getWidth () - m_pcEncTop->getPad(0);
2496  iHeight = pcPicD->getHeight() - m_pcEncTop->getPad(1);
2497 
2498  Int   iSize   = iWidth*iHeight;
2499 
2500  for( y = 0; y < iHeight; y++ )
2501  {
2502    for( x = 0; x < iWidth; x++ )
2503    {
2504      Int iDiff = (Int)( pOrg[x] - pRec[x] );
2505      uiSSDY   += iDiff * iDiff;
2506    }
2507    pOrg += iStride;
2508    pRec += iStride;
2509  }
2510 
2511  iHeight >>= 1;
2512  iWidth  >>= 1;
2513  iStride >>= 1;
2514  pOrg  = pcPic ->getPicYuvOrg()->getCbAddr();
2515  pRec  = pcPicD->getCbAddr();
2516 
2517  for( y = 0; y < iHeight; y++ )
2518  {
2519    for( x = 0; x < iWidth; x++ )
2520    {
2521      Int iDiff = (Int)( pOrg[x] - pRec[x] );
2522      uiSSDU   += iDiff * iDiff;
2523    }
2524    pOrg += iStride;
2525    pRec += iStride;
2526  }
2527 
2528  pOrg  = pcPic ->getPicYuvOrg()->getCrAddr();
2529  pRec  = pcPicD->getCrAddr();
2530 
2531  for( y = 0; y < iHeight; y++ )
2532  {
2533    for( x = 0; x < iWidth; x++ )
2534    {
2535      Int iDiff = (Int)( pOrg[x] - pRec[x] );
2536      uiSSDV   += iDiff * iDiff;
2537    }
2538    pOrg += iStride;
2539    pRec += iStride;
2540  }
2541 
2542  Int maxvalY = 255 << (g_bitDepthY-8);
2543  Int maxvalC = 255 << (g_bitDepthC-8);
2544  Double fRefValueY = (Double) maxvalY * maxvalY * iSize;
2545  Double fRefValueC = (Double) maxvalC * maxvalC * iSize / 4.0;
2546  dYPSNR            = ( uiSSDY ? 10.0 * log10( fRefValueY / (Double)uiSSDY ) : 99.99 );
2547  dUPSNR            = ( uiSSDU ? 10.0 * log10( fRefValueC / (Double)uiSSDU ) : 99.99 );
2548  dVPSNR            = ( uiSSDV ? 10.0 * log10( fRefValueC / (Double)uiSSDV ) : 99.99 );
2549
2550  /* calculate the size of the access unit, excluding:
2551   *  - any AnnexB contributions (start_code_prefix, zero_byte, etc.,)
2552   *  - SEI NAL units
2553   */
2554  UInt numRBSPBytes = 0;
2555  for (AccessUnit::const_iterator it = accessUnit.begin(); it != accessUnit.end(); it++)
2556  {
2557    UInt numRBSPBytes_nal = UInt((*it)->m_nalUnitData.str().size());
2558#if VERBOSE_RATE
2559    printf("*** %6s numBytesInNALunit: %u\n", nalUnitTypeToString((*it)->m_nalUnitType), numRBSPBytes_nal);
2560#endif
2561    if ((*it)->m_nalUnitType != NAL_UNIT_PREFIX_SEI && (*it)->m_nalUnitType != NAL_UNIT_SUFFIX_SEI)
2562    {
2563      numRBSPBytes += numRBSPBytes_nal;
2564    }
2565  }
2566
2567  UInt uibits = numRBSPBytes * 8;
2568  m_vRVM_RP.push_back( uibits );
2569
2570  //===== add PSNR =====
2571#if SVC_EXTENSION
2572  m_gcAnalyzeAll[m_layerId].addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
2573  TComSlice*  pcSlice = pcPic->getSlice(0);
2574  if (pcSlice->isIntra())
2575  {
2576    m_gcAnalyzeI[m_layerId].addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
2577  }
2578  if (pcSlice->isInterP())
2579  {
2580    m_gcAnalyzeP[m_layerId].addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
2581  }
2582  if (pcSlice->isInterB())
2583  {
2584    m_gcAnalyzeB[m_layerId].addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
2585  }
2586#else
2587  m_gcAnalyzeAll.addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
2588  TComSlice*  pcSlice = pcPic->getSlice(0);
2589  if (pcSlice->isIntra())
2590  {
2591    m_gcAnalyzeI.addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
2592  }
2593  if (pcSlice->isInterP())
2594  {
2595    m_gcAnalyzeP.addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
2596  }
2597  if (pcSlice->isInterB())
2598  {
2599    m_gcAnalyzeB.addResult (dYPSNR, dUPSNR, dVPSNR, (Double)uibits);
2600  }
2601#endif
2602
2603  Char c = (pcSlice->isIntra() ? 'I' : pcSlice->isInterP() ? 'P' : 'B');
2604  if (!pcSlice->isReferenced()) c += 32;
2605
2606#if SVC_EXTENSION
2607#if ADAPTIVE_QP_SELECTION
2608  printf("POC %4d LId: %1d TId: %1d ( %c-SLICE, nQP %d QP %d ) %10d bits",
2609         pcSlice->getPOC(),
2610         pcSlice->getLayerId(),
2611         pcSlice->getTLayer(),
2612         c,
2613         pcSlice->getSliceQpBase(),
2614         pcSlice->getSliceQp(),
2615         uibits );
2616#else
2617  printf("POC %4d LId: %1d TId: %1d ( %c-SLICE, QP %d ) %10d bits",
2618         pcSlice->getPOC()-pcSlice->getLastIDR(),
2619         pcSlice->getLayerId(),
2620         pcSlice->getTLayer(),
2621         c,
2622         pcSlice->getSliceQp(),
2623         uibits );
2624#endif
2625#else
2626#if ADAPTIVE_QP_SELECTION
2627  printf("POC %4d TId: %1d ( %c-SLICE, nQP %d QP %d ) %10d bits",
2628         pcSlice->getPOC(),
2629         pcSlice->getTLayer(),
2630         c,
2631         pcSlice->getSliceQpBase(),
2632         pcSlice->getSliceQp(),
2633         uibits );
2634#else
2635  printf("POC %4d TId: %1d ( %c-SLICE, QP %d ) %10d bits",
2636         pcSlice->getPOC()-pcSlice->getLastIDR(),
2637         pcSlice->getTLayer(),
2638         c,
2639         pcSlice->getSliceQp(),
2640         uibits );
2641#endif
2642#endif
2643
2644  printf(" [Y %6.4lf dB    U %6.4lf dB    V %6.4lf dB]", dYPSNR, dUPSNR, dVPSNR );
2645  printf(" [ET %5.0f ]", dEncTime );
2646 
2647  for (Int iRefList = 0; iRefList < 2; iRefList++)
2648  {
2649    printf(" [L%d ", iRefList);
2650    for (Int iRefIndex = 0; iRefIndex < pcSlice->getNumRefIdx(RefPicList(iRefList)); iRefIndex++)
2651    {
2652#if REF_IDX_FRAMEWORK && VPS_EXTN_DIRECT_REF_LAYERS
2653      if( pcSlice->getRefPic(RefPicList(iRefList), iRefIndex)->isILR(m_layerId) )
2654      {
2655        printf( "%d(%d) ", pcSlice->getRefPOC(RefPicList(iRefList), iRefIndex)-pcSlice->getLastIDR(), pcSlice->getRefPic(RefPicList(iRefList), iRefIndex)->getLayerId() );
2656      }
2657      else
2658#endif
2659      printf ("%d ", pcSlice->getRefPOC(RefPicList(iRefList), iRefIndex)-pcSlice->getLastIDR());
2660    }
2661    printf("]");
2662  }
2663}
2664
2665/** Function for deciding the nal_unit_type.
2666 * \param pocCurr POC of the current picture
2667 * \returns the nal unit type of the picture
2668 * This function checks the configuration and returns the appropriate nal_unit_type for the picture.
2669 */
2670NalUnitType TEncGOP::getNalUnitType(Int pocCurr, Int lastIDR)
2671{
2672  if (pocCurr == 0)
2673  {
2674    return NAL_UNIT_CODED_SLICE_IDR_W_RADL;
2675  }
2676  if (pocCurr % m_pcCfg->getIntraPeriod() == 0)
2677  {
2678    if (m_pcCfg->getDecodingRefreshType() == 1)
2679    {
2680      return NAL_UNIT_CODED_SLICE_CRA;
2681    }
2682    else if (m_pcCfg->getDecodingRefreshType() == 2)
2683    {
2684      return NAL_UNIT_CODED_SLICE_IDR_W_RADL;
2685    }
2686  }
2687  if(m_pocCRA>0)
2688  {
2689    if(pocCurr<m_pocCRA)
2690    {
2691      // All leading pictures are being marked as TFD pictures here since current encoder uses all
2692      // reference pictures while encoding leading pictures. An encoder can ensure that a leading
2693      // picture can be still decodable when random accessing to a CRA/CRANT/BLA/BLANT picture by
2694      // controlling the reference pictures used for encoding that leading picture. Such a leading
2695      // picture need not be marked as a TFD picture.
2696      return NAL_UNIT_CODED_SLICE_RASL_R;
2697    }
2698  }
2699  if (lastIDR>0)
2700  {
2701    if (pocCurr < lastIDR)
2702    {
2703      return NAL_UNIT_CODED_SLICE_RADL_R;
2704    }
2705  }
2706  return NAL_UNIT_CODED_SLICE_TRAIL_R;
2707}
2708
2709Double TEncGOP::xCalculateRVM()
2710{
2711  Double dRVM = 0;
2712 
2713  if( m_pcCfg->getGOPSize() == 1 && m_pcCfg->getIntraPeriod() != 1 && m_pcCfg->getFramesToBeEncoded() > RVM_VCEGAM10_M * 2 )
2714  {
2715    // calculate RVM only for lowdelay configurations
2716    std::vector<Double> vRL , vB;
2717    size_t N = m_vRVM_RP.size();
2718    vRL.resize( N );
2719    vB.resize( N );
2720   
2721    Int i;
2722    Double dRavg = 0 , dBavg = 0;
2723    vB[RVM_VCEGAM10_M] = 0;
2724    for( i = RVM_VCEGAM10_M + 1 ; i < N - RVM_VCEGAM10_M + 1 ; i++ )
2725    {
2726      vRL[i] = 0;
2727      for( Int j = i - RVM_VCEGAM10_M ; j <= i + RVM_VCEGAM10_M - 1 ; j++ )
2728        vRL[i] += m_vRVM_RP[j];
2729      vRL[i] /= ( 2 * RVM_VCEGAM10_M );
2730      vB[i] = vB[i-1] + m_vRVM_RP[i] - vRL[i];
2731      dRavg += m_vRVM_RP[i];
2732      dBavg += vB[i];
2733    }
2734   
2735    dRavg /= ( N - 2 * RVM_VCEGAM10_M );
2736    dBavg /= ( N - 2 * RVM_VCEGAM10_M );
2737   
2738    Double dSigamB = 0;
2739    for( i = RVM_VCEGAM10_M + 1 ; i < N - RVM_VCEGAM10_M + 1 ; i++ )
2740    {
2741      Double tmp = vB[i] - dBavg;
2742      dSigamB += tmp * tmp;
2743    }
2744    dSigamB = sqrt( dSigamB / ( N - 2 * RVM_VCEGAM10_M ) );
2745   
2746    Double f = sqrt( 12.0 * ( RVM_VCEGAM10_M - 1 ) / ( RVM_VCEGAM10_M + 1 ) );
2747   
2748    dRVM = dSigamB / dRavg * f;
2749  }
2750 
2751  return( dRVM );
2752}
2753
2754/** Attaches the input bitstream to the stream in the output NAL unit
2755    Updates rNalu to contain concatenated bitstream. rpcBitstreamRedirect is cleared at the end of this function call.
2756 *  \param codedSliceData contains the coded slice data (bitstream) to be concatenated to rNalu
2757 *  \param rNalu          target NAL unit
2758 */
2759Void TEncGOP::xAttachSliceDataToNalUnit (OutputNALUnit& rNalu, TComOutputBitstream*& codedSliceData)
2760{
2761  // Byte-align
2762  rNalu.m_Bitstream.writeByteAlignment();   // Slice header byte-alignment
2763
2764  // Perform bitstream concatenation
2765  if (codedSliceData->getNumberOfWrittenBits() > 0)
2766    {
2767    rNalu.m_Bitstream.addSubstream(codedSliceData);
2768  }
2769
2770  m_pcEntropyCoder->setBitstream(&rNalu.m_Bitstream);
2771
2772  codedSliceData->clear();
2773}
2774
2775// Function will arrange the long-term pictures in the decreasing order of poc_lsb_lt,
2776// and among the pictures with the same lsb, it arranges them in increasing delta_poc_msb_cycle_lt value
2777Void TEncGOP::arrangeLongtermPicturesInRPS(TComSlice *pcSlice, TComList<TComPic*>& rcListPic)
2778{
2779  TComReferencePictureSet *rps = pcSlice->getRPS();
2780  if(!rps->getNumberOfLongtermPictures())
2781  {
2782    return;
2783  }
2784
2785  // Arrange long-term reference pictures in the correct order of LSB and MSB,
2786  // and assign values for pocLSBLT and MSB present flag
2787  Int longtermPicsPoc[MAX_NUM_REF_PICS], longtermPicsLSB[MAX_NUM_REF_PICS], indices[MAX_NUM_REF_PICS];
2788  Int longtermPicsMSB[MAX_NUM_REF_PICS];
2789  Bool mSBPresentFlag[MAX_NUM_REF_PICS];
2790  ::memset(longtermPicsPoc, 0, sizeof(longtermPicsPoc));    // Store POC values of LTRP
2791  ::memset(longtermPicsLSB, 0, sizeof(longtermPicsLSB));    // Store POC LSB values of LTRP
2792  ::memset(longtermPicsMSB, 0, sizeof(longtermPicsMSB));    // Store POC LSB values of LTRP
2793  ::memset(indices        , 0, sizeof(indices));            // Indices to aid in tracking sorted LTRPs
2794  ::memset(mSBPresentFlag , 0, sizeof(mSBPresentFlag));     // Indicate if MSB needs to be present
2795
2796  // Get the long-term reference pictures
2797  Int offset = rps->getNumberOfNegativePictures() + rps->getNumberOfPositivePictures();
2798  Int i, ctr = 0;
2799  Int maxPicOrderCntLSB = 1 << pcSlice->getSPS()->getBitsForPOC();
2800  for(i = rps->getNumberOfPictures() - 1; i >= offset; i--, ctr++)
2801  {
2802    longtermPicsPoc[ctr] = rps->getPOC(i);                                  // LTRP POC
2803    longtermPicsLSB[ctr] = getLSB(longtermPicsPoc[ctr], maxPicOrderCntLSB); // LTRP POC LSB
2804    indices[ctr]      = i; 
2805    longtermPicsMSB[ctr] = longtermPicsPoc[ctr] - longtermPicsLSB[ctr];
2806  }
2807  Int numLongPics = rps->getNumberOfLongtermPictures();
2808  assert(ctr == numLongPics);
2809
2810  // Arrange pictures in decreasing order of MSB;
2811  for(i = 0; i < numLongPics; i++)
2812  {
2813    for(Int j = 0; j < numLongPics - 1; j++)
2814    {
2815      if(longtermPicsMSB[j] < longtermPicsMSB[j+1])
2816      {
2817        std::swap(longtermPicsPoc[j], longtermPicsPoc[j+1]);
2818        std::swap(longtermPicsLSB[j], longtermPicsLSB[j+1]);
2819        std::swap(longtermPicsMSB[j], longtermPicsMSB[j+1]);
2820        std::swap(indices[j]        , indices[j+1]        );
2821      }
2822    }
2823  }
2824
2825  for(i = 0; i < numLongPics; i++)
2826  {
2827    // Check if MSB present flag should be enabled.
2828    // Check if the buffer contains any pictures that have the same LSB.
2829    TComList<TComPic*>::iterator  iterPic = rcListPic.begin(); 
2830    TComPic*                      pcPic;
2831    while ( iterPic != rcListPic.end() )
2832    {
2833      pcPic = *iterPic;
2834      if( (getLSB(pcPic->getPOC(), maxPicOrderCntLSB) == longtermPicsLSB[i])   &&     // Same LSB
2835                                      (pcPic->getSlice(0)->isReferenced())     &&    // Reference picture
2836                                        (pcPic->getPOC() != longtermPicsPoc[i])    )  // Not the LTRP itself
2837      {
2838        mSBPresentFlag[i] = true;
2839        break;
2840      }
2841      iterPic++;     
2842    }
2843  }
2844
2845  // tempArray for usedByCurr flag
2846  Bool tempArray[MAX_NUM_REF_PICS]; ::memset(tempArray, 0, sizeof(tempArray));
2847  for(i = 0; i < numLongPics; i++)
2848  {
2849    tempArray[i] = rps->getUsed(indices[i]);
2850  }
2851  // Now write the final values;
2852  ctr = 0;
2853  Int currMSB = 0, currLSB = 0;
2854  // currPicPoc = currMSB + currLSB
2855  currLSB = getLSB(pcSlice->getPOC(), maxPicOrderCntLSB); 
2856  currMSB = pcSlice->getPOC() - currLSB;
2857
2858  for(i = rps->getNumberOfPictures() - 1; i >= offset; i--, ctr++)
2859  {
2860    rps->setPOC                   (i, longtermPicsPoc[ctr]);
2861    rps->setDeltaPOC              (i, - pcSlice->getPOC() + longtermPicsPoc[ctr]);
2862    rps->setUsed                  (i, tempArray[ctr]);
2863    rps->setPocLSBLT              (i, longtermPicsLSB[ctr]);
2864    rps->setDeltaPocMSBCycleLT    (i, (currMSB - (longtermPicsPoc[ctr] - longtermPicsLSB[ctr])) / maxPicOrderCntLSB);
2865    rps->setDeltaPocMSBPresentFlag(i, mSBPresentFlag[ctr]);     
2866
2867    assert(rps->getDeltaPocMSBCycleLT(i) >= 0);   // Non-negative value
2868  }
2869  for(i = rps->getNumberOfPictures() - 1, ctr = 1; i >= offset; i--, ctr++)
2870  {
2871    for(Int j = rps->getNumberOfPictures() - 1 - ctr; j >= offset; j--)
2872    {
2873      // Here at the encoder we know that we have set the full POC value for the LTRPs, hence we
2874      // don't have to check the MSB present flag values for this constraint.
2875      assert( rps->getPOC(i) != rps->getPOC(j) ); // If assert fails, LTRP entry repeated in RPS!!!
2876    }
2877  }
2878}
2879
2880#if L0045_NON_NESTED_SEI_RESTRICTIONS
2881/** Function for finding the position to insert the first of APS and non-nested BP, PT, DU info SEI messages.
2882 * \param accessUnit Access Unit of the current picture
2883 * This function finds the position to insert the first of APS and non-nested BP, PT, DU info SEI messages.
2884 */
2885Int TEncGOP::xGetFirstSeiLocation(AccessUnit &accessUnit)
2886{
2887  // Find the location of the first SEI message
2888  AccessUnit::iterator it;
2889  Int seiStartPos = 0;
2890  for(it = accessUnit.begin(); it != accessUnit.end(); it++, seiStartPos++)
2891  {
2892     if ((*it)->isSei() || (*it)->isVcl())
2893     {
2894       break;
2895     }               
2896  }
2897//  assert(it != accessUnit.end());  // Triggers with some legit configurations
2898  return seiStartPos;
2899}
2900#endif
2901
2902#if L0386_DB_METRIC
2903Void TEncGOP::dblMetric( TComPic* pcPic, UInt uiNumSlices )
2904{
2905  TComPicYuv* pcPicYuvRec = pcPic->getPicYuvRec();
2906  Pel* Rec    = pcPicYuvRec->getLumaAddr( 0 );
2907  Pel* tempRec = Rec;
2908  Int  stride = pcPicYuvRec->getStride();
2909  UInt log2maxTB = pcPic->getSlice(0)->getSPS()->getQuadtreeTULog2MaxSize();
2910  UInt maxTBsize = (1<<log2maxTB);
2911  const UInt minBlockArtSize = 8;
2912  const UInt picWidth = pcPicYuvRec->getWidth();
2913  const UInt picHeight = pcPicYuvRec->getHeight();
2914  const UInt noCol = (picWidth>>log2maxTB);
2915  const UInt noRows = (picHeight>>log2maxTB);
2916  UInt64 *colSAD = (UInt64*)malloc(noCol*sizeof(UInt64));
2917  UInt64 *rowSAD = (UInt64*)malloc(noRows*sizeof(UInt64));
2918  UInt colIdx = 0;
2919  UInt rowIdx = 0;
2920  Pel p0, p1, p2, q0, q1, q2;
2921 
2922  Int qp = pcPic->getSlice(0)->getSliceQp();
2923  Int bitdepthScale = 1 << (g_bitDepthY-8);
2924  Int beta = TComLoopFilter::getBeta( qp ) * bitdepthScale;
2925  const Int thr2 = (beta>>2);
2926  const Int thr1 = 2*bitdepthScale;
2927  UInt a = 0;
2928 
2929  memset(colSAD, 0, noCol*sizeof(UInt64));
2930  memset(rowSAD, 0, noRows*sizeof(UInt64));
2931 
2932  if (maxTBsize > minBlockArtSize)
2933  {
2934    // Analyze vertical artifact edges
2935    for(Int c = maxTBsize; c < picWidth; c += maxTBsize)
2936    {
2937      for(Int r = 0; r < picHeight; r++)
2938      {
2939        p2 = Rec[c-3];
2940        p1 = Rec[c-2];
2941        p0 = Rec[c-1];
2942        q0 = Rec[c];
2943        q1 = Rec[c+1];
2944        q2 = Rec[c+2];
2945        a = ((abs(p2-(p1<<1)+p0)+abs(q0-(q1<<1)+q2))<<1);
2946        if ( thr1 < a && a < thr2)
2947        {
2948          colSAD[colIdx] += abs(p0 - q0);
2949        }
2950        Rec += stride;
2951      }
2952      colIdx++;
2953      Rec = tempRec;
2954    }
2955   
2956    // Analyze horizontal artifact edges
2957    for(Int r = maxTBsize; r < picHeight; r += maxTBsize)
2958    {
2959      for(Int c = 0; c < picWidth; c++)
2960      {
2961        p2 = Rec[c + (r-3)*stride];
2962        p1 = Rec[c + (r-2)*stride];
2963        p0 = Rec[c + (r-1)*stride];
2964        q0 = Rec[c + r*stride];
2965        q1 = Rec[c + (r+1)*stride];
2966        q2 = Rec[c + (r+2)*stride];
2967        a = ((abs(p2-(p1<<1)+p0)+abs(q0-(q1<<1)+q2))<<1);
2968        if (thr1 < a && a < thr2)
2969        {
2970          rowSAD[rowIdx] += abs(p0 - q0);
2971        }
2972      }
2973      rowIdx++;
2974    }
2975  }
2976 
2977  UInt64 colSADsum = 0;
2978  UInt64 rowSADsum = 0;
2979  for(Int c = 0; c < noCol-1; c++)
2980  {
2981    colSADsum += colSAD[c];
2982  }
2983  for(Int r = 0; r < noRows-1; r++)
2984  {
2985    rowSADsum += rowSAD[r];
2986  }
2987 
2988  colSADsum <<= 10;
2989  rowSADsum <<= 10;
2990  colSADsum /= (noCol-1);
2991  colSADsum /= picHeight;
2992  rowSADsum /= (noRows-1);
2993  rowSADsum /= picWidth;
2994 
2995  UInt64 avgSAD = ((colSADsum + rowSADsum)>>1);
2996  avgSAD >>= (g_bitDepthY-8);
2997 
2998  if ( avgSAD > 2048 )
2999  {
3000    avgSAD >>= 9;
3001    Int offset = Clip3(2,6,(Int)avgSAD);
3002    for (Int i=0; i<uiNumSlices; i++)
3003    {
3004      pcPic->getSlice(i)->setDeblockingFilterOverrideFlag(true);
3005      pcPic->getSlice(i)->setDeblockingFilterDisable(false);
3006      pcPic->getSlice(i)->setDeblockingFilterBetaOffsetDiv2( offset );
3007      pcPic->getSlice(i)->setDeblockingFilterTcOffsetDiv2( offset );
3008    }
3009  }
3010  else
3011  {
3012    for (Int i=0; i<uiNumSlices; i++)
3013    {
3014      pcPic->getSlice(i)->setDeblockingFilterOverrideFlag(false);
3015      pcPic->getSlice(i)->setDeblockingFilterDisable(        pcPic->getSlice(i)->getPPS()->getPicDisableDeblockingFilterFlag() );
3016      pcPic->getSlice(i)->setDeblockingFilterBetaOffsetDiv2( pcPic->getSlice(i)->getPPS()->getDeblockingFilterBetaOffsetDiv2() );
3017      pcPic->getSlice(i)->setDeblockingFilterTcOffsetDiv2(   pcPic->getSlice(i)->getPPS()->getDeblockingFilterTcOffsetDiv2()   );
3018    }
3019  }
3020 
3021  free(colSAD);
3022  free(rowSAD);
3023}
3024#endif
3025
3026#if M0457_COL_PICTURE_SIGNALING
3027TComPic* TEncGOP::getMotionPredIlp(TComSlice* pcSlice)
3028{
3029  TComPic* ilpPic = NULL;
3030  Int activeMotionPredReflayerIdx = 0;
3031
3032  for( Int i = 0; i < pcSlice->getActiveNumILRRefIdx(); i++ )
3033  {
3034    UInt refLayerIdc = pcSlice->getInterLayerPredLayerIdc(i);
3035    if (m_pcEncTop->getMotionPredEnabledFlag(refLayerIdc))
3036    {
3037      if (activeMotionPredReflayerIdx == pcSlice->getColRefLayerIdx())
3038      {
3039        ilpPic = m_pcEncTop->getIlpList()[i];
3040        break;
3041      }
3042      else
3043      {
3044        activeMotionPredReflayerIdx++;
3045      }
3046    }
3047  }
3048
3049  assert(ilpPic != NULL);
3050
3051  return ilpPic;
3052}
3053#endif
3054
3055//! \}
Note: See TracBrowser for help on using the repository browser.