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

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

Implementation of JCTVC-M0040 "Adaptive Resolution Change"

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