cqa-ff-mac-scheduler.cc

Go to the documentation of this file.

/*
 * Copyright (c) 2012 Centre Tecnologic de Telecomunicacions de Catalunya (CTTC)
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License version 2 as
 * published by the Free Software Foundation;
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 * Authors: Biljana Bojovic <bbojovic@cttc.es>, Nicola Baldo<nbaldo@cttc.es>.
 *
 * Note:
 * Implementation is using many common scheduler functionalities in its
 * original version implemented by Marco Miozzo<mmiozzo@cttc.es> in
 * Proportional Fair and Round Robin schedulers implementations.
 */
 
#include "cqa-ff-mac-scheduler.h"
 
#include "ff-mac-common.h"
#include "lte-amc.h"
#include "lte-vendor-specific-parameters.h"
 
#include <ns3/boolean.h>
#include <ns3/integer.h>
#include <ns3/log.h>
#include <ns3/math.h>
#include <ns3/pointer.h>
#include <ns3/simulator.h>
#include <ns3/string.h>
 
#include <cfloat>
#include <set>
#include <stdexcept>
 
namespace ns3
{
 
NS_LOG_COMPONENT_DEFINE("CqaFfMacScheduler");
 
static const int CqaType0AllocationRbg[4] = {
    10,  // RGB size 1
    26,  // RGB size 2
    63,  // RGB size 3
    110, // RGB size 4
};       // see table 7.1.6.1-1 of 36.213
 
NS_OBJECT_ENSURE_REGISTERED(CqaFfMacScheduler);
 
struct qos_rb_and_CQI_assigned_to_lc
{
    uint16_t resource_block_index; 
    uint8_t cqi_value_for_lc;      
};
 
bool
CQIValueDescComparator(uint8_t key1, uint8_t key2)
{
    return key1 > key2;
}
 
bool
CqaGroupDescComparator(int key1, int key2)
{
    return key1 > key2;
}
 
typedef uint8_t CQI_value;
typedef int RBG_index;
typedef int HOL_group;
 
typedef std::map<CQI_value, LteFlowId_t, bool (*)(uint8_t, uint8_t)> t_map_CQIToUE; // sorted
typedef std::map<RBG_index, t_map_CQIToUE> t_map_RBGToCQIsSorted;
typedef std::map<HOL_group, t_map_RBGToCQIsSorted> t_map_HOLGroupToRBGs;
 
typedef std::map<CQI_value, LteFlowId_t, bool (*)(uint8_t, uint8_t)>::iterator
    t_it_CQIToUE; // sorted
typedef std::map<RBG_index, t_map_CQIToUE>::iterator t_it_RBGToCQIsSorted;
typedef std::map<HOL_group, t_map_RBGToCQIsSorted>::iterator t_it_HOLGroupToRBGs;
 
typedef std::multimap<HOL_group, std::set<LteFlowId_t>, bool (*)(int, int)> t_map_HOLgroupToUEs;
typedef std::map<HOL_group, std::set<LteFlowId_t>>::iterator t_it_HOLgroupToUEs;
 
// typedef std::map<RBG_index,CQI_value>  map_RBG_to_CQI;
// typedef std::map<LteFlowId_t,map_RBG_to_CQI> map_flowId_to_CQI_map;
 
bool
CqaKeyDescComparator(uint16_t key1, uint16_t key2)
{
    return key1 > key2;
}
 
CqaFfMacScheduler::CqaFfMacScheduler()
    : m_cschedSapUser(nullptr),
      m_schedSapUser(nullptr),
      m_timeWindow(99.0),
      m_nextRntiUl(0)
{
    m_amc = CreateObject<LteAmc>();
    m_cschedSapProvider = new MemberCschedSapProvider<CqaFfMacScheduler>(this);
    m_schedSapProvider = new MemberSchedSapProvider<CqaFfMacScheduler>(this);
    m_ffrSapProvider = nullptr;
    m_ffrSapUser = new MemberLteFfrSapUser<CqaFfMacScheduler>(this);
}
 
CqaFfMacScheduler::~CqaFfMacScheduler()
{
    NS_LOG_FUNCTION(this);
}
 
void
CqaFfMacScheduler::DoDispose()
{
    NS_LOG_FUNCTION(this);
    m_dlHarqProcessesDciBuffer.clear();
    m_dlHarqProcessesTimer.clear();
    m_dlHarqProcessesRlcPduListBuffer.clear();
    m_dlInfoListBuffered.clear();
    m_ulHarqCurrentProcessId.clear();
    m_ulHarqProcessesStatus.clear();
    m_ulHarqProcessesDciBuffer.clear();
    delete m_cschedSapProvider;
    delete m_schedSapProvider;
    delete m_ffrSapUser;
}
 
TypeId
CqaFfMacScheduler::GetTypeId()
{
    static TypeId tid =
        TypeId("ns3::CqaFfMacScheduler")
            .SetParent<FfMacScheduler>()
            .SetGroupName("Lte")
            .AddConstructor<CqaFfMacScheduler>()
            .AddAttribute("CqiTimerThreshold",
                          "The number of TTIs a CQI is valid (default 1000 - 1 sec.)",
                          UintegerValue(1000),
                          MakeUintegerAccessor(&CqaFfMacScheduler::m_cqiTimersThreshold),
                          MakeUintegerChecker<uint32_t>())
            .AddAttribute("CqaMetric",
                          "CqaFfMacScheduler metric type that can be: CqaFf, CqaPf",
                          StringValue("CqaFf"),
                          MakeStringAccessor(&CqaFfMacScheduler::m_CqaMetric),
                          MakeStringChecker())
            .AddAttribute("HarqEnabled",
                          "Activate/Deactivate the HARQ [by default is active].",
                          BooleanValue(true),
                          MakeBooleanAccessor(&CqaFfMacScheduler::m_harqOn),
                          MakeBooleanChecker())
            .AddAttribute("UlGrantMcs",
                          "The MCS of the UL grant, must be [0..15] (default 0)",
                          UintegerValue(0),
                          MakeUintegerAccessor(&CqaFfMacScheduler::m_ulGrantMcs),
                          MakeUintegerChecker<uint8_t>());
    return tid;
}
 
void
CqaFfMacScheduler::SetFfMacCschedSapUser(FfMacCschedSapUser* s)
{
    m_cschedSapUser = s;
}
 
void
CqaFfMacScheduler::SetFfMacSchedSapUser(FfMacSchedSapUser* s)
{
    m_schedSapUser = s;
}
 
FfMacCschedSapProvider*
CqaFfMacScheduler::GetFfMacCschedSapProvider()
{
    return m_cschedSapProvider;
}
 
FfMacSchedSapProvider*
CqaFfMacScheduler::GetFfMacSchedSapProvider()
{
    return m_schedSapProvider;
}
 
void
CqaFfMacScheduler::SetLteFfrSapProvider(LteFfrSapProvider* s)
{
    m_ffrSapProvider = s;
}
 
LteFfrSapUser*
CqaFfMacScheduler::GetLteFfrSapUser()
{
    return m_ffrSapUser;
}
 
void
CqaFfMacScheduler::DoCschedCellConfigReq(
    const FfMacCschedSapProvider::CschedCellConfigReqParameters& params)
{
    NS_LOG_FUNCTION(this);
    // Read the subset of parameters used
    m_cschedCellConfig = params;
    m_rachAllocationMap.resize(m_cschedCellConfig.m_ulBandwidth, 0);
    FfMacCschedSapUser::CschedUeConfigCnfParameters cnf;
    cnf.m_result = SUCCESS;
    m_cschedSapUser->CschedUeConfigCnf(cnf);
}
 
void
CqaFfMacScheduler::DoCschedUeConfigReq(
    const FfMacCschedSapProvider::CschedUeConfigReqParameters& params)
{
    NS_LOG_FUNCTION(this << " RNTI " << params.m_rnti << " txMode "
                         << (uint16_t)params.m_transmissionMode);
    std::map<uint16_t, uint8_t>::iterator it = m_uesTxMode.find(params.m_rnti);
    if (it == m_uesTxMode.end())
    {
        m_uesTxMode.insert(std::pair<uint16_t, uint8_t>(params.m_rnti, params.m_transmissionMode));
        // generate HARQ buffers
        m_dlHarqCurrentProcessId.insert(std::pair<uint16_t, uint8_t>(params.m_rnti, 0));
        DlHarqProcessesStatus_t dlHarqPrcStatus;
        dlHarqPrcStatus.resize(8, 0);
        m_dlHarqProcessesStatus.insert(
            std::pair<uint16_t, DlHarqProcessesStatus_t>(params.m_rnti, dlHarqPrcStatus));
        DlHarqProcessesTimer_t dlHarqProcessesTimer;
        dlHarqProcessesTimer.resize(8, 0);
        m_dlHarqProcessesTimer.insert(
            std::pair<uint16_t, DlHarqProcessesTimer_t>(params.m_rnti, dlHarqProcessesTimer));
        DlHarqProcessesDciBuffer_t dlHarqdci;
        dlHarqdci.resize(8);
        m_dlHarqProcessesDciBuffer.insert(
            std::pair<uint16_t, DlHarqProcessesDciBuffer_t>(params.m_rnti, dlHarqdci));
        DlHarqRlcPduListBuffer_t dlHarqRlcPdu;
        dlHarqRlcPdu.resize(2);
        dlHarqRlcPdu.at(0).resize(8);
        dlHarqRlcPdu.at(1).resize(8);
        m_dlHarqProcessesRlcPduListBuffer.insert(
            std::pair<uint16_t, DlHarqRlcPduListBuffer_t>(params.m_rnti, dlHarqRlcPdu));
        m_ulHarqCurrentProcessId.insert(std::pair<uint16_t, uint8_t>(params.m_rnti, 0));
        UlHarqProcessesStatus_t ulHarqPrcStatus;
        ulHarqPrcStatus.resize(8, 0);
        m_ulHarqProcessesStatus.insert(
            std::pair<uint16_t, UlHarqProcessesStatus_t>(params.m_rnti, ulHarqPrcStatus));
        UlHarqProcessesDciBuffer_t ulHarqdci;
        ulHarqdci.resize(8);
        m_ulHarqProcessesDciBuffer.insert(
            std::pair<uint16_t, UlHarqProcessesDciBuffer_t>(params.m_rnti, ulHarqdci));
    }
    else
    {
        (*it).second = params.m_transmissionMode;
    }
}
 
void
CqaFfMacScheduler::DoCschedLcConfigReq(
    const FfMacCschedSapProvider::CschedLcConfigReqParameters& params)
{
    NS_LOG_FUNCTION(this << " New LC, rnti: " << params.m_rnti);
 
    NS_LOG_FUNCTION("LC configuration. Number of LCs:" << params.m_logicalChannelConfigList.size());
 
    // m_reconfigureFlat indicates if this is a reconfiguration or new UE is added, table  4.1.5 in
    // LTE MAC scheduler specification
    if (params.m_reconfigureFlag)
    {
        std::vector<LogicalChannelConfigListElement_s>::const_iterator lcit;
 
        for (lcit = params.m_logicalChannelConfigList.begin();
             lcit != params.m_logicalChannelConfigList.end();
             lcit++)
        {
            LteFlowId_t flowid = LteFlowId_t(params.m_rnti, lcit->m_logicalChannelIdentity);
 
            if (m_ueLogicalChannelsConfigList.find(flowid) == m_ueLogicalChannelsConfigList.end())
            {
                NS_LOG_ERROR("UE logical channels can not be reconfigured because it was not "
                             "configured before.");
            }
            else
            {
                m_ueLogicalChannelsConfigList.find(flowid)->second = *lcit;
            }
        }
 
    } // else new UE is added
    else
    {
        std::vector<LogicalChannelConfigListElement_s>::const_iterator lcit;
 
        for (lcit = params.m_logicalChannelConfigList.begin();
             lcit != params.m_logicalChannelConfigList.end();
             lcit++)
        {
            LteFlowId_t flowId = LteFlowId_t(params.m_rnti, lcit->m_logicalChannelIdentity);
            m_ueLogicalChannelsConfigList.insert(
                std::pair<LteFlowId_t, LogicalChannelConfigListElement_s>(flowId, *lcit));
        }
    }
 
    std::map<uint16_t, CqasFlowPerf_t>::iterator it;
 
    for (std::size_t i = 0; i < params.m_logicalChannelConfigList.size(); i++)
    {
        it = m_flowStatsDl.find(params.m_rnti);
 
        if (it == m_flowStatsDl.end())
        {
            double tbrDlInBytes =
                params.m_logicalChannelConfigList.at(i).m_eRabGuaranteedBitrateDl / 8; // byte/s
            double tbrUlInBytes =
                params.m_logicalChannelConfigList.at(i).m_eRabGuaranteedBitrateUl / 8; // byte/s
 
            CqasFlowPerf_t flowStatsDl;
            flowStatsDl.flowStart = Simulator::Now();
            flowStatsDl.totalBytesTransmitted = 0;
            flowStatsDl.lastTtiBytesTransmitted = 0;
            flowStatsDl.lastAveragedThroughput = 1;
            flowStatsDl.secondLastAveragedThroughput = 1;
            flowStatsDl.targetThroughput = tbrDlInBytes;
            m_flowStatsDl.insert(std::pair<uint16_t, CqasFlowPerf_t>(params.m_rnti, flowStatsDl));
            CqasFlowPerf_t flowStatsUl;
            flowStatsUl.flowStart = Simulator::Now();
            flowStatsUl.totalBytesTransmitted = 0;
            flowStatsUl.lastTtiBytesTransmitted = 0;
            flowStatsUl.lastAveragedThroughput = 1;
            flowStatsUl.secondLastAveragedThroughput = 1;
            flowStatsUl.targetThroughput = tbrUlInBytes;
            m_flowStatsUl.insert(std::pair<uint16_t, CqasFlowPerf_t>(params.m_rnti, flowStatsUl));
        }
        else
        {
            // update GBR from UeManager::SetupDataRadioBearer ()
            double tbrDlInBytes =
                params.m_logicalChannelConfigList.at(i).m_eRabGuaranteedBitrateDl / 8; // byte/s
            double tbrUlInBytes =
                params.m_logicalChannelConfigList.at(i).m_eRabGuaranteedBitrateUl / 8; // byte/s
            m_flowStatsDl[(*it).first].targetThroughput = tbrDlInBytes;
            m_flowStatsUl[(*it).first].targetThroughput = tbrUlInBytes;
        }
    }
}
 
void
CqaFfMacScheduler::DoCschedLcReleaseReq(
    const FfMacCschedSapProvider::CschedLcReleaseReqParameters& params)
{
    NS_LOG_FUNCTION(this);
    std::vector<uint8_t>::const_iterator it;
 
    for (it = params.m_logicalChannelIdentity.begin(); it != params.m_logicalChannelIdentity.end();
         it++)
    {
        LteFlowId_t flowId = LteFlowId_t(params.m_rnti, *it);
 
        // find the logical channel with the same Logical Channel Identity in the current list,
        // release it
        if (m_ueLogicalChannelsConfigList.find(flowId) != m_ueLogicalChannelsConfigList.end())
        {
            m_ueLogicalChannelsConfigList.erase(flowId);
        }
        else
        {
            NS_FATAL_ERROR("Logical channels cannot be released because it can not be found in the "
                           "list of active LCs");
        }
    }
 
    for (std::size_t i = 0; i < params.m_logicalChannelIdentity.size(); i++)
    {
        std::map<LteFlowId_t, FfMacSchedSapProvider::SchedDlRlcBufferReqParameters>::iterator it =
            m_rlcBufferReq.begin();
        std::map<LteFlowId_t, FfMacSchedSapProvider::SchedDlRlcBufferReqParameters>::iterator temp;
        while (it != m_rlcBufferReq.end())
        {
            if (((*it).first.m_rnti == params.m_rnti) &&
                ((*it).first.m_lcId == params.m_logicalChannelIdentity.at(i)))
            {
                temp = it;
                it++;
                m_rlcBufferReq.erase(temp);
            }
            else
            {
                it++;
            }
        }
    }
}
 
void
CqaFfMacScheduler::DoCschedUeReleaseReq(
    const FfMacCschedSapProvider::CschedUeReleaseReqParameters& params)
{
    NS_LOG_FUNCTION(this);
 
    for (int i = 0; i < MAX_LC_LIST; i++)
    {
        LteFlowId_t flowId = LteFlowId_t(params.m_rnti, i);
        // find the logical channel with the same Logical Channel Identity in the current list,
        // release it
        if (m_ueLogicalChannelsConfigList.find(flowId) != m_ueLogicalChannelsConfigList.end())
        {
            m_ueLogicalChannelsConfigList.erase(flowId);
        }
    }
 
    m_uesTxMode.erase(params.m_rnti);
    m_dlHarqCurrentProcessId.erase(params.m_rnti);
    m_dlHarqProcessesStatus.erase(params.m_rnti);
    m_dlHarqProcessesTimer.erase(params.m_rnti);
    m_dlHarqProcessesDciBuffer.erase(params.m_rnti);
    m_dlHarqProcessesRlcPduListBuffer.erase(params.m_rnti);
    m_ulHarqCurrentProcessId.erase(params.m_rnti);
    m_ulHarqProcessesStatus.erase(params.m_rnti);
    m_ulHarqProcessesDciBuffer.erase(params.m_rnti);
    m_flowStatsDl.erase(params.m_rnti);
    m_flowStatsUl.erase(params.m_rnti);
    m_ceBsrRxed.erase(params.m_rnti);
    std::map<LteFlowId_t, FfMacSchedSapProvider::SchedDlRlcBufferReqParameters>::iterator it =
        m_rlcBufferReq.begin();
    std::map<LteFlowId_t, FfMacSchedSapProvider::SchedDlRlcBufferReqParameters>::iterator temp;
    while (it != m_rlcBufferReq.end())
    {
        if ((*it).first.m_rnti == params.m_rnti)
        {
            temp = it;
            it++;
            m_rlcBufferReq.erase(temp);
        }
        else
        {
            it++;
        }
    }
    if (m_nextRntiUl == params.m_rnti)
    {
        m_nextRntiUl = 0;
    }
}
 
void
CqaFfMacScheduler::DoSchedDlRlcBufferReq(
    const FfMacSchedSapProvider::SchedDlRlcBufferReqParameters& params)
{
    NS_LOG_FUNCTION(this << params.m_rnti << (uint32_t)params.m_logicalChannelIdentity);
    // API generated by RLC for updating RLC parameters on a LC (tx and retx queues)
 
    std::map<LteFlowId_t, FfMacSchedSapProvider::SchedDlRlcBufferReqParameters>::iterator it;
 
    LteFlowId_t flow(params.m_rnti, params.m_logicalChannelIdentity);
 
    it = m_rlcBufferReq.find(flow);
 
    if (it == m_rlcBufferReq.end())
    {
        m_rlcBufferReq.insert(
            std::pair<LteFlowId_t, FfMacSchedSapProvider::SchedDlRlcBufferReqParameters>(flow,
                                                                                         params));
    }
    else
    {
        (*it).second = params;
    }
}
 
void
CqaFfMacScheduler::DoSchedDlPagingBufferReq(
    const FfMacSchedSapProvider::SchedDlPagingBufferReqParameters& params)
{
    NS_LOG_FUNCTION(this);
    NS_FATAL_ERROR("method not implemented");
}
 
void
CqaFfMacScheduler::DoSchedDlMacBufferReq(
    const FfMacSchedSapProvider::SchedDlMacBufferReqParameters& params)
{
    NS_LOG_FUNCTION(this);
    NS_FATAL_ERROR("method not implemented");
}
 
int
CqaFfMacScheduler::GetRbgSize(int dlbandwidth)
{
    for (int i = 0; i < 4; i++)
    {
        if (dlbandwidth < CqaType0AllocationRbg[i])
        {
            return (i + 1);
        }
    }
 
    return (-1);
}
 
unsigned int
CqaFfMacScheduler::LcActivePerFlow(uint16_t rnti)
{
    std::map<LteFlowId_t, FfMacSchedSapProvider::SchedDlRlcBufferReqParameters>::iterator it;
    unsigned int lcActive = 0;
    for (it = m_rlcBufferReq.begin(); it != m_rlcBufferReq.end(); it++)
    {
        if (((*it).first.m_rnti == rnti) && (((*it).second.m_rlcTransmissionQueueSize > 0) ||
                                             ((*it).second.m_rlcRetransmissionQueueSize > 0) ||
                                             ((*it).second.m_rlcStatusPduSize > 0)))
        {
            lcActive++;
        }
        if ((*it).first.m_rnti > rnti)
        {
            break;
        }
    }
    return (lcActive);
}
 
bool
CqaFfMacScheduler::HarqProcessAvailability(uint16_t rnti)
{
    NS_LOG_FUNCTION(this << rnti);
 
    std::map<uint16_t, uint8_t>::iterator it = m_dlHarqCurrentProcessId.find(rnti);
    if (it == m_dlHarqCurrentProcessId.end())
    {
        NS_FATAL_ERROR("No Process Id found for this RNTI " << rnti);
    }
    std::map<uint16_t, DlHarqProcessesStatus_t>::iterator itStat =
        m_dlHarqProcessesStatus.find(rnti);
    if (itStat == m_dlHarqProcessesStatus.end())
    {
        NS_FATAL_ERROR("No Process Id Statusfound for this RNTI " << rnti);
    }
    uint8_t i = (*it).second;
    do
    {
        i = (i + 1) % HARQ_PROC_NUM;
    } while (((*itStat).second.at(i) != 0) && (i != (*it).second));
 
    return (*itStat).second.at(i) == 0;
}
 
uint8_t
CqaFfMacScheduler::UpdateHarqProcessId(uint16_t rnti)
{
    NS_LOG_FUNCTION(this << rnti);
 
    if (!m_harqOn)
    {
        return (0);
    }
 
    std::map<uint16_t, uint8_t>::iterator it = m_dlHarqCurrentProcessId.find(rnti);
    if (it == m_dlHarqCurrentProcessId.end())
    {
        NS_FATAL_ERROR("No Process Id found for this RNTI " << rnti);
    }
    std::map<uint16_t, DlHarqProcessesStatus_t>::iterator itStat =
        m_dlHarqProcessesStatus.find(rnti);
    if (itStat == m_dlHarqProcessesStatus.end())
    {
        NS_FATAL_ERROR("No Process Id Statusfound for this RNTI " << rnti);
    }
    uint8_t i = (*it).second;
    do
    {
        i = (i + 1) % HARQ_PROC_NUM;
    } while (((*itStat).second.at(i) != 0) && (i != (*it).second));
    if ((*itStat).second.at(i) == 0)
    {
        (*it).second = i;
        (*itStat).second.at(i) = 1;
    }
    else
    {
        NS_FATAL_ERROR("No HARQ process available for RNTI "
                       << rnti << " check before update with HarqProcessAvailability");
    }
 
    return ((*it).second);
}
 
void
CqaFfMacScheduler::RefreshHarqProcesses()
{
    NS_LOG_FUNCTION(this);
 
    std::map<uint16_t, DlHarqProcessesTimer_t>::iterator itTimers;
    for (itTimers = m_dlHarqProcessesTimer.begin(); itTimers != m_dlHarqProcessesTimer.end();
         itTimers++)
    {
        for (uint16_t i = 0; i < HARQ_PROC_NUM; i++)
        {
            if ((*itTimers).second.at(i) == HARQ_DL_TIMEOUT)
            {
                // reset HARQ process
 
                NS_LOG_DEBUG(this << " Reset HARQ proc " << i << " for RNTI " << (*itTimers).first);
                std::map<uint16_t, DlHarqProcessesStatus_t>::iterator itStat =
                    m_dlHarqProcessesStatus.find((*itTimers).first);
                if (itStat == m_dlHarqProcessesStatus.end())
                {
                    NS_FATAL_ERROR("No Process Id Status found for this RNTI "
                                   << (*itTimers).first);
                }
                (*itStat).second.at(i) = 0;
                (*itTimers).second.at(i) = 0;
            }
            else
            {
                (*itTimers).second.at(i)++;
            }
        }
    }
}
 
void
CqaFfMacScheduler::DoSchedDlTriggerReq(
    const FfMacSchedSapProvider::SchedDlTriggerReqParameters& params)
{
    NS_LOG_FUNCTION(this << " Frame no. " << (params.m_sfnSf >> 4) << " subframe no. "
                         << (0xF & params.m_sfnSf));
    // API generated by RLC for triggering the scheduling of a DL subframe
    // evaluate the relative channel quality indicator for each UE per each RBG
    // (since we are using allocation type 0 the small unit of allocation is RBG)
    // Resource allocation type 0 (see sec 7.1.6.1 of 36.213)
 
    RefreshDlCqiMaps();
 
    int rbgSize = GetRbgSize(m_cschedCellConfig.m_dlBandwidth);
    int numberOfRBGs = m_cschedCellConfig.m_dlBandwidth / rbgSize;
    std::map<uint16_t, std::multimap<uint8_t, qos_rb_and_CQI_assigned_to_lc>>
        allocationMapPerRntiPerLCId;
    std::map<uint16_t, std::multimap<uint8_t, qos_rb_and_CQI_assigned_to_lc>>::iterator itMap;
    allocationMapPerRntiPerLCId.clear();
    bool (*key_function_pointer_groups)(int, int) = CqaGroupDescComparator;
    t_map_HOLgroupToUEs map_GBRHOLgroupToUE(key_function_pointer_groups);
    t_map_HOLgroupToUEs map_nonGBRHOLgroupToUE(key_function_pointer_groups);
    int grouping_parameter = 1000;
    double tolerance = 1.1;
    std::map<LteFlowId_t, int> UEtoHOL;
    std::vector<bool> rbgMap; // global RBGs map
    uint16_t rbgAllocatedNum = 0;
    std::set<uint16_t> rntiAllocated;
    rbgMap.resize(m_cschedCellConfig.m_dlBandwidth / rbgSize, false);
 
    rbgMap = m_ffrSapProvider->GetAvailableDlRbg();
    for (std::vector<bool>::iterator it = rbgMap.begin(); it != rbgMap.end(); it++)
    {
        if (*it)
        {
            rbgAllocatedNum++;
        }
    }
 
    FfMacSchedSapUser::SchedDlConfigIndParameters ret;
 
    //   update UL HARQ proc id
    std::map<uint16_t, uint8_t>::iterator itProcId;
    for (itProcId = m_ulHarqCurrentProcessId.begin(); itProcId != m_ulHarqCurrentProcessId.end();
         itProcId++)
    {
        (*itProcId).second = ((*itProcId).second + 1) % HARQ_PROC_NUM;
    }
 
    // RACH Allocation
    std::vector<bool> ulRbMap;
    ulRbMap.resize(m_cschedCellConfig.m_ulBandwidth, false);
    ulRbMap = m_ffrSapProvider->GetAvailableUlRbg();
    uint8_t maxContinuousUlBandwidth = 0;
    uint8_t tmpMinBandwidth = 0;
    uint16_t ffrRbStartOffset = 0;
    uint16_t tmpFfrRbStartOffset = 0;
    uint16_t index = 0;
 
    for (std::vector<bool>::iterator it = ulRbMap.begin(); it != ulRbMap.end(); it++)
    {
        if (*it)
        {
            if (tmpMinBandwidth > maxContinuousUlBandwidth)
            {
                maxContinuousUlBandwidth = tmpMinBandwidth;
                ffrRbStartOffset = tmpFfrRbStartOffset;
            }
            tmpMinBandwidth = 0;
        }
        else
        {
            if (tmpMinBandwidth == 0)
            {
                tmpFfrRbStartOffset = index;
            }
            tmpMinBandwidth++;
        }
        index++;
    }
 
    if (tmpMinBandwidth > maxContinuousUlBandwidth)
    {
        maxContinuousUlBandwidth = tmpMinBandwidth;
        ffrRbStartOffset = tmpFfrRbStartOffset;
    }
 
    m_rachAllocationMap.resize(m_cschedCellConfig.m_ulBandwidth, 0);
    uint16_t rbStart = 0;
    rbStart = ffrRbStartOffset;
    std::vector<RachListElement_s>::iterator itRach;
    for (itRach = m_rachList.begin(); itRach != m_rachList.end(); itRach++)
    {
        NS_ASSERT_MSG(m_amc->GetUlTbSizeFromMcs(m_ulGrantMcs, m_cschedCellConfig.m_ulBandwidth) >
                          (*itRach).m_estimatedSize,
                      " Default UL Grant MCS does not allow to send RACH messages");
        BuildRarListElement_s newRar;
        newRar.m_rnti = (*itRach).m_rnti;
        // DL-RACH Allocation
        // Ideal: no needs of configuring m_dci
        // UL-RACH Allocation
        newRar.m_grant.m_rnti = newRar.m_rnti;
        newRar.m_grant.m_mcs = m_ulGrantMcs;
        uint16_t rbLen = 1;
        uint16_t tbSizeBits = 0;
        // find lowest TB size that fits UL grant estimated size
        while ((tbSizeBits < (*itRach).m_estimatedSize) &&
               (rbStart + rbLen < (ffrRbStartOffset + maxContinuousUlBandwidth)))
        {
            rbLen++;
            tbSizeBits = m_amc->GetUlTbSizeFromMcs(m_ulGrantMcs, rbLen);
        }
        if (tbSizeBits < (*itRach).m_estimatedSize)
        {
            // no more allocation space: finish allocation
            break;
        }
        newRar.m_grant.m_rbStart = rbStart;
        newRar.m_grant.m_rbLen = rbLen;
        newRar.m_grant.m_tbSize = tbSizeBits / 8;
        newRar.m_grant.m_hopping = false;
        newRar.m_grant.m_tpc = 0;
        newRar.m_grant.m_cqiRequest = false;
        newRar.m_grant.m_ulDelay = false;
        NS_LOG_INFO(this << " UL grant allocated to RNTI " << (*itRach).m_rnti << " rbStart "
                         << rbStart << " rbLen " << rbLen << " MCS " << m_ulGrantMcs << " tbSize "
                         << newRar.m_grant.m_tbSize);
        for (uint16_t i = rbStart; i < rbStart + rbLen; i++)
        {
            m_rachAllocationMap.at(i) = (*itRach).m_rnti;
        }
 
        if (m_harqOn)
        {
            // generate UL-DCI for HARQ retransmissions
            UlDciListElement_s uldci;
            uldci.m_rnti = newRar.m_rnti;
            uldci.m_rbLen = rbLen;
            uldci.m_rbStart = rbStart;
            uldci.m_mcs = m_ulGrantMcs;
            uldci.m_tbSize = tbSizeBits / 8;
            uldci.m_ndi = 1;
            uldci.m_cceIndex = 0;
            uldci.m_aggrLevel = 1;
            uldci.m_ueTxAntennaSelection = 3; // antenna selection OFF
            uldci.m_hopping = false;
            uldci.m_n2Dmrs = 0;
            uldci.m_tpc = 0;            // no power control
            uldci.m_cqiRequest = false; // only period CQI at this stage
            uldci.m_ulIndex = 0;        // TDD parameter
            uldci.m_dai = 1;            // TDD parameter
            uldci.m_freqHopping = 0;
            uldci.m_pdcchPowerOffset = 0; // not used
 
            uint8_t harqId = 0;
            std::map<uint16_t, uint8_t>::iterator itProcId;
            itProcId = m_ulHarqCurrentProcessId.find(uldci.m_rnti);
            if (itProcId == m_ulHarqCurrentProcessId.end())
            {
                NS_FATAL_ERROR("No info find in HARQ buffer for UE " << uldci.m_rnti);
            }
            harqId = (*itProcId).second;
            std::map<uint16_t, UlHarqProcessesDciBuffer_t>::iterator itDci =
                m_ulHarqProcessesDciBuffer.find(uldci.m_rnti);
            if (itDci == m_ulHarqProcessesDciBuffer.end())
            {
                NS_FATAL_ERROR("Unable to find RNTI entry in UL DCI HARQ buffer for RNTI "
                               << uldci.m_rnti);
            }
            (*itDci).second.at(harqId) = uldci;
        }
 
        rbStart = rbStart + rbLen;
        ret.m_buildRarList.push_back(newRar);
    }
    m_rachList.clear();
 
    // Process DL HARQ feedback
    RefreshHarqProcesses();
    // retrieve past HARQ retx buffered
    if (!m_dlInfoListBuffered.empty())
    {
        if (!params.m_dlInfoList.empty())
        {
            NS_LOG_INFO(this << " Received DL-HARQ feedback");
            m_dlInfoListBuffered.insert(m_dlInfoListBuffered.end(),
                                        params.m_dlInfoList.begin(),
                                        params.m_dlInfoList.end());
        }
    }
    else
    {
        if (!params.m_dlInfoList.empty())
        {
            m_dlInfoListBuffered = params.m_dlInfoList;
        }
    }
    if (!m_harqOn)
    {
        // Ignore HARQ feedback
        m_dlInfoListBuffered.clear();
    }
    std::vector<DlInfoListElement_s> dlInfoListUntxed;
    for (std::size_t i = 0; i < m_dlInfoListBuffered.size(); i++)
    {
        std::set<uint16_t>::iterator itRnti = rntiAllocated.find(m_dlInfoListBuffered.at(i).m_rnti);
        if (itRnti != rntiAllocated.end())
        {
            // RNTI already allocated for retx
            continue;
        }
        auto nLayers = m_dlInfoListBuffered.at(i).m_harqStatus.size();
        std::vector<bool> retx;
        NS_LOG_INFO(this << " Processing DLHARQ feedback");
        if (nLayers == 1)
        {
            retx.push_back(m_dlInfoListBuffered.at(i).m_harqStatus.at(0) ==
                           DlInfoListElement_s::NACK);
            retx.push_back(false);
        }
        else
        {
            retx.push_back(m_dlInfoListBuffered.at(i).m_harqStatus.at(0) ==
                           DlInfoListElement_s::NACK);
            retx.push_back(m_dlInfoListBuffered.at(i).m_harqStatus.at(1) ==
                           DlInfoListElement_s::NACK);
        }
        if (retx.at(0) || retx.at(1))
        {
            // retrieve HARQ process information
            uint16_t rnti = m_dlInfoListBuffered.at(i).m_rnti;
            uint8_t harqId = m_dlInfoListBuffered.at(i).m_harqProcessId;
            NS_LOG_INFO(this << " HARQ retx RNTI " << rnti << " harqId " << (uint16_t)harqId);
            std::map<uint16_t, DlHarqProcessesDciBuffer_t>::iterator itHarq =
                m_dlHarqProcessesDciBuffer.find(rnti);
            if (itHarq == m_dlHarqProcessesDciBuffer.end())
            {
                NS_FATAL_ERROR("No info find in HARQ buffer for UE " << rnti);
            }
 
            DlDciListElement_s dci = (*itHarq).second.at(harqId);
            int rv = 0;
            if (dci.m_rv.size() == 1)
            {
                rv = dci.m_rv.at(0);
            }
            else
            {
                rv = (dci.m_rv.at(0) > dci.m_rv.at(1) ? dci.m_rv.at(0) : dci.m_rv.at(1));
            }
 
            if (rv == 3)
            {
                // maximum number of retx reached -> drop process
                NS_LOG_INFO("Maximum number of retransmissions reached -> drop process");
                std::map<uint16_t, DlHarqProcessesStatus_t>::iterator it =
                    m_dlHarqProcessesStatus.find(rnti);
                if (it == m_dlHarqProcessesStatus.end())
                {
                    NS_LOG_ERROR("No info find in HARQ buffer for UE (might change eNB) "
                                 << m_dlInfoListBuffered.at(i).m_rnti);
                }
                (*it).second.at(harqId) = 0;
                std::map<uint16_t, DlHarqRlcPduListBuffer_t>::iterator itRlcPdu =
                    m_dlHarqProcessesRlcPduListBuffer.find(rnti);
                if (itRlcPdu == m_dlHarqProcessesRlcPduListBuffer.end())
                {
                    NS_FATAL_ERROR("Unable to find RlcPdcList in HARQ buffer for RNTI "
                                   << m_dlInfoListBuffered.at(i).m_rnti);
                }
                for (std::size_t k = 0; k < (*itRlcPdu).second.size(); k++)
                {
                    (*itRlcPdu).second.at(k).at(harqId).clear();
                }
                continue;
            }
            // check the feasibility of retransmitting on the same RBGs
            // translate the DCI to Spectrum framework
            std::vector<int> dciRbg;
            uint32_t mask = 0x1;
            NS_LOG_INFO("Original RBGs " << dci.m_rbBitmap << " rnti " << dci.m_rnti);
            for (int j = 0; j < 32; j++)
            {
                if (((dci.m_rbBitmap & mask) >> j) == 1)
                {
                    dciRbg.push_back(j);
                    NS_LOG_INFO("\t" << j);
                }
                mask = (mask << 1);
            }
            bool free = true;
            for (std::size_t j = 0; j < dciRbg.size(); j++)
            {
                if (rbgMap.at(dciRbg.at(j)))
                {
                    free = false;
                    break;
                }
            }
            if (free)
            {
                // use the same RBGs for the retx
                // reserve RBGs
                for (std::size_t j = 0; j < dciRbg.size(); j++)
                {
                    rbgMap.at(dciRbg.at(j)) = true;
                    NS_LOG_INFO("RBG " << dciRbg.at(j) << " assigned");
                    rbgAllocatedNum++;
                }
 
                NS_LOG_INFO(this << " Send retx in the same RBGs");
            }
            else
            {
                // find RBGs for sending HARQ retx
                uint8_t j = 0;
                uint8_t rbgId = (dciRbg.at(dciRbg.size() - 1) + 1) % numberOfRBGs;
                uint8_t startRbg = dciRbg.at(dciRbg.size() - 1);
                std::vector<bool> rbgMapCopy = rbgMap;
                while ((j < dciRbg.size()) && (startRbg != rbgId))
                {
                    if (!rbgMapCopy.at(rbgId))
                    {
                        rbgMapCopy.at(rbgId) = true;
                        dciRbg.at(j) = rbgId;
                        j++;
                    }
                    rbgId = (rbgId + 1) % numberOfRBGs;
                }
                if (j == dciRbg.size())
                {
                    // find new RBGs -> update DCI map
                    uint32_t rbgMask = 0;
                    for (std::size_t k = 0; k < dciRbg.size(); k++)
                    {
                        rbgMask = rbgMask + (0x1 << dciRbg.at(k));
                        rbgAllocatedNum++;
                    }
                    dci.m_rbBitmap = rbgMask;
                    rbgMap = rbgMapCopy;
                    NS_LOG_INFO(this << " Move retx in RBGs " << dciRbg.size());
                }
                else
                {
                    // HARQ retx cannot be performed on this TTI -> store it
                    dlInfoListUntxed.push_back(m_dlInfoListBuffered.at(i));
                    NS_LOG_INFO(this << " No resource for this retx -> buffer it");
                }
            }
            // retrieve RLC PDU list for retx TBsize and update DCI
            BuildDataListElement_s newEl;
            std::map<uint16_t, DlHarqRlcPduListBuffer_t>::iterator itRlcPdu =
                m_dlHarqProcessesRlcPduListBuffer.find(rnti);
            if (itRlcPdu == m_dlHarqProcessesRlcPduListBuffer.end())
            {
                NS_FATAL_ERROR("Unable to find RlcPdcList in HARQ buffer for RNTI " << rnti);
            }
            for (std::size_t j = 0; j < nLayers; j++)
            {
                if (retx.at(j))
                {
                    if (j >= dci.m_ndi.size())
                    {
                        // for avoiding errors in MIMO transient phases
                        dci.m_ndi.push_back(0);
                        dci.m_rv.push_back(0);
                        dci.m_mcs.push_back(0);
                        dci.m_tbsSize.push_back(0);
                        NS_LOG_INFO(this << " layer " << (uint16_t)j
                                         << " no txed (MIMO transition)");
                    }
                    else
                    {
                        dci.m_ndi.at(j) = 0;
                        dci.m_rv.at(j)++;
                        (*itHarq).second.at(harqId).m_rv.at(j)++;
                        NS_LOG_INFO(this << " layer " << (uint16_t)j << " RV "
                                         << (uint16_t)dci.m_rv.at(j));
                    }
                }
                else
                {
                    // empty TB of layer j
                    dci.m_ndi.at(j) = 0;
                    dci.m_rv.at(j) = 0;
                    dci.m_mcs.at(j) = 0;
                    dci.m_tbsSize.at(j) = 0;
                    NS_LOG_INFO(this << " layer " << (uint16_t)j << " no retx");
                }
            }
            for (std::size_t k = 0; k < (*itRlcPdu).second.at(0).at(dci.m_harqProcess).size(); k++)
            {
                std::vector<RlcPduListElement_s> rlcPduListPerLc;
                for (std::size_t j = 0; j < nLayers; j++)
                {
                    if (retx.at(j))
                    {
                        if (j < dci.m_ndi.size())
                        {
                            NS_LOG_INFO(" layer " << (uint16_t)j << " tb size "
                                                  << dci.m_tbsSize.at(j));
                            rlcPduListPerLc.push_back(
                                (*itRlcPdu).second.at(j).at(dci.m_harqProcess).at(k));
                        }
                    }
                    else
                    { // if no retx needed on layer j, push an RlcPduListElement_s object with
                      // m_size=0 to keep the size of rlcPduListPerLc vector = 2 in case of MIMO
                        NS_LOG_INFO(" layer " << (uint16_t)j << " tb size " << dci.m_tbsSize.at(j));
                        RlcPduListElement_s emptyElement;
                        emptyElement.m_logicalChannelIdentity = (*itRlcPdu)
                                                                    .second.at(j)
                                                                    .at(dci.m_harqProcess)
                                                                    .at(k)
                                                                    .m_logicalChannelIdentity;
                        emptyElement.m_size = 0;
                        rlcPduListPerLc.push_back(emptyElement);
                    }
                }
 
                if (!rlcPduListPerLc.empty())
                {
                    newEl.m_rlcPduList.push_back(rlcPduListPerLc);
                }
            }
            newEl.m_rnti = rnti;
            newEl.m_dci = dci;
            (*itHarq).second.at(harqId).m_rv = dci.m_rv;
            // refresh timer
            std::map<uint16_t, DlHarqProcessesTimer_t>::iterator itHarqTimer =
                m_dlHarqProcessesTimer.find(rnti);
            if (itHarqTimer == m_dlHarqProcessesTimer.end())
            {
                NS_FATAL_ERROR("Unable to find HARQ timer for RNTI " << (uint16_t)rnti);
            }
            (*itHarqTimer).second.at(harqId) = 0;
            ret.m_buildDataList.push_back(newEl);
            rntiAllocated.insert(rnti);
        }
        else
        {
            // update HARQ process status
            NS_LOG_INFO(this << " HARQ received ACK for UE " << m_dlInfoListBuffered.at(i).m_rnti);
            std::map<uint16_t, DlHarqProcessesStatus_t>::iterator it =
                m_dlHarqProcessesStatus.find(m_dlInfoListBuffered.at(i).m_rnti);
            if (it == m_dlHarqProcessesStatus.end())
            {
                NS_FATAL_ERROR("No info find in HARQ buffer for UE "
                               << m_dlInfoListBuffered.at(i).m_rnti);
            }
            (*it).second.at(m_dlInfoListBuffered.at(i).m_harqProcessId) = 0;
            std::map<uint16_t, DlHarqRlcPduListBuffer_t>::iterator itRlcPdu =
                m_dlHarqProcessesRlcPduListBuffer.find(m_dlInfoListBuffered.at(i).m_rnti);
            if (itRlcPdu == m_dlHarqProcessesRlcPduListBuffer.end())
            {
                NS_FATAL_ERROR("Unable to find RlcPdcList in HARQ buffer for RNTI "
                               << m_dlInfoListBuffered.at(i).m_rnti);
            }
            for (std::size_t k = 0; k < (*itRlcPdu).second.size(); k++)
            {
                (*itRlcPdu).second.at(k).at(m_dlInfoListBuffered.at(i).m_harqProcessId).clear();
            }
        }
    }
    m_dlInfoListBuffered.clear();
    m_dlInfoListBuffered = dlInfoListUntxed;
 
    if (rbgAllocatedNum == numberOfRBGs)
    {
        // all the RBGs are already allocated -> exit
        if (!ret.m_buildDataList.empty() || !ret.m_buildRarList.empty())
        {
            m_schedSapUser->SchedDlConfigInd(ret);
        }
        return;
    }
 
    std::map<LteFlowId_t, LogicalChannelConfigListElement_s>::iterator itLogicalChannels;
 
    for (itLogicalChannels = m_ueLogicalChannelsConfigList.begin();
         itLogicalChannels != m_ueLogicalChannelsConfigList.end();
         itLogicalChannels++)
    {
        std::set<uint16_t>::iterator itRnti = rntiAllocated.find(itLogicalChannels->first.m_rnti);
        if ((itRnti != rntiAllocated.end()) ||
            (!HarqProcessAvailability(itLogicalChannels->first.m_rnti)))
        {
            // UE already allocated for HARQ or without HARQ process available -> drop it
            if (itRnti != rntiAllocated.end())
            {
                NS_LOG_DEBUG(this << " RNTI discarded for HARQ tx"
                                  << (uint16_t)(itLogicalChannels->first.m_rnti));
            }
            if (!HarqProcessAvailability(itLogicalChannels->first.m_rnti))
            {
                NS_LOG_DEBUG(this << " RNTI discarded for HARQ id"
                                  << (uint16_t)(itLogicalChannels->first.m_rnti));
            }
            continue;
        }
 
        std::map<LteFlowId_t, FfMacSchedSapProvider::SchedDlRlcBufferReqParameters>::iterator
            itRlcBufferReq = m_rlcBufferReq.find(itLogicalChannels->first);
        if (itRlcBufferReq == m_rlcBufferReq.end())
        {
            continue;
        }
 
        int group = -1;
        int delay = 0;
 
        if (itRlcBufferReq->second.m_rlcRetransmissionQueueSize > 0)
        {
            delay = itRlcBufferReq->second.m_rlcRetransmissionHolDelay;
            group = delay / grouping_parameter;
        }
        else if (itRlcBufferReq->second.m_rlcTransmissionQueueSize > 0)
        {
            delay = itRlcBufferReq->second.m_rlcTransmissionQueueHolDelay;
            group = delay / grouping_parameter;
        }
        else
        {
            continue;
        }
 
        UEtoHOL.insert(std::pair<LteFlowId_t, int>(itLogicalChannels->first, delay));
 
        if (itLogicalChannels->second.m_qosBearerType ==
            LogicalChannelConfigListElement_s::QBT_NON_GBR)
        {
            if (map_nonGBRHOLgroupToUE.count(group) == 0)
            {
                std::set<LteFlowId_t> v;
                v.insert(itRlcBufferReq->first);
                map_nonGBRHOLgroupToUE.insert(std::pair<int, std::set<LteFlowId_t>>(group, v));
            }
            else
            {
                map_nonGBRHOLgroupToUE.find(group)->second.insert(itRlcBufferReq->first);
            }
        }
        else if (itLogicalChannels->second.m_qosBearerType ==
                     LogicalChannelConfigListElement_s::QBT_GBR ||
                 itLogicalChannels->second.m_qosBearerType ==
                     LogicalChannelConfigListElement_s::QBT_DGBR)
        {
            if (map_GBRHOLgroupToUE.count(group) == 0)
            {
                std::set<LteFlowId_t> v;
                v.insert(itRlcBufferReq->first);
                map_GBRHOLgroupToUE.insert(std::pair<int, std::set<LteFlowId_t>>(group, v));
            }
            else
            {
                map_GBRHOLgroupToUE.find(group)->second.insert(itRlcBufferReq->first);
            }
        }
    };
 
    // Prepare data for the scheduling mechanism
    // map: UE, to the amount of traffic they have to transfer
    std::map<LteFlowId_t, int> UeToAmountOfDataToTransfer;
    // Initialize the map per UE, how much resources is already assigned to the user
    std::map<LteFlowId_t, int> UeToAmountOfAssignedResources;
    // prepare values to calculate FF metric, this metric will be the same for all flows(logical
    // channels) that belong to the same RNTI
    std::map<uint16_t, uint8_t> sbCqiSum;
 
    for (std::map<LteFlowId_t, FfMacSchedSapProvider::SchedDlRlcBufferReqParameters>::iterator
             itrbr = m_rlcBufferReq.begin();
         itrbr != m_rlcBufferReq.end();
         itrbr++)
    {
        LteFlowId_t flowId = itrbr->first; // Prepare data for the scheduling mechanism
        // check first the channel conditions for this UE, if CQI!=0
        std::map<uint16_t, SbMeasResult_s>::iterator itCqi;
        itCqi = m_a30CqiRxed.find((*itrbr).first.m_rnti);
        std::map<uint16_t, uint8_t>::iterator itTxMode;
        itTxMode = m_uesTxMode.find((*itrbr).first.m_rnti);
        if (itTxMode == m_uesTxMode.end())
        {
            NS_FATAL_ERROR("No Transmission Mode info on user " << (*itrbr).first.m_rnti);
        }
        auto nLayer = TransmissionModesLayers::TxMode2LayerNum((*itTxMode).second);
 
        uint8_t cqiSum = 0;
        for (int k = 0; k < numberOfRBGs; k++)
        {
            for (uint8_t j = 0; j < nLayer; j++)
            {
                if (itCqi == m_a30CqiRxed.end())
                {
                    cqiSum += 1; // no info on this user -> lowest MCS
                }
                else
                {
                    cqiSum += (*itCqi).second.m_higherLayerSelected.at(k).m_sbCqi.at(j);
                }
            }
        }
 
        if (cqiSum == 0)
        {
            NS_LOG_INFO("Skip this flow, CQI==0, rnti:" << (*itrbr).first.m_rnti);
            continue;
        }
 
        // map: UE, to the amount of traffic they have to transfer
        int amountOfDataToTransfer =
            8 * ((int)m_rlcBufferReq.find(flowId)->second.m_rlcRetransmissionQueueSize +
                 (int)m_rlcBufferReq.find(flowId)->second.m_rlcTransmissionQueueSize);
 
        UeToAmountOfDataToTransfer.insert(
            std::pair<LteFlowId_t, int>(flowId, amountOfDataToTransfer));
        UeToAmountOfAssignedResources.insert(std::pair<LteFlowId_t, int>(flowId, 0));
 
        uint8_t sum = 0;
        for (int i = 0; i < numberOfRBGs; i++)
        {
            std::map<uint16_t, SbMeasResult_s>::iterator itCqi;
            itCqi = m_a30CqiRxed.find((*itrbr).first.m_rnti);
            std::map<uint16_t, uint8_t>::iterator itTxMode;
            itTxMode = m_uesTxMode.find((*itrbr).first.m_rnti);
            if (itTxMode == m_uesTxMode.end())
            {
                NS_FATAL_ERROR("No Transmission Mode info on user " << (*itrbr).first.m_rnti);
            }
            auto nLayer = TransmissionModesLayers::TxMode2LayerNum((*itTxMode).second);
            std::vector<uint8_t> sbCqis;
            if (itCqi == m_a30CqiRxed.end())
            {
                for (uint8_t k = 0; k < nLayer; k++)
                {
                    sbCqis.push_back(1); // start with lowest value
                }
            }
            else
            {
                sbCqis = (*itCqi).second.m_higherLayerSelected.at(i).m_sbCqi;
            }
 
            uint8_t cqi1 = sbCqis.at(0);
            uint8_t cqi2 = 0;
            if (sbCqis.size() > 1)
            {
                cqi2 = sbCqis.at(1);
            }
 
            uint8_t sbCqi = 0;
            if ((cqi1 > 0) ||
                (cqi2 > 0)) // CQI == 0 means "out of range" (see table 7.2.3-1 of 36.213)
            {
                for (uint8_t k = 0; k < nLayer; k++)
                {
                    if (sbCqis.size() > k)
                    {
                        sbCqi = sbCqis.at(k);
                    }
                    else
                    {
                        // no info on this subband
                        sbCqi = 0;
                    }
                    sum += sbCqi;
                }
            } // end if cqi
        }     // end of rbgNum
 
        sbCqiSum.insert(std::pair<uint16_t, uint8_t>((*itrbr).first.m_rnti, sum));
    }
 
    // availableRBGs - set that contains indexes of available resource block groups
    std::set<int> availableRBGs;
    for (int i = 0; i < numberOfRBGs; i++)
    {
        if (!rbgMap.at(i))
        {
            availableRBGs.insert(i);
        }
    }
 
    t_it_HOLgroupToUEs itGBRgroups = map_GBRHOLgroupToUE.begin();
    t_it_HOLgroupToUEs itnonGBRgroups = map_nonGBRHOLgroupToUE.begin();
 
    // while there are more resources available, loop through the users that are grouped by HOL
    // value
    while (!availableRBGs.empty())
    {
        if (UeToAmountOfDataToTransfer.empty())
        {
            NS_LOG_INFO("No UEs to be scheduled (no data or CQI==0),");
            break;
        }
        std::set<LteFlowId_t> vUEs;
        t_it_HOLgroupToUEs itCurrentGroup;
 
        if (itGBRgroups != map_GBRHOLgroupToUE.end())
        {
            itCurrentGroup = itGBRgroups;
            itGBRgroups++;
        }
        else if (itnonGBRgroups !=
                 map_nonGBRHOLgroupToUE
                     .end()) // if there are no more flows with retransmission queue start to
                             // scheduler flows with transmission queue
        {
            itCurrentGroup = itnonGBRgroups;
            itnonGBRgroups++;
        }
        else
        {
            NS_LOG_INFO("Available RBGs:" << availableRBGs.size() << "but no users");
            break;
        }
 
        while (!availableRBGs.empty() && !itCurrentGroup->second.empty())
        {
            bool currentRBchecked = false;
            int currentRB = *(availableRBGs.begin());
            std::map<LteFlowId_t, CQI_value> UeToCQIValue;
            std::map<LteFlowId_t, double> UeToCoitaMetric;
            std::map<LteFlowId_t, bool> UeHasReachedGBR;
            double maximumValueMetric = 0;
            LteFlowId_t userWithMaximumMetric;
            UeToCQIValue.clear();
            UeToCoitaMetric.clear();
 
            // Iterate through the users and calculate which user will use the best of the current
            // resource block.end() and assign to that user.
            for (std::set<LteFlowId_t>::iterator it = itCurrentGroup->second.begin();
                 it != itCurrentGroup->second.end();
                 it++)
            {
                LteFlowId_t flowId = *it;
                uint8_t cqi_value = 1; // higher better, maximum is 15
                double coita_metric = 1;
                double coita_sum = 0;
                double metric = 0;
                uint8_t worstCQIAmongRBGsAllocatedForThisUser = 15;
                int numberOfRBGAllocatedForThisUser = 0;
                LogicalChannelConfigListElement_s lc =
                    m_ueLogicalChannelsConfigList.find(flowId)->second;
                std::map<uint16_t, SbMeasResult_s>::iterator itRntiCQIsMap =
                    m_a30CqiRxed.find(flowId.m_rnti);
 
                std::map<uint16_t, CqasFlowPerf_t>::iterator itStats;
 
                if (!m_ffrSapProvider->IsDlRbgAvailableForUe(currentRB, flowId.m_rnti))
                {
                    continue;
                }
 
                if (m_flowStatsDl.find(flowId.m_rnti) == m_flowStatsDl.end())
                {
                    continue; // TO DO:  check if this should be logged and how.
                }
                currentRBchecked = true;
 
                itStats = m_flowStatsDl.find(flowId.m_rnti);
                double tbr_weight =
                    (*itStats).second.targetThroughput / (*itStats).second.lastAveragedThroughput;
                if (tbr_weight < 1.0)
                {
                    tbr_weight = 1.0;
                }
 
                if (itRntiCQIsMap != m_a30CqiRxed.end())
                {
                    for (std::set<int>::iterator it = availableRBGs.begin();
                         it != availableRBGs.end();
                         it++)
                    {
                        try
                        {
                            int val =
                                (itRntiCQIsMap->second.m_higherLayerSelected.at(*it).m_sbCqi.at(0));
                            if (val == 0)
                            {
                                val = 1; // if no info, use minimum
                            }
                            if (*it == currentRB)
                            {
                                cqi_value = val;
                            }
                            coita_sum += val;
                        }
                        catch (std::out_of_range&)
                        {
                            coita_sum += 1; // if no info on channel use the worst cqi
                            NS_LOG_INFO("No CQI for lcId:" << flowId.m_lcId
                                                           << " rnti:" << flowId.m_rnti
                                                           << " at subband:" << currentRB);
                            // std::cout<<"\n No CQI for
                            // lcId:.....................................";
                        }
                    }
                    coita_metric = cqi_value / coita_sum;
                    UeToCQIValue.insert(std::pair<LteFlowId_t, CQI_value>(flowId, cqi_value));
                    UeToCoitaMetric.insert(std::pair<LteFlowId_t, double>(flowId, coita_metric));
                }
 
                if (allocationMapPerRntiPerLCId.find(flowId.m_rnti) ==
                    allocationMapPerRntiPerLCId.end())
                {
                    worstCQIAmongRBGsAllocatedForThisUser = cqi_value;
                }
                else
                {
                    numberOfRBGAllocatedForThisUser =
                        (allocationMapPerRntiPerLCId.find(flowId.m_rnti)->second.size());
 
                    for (std::multimap<uint8_t, qos_rb_and_CQI_assigned_to_lc>::iterator itRBG =
                             allocationMapPerRntiPerLCId.find(flowId.m_rnti)->second.begin();
                         itRBG != allocationMapPerRntiPerLCId.find(flowId.m_rnti)->second.end();
                         itRBG++)
                    {
                        qos_rb_and_CQI_assigned_to_lc e = itRBG->second;
                        if (e.cqi_value_for_lc < worstCQIAmongRBGsAllocatedForThisUser)
                        {
                            worstCQIAmongRBGsAllocatedForThisUser = e.cqi_value_for_lc;
                        }
                    }
 
                    if (cqi_value < worstCQIAmongRBGsAllocatedForThisUser)
                    {
                        worstCQIAmongRBGsAllocatedForThisUser = cqi_value;
                    }
                }
 
                int mcsForThisUser = m_amc->GetMcsFromCqi(worstCQIAmongRBGsAllocatedForThisUser);
                int tbSize =
                    m_amc->GetDlTbSizeFromMcs(mcsForThisUser,
                                              (numberOfRBGAllocatedForThisUser + 1) * rbgSize) /
                    8; // similar to calculation of TB size (size of TB in bytes according to
                       // table 7.1.7.2.1-1 of 36.213)
 
                double achievableRate =
                    ((m_amc->GetDlTbSizeFromMcs(mcsForThisUser, rbgSize) / 8) / 0.001);
                double pf_weight = achievableRate / (*itStats).second.secondLastAveragedThroughput;
 
                UeToAmountOfAssignedResources.find(flowId)->second = 8 * tbSize;
                FfMacSchedSapProvider::SchedDlRlcBufferReqParameters lcBufferInfo =
                    m_rlcBufferReq.find(flowId)->second;
 
                if (UeToAmountOfDataToTransfer.find(flowId)->second -
                        UeToAmountOfAssignedResources.find(flowId)->second <
                    0)
                {
                    UeHasReachedGBR.insert(std::pair<LteFlowId_t, bool>(flowId, false));
                }
 
                double bitRateWithNewRBG = 0;
 
                if (m_flowStatsDl.find(flowId.m_rnti) !=
                    m_flowStatsDl.end()) // there are some statistics{
                {
                    bitRateWithNewRBG =
                        (1.0 - (1.0 / m_timeWindow)) *
                            (m_flowStatsDl.find(flowId.m_rnti)->second.lastAveragedThroughput) +
                        ((1.0 / m_timeWindow) * (double)(tbSize * 1000));
                }
                else
                {
                    bitRateWithNewRBG = (1.0 / m_timeWindow) * (double)(tbSize * 1000);
                }
 
                if (bitRateWithNewRBG > lc.m_eRabGuaranteedBitrateDl)
                {
                    UeHasReachedGBR.insert(std::pair<LteFlowId_t, bool>(flowId, true));
                }
                else
                {
                    UeHasReachedGBR.insert(std::pair<LteFlowId_t, bool>(flowId, false));
                }
 
                int hol = UEtoHOL.find(flowId)->second;
 
                if (hol == 0)
                {
                    hol = 1;
                }
 
                if (m_CqaMetric == "CqaFf")
                {
                    metric = coita_metric * tbr_weight * hol;
                }
                else if (m_CqaMetric == "CqaPf")
                {
                    metric = tbr_weight * pf_weight * hol;
                }
                else
                {
                    metric = 1;
                }
 
                if (metric >= maximumValueMetric)
                {
                    maximumValueMetric = metric;
                    userWithMaximumMetric = flowId;
                }
            }
 
            if (!currentRBchecked)
            {
                // erase current RBG from the list of available RBG
                availableRBGs.erase(currentRB);
                continue;
            }
 
            qos_rb_and_CQI_assigned_to_lc s;
            s.cqi_value_for_lc = UeToCQIValue.find(userWithMaximumMetric)->second;
            s.resource_block_index = currentRB;
 
            itMap = allocationMapPerRntiPerLCId.find(userWithMaximumMetric.m_rnti);
 
            if (itMap == allocationMapPerRntiPerLCId.end())
            {
                std::multimap<uint8_t, qos_rb_and_CQI_assigned_to_lc> tempMap;
                tempMap.insert(
                    std::pair<uint8_t, qos_rb_and_CQI_assigned_to_lc>(userWithMaximumMetric.m_lcId,
                                                                      s));
                allocationMapPerRntiPerLCId.insert(
                    std::pair<uint16_t, std::multimap<uint8_t, qos_rb_and_CQI_assigned_to_lc>>(
                        userWithMaximumMetric.m_rnti,
                        tempMap));
            }
            else
            {
                itMap->second.insert(
                    std::pair<uint8_t, qos_rb_and_CQI_assigned_to_lc>(userWithMaximumMetric.m_lcId,
                                                                      s));
            }
 
            // erase current RBG from the list of available RBG
            availableRBGs.erase(currentRB);
 
            if (UeToAmountOfDataToTransfer.find(userWithMaximumMetric)->second <=
                UeToAmountOfAssignedResources.find(userWithMaximumMetric)->second * tolerance)
            //||(UeHasReachedGBR.find(userWithMaximumMetric)->second == true))
            {
                itCurrentGroup->second.erase(userWithMaximumMetric);
            }
 
        } // while there are more users in current group
    }     // while there are more groups of users
 
    // reset TTI stats of users
    std::map<uint16_t, CqasFlowPerf_t>::iterator itStats;
    for (itStats = m_flowStatsDl.begin(); itStats != m_flowStatsDl.end(); itStats++)
    {
        (*itStats).second.lastTtiBytesTransmitted = 0;
    }
 
    // 3) Creating the correspondent DCIs (Generate the transmission opportunities by grouping the
    // RBGs of the same RNTI)
    // FfMacSchedSapUser::SchedDlConfigIndParameters ret;
    itMap = allocationMapPerRntiPerLCId.begin();
    std::map<uint16_t, double> m_rnti_per_ratio;
 
    while (itMap != allocationMapPerRntiPerLCId.end())
    {
        // create new BuildDataListElement_s for this LC
        BuildDataListElement_s newEl;
        newEl.m_rnti = (*itMap).first;
        NS_LOG_INFO("Scheduled RNTI:" << newEl.m_rnti);
        // create the DlDciListElement_s
        DlDciListElement_s newDci;
        std::vector<RlcPduListElement_s> newRlcPduLe;
        newDci.m_rnti = (*itMap).first;
        newDci.m_harqProcess = UpdateHarqProcessId((*itMap).first);
        uint16_t lcActives = LcActivePerFlow(itMap->first);
        if (lcActives == 0)
        { // if there is still no buffer report information on any flow
            lcActives = 1;
        }
        // NS_LOG_DEBUG (this << "Allocate user " << newEl.m_rnti << " rbg " << lcActives);
        uint16_t RgbPerRnti = (*itMap).second.size();
        double doubleRBgPerRnti = RgbPerRnti;
        double doubleRbgNum = numberOfRBGs;
        double rrRatio = doubleRBgPerRnti / doubleRbgNum;
        m_rnti_per_ratio.insert(std::pair<uint16_t, double>((*itMap).first, rrRatio));
        std::map<uint16_t, SbMeasResult_s>::iterator itCqi;
        itCqi = m_a30CqiRxed.find((*itMap).first);
        uint8_t worstCqi = 15;
 
        // assign the worst value of CQI that user experienced on any of its subbands
        for (std::multimap<uint8_t, qos_rb_and_CQI_assigned_to_lc>::iterator it =
                 (*itMap).second.begin();
             it != (*itMap).second.end();
             it++)
        {
            if (it->second.cqi_value_for_lc < worstCqi)
            {
                worstCqi = it->second.cqi_value_for_lc;
            }
        }
 
        newDci.m_mcs.push_back(m_amc->GetMcsFromCqi(worstCqi));
        int tbSize = (m_amc->GetDlTbSizeFromMcs(newDci.m_mcs.at(0), RgbPerRnti * rbgSize) /
                      8); // (size of TB in bytes according to table 7.1.7.2.1-1 of 36.213)
        newDci.m_tbsSize.push_back(tbSize);
        newDci.m_resAlloc = 0; // only allocation type 0 at this stage
        newDci.m_rbBitmap = 0; // TBD (32 bit bitmap see 7.1.6 of 36.213)
        uint32_t rbgMask = 0;
        std::multimap<uint8_t, qos_rb_and_CQI_assigned_to_lc>::iterator itRBGsPerRNTI;
        for (itRBGsPerRNTI = (*itMap).second.begin(); itRBGsPerRNTI != (*itMap).second.end();
             itRBGsPerRNTI++)
        {
            rbgMask = rbgMask + (0x1 << itRBGsPerRNTI->second.resource_block_index);
        }
        newDci.m_rbBitmap = rbgMask; // (32 bit bitmap see 7.1.6 of 36.213)
        // NOTE: In this first version of CqaFfMacScheduler, it is assumed one flow per user.
        // create the rlc PDUs -> equally divide resources among active LCs
        std::map<LteFlowId_t, FfMacSchedSapProvider::SchedDlRlcBufferReqParameters>::iterator
            itBufReq;
        for (itBufReq = m_rlcBufferReq.begin(); itBufReq != m_rlcBufferReq.end(); itBufReq++)
        {
            if (((*itBufReq).first.m_rnti == (*itMap).first) &&
                (((*itBufReq).second.m_rlcTransmissionQueueSize > 0) ||
                 ((*itBufReq).second.m_rlcRetransmissionQueueSize > 0) ||
                 ((*itBufReq).second.m_rlcStatusPduSize > 0)))
            {
                std::vector<RlcPduListElement_s> newRlcPduLe;
                // for (uint8_t j = 0; j < nLayer; j++)
                //{
                RlcPduListElement_s newRlcEl;
                newRlcEl.m_logicalChannelIdentity = (*itBufReq).first.m_lcId;
                // newRlcEl.m_size = newDci.m_tbsSize.at (j) / lcActives;
                newRlcEl.m_size = tbSize / lcActives;
                // NS_LOG_INFO (this << " LCID " << (uint32_t) newRlcEl.m_logicalChannelIdentity <<
                // " size " << newRlcEl.m_size << " layer " << (uint16_t)j);
                newRlcPduLe.push_back(newRlcEl);
                UpdateDlRlcBufferInfo(newDci.m_rnti,
                                      newRlcEl.m_logicalChannelIdentity,
                                      newRlcEl.m_size);
                if (m_harqOn)
                {
                    // store RLC PDU list for HARQ
                    std::map<uint16_t, DlHarqRlcPduListBuffer_t>::iterator itRlcPdu =
                        m_dlHarqProcessesRlcPduListBuffer.find((*itMap).first);
                    if (itRlcPdu == m_dlHarqProcessesRlcPduListBuffer.end())
                    {
                        NS_FATAL_ERROR("Unable to find RlcPdcList in HARQ buffer for RNTI "
                                       << (*itMap).first);
                    }
                    int j = 0;
                    (*itRlcPdu).second.at(j).at(newDci.m_harqProcess).push_back(newRlcEl);
                }
                // }
                newEl.m_rlcPduList.push_back(newRlcPduLe);
            }
            if ((*itBufReq).first.m_rnti > (*itMap).first)
            {
                break;
            }
        }
        // for (uint8_t j = 0; j < nLayer; j++)
        // {
        newDci.m_ndi.push_back(1);
        newDci.m_rv.push_back(0);
        //}
 
        newDci.m_tpc = m_ffrSapProvider->GetTpc((*itMap).first);
 
        newEl.m_dci = newDci;
 
        if (m_harqOn)
        {
            // store DCI for HARQ
            std::map<uint16_t, DlHarqProcessesDciBuffer_t>::iterator itDci =
                m_dlHarqProcessesDciBuffer.find(newEl.m_rnti);
            if (itDci == m_dlHarqProcessesDciBuffer.end())
            {
                NS_FATAL_ERROR("Unable to find RNTI entry in DCI HARQ buffer for RNTI "
                               << newEl.m_rnti);
            }
            (*itDci).second.at(newDci.m_harqProcess) = newDci;
            // refresh timer
            std::map<uint16_t, DlHarqProcessesTimer_t>::iterator itHarqTimer =
                m_dlHarqProcessesTimer.find(newEl.m_rnti);
            if (itHarqTimer == m_dlHarqProcessesTimer.end())
            {
                NS_FATAL_ERROR("Unable to find HARQ timer for RNTI " << (uint16_t)newEl.m_rnti);
            }
            (*itHarqTimer).second.at(newDci.m_harqProcess) = 0;
        }
 
        // ...more parameters -> ignored in this version
 
        ret.m_buildDataList.push_back(newEl);
        // update UE stats
        std::map<uint16_t, CqasFlowPerf_t>::iterator it;
        it = m_flowStatsDl.find((*itMap).first);
        if (it != m_flowStatsDl.end())
        {
            (*it).second.lastTtiBytesTransmitted = tbSize;
        }
        else
        {
            NS_FATAL_ERROR(this << " No Stats for this allocated UE");
        }
 
        itMap++;
    }                               // end while allocation
    ret.m_nrOfPdcchOfdmSymbols = 1; // TODO: check correct value according the DCIs txed
 
    // update UEs stats
    NS_LOG_INFO(this << " Update UEs statistics");
    for (itStats = m_flowStatsDl.begin(); itStats != m_flowStatsDl.end(); itStats++)
    {
        if (allocationMapPerRntiPerLCId.find(itStats->first) != allocationMapPerRntiPerLCId.end())
        {
            (*itStats).second.secondLastAveragedThroughput =
                ((1.0 - (1 / m_timeWindow)) * (*itStats).second.secondLastAveragedThroughput) +
                ((1 / m_timeWindow) * (double)((*itStats).second.lastTtiBytesTransmitted / 0.001));
        }
 
        (*itStats).second.totalBytesTransmitted += (*itStats).second.lastTtiBytesTransmitted;
        // update average throughput (see eq. 12.3 of Sec 12.3.1.2 of LTE – The UMTS Long Term
        // Evolution, Ed Wiley)
        (*itStats).second.lastAveragedThroughput =
            ((1.0 - (1.0 / m_timeWindow)) * (*itStats).second.lastAveragedThroughput) +
            ((1.0 / m_timeWindow) * (double)((*itStats).second.lastTtiBytesTransmitted / 0.001));
        NS_LOG_INFO(this << " UE total bytes " << (*itStats).second.totalBytesTransmitted);
        NS_LOG_INFO(this << " UE average throughput " << (*itStats).second.lastAveragedThroughput);
        (*itStats).second.lastTtiBytesTransmitted = 0;
    }
 
    m_schedSapUser->SchedDlConfigInd(ret);
 
    int count_allocated_resource_blocks = 0;
    for (std::map<uint16_t, std::multimap<uint8_t, qos_rb_and_CQI_assigned_to_lc>>::iterator itMap =
             allocationMapPerRntiPerLCId.begin();
         itMap != allocationMapPerRntiPerLCId.end();
         itMap++)
    {
        count_allocated_resource_blocks += itMap->second.size();
    }
    NS_LOG_INFO(this << " Allocated RBs:" << count_allocated_resource_blocks);
}
 
void
CqaFfMacScheduler::DoSchedDlRachInfoReq(
    const FfMacSchedSapProvider::SchedDlRachInfoReqParameters& params)
{
    NS_LOG_FUNCTION(this);
 
    m_rachList = params.m_rachList;
}
 
void
CqaFfMacScheduler::DoSchedDlCqiInfoReq(
    const FfMacSchedSapProvider::SchedDlCqiInfoReqParameters& params)
{
    NS_LOG_FUNCTION(this);
    m_ffrSapProvider->ReportDlCqiInfo(params);
 
    for (unsigned int i = 0; i < params.m_cqiList.size(); i++)
    {
        if (params.m_cqiList.at(i).m_cqiType == CqiListElement_s::P10)
        {
            NS_LOG_LOGIC("wideband CQI " << (uint32_t)params.m_cqiList.at(i).m_wbCqi.at(0)
                                         << " reported");
            std::map<uint16_t, uint8_t>::iterator it;
            uint16_t rnti = params.m_cqiList.at(i).m_rnti;
            it = m_p10CqiRxed.find(rnti);
            if (it == m_p10CqiRxed.end())
            {
                // create the new entry
                m_p10CqiRxed.insert(std::pair<uint16_t, uint8_t>(
                    rnti,
                    params.m_cqiList.at(i).m_wbCqi.at(0))); // only codeword 0 at this stage (SISO)
                // generate correspondent timer
                m_p10CqiTimers.insert(std::pair<uint16_t, uint32_t>(rnti, m_cqiTimersThreshold));
            }
            else
            {
                // update the CQI value and refresh correspondent timer
                (*it).second = params.m_cqiList.at(i).m_wbCqi.at(0);
                // update correspondent timer
                std::map<uint16_t, uint32_t>::iterator itTimers;
                itTimers = m_p10CqiTimers.find(rnti);
                (*itTimers).second = m_cqiTimersThreshold;
            }
        }
        else if (params.m_cqiList.at(i).m_cqiType == CqiListElement_s::A30)
        {
            // subband CQI reporting high layer configured
            std::map<uint16_t, SbMeasResult_s>::iterator it;
            uint16_t rnti = params.m_cqiList.at(i).m_rnti;
            it = m_a30CqiRxed.find(rnti);
            if (it == m_a30CqiRxed.end())
            {
                // create the new entry
                m_a30CqiRxed.insert(
                    std::pair<uint16_t, SbMeasResult_s>(rnti,
                                                        params.m_cqiList.at(i).m_sbMeasResult));
                m_a30CqiTimers.insert(std::pair<uint16_t, uint32_t>(rnti, m_cqiTimersThreshold));
            }
            else
            {
                // update the CQI value and refresh correspondent timer
                (*it).second = params.m_cqiList.at(i).m_sbMeasResult;
                std::map<uint16_t, uint32_t>::iterator itTimers;
                itTimers = m_a30CqiTimers.find(rnti);
                (*itTimers).second = m_cqiTimersThreshold;
            }
        }
        else
        {
            NS_LOG_ERROR(this << " CQI type unknown");
        }
    }
}
 
double
CqaFfMacScheduler::EstimateUlSinr(uint16_t rnti, uint16_t rb)
{
    std::map<uint16_t, std::vector<double>>::iterator itCqi = m_ueCqi.find(rnti);
    if (itCqi == m_ueCqi.end())
    {
        // no cqi info about this UE
        return (NO_SINR);
    }
    else
    {
        // take the average SINR value among the available
        double sinrSum = 0;
        unsigned int sinrNum = 0;
        for (uint32_t i = 0; i < m_cschedCellConfig.m_ulBandwidth; i++)
        {
            double sinr = (*itCqi).second.at(i);
            if (sinr != NO_SINR)
            {
                sinrSum += sinr;
                sinrNum++;
            }
        }
        double estimatedSinr = (sinrNum > 0) ? (sinrSum / sinrNum) : DBL_MAX;
        // store the value
        (*itCqi).second.at(rb) = estimatedSinr;
        return (estimatedSinr);
    }
}
 
void
CqaFfMacScheduler::DoSchedUlTriggerReq(
    const FfMacSchedSapProvider::SchedUlTriggerReqParameters& params)
{
    NS_LOG_FUNCTION(this << " UL - Frame no. " << (params.m_sfnSf >> 4) << " subframe no. "
                         << (0xF & params.m_sfnSf) << " size " << params.m_ulInfoList.size());
 
    RefreshUlCqiMaps();
    m_ffrSapProvider->ReportUlCqiInfo(m_ueCqi);
 
    // Generate RBs map
    FfMacSchedSapUser::SchedUlConfigIndParameters ret;
    std::vector<bool> rbMap;
    uint16_t rbAllocatedNum = 0;
    std::set<uint16_t> rntiAllocated;
    std::vector<uint16_t> rbgAllocationMap;
    // update with RACH allocation map
    rbgAllocationMap = m_rachAllocationMap;
    // rbgAllocationMap.resize (m_cschedCellConfig.m_ulBandwidth, 0);
    m_rachAllocationMap.clear();
    m_rachAllocationMap.resize(m_cschedCellConfig.m_ulBandwidth, 0);
 
    rbMap.resize(m_cschedCellConfig.m_ulBandwidth, false);
 
    rbMap = m_ffrSapProvider->GetAvailableUlRbg();
 
    for (std::vector<bool>::iterator it = rbMap.begin(); it != rbMap.end(); it++)
    {
        if (*it)
        {
            rbAllocatedNum++;
        }
    }
 
    uint8_t minContinuousUlBandwidth = m_ffrSapProvider->GetMinContinuousUlBandwidth();
    uint8_t ffrUlBandwidth = m_cschedCellConfig.m_ulBandwidth - rbAllocatedNum;
 
    // remove RACH allocation
    for (uint16_t i = 0; i < m_cschedCellConfig.m_ulBandwidth; i++)
    {
        if (rbgAllocationMap.at(i) != 0)
        {
            rbMap.at(i) = true;
            NS_LOG_DEBUG(this << " Allocated for RACH " << i);
        }
    }
 
    if (m_harqOn)
    {
        //   Process UL HARQ feedback
        for (std::size_t i = 0; i < params.m_ulInfoList.size(); i++)
        {
            if (params.m_ulInfoList.at(i).m_receptionStatus == UlInfoListElement_s::NotOk)
            {
                // retx correspondent block: retrieve the UL-DCI
                uint16_t rnti = params.m_ulInfoList.at(i).m_rnti;
                std::map<uint16_t, uint8_t>::iterator itProcId =
                    m_ulHarqCurrentProcessId.find(rnti);
                if (itProcId == m_ulHarqCurrentProcessId.end())
                {
                    NS_LOG_ERROR("No info find in HARQ buffer for UE (might change eNB) " << rnti);
                }
                uint8_t harqId = (uint8_t)((*itProcId).second - HARQ_PERIOD) % HARQ_PROC_NUM;
                NS_LOG_INFO(this << " UL-HARQ retx RNTI " << rnti << " harqId " << (uint16_t)harqId
                                 << " i " << i << " size " << params.m_ulInfoList.size());
                std::map<uint16_t, UlHarqProcessesDciBuffer_t>::iterator itHarq =
                    m_ulHarqProcessesDciBuffer.find(rnti);
                if (itHarq == m_ulHarqProcessesDciBuffer.end())
                {
                    NS_LOG_ERROR("No info find in HARQ buffer for UE (might change eNB) " << rnti);
                    continue;
                }
                UlDciListElement_s dci = (*itHarq).second.at(harqId);
                std::map<uint16_t, UlHarqProcessesStatus_t>::iterator itStat =
                    m_ulHarqProcessesStatus.find(rnti);
                if (itStat == m_ulHarqProcessesStatus.end())
                {
                    NS_LOG_ERROR("No info find in HARQ buffer for UE (might change eNB) " << rnti);
                }
                if ((*itStat).second.at(harqId) >= 3)
                {
                    NS_LOG_INFO("Max number of retransmissions reached (UL)-> drop process");
                    continue;
                }
                bool free = true;
 
                for (int j = dci.m_rbStart; j < dci.m_rbStart + dci.m_rbLen; j++)
                {
                    if (rbMap.at(j))
                    {
                        free = false;
                        NS_LOG_INFO(this << " BUSY " << j);
                    }
                }
                if (free)
                {
                    // retx on the same RBs
                    for (int j = dci.m_rbStart; j < dci.m_rbStart + dci.m_rbLen; j++)
                    {
                        rbMap.at(j) = true;
                        rbgAllocationMap.at(j) = dci.m_rnti;
                        NS_LOG_INFO("\tRB " << j);
                        rbAllocatedNum++;
                    }
                    NS_LOG_INFO(this << " Send retx in the same RBs " << (uint16_t)dci.m_rbStart
                                     << " to " << dci.m_rbStart + dci.m_rbLen << " RV "
                                     << (*itStat).second.at(harqId) + 1);
                }
                else
                {
                    NS_LOG_INFO("Cannot allocate retx due to RACH allocations for UE " << rnti);
                    continue;
                }
                dci.m_ndi = 0;
                // Update HARQ buffers with new HarqId
                (*itStat).second.at((*itProcId).second) = (*itStat).second.at(harqId) + 1;
                (*itStat).second.at(harqId) = 0;
                (*itHarq).second.at((*itProcId).second) = dci;
                ret.m_dciList.push_back(dci);
                rntiAllocated.insert(dci.m_rnti);
            }
            else
            {
                NS_LOG_INFO(this << " HARQ-ACK feedback from RNTI "
                                 << params.m_ulInfoList.at(i).m_rnti);
            }
        }
    }
 
    std::map<uint16_t, uint32_t>::iterator it;
    int nflows = 0;
 
    for (it = m_ceBsrRxed.begin(); it != m_ceBsrRxed.end(); it++)
    {
        std::set<uint16_t>::iterator itRnti = rntiAllocated.find((*it).first);
        // select UEs with queues not empty and not yet allocated for HARQ
        if (((*it).second > 0) && (itRnti == rntiAllocated.end()))
        {
            nflows++;
        }
    }
 
    if (nflows == 0)
    {
        if (!ret.m_dciList.empty())
        {
            m_allocationMaps.insert(
                std::pair<uint16_t, std::vector<uint16_t>>(params.m_sfnSf, rbgAllocationMap));
            m_schedSapUser->SchedUlConfigInd(ret);
        }
 
        return; // no flows to be scheduled
    }
 
    // Divide the remaining resources equally among the active users starting from the subsequent
    // one served last scheduling trigger
    uint16_t tempRbPerFlow = (ffrUlBandwidth) / (nflows + rntiAllocated.size());
    uint16_t rbPerFlow =
        (minContinuousUlBandwidth < tempRbPerFlow) ? minContinuousUlBandwidth : tempRbPerFlow;
 
    if (rbPerFlow < 3)
    {
        rbPerFlow = 3; // at least 3 rbg per flow (till available resource) to ensure TxOpportunity
                       // >= 7 bytes
    }
    int rbAllocated = 0;
 
    std::map<uint16_t, CqasFlowPerf_t>::iterator itStats;
    if (m_nextRntiUl != 0)
    {
        for (it = m_ceBsrRxed.begin(); it != m_ceBsrRxed.end(); it++)
        {
            if ((*it).first == m_nextRntiUl)
            {
                break;
            }
        }
        if (it == m_ceBsrRxed.end())
        {
            NS_LOG_ERROR(this << " no user found");
        }
    }
    else
    {
        it = m_ceBsrRxed.begin();
        m_nextRntiUl = (*it).first;
    }
    do
    {
        std::set<uint16_t>::iterator itRnti = rntiAllocated.find((*it).first);
        if ((itRnti != rntiAllocated.end()) || ((*it).second == 0))
        {
            // UE already allocated for UL-HARQ -> skip it
            NS_LOG_DEBUG(this << " UE already allocated in HARQ -> discarded, RNTI "
                              << (*it).first);
            it++;
            if (it == m_ceBsrRxed.end())
            {
                // restart from the first
                it = m_ceBsrRxed.begin();
            }
            continue;
        }
        if (rbAllocated + rbPerFlow - 1 > m_cschedCellConfig.m_ulBandwidth)
        {
            // limit to physical resources last resource assignment
            rbPerFlow = m_cschedCellConfig.m_ulBandwidth - rbAllocated;
            // at least 3 rbg per flow to ensure TxOpportunity >= 7 bytes
            if (rbPerFlow < 3)
            {
                // terminate allocation
                rbPerFlow = 0;
            }
        }
 
        rbAllocated = 0;
        UlDciListElement_s uldci;
        uldci.m_rnti = (*it).first;
        uldci.m_rbLen = rbPerFlow;
        bool allocated = false;
        NS_LOG_INFO(this << " RB Allocated " << rbAllocated << " rbPerFlow " << rbPerFlow
                         << " flows " << nflows);
        while ((!allocated) && ((rbAllocated + rbPerFlow - m_cschedCellConfig.m_ulBandwidth) < 1) &&
               (rbPerFlow != 0))
        {
            // check availability
            bool free = true;
            for (int j = rbAllocated; j < rbAllocated + rbPerFlow; j++)
            {
                if (rbMap.at(j))
                {
                    free = false;
                    break;
                }
                if (!m_ffrSapProvider->IsUlRbgAvailableForUe(j, (*it).first))
                {
                    free = false;
                    break;
                }
            }
            if (free)
            {
                NS_LOG_INFO(this << "RNTI: " << (*it).first << " RB Allocated " << rbAllocated
                                 << " rbPerFlow " << rbPerFlow << " flows " << nflows);
                uldci.m_rbStart = rbAllocated;
 
                for (int j = rbAllocated; j < rbAllocated + rbPerFlow; j++)
                {
                    rbMap.at(j) = true;
                    // store info on allocation for managing ul-cqi interpretation
                    rbgAllocationMap.at(j) = (*it).first;
                }
                rbAllocated += rbPerFlow;
                allocated = true;
                break;
            }
            rbAllocated++;
            if (rbAllocated + rbPerFlow - 1 > m_cschedCellConfig.m_ulBandwidth)
            {
                // limit to physical resources last resource assignment
                rbPerFlow = m_cschedCellConfig.m_ulBandwidth - rbAllocated;
                // at least 3 rbg per flow to ensure TxOpportunity >= 7 bytes
                if (rbPerFlow < 3)
                {
                    // terminate allocation
                    rbPerFlow = 0;
                }
            }
        }
        if (!allocated)
        {
            // unable to allocate new resource: finish scheduling
            //          m_nextRntiUl = (*it).first;
            //          if (ret.m_dciList.size () > 0)
            //            {
            //              m_schedSapUser->SchedUlConfigInd (ret);
            //            }
            //          m_allocationMaps.insert (std::pair <uint16_t, std::vector <uint16_t> >
            //          (params.m_sfnSf, rbgAllocationMap)); return;
            break;
        }
 
        std::map<uint16_t, std::vector<double>>::iterator itCqi = m_ueCqi.find((*it).first);
        int cqi = 0;
        if (itCqi == m_ueCqi.end())
        {
            // no cqi info about this UE
            uldci.m_mcs = 0; // MCS 0 -> UL-AMC TBD
        }
        else
        {
            // take the lowest CQI value (worst RB)
            NS_ABORT_MSG_IF((*itCqi).second.empty(),
                            "CQI of RNTI = " << (*it).first << " has expired");
            double minSinr = (*itCqi).second.at(uldci.m_rbStart);
            if (minSinr == NO_SINR)
            {
                minSinr = EstimateUlSinr((*it).first, uldci.m_rbStart);
            }
            for (uint16_t i = uldci.m_rbStart; i < uldci.m_rbStart + uldci.m_rbLen; i++)
            {
                double sinr = (*itCqi).second.at(i);
                if (sinr == NO_SINR)
                {
                    sinr = EstimateUlSinr((*it).first, i);
                }
                if (sinr < minSinr)
                {
                    minSinr = sinr;
                }
            }
 
            // translate SINR -> cqi: WILD ACK: same as DL
            double s = log2(1 + (std::pow(10, minSinr / 10) / ((-std::log(5.0 * 0.00005)) / 1.5)));
            cqi = m_amc->GetCqiFromSpectralEfficiency(s);
            if (cqi == 0)
            {
                it++;
                if (it == m_ceBsrRxed.end())
                {
                    // restart from the first
                    it = m_ceBsrRxed.begin();
                }
                NS_LOG_DEBUG(this << " UE discarded for CQI = 0, RNTI " << uldci.m_rnti);
                // remove UE from allocation map
                for (uint16_t i = uldci.m_rbStart; i < uldci.m_rbStart + uldci.m_rbLen; i++)
                {
                    rbgAllocationMap.at(i) = 0;
                }
                continue; // CQI == 0 means "out of range" (see table 7.2.3-1 of 36.213)
            }
            uldci.m_mcs = m_amc->GetMcsFromCqi(cqi);
        }
 
        uldci.m_tbSize = (m_amc->GetUlTbSizeFromMcs(uldci.m_mcs, rbPerFlow) / 8);
        UpdateUlRlcBufferInfo(uldci.m_rnti, uldci.m_tbSize);
        uldci.m_ndi = 1;
        uldci.m_cceIndex = 0;
        uldci.m_aggrLevel = 1;
        uldci.m_ueTxAntennaSelection = 3; // antenna selection OFF
        uldci.m_hopping = false;
        uldci.m_n2Dmrs = 0;
        uldci.m_tpc = 0;            // no power control
        uldci.m_cqiRequest = false; // only period CQI at this stage
        uldci.m_ulIndex = 0;        // TDD parameter
        uldci.m_dai = 1;            // TDD parameter
        uldci.m_freqHopping = 0;
        uldci.m_pdcchPowerOffset = 0; // not used
        ret.m_dciList.push_back(uldci);
        // store DCI for HARQ_PERIOD
        uint8_t harqId = 0;
        if (m_harqOn)
        {
            std::map<uint16_t, uint8_t>::iterator itProcId;
            itProcId = m_ulHarqCurrentProcessId.find(uldci.m_rnti);
            if (itProcId == m_ulHarqCurrentProcessId.end())
            {
                NS_FATAL_ERROR("No info find in HARQ buffer for UE " << uldci.m_rnti);
            }
            harqId = (*itProcId).second;
            std::map<uint16_t, UlHarqProcessesDciBuffer_t>::iterator itDci =
                m_ulHarqProcessesDciBuffer.find(uldci.m_rnti);
            if (itDci == m_ulHarqProcessesDciBuffer.end())
            {
                NS_FATAL_ERROR("Unable to find RNTI entry in UL DCI HARQ buffer for RNTI "
                               << uldci.m_rnti);
            }
            (*itDci).second.at(harqId) = uldci;
            // Update HARQ process status (RV 0)
            std::map<uint16_t, UlHarqProcessesStatus_t>::iterator itStat =
                m_ulHarqProcessesStatus.find(uldci.m_rnti);
            if (itStat == m_ulHarqProcessesStatus.end())
            {
                NS_LOG_ERROR("No info find in HARQ buffer for UE (might change eNB) "
                             << uldci.m_rnti);
            }
            (*itStat).second.at(harqId) = 0;
        }
 
        NS_LOG_INFO(this << " UE Allocation RNTI " << (*it).first << " startPRB "
                         << (uint32_t)uldci.m_rbStart << " nPRB " << (uint32_t)uldci.m_rbLen
                         << " CQI " << cqi << " MCS " << (uint32_t)uldci.m_mcs << " TBsize "
                         << uldci.m_tbSize << " RbAlloc " << rbAllocated << " harqId "
                         << (uint16_t)harqId);
 
        // update TTI  UE stats
        itStats = m_flowStatsUl.find((*it).first);
        if (itStats != m_flowStatsUl.end())
        {
            (*itStats).second.lastTtiBytesTransmitted = uldci.m_tbSize;
        }
        else
        {
            NS_LOG_DEBUG(this << " No Stats for this allocated UE");
        }
 
        it++;
        if (it == m_ceBsrRxed.end())
        {
            // restart from the first
            it = m_ceBsrRxed.begin();
        }
        if ((rbAllocated == m_cschedCellConfig.m_ulBandwidth) || (rbPerFlow == 0))
        {
            // Stop allocation: no more PRBs
            m_nextRntiUl = (*it).first;
            break;
        }
    } while (((*it).first != m_nextRntiUl) && (rbPerFlow != 0));
 
    // Update global UE stats
    // update UEs stats
    for (itStats = m_flowStatsUl.begin(); itStats != m_flowStatsUl.end(); itStats++)
    {
        (*itStats).second.totalBytesTransmitted += (*itStats).second.lastTtiBytesTransmitted;
        // update average throughput (see eq. 12.3 of Sec 12.3.1.2 of LTE – The UMTS Long Term
        // Evolution, Ed Wiley)
        (*itStats).second.lastAveragedThroughput =
            ((1.0 - (1.0 / m_timeWindow)) * (*itStats).second.lastAveragedThroughput) +
            ((1.0 / m_timeWindow) * (double)((*itStats).second.lastTtiBytesTransmitted / 0.001));
        NS_LOG_INFO(this << " UE total bytes " << (*itStats).second.totalBytesTransmitted);
        NS_LOG_INFO(this << " UE average throughput " << (*itStats).second.lastAveragedThroughput);
        (*itStats).second.lastTtiBytesTransmitted = 0;
    }
    m_allocationMaps.insert(
        std::pair<uint16_t, std::vector<uint16_t>>(params.m_sfnSf, rbgAllocationMap));
    m_schedSapUser->SchedUlConfigInd(ret);
}
 
void
CqaFfMacScheduler::DoSchedUlNoiseInterferenceReq(
    const FfMacSchedSapProvider::SchedUlNoiseInterferenceReqParameters& params)
{
    NS_LOG_FUNCTION(this);
}
 
void
CqaFfMacScheduler::DoSchedUlSrInfoReq(
    const FfMacSchedSapProvider::SchedUlSrInfoReqParameters& params)
{
    NS_LOG_FUNCTION(this);
}
 
void
CqaFfMacScheduler::DoSchedUlMacCtrlInfoReq(
    const FfMacSchedSapProvider::SchedUlMacCtrlInfoReqParameters& params)
{
    NS_LOG_FUNCTION(this);
 
    std::map<uint16_t, uint32_t>::iterator it;
 
    for (unsigned int i = 0; i < params.m_macCeList.size(); i++)
    {
        if (params.m_macCeList.at(i).m_macCeType == MacCeListElement_s::BSR)
        {
            // buffer status report
            // note that this scheduler does not differentiate the
            // allocation according to which LCGs have more/less bytes
            // to send.
            // Hence the BSR of different LCGs are just summed up to get
            // a total queue size that is used for allocation purposes.
 
            uint32_t buffer = 0;
            for (uint8_t lcg = 0; lcg < 4; ++lcg)
            {
                uint8_t bsrId = params.m_macCeList.at(i).m_macCeValue.m_bufferStatus.at(lcg);
                buffer += BufferSizeLevelBsr::BsrId2BufferSize(bsrId);
            }
 
            uint16_t rnti = params.m_macCeList.at(i).m_rnti;
            NS_LOG_LOGIC(this << "RNTI=" << rnti << " buffer=" << buffer);
            it = m_ceBsrRxed.find(rnti);
            if (it == m_ceBsrRxed.end())
            {
                // create the new entry
                m_ceBsrRxed.insert(std::pair<uint16_t, uint32_t>(rnti, buffer));
            }
            else
            {
                // update the buffer size value
                (*it).second = buffer;
            }
        }
    }
}
 
void
CqaFfMacScheduler::DoSchedUlCqiInfoReq(
    const FfMacSchedSapProvider::SchedUlCqiInfoReqParameters& params)
{
    NS_LOG_FUNCTION(this);
    // retrieve the allocation for this subframe
    switch (m_ulCqiFilter)
    {
    case FfMacScheduler::SRS_UL_CQI: {
        // filter all the CQIs that are not SRS based
        if (params.m_ulCqi.m_type != UlCqi_s::SRS)
        {
            return;
        }
    }
    break;
    case FfMacScheduler::PUSCH_UL_CQI: {
        // filter all the CQIs that are not SRS based
        if (params.m_ulCqi.m_type != UlCqi_s::PUSCH)
        {
            return;
        }
    }
    break;
    default:
        NS_FATAL_ERROR("Unknown UL CQI type");
    }
 
    switch (params.m_ulCqi.m_type)
    {
    case UlCqi_s::PUSCH: {
        std::map<uint16_t, std::vector<uint16_t>>::iterator itMap;
        std::map<uint16_t, std::vector<double>>::iterator itCqi;
        NS_LOG_DEBUG(this << " Collect PUSCH CQIs of Frame no. " << (params.m_sfnSf >> 4)
                          << " subframe no. " << (0xF & params.m_sfnSf));
        itMap = m_allocationMaps.find(params.m_sfnSf);
        if (itMap == m_allocationMaps.end())
        {
            return;
        }
        for (uint32_t i = 0; i < (*itMap).second.size(); i++)
        {
            // convert from fixed point notation Sxxxxxxxxxxx.xxx to double
            double sinr = LteFfConverter::fpS11dot3toDouble(params.m_ulCqi.m_sinr.at(i));
            itCqi = m_ueCqi.find((*itMap).second.at(i));
            if (itCqi == m_ueCqi.end())
            {
                // create a new entry
                std::vector<double> newCqi;
                for (uint32_t j = 0; j < m_cschedCellConfig.m_ulBandwidth; j++)
                {
                    if (i == j)
                    {
                        newCqi.push_back(sinr);
                    }
                    else
                    {
                        // initialize with NO_SINR value.
                        newCqi.push_back(NO_SINR);
                    }
                }
                m_ueCqi.insert(
                    std::pair<uint16_t, std::vector<double>>((*itMap).second.at(i), newCqi));
                // generate correspondent timer
                m_ueCqiTimers.insert(
                    std::pair<uint16_t, uint32_t>((*itMap).second.at(i), m_cqiTimersThreshold));
            }
            else
            {
                // update the value
                (*itCqi).second.at(i) = sinr;
                NS_LOG_DEBUG(this << " RNTI " << (*itMap).second.at(i) << " RB " << i << " SINR "
                                  << sinr);
                // update correspondent timer
                std::map<uint16_t, uint32_t>::iterator itTimers;
                itTimers = m_ueCqiTimers.find((*itMap).second.at(i));
                (*itTimers).second = m_cqiTimersThreshold;
            }
        }
        // remove obsolete info on allocation
        m_allocationMaps.erase(itMap);
    }
    break;
    case UlCqi_s::SRS: {
        NS_LOG_DEBUG(this << " Collect SRS CQIs of Frame no. " << (params.m_sfnSf >> 4)
                          << " subframe no. " << (0xF & params.m_sfnSf));
        // get the RNTI from vendor specific parameters
        uint16_t rnti = 0;
        NS_ASSERT(!params.m_vendorSpecificList.empty());
        for (std::size_t i = 0; i < params.m_vendorSpecificList.size(); i++)
        {
            if (params.m_vendorSpecificList.at(i).m_type == SRS_CQI_RNTI_VSP)
            {
                Ptr<SrsCqiRntiVsp> vsp =
                    DynamicCast<SrsCqiRntiVsp>(params.m_vendorSpecificList.at(i).m_value);
                rnti = vsp->GetRnti();
            }
        }
        std::map<uint16_t, std::vector<double>>::iterator itCqi;
        itCqi = m_ueCqi.find(rnti);
        if (itCqi == m_ueCqi.end())
        {
            // create a new entry
            std::vector<double> newCqi;
            for (uint32_t j = 0; j < m_cschedCellConfig.m_ulBandwidth; j++)
            {
                double sinr = LteFfConverter::fpS11dot3toDouble(params.m_ulCqi.m_sinr.at(j));
                newCqi.push_back(sinr);
                NS_LOG_INFO(this << " RNTI " << rnti << " new SRS-CQI for RB  " << j << " value "
                                 << sinr);
            }
            m_ueCqi.insert(std::pair<uint16_t, std::vector<double>>(rnti, newCqi));
            // generate correspondent timer
            m_ueCqiTimers.insert(std::pair<uint16_t, uint32_t>(rnti, m_cqiTimersThreshold));
        }
        else
        {
            // update the values
            for (uint32_t j = 0; j < m_cschedCellConfig.m_ulBandwidth; j++)
            {
                double sinr = LteFfConverter::fpS11dot3toDouble(params.m_ulCqi.m_sinr.at(j));
                (*itCqi).second.at(j) = sinr;
                NS_LOG_INFO(this << " RNTI " << rnti << " update SRS-CQI for RB  " << j << " value "
                                 << sinr);
            }
            // update correspondent timer
            std::map<uint16_t, uint32_t>::iterator itTimers;
            itTimers = m_ueCqiTimers.find(rnti);
            (*itTimers).second = m_cqiTimersThreshold;
        }
    }
    break;
    case UlCqi_s::PUCCH_1:
    case UlCqi_s::PUCCH_2:
    case UlCqi_s::PRACH: {
        NS_FATAL_ERROR("PfFfMacScheduler supports only PUSCH and SRS UL-CQIs");
    }
    break;
    default:
        NS_FATAL_ERROR("Unknown type of UL-CQI");
    }
}
 
void
CqaFfMacScheduler::RefreshDlCqiMaps()
{
    // refresh DL CQI P01 Map
    std::map<uint16_t, uint32_t>::iterator itP10 = m_p10CqiTimers.begin();
    while (itP10 != m_p10CqiTimers.end())
    {
        NS_LOG_INFO(this << " P10-CQI for user " << (*itP10).first << " is "
                         << (uint32_t)(*itP10).second << " thr " << (uint32_t)m_cqiTimersThreshold);
        if ((*itP10).second == 0)
        {
            // delete correspondent entries
            std::map<uint16_t, uint8_t>::iterator itMap = m_p10CqiRxed.find((*itP10).first);
            NS_ASSERT_MSG(itMap != m_p10CqiRxed.end(),
                          " Does not find CQI report for user " << (*itP10).first);
            NS_LOG_INFO(this << " P10-CQI expired for user " << (*itP10).first);
            m_p10CqiRxed.erase(itMap);
            std::map<uint16_t, uint32_t>::iterator temp = itP10;
            itP10++;
            m_p10CqiTimers.erase(temp);
        }
        else
        {
            (*itP10).second--;
            itP10++;
        }
    }
 
    // refresh DL CQI A30 Map
    std::map<uint16_t, uint32_t>::iterator itA30 = m_a30CqiTimers.begin();
    while (itA30 != m_a30CqiTimers.end())
    {
        NS_LOG_INFO(this << " A30-CQI for user " << (*itA30).first << " is "
                         << (uint32_t)(*itA30).second << " thr " << (uint32_t)m_cqiTimersThreshold);
        if ((*itA30).second == 0)
        {
            // delete correspondent entries
            std::map<uint16_t, SbMeasResult_s>::iterator itMap = m_a30CqiRxed.find((*itA30).first);
            NS_ASSERT_MSG(itMap != m_a30CqiRxed.end(),
                          " Does not find CQI report for user " << (*itA30).first);
            NS_LOG_INFO(this << " A30-CQI expired for user " << (*itA30).first);
            m_a30CqiRxed.erase(itMap);
            std::map<uint16_t, uint32_t>::iterator temp = itA30;
            itA30++;
            m_a30CqiTimers.erase(temp);
        }
        else
        {
            (*itA30).second--;
            itA30++;
        }
    }
}
 
void
CqaFfMacScheduler::RefreshUlCqiMaps()
{
    // refresh UL CQI  Map
    std::map<uint16_t, uint32_t>::iterator itUl = m_ueCqiTimers.begin();
    while (itUl != m_ueCqiTimers.end())
    {
        NS_LOG_INFO(this << " UL-CQI for user " << (*itUl).first << " is "
                         << (uint32_t)(*itUl).second << " thr " << (uint32_t)m_cqiTimersThreshold);
        if ((*itUl).second == 0)
        {
            // delete correspondent entries
            std::map<uint16_t, std::vector<double>>::iterator itMap = m_ueCqi.find((*itUl).first);
            NS_ASSERT_MSG(itMap != m_ueCqi.end(),
                          " Does not find CQI report for user " << (*itUl).first);
            NS_LOG_INFO(this << " UL-CQI exired for user " << (*itUl).first);
            (*itMap).second.clear();
            m_ueCqi.erase(itMap);
            std::map<uint16_t, uint32_t>::iterator temp = itUl;
            itUl++;
            m_ueCqiTimers.erase(temp);
        }
        else
        {
            (*itUl).second--;
            itUl++;
        }
    }
}
 
void
CqaFfMacScheduler::UpdateDlRlcBufferInfo(uint16_t rnti, uint8_t lcid, uint16_t size)
{
    std::map<LteFlowId_t, FfMacSchedSapProvider::SchedDlRlcBufferReqParameters>::iterator it;
    LteFlowId_t flow(rnti, lcid);
    it = m_rlcBufferReq.find(flow);
    if (it != m_rlcBufferReq.end())
    {
        NS_LOG_INFO(this << " UE " << rnti << " LC " << (uint16_t)lcid << " txqueue "
                         << (*it).second.m_rlcTransmissionQueueSize << " retxqueue "
                         << (*it).second.m_rlcRetransmissionQueueSize << " status "
                         << (*it).second.m_rlcStatusPduSize << " decrease " << size);
        // Update queues: RLC tx order Status, ReTx, Tx
        // Update status queue
        if (((*it).second.m_rlcStatusPduSize > 0) && (size >= (*it).second.m_rlcStatusPduSize))
        {
            (*it).second.m_rlcStatusPduSize = 0;
        }
        else if (((*it).second.m_rlcRetransmissionQueueSize > 0) &&
                 (size >= (*it).second.m_rlcRetransmissionQueueSize))
        {
            (*it).second.m_rlcRetransmissionQueueSize = 0;
        }
        else if ((*it).second.m_rlcTransmissionQueueSize > 0)
        {
            uint32_t rlcOverhead;
            if (lcid == 1)
            {
                // for SRB1 (using RLC AM) it's better to
                // overestimate RLC overhead rather than
                // underestimate it and risk unneeded
                // segmentation which increases delay
                rlcOverhead = 4;
            }
            else
            {
                // minimum RLC overhead due to header
                rlcOverhead = 2;
            }
            // update transmission queue
            if ((*it).second.m_rlcTransmissionQueueSize <= size - rlcOverhead)
            {
                (*it).second.m_rlcTransmissionQueueSize = 0;
            }
            else
            {
                (*it).second.m_rlcTransmissionQueueSize -= size - rlcOverhead;
            }
        }
    }
    else
    {
        NS_LOG_ERROR(this << " Does not find DL RLC Buffer Report of UE " << rnti);
    }
}
 
void
CqaFfMacScheduler::UpdateUlRlcBufferInfo(uint16_t rnti, uint16_t size)
{
    size = size - 2; // remove the minimum RLC overhead
    std::map<uint16_t, uint32_t>::iterator it = m_ceBsrRxed.find(rnti);
    if (it != m_ceBsrRxed.end())
    {
        NS_LOG_INFO(this << " UE " << rnti << " size " << size << " BSR " << (*it).second);
        if ((*it).second >= size)
        {
            (*it).second -= size;
        }
        else
        {
            (*it).second = 0;
        }
    }
    else
    {
        NS_LOG_ERROR(this << " Does not find BSR report info of UE " << rnti);
    }
}
 
void
CqaFfMacScheduler::TransmissionModeConfigurationUpdate(uint16_t rnti, uint8_t txMode)
{
    NS_LOG_FUNCTION(this << " RNTI " << rnti << " txMode " << (uint16_t)txMode);
    FfMacCschedSapUser::CschedUeConfigUpdateIndParameters params;
    params.m_rnti = rnti;
    params.m_transmissionMode = txMode;
    m_cschedSapUser->CschedUeConfigUpdateInd(params);
}
 
} // namespace ns3