ulsim.c 69.1 KB
Newer Older
1
/*
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
 * Licensed to the OpenAirInterface (OAI) Software Alliance under one or more
 * contributor license agreements.  See the NOTICE file distributed with
 * this work for additional information regarding copyright ownership.
 * The OpenAirInterface Software Alliance licenses this file to You under
 * the OAI Public License, Version 1.1  (the "License"); you may not use this file
 * except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *      http://www.openairinterface.org/?page_id=698
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 *-------------------------------------------------------------------------------
 * For more information about the OpenAirInterface (OAI) Software Alliance:
 *      contact@openairinterface.org
 */
21 22

/*! \file ulsim.c
Thomas Laurent's avatar
Thomas Laurent committed
23 24 25 26 27 28 29 30
  \brief Top-level UL simulator
  \author R. Knopp
  \date 2011 - 2014
  \version 0.1
  \company Eurecom
  \email: knopp@eurecom.fr
  \note
  \warning
31 32
*/

33 34 35
#include <string.h>
#include <math.h>
#include <unistd.h>
36
#include "SIMULATION/TOOLS/sim.h"
37
#include "PHY/types.h"
38 39 40 41
#include "PHY/defs_common.h"
#include "PHY/defs_eNB.h"
#include "PHY/defs_UE.h"
#include "PHY/phy_vars.h"
42

43 44 45 46 47

#include "SCHED/sched_common_vars.h"
#include "SCHED/sched_eNB.h"
#include "SCHED_UE/sched_UE.h"
#include "LAYER2/MAC/mac_vars.h"
48
#include "OCG_vars.h"
49 50 51

#include "PHY/LTE_TRANSPORT/transport_proto.h"
#include "PHY/LTE_UE_TRANSPORT/transport_proto_ue.h"
52
#include "PHY/INIT/phy_init.h"
53

54 55
#include "unitary_defs.h"

56
#include "PHY/TOOLS/lte_phy_scope.h"
57
#include "dummy_functions.c"
58

59
#include "common/config/config_load_configmodule.h"
60
double cpuf;
Thomas Laurent's avatar
Thomas Laurent committed
61 62
#define inMicroS(a) (((double)(a))/(cpu_freq_GHz*1000.0))
//#define MCS_COUNT 23//added for PHY abstraction
63
#include <openair1/SIMULATION/LTE_PHY/common_sim.h>
64 65
channel_desc_t *eNB2UE[NUMBER_OF_eNB_MAX][NUMBER_OF_UE_MAX];
channel_desc_t *UE2eNB[NUMBER_OF_UE_MAX][NUMBER_OF_eNB_MAX];
Thomas Laurent's avatar
Thomas Laurent committed
66
//Added for PHY abstractionopenair1/PHY/TOOLS/lte_phy_scope.h
67
node_desc_t *enb_data[NUMBER_OF_eNB_MAX];
68 69 70
node_desc_t *ue_data[NUMBER_OF_UE_MAX];
//double sinr_bler_map[MCS_COUNT][2][16];

71
extern uint16_t beta_ack[16],beta_ri[16],beta_cqi[16];
72
//extern  char* namepointer_chMag ;
73
int xforms=0;
74 75 76
FD_lte_phy_scope_enb *form_enb;
char title[255];

77 78 79 80 81 82 83
/*the following parameters are used to control the processing times*/
double t_tx_max = -1000000000; /*!< \brief initial max process time for tx */
double t_rx_max = -1000000000; /*!< \brief initial max process time for rx */
double t_tx_min = 1000000000; /*!< \brief initial min process time for tx */
double t_rx_min = 1000000000; /*!< \brief initial min process time for tx */
int n_tx_dropped = 0; /*!< \brief initial max process time for tx */
int n_rx_dropped = 0; /*!< \brief initial max process time for rx */
84 85 86 87
int nfapi_mode = 0;

extern void fep_full(RU_t *ru);
extern void ru_fep_full_2thread(RU_t *ru);
Thomas Laurent's avatar
Thomas Laurent committed
88 89
extern void eNB_fep_full(PHY_VARS_eNB *eNB,eNB_rxtx_proc_t *proc);
extern void eNB_fep_full_2thread(PHY_VARS_eNB *eNB,eNB_rxtx_proc_t *proc);
90 91 92 93

nfapi_dl_config_request_t DL_req;
nfapi_ul_config_request_t UL_req;
nfapi_hi_dci0_request_t HI_DCI0_req;
Cedric Roux's avatar
Cedric Roux committed
94
nfapi_ul_config_request_pdu_t ul_config_pdu_list[MAX_NUM_DL_PDU];
95 96 97 98
nfapi_tx_request_pdu_t tx_pdu_list[MAX_NUM_TX_REQUEST_PDU];
nfapi_tx_request_t TX_req;
Sched_Rsp_t sched_resp;

99 100
char *parallel_config = NULL;
char *worker_config = NULL;
101
static THREAD_STRUCT thread_struct;
102
void set_parallel_conf(char *parallel_conf)
103
{
104 105 106 107
  if(strcmp(parallel_conf,"PARALLEL_SINGLE_THREAD")==0)           thread_struct.parallel_conf = PARALLEL_SINGLE_THREAD;
  else if(strcmp(parallel_conf,"PARALLEL_RU_L1_SPLIT")==0)        thread_struct.parallel_conf = PARALLEL_RU_L1_SPLIT;
  else if(strcmp(parallel_conf,"PARALLEL_RU_L1_TRX_SPLIT")==0)    thread_struct.parallel_conf = PARALLEL_RU_L1_TRX_SPLIT;
  printf("[CONFIG] parallel conf is set to %d\n",thread_struct.parallel_conf);
108
} 
109
void set_worker_conf(char *worker_conf)
110
{
111 112 113
  if(strcmp(worker_conf,"WORKER_DISABLE")==0)                     thread_struct.worker_conf = WORKER_DISABLE;
  else if(strcmp(worker_conf,"WORKER_ENABLE")==0)                 thread_struct.worker_conf = WORKER_ENABLE;
  printf("[CONFIG] worker conf is set to %d\n",thread_struct.worker_conf);
114 115 116 117 118 119 120 121 122
} 
PARALLEL_CONF_t get_thread_parallel_conf(void)
{
  return thread_struct.parallel_conf;
} 
WORKER_CONF_t get_thread_worker_conf(void)
{
  return thread_struct.worker_conf;
} 
123

124
void
Cedric Roux's avatar
Cedric Roux committed
125
fill_nfapi_ulsch_config_request(nfapi_ul_config_request_pdu_t *ul_config_pdu,
Thomas Laurent's avatar
Thomas Laurent committed
126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147
                                uint8_t                        cqi_req,
                                uint8_t                        p_eNB,
                                uint8_t                        cqi_ReportModeAperiodic,
                                uint8_t                        betaOffset_CQI_Index,
                                uint8_t                        betaOffset_RI_Index,
                                uint8_t                        dl_cqi_pmi_size,
                                uint8_t                        tmode,
                                uint32_t                       handle,
                                uint16_t                       rnti,
                                uint8_t                        resource_block_start,
                                uint8_t                        number_of_resource_blocks,
                                uint8_t                        modulation_type,
                                uint8_t                        cyclic_shift_2_for_drms,
                                uint8_t                        frequency_hopping_enabled_flag,
                                uint8_t                        frequency_hopping_bits,
                                uint8_t                        new_data_indication,
                                uint8_t                        redundancy_version,
                                uint8_t                        harq_process_number,
                                uint8_t                        ul_tx_mode,
                                uint8_t                        current_tx_nb,
                                uint8_t                        n_srs,
                                uint16_t                       size) {
148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167
  memset((void *) ul_config_pdu, 0, sizeof(nfapi_ul_config_request_pdu_t));
  ul_config_pdu->pdu_type                                                    = NFAPI_UL_CONFIG_ULSCH_PDU_TYPE;
  ul_config_pdu->pdu_size                                                    = (uint8_t) (2 + sizeof(nfapi_ul_config_ulsch_pdu));
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.tl.tag                             = NFAPI_UL_CONFIG_REQUEST_ULSCH_PDU_REL8_TAG;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.handle                             = handle;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.rnti                               = rnti;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.resource_block_start               = resource_block_start;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.number_of_resource_blocks          = number_of_resource_blocks;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.modulation_type                    = modulation_type;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.cyclic_shift_2_for_drms            = cyclic_shift_2_for_drms;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.frequency_hopping_enabled_flag     = frequency_hopping_enabled_flag;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.frequency_hopping_bits             = frequency_hopping_bits;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.new_data_indication                = new_data_indication;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.redundancy_version                 = redundancy_version;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.harq_process_number                = harq_process_number;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.ul_tx_mode                         = ul_tx_mode;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.current_tx_nb                      = current_tx_nb;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.n_srs                              = n_srs;
  ul_config_pdu->ulsch_pdu.ulsch_pdu_rel8.size                               = size;

knopp's avatar
knopp committed
168 169
  //printf("Filling ul_config_pdu : Q %d, TBS %d, rv %d, ndi %d\n", modulation_type,size,redundancy_version,new_data_indication);

170 171 172 173 174 175 176 177 178 179
  if (cqi_req == 1) {
    // Add CQI portion
    ul_config_pdu->pdu_type = NFAPI_UL_CONFIG_ULSCH_CQI_RI_PDU_TYPE;
    ul_config_pdu->pdu_size = (uint8_t) (2 + sizeof(nfapi_ul_config_ulsch_cqi_ri_pdu));
    ul_config_pdu->ulsch_cqi_ri_pdu.cqi_ri_information.cqi_ri_information_rel9.tl.tag = NFAPI_UL_CONFIG_REQUEST_CQI_RI_INFORMATION_REL9_TAG;
    ul_config_pdu->ulsch_cqi_ri_pdu.cqi_ri_information.cqi_ri_information_rel9.report_type = 1;
    ul_config_pdu->ulsch_cqi_ri_pdu.cqi_ri_information.cqi_ri_information_rel9.aperiodic_cqi_pmi_ri_report.number_of_cc = 1;
    LOG_D(MAC, "report_type %d\n",ul_config_pdu->ulsch_cqi_ri_pdu.cqi_ri_information.cqi_ri_information_rel9.report_type);

    if (p_eNB <= 2
Thomas Laurent's avatar
Thomas Laurent committed
180
        && (tmode == 3 || tmode == 4 || tmode == 8 || tmode == 9 || tmode == 10))
181 182 183 184
      ul_config_pdu->ulsch_cqi_ri_pdu.cqi_ri_information.cqi_ri_information_rel9.aperiodic_cqi_pmi_ri_report.cc[0].ri_size = 1;
    else if (p_eNB <= 2) ul_config_pdu->ulsch_cqi_ri_pdu.cqi_ri_information.cqi_ri_information_rel9.aperiodic_cqi_pmi_ri_report.cc[0].ri_size = 0;
    else if (p_eNB == 4) ul_config_pdu->ulsch_cqi_ri_pdu.cqi_ri_information.cqi_ri_information_rel9.aperiodic_cqi_pmi_ri_report.cc[0].ri_size = 2;

Cedric Roux's avatar
Cedric Roux committed
185
    for (int ri = 0;
Thomas Laurent's avatar
Thomas Laurent committed
186 187 188
         ri <  (1 << ul_config_pdu->ulsch_cqi_ri_pdu.cqi_ri_information.cqi_ri_information_rel9.aperiodic_cqi_pmi_ri_report.cc[0].ri_size);
         ri++)
      ul_config_pdu->ulsch_cqi_ri_pdu.cqi_ri_information.cqi_ri_information_rel9.aperiodic_cqi_pmi_ri_report.cc[0].dl_cqi_pmi_size[ri] =  dl_cqi_pmi_size;
Cedric Roux's avatar
Cedric Roux committed
189

190 191 192 193
    ul_config_pdu->ulsch_cqi_ri_pdu.cqi_ri_information.cqi_ri_information_rel9.delta_offset_cqi = betaOffset_CQI_Index;
    ul_config_pdu->ulsch_cqi_ri_pdu.cqi_ri_information.cqi_ri_information_rel9.delta_offset_ri  = betaOffset_RI_Index;
  }
}
194

195
void fill_ulsch_dci(PHY_VARS_eNB *eNB,
Thomas Laurent's avatar
Thomas Laurent committed
196 197 198 199 200 201 202 203 204 205
                    int frame,
                    int subframe,
                    Sched_Rsp_t *sched_resp,
                    uint16_t rnti,
                    void *UL_dci,
                    int first_rb,
                    int nb_rb,
                    int mcs,
                    int modulation_type,
                    int ndi,
knopp's avatar
knopp committed
206
                    int TBS,
Thomas Laurent's avatar
Thomas Laurent committed
207 208 209 210
                    int cqi_flag,
                    uint8_t beta_CQI,
                    uint8_t beta_RI,
                    uint8_t cqi_size) {
211 212 213
  nfapi_ul_config_request_body_t *ul_req=&sched_resp->UL_req->ul_config_request_body;
  int harq_pid = ((frame*10)+subframe)&7;

knopp's avatar
knopp committed
214
  //printf("ulsch in frame %d, subframe %d => harq_pid %d, mcs %d, ndi %d\n",frame,subframe,harq_pid,mcs,ndi);
215 216

  switch (eNB->frame_parms.N_RB_UL) {
Thomas Laurent's avatar
Thomas Laurent committed
217 218
    case 6:
      break;
219

Thomas Laurent's avatar
Thomas Laurent committed
220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240
    case 25:
      if (eNB->frame_parms.frame_type == TDD) {
        ((DCI0_5MHz_TDD_1_6_t *)UL_dci)->type    = 0;
        ((DCI0_5MHz_TDD_1_6_t *)UL_dci)->rballoc = computeRIV(eNB->frame_parms.N_RB_UL,first_rb,nb_rb); // 12 RBs from position 8
        //printf("nb_rb %d/%d, rballoc %d (dci %x)\n",nb_rb,eNB->frame_parms.N_RB_UL,((DCI0_5MHz_TDD_1_6_t*)UL_dci)->rballoc,*(uint32_t *)UL_dci);
        ((DCI0_5MHz_TDD_1_6_t *)UL_dci)->mcs     = mcs;
        ((DCI0_5MHz_TDD_1_6_t *)UL_dci)->ndi     = ndi;
        ((DCI0_5MHz_TDD_1_6_t *)UL_dci)->TPC     = 0;
        ((DCI0_5MHz_TDD_1_6_t *)UL_dci)->cqi_req = cqi_flag&1;
        ((DCI0_5MHz_TDD_1_6_t *)UL_dci)->cshift  = 0;
        ((DCI0_5MHz_TDD_1_6_t *)UL_dci)->dai     = 1;
      } else {
        ((DCI0_5MHz_FDD_t *)UL_dci)->type    = 0;
        ((DCI0_5MHz_FDD_t *)UL_dci)->rballoc = computeRIV(eNB->frame_parms.N_RB_UL,first_rb,nb_rb); // 12 RBs from position 8
        //      printf("nb_rb %d/%d, rballoc %d (dci %x) (dcip %p)\n",nb_rb,eNB->frame_parms.N_RB_UL,((DCI0_5MHz_FDD_t*)UL_dci)->rballoc,*(uint32_t *)UL_dci,UL_dci);
        ((DCI0_5MHz_FDD_t *)UL_dci)->mcs     = mcs;
        ((DCI0_5MHz_FDD_t *)UL_dci)->ndi     = ndi;
        ((DCI0_5MHz_FDD_t *)UL_dci)->TPC     = 0;
        ((DCI0_5MHz_FDD_t *)UL_dci)->cqi_req = cqi_flag&1;
        ((DCI0_5MHz_FDD_t *)UL_dci)->cshift  = 0;
      }
241

Thomas Laurent's avatar
Thomas Laurent committed
242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288
      break;

    case 50:
      if (eNB->frame_parms.frame_type == TDD) {
        ((DCI0_10MHz_TDD_1_6_t *)UL_dci)->type    = 0;
        ((DCI0_10MHz_TDD_1_6_t *)UL_dci)->rballoc = computeRIV(eNB->frame_parms.N_RB_UL,first_rb,nb_rb); // 12 RBs from position 8
        //      printf("nb_rb %d/%d, rballoc %d (dci %x)\n",nb_rb,eNB->frame_parms.N_RB_UL,((DCI0_10MHz_TDD_1_6_t*)UL_dci)->rballoc,*(uint32_t *)UL_dci);
        ((DCI0_10MHz_TDD_1_6_t *)UL_dci)->mcs     = mcs;
        ((DCI0_10MHz_TDD_1_6_t *)UL_dci)->ndi     = ndi;
        ((DCI0_10MHz_TDD_1_6_t *)UL_dci)->TPC     = 0;
        ((DCI0_10MHz_TDD_1_6_t *)UL_dci)->cqi_req = cqi_flag&1;
        ((DCI0_10MHz_TDD_1_6_t *)UL_dci)->cshift  = 0;
        ((DCI0_10MHz_TDD_1_6_t *)UL_dci)->dai     = 1;
      } else {
        ((DCI0_10MHz_FDD_t *)UL_dci)->type    = 0;
        ((DCI0_10MHz_FDD_t *)UL_dci)->rballoc = computeRIV(eNB->frame_parms.N_RB_UL,first_rb,nb_rb); // 12 RBs from position 8
        //printf("nb_rb %d/%d, rballoc %d (dci %x)\n",nb_rb,eNB->frame_parms.N_RB_UL,((DCI0_10MHz_FDD_t*)UL_dci)->rballoc,*(uint32_t *)UL_dci);
        ((DCI0_10MHz_FDD_t *)UL_dci)->mcs     = mcs;
        ((DCI0_10MHz_FDD_t *)UL_dci)->ndi     = ndi;
        ((DCI0_10MHz_FDD_t *)UL_dci)->TPC     = 0;
        ((DCI0_10MHz_FDD_t *)UL_dci)->cqi_req = cqi_flag&1;
        ((DCI0_10MHz_FDD_t *)UL_dci)->cshift  = 0;
      }

      break;

    case 100:
      if (eNB->frame_parms.frame_type == TDD) {
        ((DCI0_20MHz_TDD_1_6_t *)UL_dci)->type    = 0;
        ((DCI0_20MHz_TDD_1_6_t *)UL_dci)->rballoc = computeRIV(eNB->frame_parms.N_RB_UL,first_rb,nb_rb); // 12 RBs from position 8
        //      printf("nb_rb %d/%d, rballoc %d (dci %x)\n",nb_rb,eNB->frame_parms.N_RB_UL,((DCI0_20MHz_TDD_1_6_t*)UL_dci)->rballoc,*(uint32_t *)UL_dci);
        ((DCI0_20MHz_TDD_1_6_t *)UL_dci)->mcs     = mcs;
        ((DCI0_20MHz_TDD_1_6_t *)UL_dci)->ndi     = ndi;
        ((DCI0_20MHz_TDD_1_6_t *)UL_dci)->TPC     = 0;
        ((DCI0_20MHz_TDD_1_6_t *)UL_dci)->cqi_req = cqi_flag&1;
        ((DCI0_20MHz_TDD_1_6_t *)UL_dci)->cshift  = 0;
        ((DCI0_20MHz_TDD_1_6_t *)UL_dci)->dai     = 1;
      } else {
        ((DCI0_20MHz_FDD_t *)UL_dci)->type    = 0;
        ((DCI0_20MHz_FDD_t *)UL_dci)->rballoc = computeRIV(eNB->frame_parms.N_RB_UL,first_rb,nb_rb); // 12 RBs from position 8
        //   printf("nb_rb %d/%d, rballoc %d (dci %x) (UL_dci %p)\n",nb_rb,eNB->frame_parms.N_RB_UL,((DCI0_20MHz_FDD_t*)UL_dci)->rballoc,*(uint32_t *)UL_dci,(void*)UL_dci);
        ((DCI0_20MHz_FDD_t *)UL_dci)->mcs     = mcs;
        ((DCI0_20MHz_FDD_t *)UL_dci)->ndi     = ndi;
        ((DCI0_20MHz_FDD_t *)UL_dci)->TPC     = 0;
        ((DCI0_20MHz_FDD_t *)UL_dci)->cqi_req = cqi_flag&1;
        ((DCI0_20MHz_FDD_t *)UL_dci)->cshift  = 0;
      }
289

Thomas Laurent's avatar
Thomas Laurent committed
290
      break;
291

Thomas Laurent's avatar
Thomas Laurent committed
292 293
    default:
      break;
294 295
  }

Cedric Roux's avatar
Cedric Roux committed
296
  fill_nfapi_ulsch_config_request(&ul_req->ul_config_pdu_list[0],
Thomas Laurent's avatar
Thomas Laurent committed
297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320
                                  cqi_flag&1,
                                  1,  // p_eNB
                                  0,  // reportmode Aperiodic
                                  beta_CQI,
                                  beta_RI,
                                  cqi_size,
                                  //cc,
                                  //UE_template->physicalConfigDedicated,
                                  1,
                                  0,
                                  14,     // rnti
                                  first_rb, // resource_block_start
                                  nb_rb,  // number_of_resource_blocks
                                  modulation_type,
                                  0,  // cyclic_shift_2_for_drms
                                  0,  // frequency_hopping_enabled_flag
                                  0,  // frequency_hopping_bits
                                  ndi,  // new_data_indication
                                  mcs>28?(mcs-28):0,  // redundancy_version
                                  harq_pid, // harq_process_number
                                  0,  // ul_tx_mode
                                  0,  // current_tx_nb
                                  0,  // n_srs
                                  TBS);
321 322 323
  sched_resp->UL_req->header.message_id = NFAPI_UL_CONFIG_REQUEST;
  ul_req->number_of_pdus=1;
  ul_req->tl.tag = NFAPI_UL_CONFIG_REQUEST_BODY_TAG;
Thomas Laurent's avatar
Thomas Laurent committed
324
}
Cedric Roux's avatar
Cedric Roux committed
325

Thomas Laurent's avatar
Thomas Laurent committed
326
enum eTypes { eBool, eInt, eFloat, eText };
327
static int verbose,help,disable_bundling=0,cqi_flag=0, extended_prefix_flag=0, test_perf=0, subframe=3, transmission_m=1,n_rx=1;
328

Thomas Laurent's avatar
Thomas Laurent committed
329
int main(int argc, char **argv) {
330
  int i,j,aa,u;
331 332 333
  PHY_VARS_eNB *eNB;
  PHY_VARS_UE *UE;
  RU_t *ru;
334 335
  int aarx,aatx;
  double channelx,channely;
Thomas Laurent's avatar
Thomas Laurent committed
336 337
  static double sigma2, sigma2_dB=10,SNR,SNR2=0,snr0=-2.0,snr1,SNRmeas,rate,saving_bler=0;
  static double input_snr_step=.2,snr_int=30;
338
  double blerr;
Thomas Laurent's avatar
Thomas Laurent committed
339
  int rvidx[8]= {0,2,3,1,0,2,3,1};
340 341
  int **txdata;
  LTE_DL_FRAME_PARMS *frame_parms;
342 343
  double s_re0[30720],s_im0[30720],r_re0[30720],r_im0[30720];
  double s_re1[30720],s_im1[30720],r_re1[30720],r_im1[30720];
344 345
  double r_re2[30720],r_im2[30720];
  double r_re3[30720],r_im3[30720];
Thomas Laurent's avatar
Thomas Laurent committed
346 347 348 349
  double *s_re[2]= {s_re0,s_re1};
  double *s_im[2]= {s_im0,s_im1};
  double *r_re[4]= {r_re0,r_re1,r_re2,r_re3};
  double *r_im[4]= {r_im0,r_im1,r_im2,r_im3};
350 351
  double forgetting_factor=0.0; //in [0,1] 0 means a new channel every time, 1 means keep the same channel
  double iqim=0.0;
Thomas Laurent's avatar
Thomas Laurent committed
352
  int cqi_error,cqi_errors,ack_errors,cqi_crc_falsepositives,cqi_crc_falsenegatives;
353 354 355 356
  int ch_realization;
  int eNB_id = 0;
  int chMod = 0 ;
  int UE_id = 0;
Thomas Laurent's avatar
Thomas Laurent committed
357
  static int nb_rb=25,first_rb=0,mcs=0,round=0;
358
  //unsigned char l;
Thomas Laurent's avatar
Thomas Laurent committed
359
  static int awgn_flag = 0 ;
360
  SCM_t channel_model=Rice1;
Thomas Laurent's avatar
Thomas Laurent committed
361
  unsigned char *input_buffer=0,harq_pid;
362 363 364
  unsigned short input_buffer_length;
  unsigned int ret;
  unsigned int coded_bits_per_codeword,nsymb;
365
  unsigned int tx_lev=0,tx_lev_dB,trials,errs[4]= {0,0,0,0},round_trials[4]= {0,0,0,0};
366
  FILE *bler_fd=NULL;
367
  char bler_fname[512];
368 369
  FILE *time_meas_fd=NULL;
  char time_meas_fname[256];
370 371 372 373
  FILE *input_fdUL=NULL,*trch_out_fdUL=NULL;
  //  unsigned char input_file=0;
  char input_val_str[50],input_val_str2[50];
  //  FILE *rx_frame_file;
374
  FILE *csv_fdUL=NULL;
375
  /*
Thomas Laurent's avatar
Thomas Laurent committed
376 377
    FILE *fperen=NULL;
    char fperen_name[512];
378

Thomas Laurent's avatar
Thomas Laurent committed
379 380
    FILE *fmageren=NULL;
    char fmageren_name[512];
381

Thomas Laurent's avatar
Thomas Laurent committed
382 383
    FILE *flogeren=NULL;
    char flogeren_name[512];
384
  */
385 386 387 388
  /* FILE *ftxlev;
     char ftxlev_name[512];
  */
  char csv_fname[512];
Thomas Laurent's avatar
Thomas Laurent committed
389 390 391
  static int n_frames=5000;
  static int n_ch_rlz = 1;
  static int abstx = 0;
392
  int hold_channel=0;
393
  channel_desc_t *UE2eNB;
394
  //uint8_t control_only_flag = 0;
Thomas Laurent's avatar
Thomas Laurent committed
395 396 397 398
  static int delay = 0;
  static double maxDoppler = 0.0;
  static int srs_flag = 0;
  static int N_RB_DL=25,osf=1;
399
  //uint8_t cyclic_shift = 0;
Thomas Laurent's avatar
Thomas Laurent committed
400 401 402 403
  static uint8_t beta_ACK=0,beta_RI=0,beta_CQI=2,cqi_size=11;
  static uint8_t tdd_config=3,frame_type=FDD;
  static int N0=30;
  static double tx_gain=1.0;
404
  double cpu_freq_GHz;
Thomas Laurent's avatar
Thomas Laurent committed
405
  int iter_trials;
406 407 408
  uint32_t UL_alloc_pdu;
  int s,Kr,Kr_bytes;
  int dump_perf=0;
Thomas Laurent's avatar
Thomas Laurent committed
409
  static int dump_table =0;
410
  double effective_rate=0.0;
Thomas Laurent's avatar
Thomas Laurent committed
411 412 413
  char channel_model_input[10]= {0};
  static int max_turbo_iterations=4;
  static int parallel_flag=0;
414
  int nb_rb_set = 0;
415
  int sf;
Thomas Laurent's avatar
Thomas Laurent committed
416
  static int threequarter_fs=0;
417
  int ndi;
knopp's avatar
knopp committed
418
  opp_enabled=1; // to enable the time meas
419 420 421 422
  sched_resp.DL_req = &DL_req;
  sched_resp.UL_req = &UL_req;
  sched_resp.HI_DCI0_req = &HI_DCI0_req;
  sched_resp.TX_req = &TX_req;
Thomas Laurent's avatar
Thomas Laurent committed
423 424 425 426
  memset((void *)&DL_req,0,sizeof(DL_req));
  memset((void *)&UL_req,0,sizeof(UL_req));
  memset((void *)&HI_DCI0_req,0,sizeof(HI_DCI0_req));
  memset((void *)&TX_req,0,sizeof(TX_req));
427 428
  UL_req.ul_config_request_body.ul_config_pdu_list = ul_config_pdu_list;
  TX_req.tx_request_body.tx_pdu_list = tx_pdu_list;
429
  cpu_freq_GHz = (double)get_cpu_freq_GHz();
430
  cpuf = cpu_freq_GHz;
431
  printf("Detected cpu_freq %f GHz\n",cpu_freq_GHz);
432
  AssertFatal(load_configmodule(argc,argv) != NULL,
Thomas Laurent's avatar
Thomas Laurent committed
433
              "cannot load configuration module, exiting\n");
434
  logInit();
Calisson's avatar
Calisson committed
435 436 437
  // enable these lines if you need debug info
  // however itti will catch all signals, so ctrl-c won't work anymore
  // alternatively you can disable ITTI completely in CMakeLists.txt
knopp's avatar
knopp committed
438 439 440
  //itti_init(TASK_MAX, THREAD_MAX, MESSAGES_ID_MAX, tasks_info, messages_info, messages_definition_xml, NULL);
  //set_comp_log(PHY,LOG_DEBUG,LOG_MED,1);
  //set_glog(LOG_DEBUG,LOG_MED);
Thomas Laurent's avatar
Thomas Laurent committed
441 442
  //hapZEbm:n:Y:X:x:s:w:e:q:d:D:O:c:r:i:f:y:c:oA:C:R:g:N:l:S:T:QB:PI:LF
  static paramdef_t options[] = {
443
    { "awgn", "Use AWGN channel and not multipath", PARAMFLAG_BOOL, strptr:NULL, defintval:0, TYPE_INT, 0, NULL, NULL },
Thomas Laurent's avatar
Thomas Laurent committed
444 445 446 447
    { "BnbRBs", "The LTE bandwith in RBs (100 is 20MHz)",0, iptr:&N_RB_DL,  defintval:25, TYPE_INT, 0 },
    { "mcs", "The MCS to use", 0, iptr:&mcs,  defintval:10, TYPE_INT, 0 },
    { "nb_frame", "number of frame in a test",0, iptr:&n_frames,  defintval:1, TYPE_INT, 0 },
    { "snr", "starting snr", 0, dblptr:&snr0,  defdblval:-2.9, TYPE_DOUBLE, 0 },
448
    { "wsnrInterrupt", "snr int ?", 0, dblptr:&snr_int,  defdblval:30, TYPE_DOUBLE, 0 },
449
    { "e_snr_step", "step increasing snr",0, dblptr:&input_snr_step,  defdblval:0.2, TYPE_DOUBLE, 0 },
Thomas Laurent's avatar
Thomas Laurent committed
450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475
    { "rb_dynamic", "number of rb in dynamic allocation",0, iptr:NULL,  defintval:0, TYPE_INT, 0 },
    { "first_rb", "first rb used in dynamic allocation",0, iptr:&first_rb,  defintval:0, TYPE_INT, 0 },
    { "osrs", "enable srs generation",PARAMFLAG_BOOL, iptr:&srs_flag,  defintval:0, TYPE_INT, 0 },
    { "gchannel", "[A:M] Use 3GPP 25.814 SCM-A/B/C/D('A','B','C','D') or 36-101 EPA('E'), EVA ('F'),ETU('G') models (ignores delay spread and Ricean factor), Rayghleigh8 ('H'), Rayleigh1('I'), Rayleigh1_corr('J'), Rayleigh1_anticorr ('K'),  Rice8('L'), Rice1('M')",0, strptr:NULL,  defstrval:NULL, TYPE_STRING, 0 },
    { "delay_chan", "Channel delay",0, iptr:&delay,  defintval:0, TYPE_INT, 0 },
    { "Doppler", "Maximum doppler shift",0, dblptr:&maxDoppler,  defdblval:0.0, TYPE_DOUBLE, 0 },
    { "Zdump", "dump table",PARAMFLAG_BOOL,  iptr:&dump_table, defintval:0, TYPE_INT, 0 },
    { "Forms", "Display the soft scope", PARAMFLAG_BOOL, iptr:&xforms,  defintval:0, TYPE_INT, 0 },
    { "Lparallel", "Enable parallel execution", PARAMFLAG_BOOL, iptr:&parallel_flag,  defintval:0, TYPE_INT, 0 },
    { "Iterations", "Number of iterations of turbo decoder", 0, iptr:&max_turbo_iterations,  defintval:4, TYPE_INT, 0 },
    { "Performance", "Display CPU perfomance of each L1 piece", PARAMFLAG_BOOL,  iptr:NULL,  defintval:0, TYPE_INT, 0 },
    { "Q_cqi", "Enable CQI", PARAMFLAG_BOOL, iptr:&cqi_flag,  defintval:0, TYPE_INT, 0 },
    { "prefix_extended","Extended prefix", PARAMFLAG_BOOL, iptr:&extended_prefix_flag,  defintval:0, TYPE_INT, 0 },
    { "RI_beta", "TBD", 0, iptr:NULL,  defintval:0, TYPE_INT, 0 },
    { "CQI_beta", "TBD",0, iptr:NULL,  defintval:0, TYPE_INT, 0 },
    { "ACK_beta", "TBD",0, iptr:NULL,  defintval:0, TYPE_INT, 0 },
    { "input_file", "input IQ data file",0, iptr:NULL,  defintval:0, TYPE_INT, 0 },
    { "N0", "N0",0, iptr:&N0,  defintval:30, TYPE_INT, 0 },
    { "EsubSampling","three quarters sub-sampling",PARAMFLAG_BOOL, iptr:&threequarter_fs, defintval:0, TYPE_INT, 0 },
    { "TDD", "Enable TDD and set the tdd configuration mode",0, iptr:NULL,  defintval:25, TYPE_INT, 0 },
    { "Subframe", "subframe to use",0, iptr:&subframe,  defintval:3, TYPE_INT, 0 },
    { "xTransmission","transmission mode (1 or 2 are supported)",0, iptr:NULL,  defintval:25, TYPE_INT, 0 },
    { "yN_rx", "TBD: n_rx",0, iptr:&n_rx,  defintval:1, TYPE_INT, 0 },
    { "bundling_disable", "bundling disable",PARAMFLAG_BOOL,  iptr:&disable_bundling, defintval:0, TYPE_INT, 0 },
    { "Y",  "n_ch_rlz",0, iptr:&n_ch_rlz,  defintval:1, TYPE_INT, 0 },
    { "X",  "abstx", PARAMFLAG_BOOL,  iptr:&abstx, defintval:0, TYPE_INT, 0 },
476
    { "Operf", "Set the percentage of effective rate to testbench the modem performance (typically 30 and 70, range 1-100)",0, iptr:&test_perf,  defintval:0, TYPE_INT, 0 },
Thomas Laurent's avatar
Thomas Laurent committed
477
    { "verbose", "display debug text", PARAMFLAG_BOOL,  iptr:&verbose, defintval:0, TYPE_INT, 0 },
478
    { "help", "display help and exit", PARAMFLAG_BOOL,  iptr:&help, defintval:0, TYPE_INT, 0 },
Thomas Laurent's avatar
Thomas Laurent committed
479 480
    { "", "",0,  iptr:NULL, defintval:0, TYPE_INT, 0 },
  };
481
  struct option * long_options = parse_oai_options(options); 
482

Thomas Laurent's avatar
Thomas Laurent committed
483
  int option_index;
484

Thomas Laurent's avatar
Thomas Laurent committed
485
  int res;
486

Thomas Laurent's avatar
Thomas Laurent committed
487 488 489 490 491
  while ((res=getopt_long_only(argc, argv, "", long_options, &option_index)) == 0) {
    if (options[option_index].voidptr != NULL ) {
      if (long_options[option_index].has_arg==no_argument)
        *(bool *)options[option_index].iptr=1;
      else switch (options[option_index].type) {
492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510
	case TYPE_INT:
	  *(int *)options[option_index].iptr=atoi(optarg);
	  break;

	case TYPE_DOUBLE:
	  *(double *)options[option_index].dblptr=atof(optarg);
	  break;

	case TYPE_UINT8:
	  *(uint8_t *)options[option_index].dblptr=atoi(optarg);
	  break;
	  	    
	case TYPE_UINT16:
	  *(uint16_t *)options[option_index].dblptr=atoi(optarg);
	  break;
	  
	default:
	  printf("not decoded type.\n");
	  exit(1);
Thomas Laurent's avatar
Thomas Laurent committed
511
        }
512

Thomas Laurent's avatar
Thomas Laurent committed
513 514
      continue;
    }
515

Thomas Laurent's avatar
Thomas Laurent committed
516
    switch (long_options[option_index].name[0]) {
517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571
    case 'T':
      tdd_config=atoi(optarg);
      frame_type=TDD;
      break;

    case 'a':
      channel_model = AWGN;
      chMod = 1;
      break;

    case 'g':
      strncpy(channel_model_input,optarg,9);
      struct tmp {
	char opt;
	int m;
	int M;
      }
      tmp[]= {
	{'A',SCM_A,2},
	{'B',SCM_B,3},
	{'C',SCM_C,4},
	{'D',SCM_D,5},
	{'E',EPA,6},
	{'G',ETU,8},
	{'H',Rayleigh8,9},
	{'I',Rayleigh1,10},
	{'J',Rayleigh1_corr,11},
	{'K',Rayleigh1_anticorr,12},
	{'L',Rice8,13},
	{'M',Rice1,14},
	{'N',AWGN,1},
	{0,0,0}
      };
      struct tmp *ptr;

      for (ptr=tmp; ptr->opt!=0; ptr++)
	if ( ptr->opt == optarg[0] ) {
	  channel_model=ptr->m;
	  chMod=ptr->M;
	  break;
	}

      AssertFatal(ptr->opt != 0, "Unsupported channel model: %s !\n", optarg );
      break;

    case 'x':
      transmission_m=atoi(optarg);
      AssertFatal(transmission_m==1 || transmission_m==2,
		  "Unsupported transmission mode %d\n",transmission_m);
      break;

    case 'r':
      nb_rb = atoi(optarg);
      nb_rb_set = 1;
      break;
572

Thomas Laurent's avatar
Thomas Laurent committed
573 574 575
      //case 'c':
      //  cyclic_shift = atoi(optarg);
      //  break;
576

577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606
    case 'i':
      input_fdUL = fopen(optarg,"r");
      printf("Reading in %s (%p)\n",optarg,input_fdUL);
      AssertFatal(input_fdUL != (FILE *)NULL,"Unknown file %s\n",optarg);
      break;

    case 'A':
      beta_ACK = atoi(optarg);
      AssertFatal(beta_ACK>15,"beta_ack must be in (0..15)\n");
      break;

    case 'C':
      beta_CQI = atoi(optarg);
      AssertFatal((beta_CQI>15)||(beta_CQI<2),"beta_cqi must be in (2..15)\n");
      break;

    case 'R':
      beta_RI = atoi(optarg);
      AssertFatal((beta_RI>15)||(beta_RI<2),"beta_ri must be in (0..13)\n");
      break;

    case 'P':
      dump_perf=1;
      opp_enabled=1;
      break;

    default:
      printf("Wrong option: %s\n",long_options[option_index].name);
      exit(1);
      break;
Thomas Laurent's avatar
Thomas Laurent committed
607 608
    }
  }
609

Thomas Laurent's avatar
Thomas Laurent committed
610 611 612 613
  if ( res != -1 ) {
    printf("A wrong option has been found\n");
    exit(1);
  }
614 615


616 617 618 619 620
  if (help || verbose )
     display_options_values(options, true);
  if (help)
    exit(0);
  
621 622
  set_parallel_conf("PARALLEL_RU_L1_TRX_SPLIT");
  set_worker_conf("WORKER_ENABLE");
623 624 625
  RC.nb_L1_inst = 1;
  RC.nb_RU = 1;
  lte_param_init(&eNB,&UE,&ru,
626
                 1,
Thomas Laurent's avatar
Thomas Laurent committed
627 628 629 630 631 632 633 634 635 636 637 638 639
                 1,
                 n_rx,
                 1,
                 1,
                 extended_prefix_flag,
                 frame_type,
                 0,
                 tdd_config,
                 N_RB_DL,
                 4,
                 threequarter_fs,
                 osf,
                 0);
640 641 642 643 644 645
  RC.eNB = (PHY_VARS_eNB ***)malloc(sizeof(PHY_VARS_eNB **));
  RC.eNB[0] = (PHY_VARS_eNB **)malloc(sizeof(PHY_VARS_eNB *));
  RC.ru = (RU_t **)malloc(sizeof(RC.ru));
  RC.eNB[0][0] = eNB;
  RC.ru[0] = ru;

Thomas Laurent's avatar
Thomas Laurent committed
646
  for (int k=0; k<eNB->RU_list[0]->nb_rx; k++) eNB->common_vars.rxdataF[k]     =  eNB->RU_list[0]->common.rxdataF[k];
Cedric Roux's avatar
Cedric Roux committed
647

Thomas Laurent's avatar
Thomas Laurent committed
648
  memset((void *)&eNB->UL_INFO,0,sizeof(eNB->UL_INFO));
649 650 651 652 653 654 655 656
  printf("Setting indication lists\n");
  eNB->UL_INFO.rx_ind.rx_indication_body.rx_pdu_list   = eNB->rx_pdu_list;
  eNB->UL_INFO.crc_ind.crc_indication_body.crc_pdu_list = eNB->crc_pdu_list;
  eNB->UL_INFO.sr_ind.sr_indication_body.sr_pdu_list = eNB->sr_pdu_list;
  eNB->UL_INFO.harq_ind.harq_indication_body.harq_pdu_list = eNB->harq_pdu_list;
  eNB->UL_INFO.cqi_ind.cqi_pdu_list = eNB->cqi_pdu_list;
  eNB->UL_INFO.cqi_ind.cqi_raw_pdu_list = eNB->cqi_raw_pdu_list;
  printf("lte_param_init done\n");
Calisson's avatar
Calisson committed
657
  // for a call to phy_reset_ue later we need PHY_vars_UE_g allocated and pointing to UE
Thomas Laurent's avatar
Thomas Laurent committed
658 659
  PHY_vars_UE_g = (PHY_VARS_UE ***)malloc(sizeof(PHY_VARS_UE **));
  PHY_vars_UE_g[0] = (PHY_VARS_UE **) malloc(sizeof(PHY_VARS_UE *));
Calisson's avatar
Calisson committed
660 661
  PHY_vars_UE_g[0][0] = UE;

662
  if (nb_rb_set == 0)
663
    nb_rb = eNB->frame_parms.N_RB_UL;
664

665
  printf("1 . rxdataF_comp[0] %p\n",eNB->pusch_vars[0]->rxdataF_comp[0]);
666
  printf("Setting mcs = %d\n",mcs);
667
  printf("n_frames = %d\n", n_frames);
knopp's avatar
knopp committed
668
  snr1 = snr0+snr_int;
669
  printf("SNR0 %f, SNR1 %f\n",snr0,snr1);
670 671 672
  frame_parms = &eNB->frame_parms;
  txdata = UE->common_vars.txdata;
  nsymb = (eNB->frame_parms.Ncp == NORMAL) ? 14 : 12;
673
  sprintf(bler_fname,"ULbler_mcs%d_nrb%d_ChannelModel%d_nsim%d.csv",mcs,nb_rb,chMod,n_frames);
674
  bler_fd = fopen(bler_fname,"w");
Thomas Laurent's avatar
Thomas Laurent committed
675

676 677 678 679
  if (bler_fd==NULL) {
    fprintf(stderr,"Problem creating file %s\n",bler_fname);
    exit(-1);
  }
680

681
  fprintf(bler_fd,"#SNR;mcs;nb_rb;TBS;rate;errors[0];trials[0];errors[1];trials[1];errors[2];trials[2];errors[3];trials[3]\n");
682 683

  if (test_perf != 0) {
684 685 686 687
    char hostname[1024];
    hostname[1023] = '\0';
    gethostname(hostname, 1023);
    printf("Hostname: %s\n", hostname);
688 689 690 691
    //char dirname[FILENAME_MAX];
    //sprintf(dirname, "%s//SIMU/USER/pre-ci-logs-%s", getenv("OPENAIR_TARGETS"),hostname);
    //mkdir(dirname, 0777);
    sprintf(time_meas_fname,"time_meas_prb%d_mcs%d_antrx%d_channel%s_tx%d.csv",
Thomas Laurent's avatar
Thomas Laurent committed
692
            N_RB_DL,mcs,n_rx,channel_model_input,transmission_m);
693
    time_meas_fd = fopen(time_meas_fname,"w");
Thomas Laurent's avatar
Thomas Laurent committed
694

695 696 697 698
    if (time_meas_fd==NULL) {
      fprintf(stderr,"Cannot create file %s!\n",time_meas_fname);
      exit(-1);
    }
699
  }
700 701 702

  if(abstx) {
    // CSV file
Thomas Laurent's avatar
Thomas Laurent committed
703
    sprintf(csv_fname,"EULdataout_tx%d_mcs%d_nbrb%d_chan%d_nsimus%d_eren.m",transmission_m,mcs,nb_rb,chMod,n_frames);
704
    csv_fdUL = fopen(csv_fname,"w");
Thomas Laurent's avatar
Thomas Laurent committed
705

706 707 708 709
    if (csv_fdUL == NULL) {
      fprintf(stderr,"Problem opening file %s\n",csv_fname);
      exit(-1);
    }
Thomas Laurent's avatar
Thomas Laurent committed
710

711 712
    fprintf(csv_fdUL,"data_all%d=[",mcs);
  }
713

714 715 716 717 718 719
  if (xforms==1) {
    fl_initialize (&argc, argv, NULL, 0, 0);
    form_enb = create_lte_phy_scope_enb();
    sprintf (title, "LTE PHY SCOPE eNB");
    fl_show_form (form_enb->lte_phy_scope_enb, FL_PLACE_HOTSPOT, FL_FULLBORDER, title);
  }
720

Bilel's avatar
Bilel committed
721
  UE->pdcch_vars[0][0]->crnti = 14;
722
  UE->frame_parms.soundingrs_ul_config_common.enabled_flag = srs_flag;
723
  UE->frame_parms.soundingrs_ul_config_common.srs_BandwidthConfig = 2;
724 725 726 727
  UE->frame_parms.soundingrs_ul_config_common.srs_SubframeConfig = 3;
  UE->soundingrs_ul_config_dedicated[eNB_id].srsConfigDedicatedSetup = srs_flag;
  UE->soundingrs_ul_config_dedicated[eNB_id].duration = 1;
  UE->soundingrs_ul_config_dedicated[eNB_id].srs_ConfigIndex = 2;
728 729 730
  UE->soundingrs_ul_config_dedicated[eNB_id].srs_Bandwidth = 0;
  UE->soundingrs_ul_config_dedicated[eNB_id].transmissionComb = 0;
  UE->soundingrs_ul_config_dedicated[eNB_id].freqDomainPosition = 0;
731
  UE->soundingrs_ul_config_dedicated[eNB_id].cyclicShift = 0;
732
  eNB->frame_parms.soundingrs_ul_config_common.enabled_flag = srs_flag;
733
  eNB->frame_parms.soundingrs_ul_config_common.srs_BandwidthConfig = 2;
734 735 736 737
  eNB->frame_parms.soundingrs_ul_config_common.srs_SubframeConfig = 3;
  eNB->soundingrs_ul_config_dedicated[UE_id].srsConfigDedicatedSetup = srs_flag;
  eNB->soundingrs_ul_config_dedicated[UE_id].duration = 1;
  eNB->soundingrs_ul_config_dedicated[UE_id].srs_ConfigIndex = 2;
738 739 740
  eNB->soundingrs_ul_config_dedicated[UE_id].srs_Bandwidth = 0;
  eNB->soundingrs_ul_config_dedicated[UE_id].transmissionComb = 0;
  eNB->soundingrs_ul_config_dedicated[UE_id].freqDomainPosition = 0;
741
  eNB->soundingrs_ul_config_dedicated[UE_id].cyclicShift = 0;
742 743 744 745 746 747 748
  eNB->pusch_config_dedicated[UE_id].betaOffset_ACK_Index = beta_ACK;
  eNB->pusch_config_dedicated[UE_id].betaOffset_RI_Index  = beta_RI;
  eNB->pusch_config_dedicated[UE_id].betaOffset_CQI_Index = beta_CQI;
  UE->pusch_config_dedicated[eNB_id].betaOffset_ACK_Index = beta_ACK;
  UE->pusch_config_dedicated[eNB_id].betaOffset_RI_Index  = beta_RI;
  UE->pusch_config_dedicated[eNB_id].betaOffset_CQI_Index = beta_CQI;
  UE->ul_power_control_dedicated[eNB_id].deltaMCS_Enabled = 1;
749 750 751
  // disable periodic cqi/ri reporting
  UE->cqi_report_config[eNB_id].CQI_ReportPeriodic.ri_ConfigIndex = -1;
  UE->cqi_report_config[eNB_id].CQI_ReportPeriodic.cqi_PMI_ConfigIndex = -1;
752
  printf("PUSCH Beta : ACK %f, RI %f, CQI %f\n",(double)beta_ack[beta_ACK]/8,(double)beta_ri[beta_RI]/8,(double)beta_cqi[beta_CQI]/8);
753 754
  UE2eNB = new_channel_desc_scm(1,
                                n_rx,
755
                                channel_model,
Thomas Laurent's avatar
Thomas Laurent committed
756 757
                                N_RB2sampling_rate(eNB->frame_parms.N_RB_UL),
                                N_RB2channel_bandwidth(eNB->frame_parms.N_RB_UL),
758 759 760 761 762 763 764
                                forgetting_factor,
                                delay,
                                0);
  // set Doppler
  UE2eNB->max_Doppler = maxDoppler;

  // NN: N_RB_UL has to be defined in ulsim
Thomas Laurent's avatar
Thomas Laurent committed
765 766
  for (int k=0; k<NUMBER_OF_UE_MAX; k++) eNB->ulsch[k] = new_eNB_ulsch(max_turbo_iterations,N_RB_DL,0);

767
  UE->ulsch[0]   = new_ue_ulsch(N_RB_DL,0);
768
  printf("ULSCH %p\n",UE->ulsch[0]);
769

770
  if (parallel_flag == 1) {
771
    extern void init_fep_thread(PHY_VARS_eNB *, pthread_attr_t *);
772
    extern void init_td_thread(PHY_VARS_eNB *);
773
    init_fep_thread(eNB,NULL);
774
    init_td_thread(eNB);
775
  }
Thomas Laurent's avatar
Thomas Laurent committed
776

777
  // Create transport channel structures for 2 transport blocks (MIMO)
778
  for (i=0; i<2; i++) {
779
    eNB->dlsch[0][i] = new_eNB_dlsch(1,8,1827072,N_RB_DL,0,&eNB->frame_parms);
Thomas Laurent's avatar
Thomas Laurent committed
780

781
    if (!eNB->dlsch[0][i]) {
782 783 784
      printf("Can't get eNB dlsch structures\n");
      exit(-1);
    }
Thomas Laurent's avatar
Thomas Laurent committed
785

786
    eNB->dlsch[0][i]->rnti = 14;
787
  }
Thomas Laurent's avatar
Thomas Laurent committed
788

789
  /* allocate memory for both subframes (only one is really used
Thomas Laurent's avatar
Thomas Laurent committed
790 791 792 793
     but there is now "copy_harq_proc_struct" which needs both
     to be valid)
     TODO: refine this somehow (necessary?)
  */
794 795 796
  for (sf = 0; sf < 2; sf++) {
    for (i=0; i<2; i++) {
      UE->dlsch[sf][0][i]  = new_ue_dlsch(1,8,1827072,MAX_TURBO_ITERATIONS,N_RB_DL,0);
Thomas Laurent's avatar
Thomas Laurent committed
797

798 799 800 801
      if (!UE->dlsch[sf][0][i]) {
        printf("Can't get ue dlsch structures\n");
        exit(-1);
      }
Thomas Laurent's avatar
Thomas Laurent committed
802

803 804 805
      UE->dlsch[sf][0][i]->rnti   = 14;
    }
  }
806

807 808
  UE->dlsch_SI[0]  = new_ue_dlsch(1,1,1827072,MAX_TURBO_ITERATIONS,N_RB_DL,0);
  UE->dlsch_ra[0]  = new_ue_dlsch(1,1,1827072,MAX_TURBO_ITERATIONS,N_RB_DL,0);
809 810
  UE->measurements.rank[0] = 0;
  UE->transmission_mode[0] = 2;
Thomas Laurent's avatar
Thomas Laurent committed
811
  UE->pucch_config_dedicated[0].tdd_AckNackFeedbackMode = disable_bundling == 0 ? bundling : multiplexing;
812
  eNB->transmission_mode[0] = 2;
Thomas Laurent's avatar
Thomas Laurent committed
813
  eNB->pucch_config_dedicated[0].tdd_AckNackFeedbackMode = disable_bundling == 0 ? bundling : multiplexing;
814 815 816 817 818 819
  UE->frame_parms.pusch_config_common.ul_ReferenceSignalsPUSCH.groupHoppingEnabled = 1;
  eNB->frame_parms.pusch_config_common.ul_ReferenceSignalsPUSCH.groupHoppingEnabled = 1;
  UE->frame_parms.pusch_config_common.ul_ReferenceSignalsPUSCH.sequenceHoppingEnabled = 0;
  eNB->frame_parms.pusch_config_common.ul_ReferenceSignalsPUSCH.sequenceHoppingEnabled = 0;
  UE->frame_parms.pusch_config_common.ul_ReferenceSignalsPUSCH.groupAssignmentPUSCH = 0;
  eNB->frame_parms.pusch_config_common.ul_ReferenceSignalsPUSCH.groupAssignmentPUSCH = 0;
820 821 822
  UE->mac_enabled=0;
  eNB_rxtx_proc_t *proc_rxtx   = &eNB->proc.proc_rxtx[subframe&1];
  UE_rxtx_proc_t *proc_rxtx_ue = &UE->proc.proc_rxtx[subframe&1];
823 824 825 826
  proc_rxtx->frame_rx=1;
  proc_rxtx->subframe_rx=subframe;
  proc_rxtx->frame_tx=pdcch_alloc2ul_frame(&eNB->frame_parms,1,subframe);
  proc_rxtx->subframe_tx=pdcch_alloc2ul_subframe(&eNB->frame_parms,subframe);
827
  proc_rxtx_ue->frame_tx = proc_rxtx->frame_rx;
828
  proc_rxtx_ue->frame_rx = (subframe<4)?(proc_rxtx->frame_tx-1):(proc_rxtx->frame_tx);
829
  proc_rxtx_ue->subframe_tx = proc_rxtx->subframe_rx;
830
  proc_rxtx_ue->subframe_rx = (proc_rxtx->subframe_tx+6)%10;
831
  printf("Init UL hopping UE\n");
832
  init_ul_hopping(&UE->frame_parms);
833
  printf("Init UL hopping eNB\n");
834
  init_ul_hopping(&eNB->frame_parms);
835
  UE->dlsch[subframe&1][0][0]->harq_ack[ul_subframe2pdcch_alloc_subframe(&eNB->frame_parms,subframe)].send_harq_status = 1;
836 837
  UE->ulsch_Msg3_active[eNB_id] = 0;
  UE->ul_power_control_dedicated[eNB_id].accumulationEnabled=1;
838 839
  coded_bits_per_codeword = nb_rb * (12 * get_Qm_ul(mcs)) * nsymb;

840
  if (cqi_flag == 1) coded_bits_per_codeword-=UE->ulsch[0]->O;
841 842 843

  rate = (double)dlsch_tbs25[get_I_TBS(mcs)][nb_rb-1]/(coded_bits_per_codeword);
  printf("Rate = %f (mod %d), coded bits %d\n",rate,get_Qm_ul(mcs),coded_bits_per_codeword);
844

845
  for (ch_realization=0; ch_realization<n_ch_rlz; ch_realization++) {
846 847 848 849 850 851
    /*
      if(abstx){
      int ulchestim_f[300*12];
      int ulchestim_t[2*(frame_parms->ofdm_symbol_size)];
      }
    */
852
    if(abstx) {
853 854 855
      printf("**********************Channel Realization Index = %d **************************\n", ch_realization);
      saving_bler=1;
    }
856

knopp's avatar
knopp committed
857
    //    if ((subframe>5) || (subframe < 4))
858
    //      UE->frame++;
859 860

    for (SNR=snr0; SNR<snr1; SNR+=input_snr_step) {
861 862 863 864 865 866 867 868 869 870 871 872 873 874
      errs[0]=0;
      errs[1]=0;
      errs[2]=0;
      errs[3]=0;
      round_trials[0] = 0;
      round_trials[1] = 0;
      round_trials[2] = 0;
      round_trials[3] = 0;
      cqi_errors=0;
      ack_errors=0;
      cqi_crc_falsepositives=0;
      cqi_crc_falsenegatives=0;
      round=0;
      //randominit(0);
875
      harq_pid = subframe2harq_pid(&UE->frame_parms,proc_rxtx_ue->frame_tx,subframe);
876
      input_buffer_length = UE->ulsch[0]->harq_processes[harq_pid]->TBS/8;
Thomas Laurent's avatar
Thomas Laurent committed
877 878 879 880

      if ( input_buffer != NULL )
        free(input_buffer);

881
      input_buffer = (unsigned char *)memalign(32,input_buffer_length+64);
Thomas Laurent's avatar
Thomas Laurent committed
882

883
      //      printf("UL frame %d/subframe %d, harq_pid %d\n",UE->frame,subframe,harq_pid);
884
      if (input_fdUL == NULL) {
885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904
        if (n_frames == 1) {
          trch_out_fdUL= fopen("ulsch_trchUL.txt","w");

          for (i=0; i<input_buffer_length; i++) {
            input_buffer[i] = taus()&0xff;

            for (j=0; j<8; j++)
              fprintf(trch_out_fdUL,"%d\n",(input_buffer[i]>>(7-j))&1);
          }

          fclose(trch_out_fdUL);
        } else {
          for (i=0; i<input_buffer_length; i++)
            input_buffer[i] = taus()&0xff;
        }
      } else {
        n_frames=1;
        i=0;

        while (!feof(input_fdUL)) {
905
          ret=fscanf(input_fdUL,"%s %s",input_val_str,input_val_str2);//&input_val1,&input_val2);
Thomas Laurent's avatar
Thomas Laurent committed
906

907
          if (ret != 2) printf("ERROR: error reading file\n");
908 909

          if ((i%4)==0) {
Thomas Laurent's avatar
Thomas Laurent committed
910 911
            ((short *)txdata[0])[i/2] = (short)((1<<15)*strtod(input_val_str,NULL));
            ((short *)txdata[0])[(i/2)+1] = (short)((1<<15)*strtod(input_val_str2,NULL));
912 913

            if ((i/4)<100)
Thomas Laurent's avatar
Thomas Laurent committed
914
              printf("sample %d => %e + j%e (%d +j%d)\n",i/4,strtod(input_val_str,NULL),strtod(input_val_str2,NULL),((short *)txdata[0])[i/4],((short *)txdata[0])[(i/4)+1]); //1,input_val2,);
915 916 917 918 919 920 921 922 923
          }

          i++;

          if (i>(FRAME_LENGTH_SAMPLES))
            break;
        }

        printf("Read in %d samples\n",i/4);
Mongazon's avatar
Mongazon committed
924 925
        //      LOG_M("txsig0UL.m","txs0", txdata[0],2*frame_parms->samples_per_tti,1,1);
        //    LOG_M("txsig1.m","txs1", txdata[1],FRAME_LENGTH_COMPLEX_SAMPLES,1,1);
926 927 928
        tx_lev = signal_energy(&txdata[0][0],
                               OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES);
        tx_lev_dB = (unsigned int) dB_fixed(tx_lev);
929 930
      }

931
      iter_trials=0;
932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949
      reset_meas(&UE->phy_proc_tx);
      reset_meas(&UE->ofdm_mod_stats);
      reset_meas(&UE->ulsch_modulation_stats);
      reset_meas(&UE->ulsch_encoding_stats);
      reset_meas(&UE->ulsch_interleaving_stats);
      reset_meas(&UE->ulsch_rate_matching_stats);
      reset_meas(&UE->ulsch_turbo_encoding_stats);
      reset_meas(&UE->ulsch_segmentation_stats);
      reset_meas(&UE->ulsch_multiplexing_stats);
      reset_meas(&eNB->phy_proc_rx);
      reset_meas(&eNB->ulsch_channel_estimation_stats);
      reset_meas(&eNB->ulsch_freq_offset_estimation_stats);
      reset_meas(&eNB->rx_dft_stats);
      reset_meas(&eNB->ulsch_decoding_stats);
      reset_meas(&eNB->ulsch_turbo_decoding_stats);
      reset_meas(&eNB->ulsch_deinterleaving_stats);
      reset_meas(&eNB->ulsch_demultiplexing_stats);
      reset_meas(&eNB->ulsch_rate_unmatching_stats);
Thomas Laurent's avatar
Thomas Laurent committed
950
      reset_meas(&eNB->ulsch_demodulation_stats);
951 952 953 954 955 956 957
      reset_meas(&eNB->ulsch_tc_init_stats);
      reset_meas(&eNB->ulsch_tc_alpha_stats);
      reset_meas(&eNB->ulsch_tc_beta_stats);
      reset_meas(&eNB->ulsch_tc_gamma_stats);
      reset_meas(&eNB->ulsch_tc_ext_stats);
      reset_meas(&eNB->ulsch_tc_intl1_stats);
      reset_meas(&eNB->ulsch_tc_intl2_stats);
Cedric Roux's avatar
Cedric Roux committed
958
      // initialization
Thomas Laurent's avatar
Thomas Laurent committed
959 960 961 962 963 964 965 966
      varArray_t *table_tx=initVarArray(1000,sizeof(double));
      varArray_t *table_tx_ifft=initVarArray(1000,sizeof(double));
      varArray_t *table_tx_mod=initVarArray(1000,sizeof(double));
      varArray_t *table_tx_enc=initVarArray(1000,sizeof(double));
      varArray_t *table_rx=initVarArray(1000,sizeof(double));
      varArray_t *table_rx_fft=initVarArray(1000,sizeof(double));
      varArray_t *table_rx_demod=initVarArray(1000,sizeof(double));
      varArray_t *table_rx_dec=initVarArray(1000,sizeof(double));
967
      ndi=0;
968
      phy_reset_ue(0,0,0);
969
      UE->UE_mode[eNB_id]=PUSCH;
Thomas Laurent's avatar
Thomas Laurent committed
970
      SET_LOG_DEBUG(UE_TIMING);
971

972 973
      for (trials = 0; trials<n_frames; trials++) {
        //      printf("*");
974 975
        //        UE->frame++;
        //        eNB->frame++;
Thomas Laurent's avatar
Thomas Laurent committed
976
        ndi = (1-ndi);
977 978 979 980
        fflush(stdout);
        round=0;

        while (round < 4) {
Thomas Laurent's avatar
Thomas Laurent committed
981 982 983 984 985 986 987 988
          proc_rxtx->frame_rx=1;
          proc_rxtx->subframe_rx=subframe;
          proc_rxtx->frame_tx=pdcch_alloc2ul_frame(&eNB->frame_parms,1,subframe);
          proc_rxtx->subframe_tx=pdcch_alloc2ul_subframe(&eNB->frame_parms,subframe);
          proc_rxtx_ue->frame_tx = proc_rxtx->frame_rx;
          proc_rxtx_ue->frame_rx = (subframe<4)?(proc_rxtx->frame_tx-1):(proc_rxtx->frame_tx);
          proc_rxtx_ue->subframe_tx = proc_rxtx->subframe_rx;
          proc_rxtx_ue->subframe_rx = (proc_rxtx->subframe_tx+6)%10;
989 990
          eNB->ulsch[0]->harq_processes[harq_pid]->round=round;
          UE->ulsch[0]->harq_processes[harq_pid]->round=round;
991

Thomas Laurent's avatar
Thomas Laurent committed
992
          if (n_frames==1) printf("filling ulsch: Trial %d : Round %d (subframe %d, frame %d)\n",trials,round,proc_rxtx_ue->subframe_tx,proc_rxtx_ue->frame_tx);
993

Thomas Laurent's avatar
Thomas Laurent committed
994 995
          round_trials[round]++;
          UL_req.sfn_sf = (1<<4)+subframe;
Cedric Roux's avatar
Cedric Roux committed
996

Thomas Laurent's avatar
Thomas Laurent committed
997
          if (n_frames==1) printf("filling ulsch: eNB prog frame %d, subframe %d (%d,%d)\n",proc_rxtx->frame_rx,subframe,sched_resp.frame,sched_resp.subframe);
998

Thomas Laurent's avatar
Thomas Laurent committed
999
          int modulation_type;
1000

Thomas Laurent's avatar
Thomas Laurent committed
1001 1002 1003 1004 1005 1006 1007
          if (mcs < 11)      modulation_type = 2;
          else if (mcs < 21) modulation_type = 4;
          else if (mcs < 29) modulation_type = 6;
          else {
            LOG_E(SIM,"mcs %i is not valid\n",mcs);
            exit(-1);
          }
1008

Thomas Laurent's avatar
Thomas Laurent committed
1009 1010 1011 1012 1013 1014
          fill_ulsch_dci(eNB, proc_rxtx->frame_rx, subframe, &sched_resp, 14,
                         (void *)&UL_alloc_pdu, first_rb,nb_rb, (round==0)?mcs:(28+rvidx[round]),
                         modulation_type, ndi, get_TBS_UL(mcs,nb_rb), cqi_flag, beta_CQI,
                         beta_RI, cqi_size);
          UE->ulsch_Msg3_active[eNB_id] = 0;
          UE->ul_power_control_dedicated[eNB_id].accumulationEnabled=1;
1015

Thomas Laurent's avatar
Thomas Laurent committed
1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033
          if (n_frames==1)
            printf("filling ulsch: ue prog SFN/SF %d/%d\n",proc_rxtx_ue->frame_rx,proc_rxtx_ue->subframe_rx);

          generate_ue_ulsch_params_from_dci((void *)&UL_alloc_pdu,
                                            14,
                                            (subframe+6)%10,
                                            format0,
                                            UE,
                                            proc_rxtx_ue,
                                            SI_RNTI,
                                            0,
                                            P_RNTI,
                                            CBA_RNTI,
                                            0,
                                            srs_flag);
          sched_resp.subframe=(subframe+6)%10;
          sched_resp.frame=(1024+eNB->proc.frame_rx+((subframe<4)?-1:0))&1023;
          schedule_response(&sched_resp);
1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051

          /////////////////////
          if (abstx) {
            if (trials==0 && round==0 && SNR==snr0) { //generate a new channel
              hold_channel = 0;
              flagMag=0;
            } else {
              hold_channel = 1;
              flagMag = 1;
            }
          } else {
            hold_channel = 0;
            flagMag=1;
          }

          ///////////////////////////////////////

          if (input_fdUL == NULL) {
Thomas Laurent's avatar
Thomas Laurent committed
1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062
            eNB->proc.frame_rx = 1;
            eNB->proc.subframe_rx = subframe;
            ru->proc.frame_rx = 1;
            ru->proc.subframe_rx = subframe;
            proc_rxtx_ue->frame_tx = proc_rxtx->frame_rx;
            proc_rxtx_ue->frame_rx = proc_rxtx->frame_tx;
            proc_rxtx_ue->subframe_tx = proc_rxtx->subframe_rx;
            proc_rxtx_ue->subframe_rx = proc_rxtx->subframe_tx;
            phy_procedures_UE_TX(UE,proc_rxtx_ue,0,0,normal_txrx);
            tx_lev = signal_energy(&UE->common_vars.txdata[0][eNB->frame_parms.samples_per_tti*subframe],
                                   eNB->frame_parms.samples_per_tti);
Cedric Roux's avatar
Cedric Roux committed
1063

1064
            if (n_frames==1) {
Mongazon's avatar
Mongazon committed
1065
              LOG_M("txsigF0UL.m","txsF0", &UE->common_vars.txdataF[0][eNB->frame_parms.ofdm_symbol_size*nsymb*subframe],eNB->frame_parms.ofdm_symbol_size*nsymb,1,
Thomas Laurent's avatar
Thomas Laurent committed
1066
                    1);
Mongazon's avatar
Mongazon committed
1067
              //LOG_M("txsigF1.m","txsF1", UE->common_vars.txdataF[0],FRAME_LENGTH_COMPLEX_SAMPLES_NO_PREFIX,1,1);
1068
            }
Thomas Laurent's avatar
Thomas Laurent committed
1069
          }  // input_fd == NULL
Cedric Roux's avatar
Cedric Roux committed
1070

1071 1072 1073
          tx_lev_dB = (unsigned int) dB_fixed_times10(tx_lev);

          if (n_frames==1) {
Mongazon's avatar
Mongazon committed
1074 1075
            LOG_M("txsig0UL.m","txs0", &txdata[0][eNB->frame_parms.samples_per_tti*subframe],2*frame_parms->samples_per_tti,1,1);
            //        LOG_M("txsig1UL.m","txs1", &txdata[1][eNB->frame_parms.samples_per_tti*subframe],2*frame_parms->samples_per_tti,1,1);
1076 1077 1078 1079
          }

          //AWGN
          //Set target wideband RX noise level to N0
1080
          sigma2_dB = N0;//-10*log10(UE->frame_parms.ofdm_symbol_size/(UE->frame_parms.N_RB_DL*12));//10*log10((double)tx_lev)  +10*log10(UE->frame_parms.ofdm_symbol_size/(UE->frame_parms.N_RB_DL*12)) - SNR;
1081 1082
          sigma2 = pow(10,sigma2_dB/10);
          // compute tx_gain to achieve target SNR (per resource element!)
1083
          tx_gain = sqrt(pow(10.0,.1*(N0+SNR))/(double)tx_lev);//*(nb_rb*12/(double)UE->frame_parms.ofdm_symbol_size)/(double)tx_lev);
1084

Thomas Laurent's avatar
Thomas Laurent committed
1085
          if (n_frames==1)
1086 1087 1088 1089
            printf("tx_lev = %d (%d.%d dB,%f), gain %f\n",tx_lev,tx_lev_dB/10,tx_lev_dB,10*log10((double)tx_lev),10*log10(tx_gain));

          // fill measurement symbol (19) with noise
          for (i=0; i<OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES; i++) {
1090
            for (aa=0; aa<eNB->frame_parms.nb_antennas_rx; aa++) {
Thomas Laurent's avatar
Thomas Laurent committed
1091 1092
              ((short *) &ru->common.rxdata[aa][(frame_parms->samples_per_tti<<1) -frame_parms->ofdm_symbol_size])[2*i] = (short) ((sqrt(sigma2/2)*gaussdouble(0.0,1.0)));
              ((short *) &ru->common.rxdata[aa][(frame_parms->samples_per_tti<<1) -frame_parms->ofdm_symbol_size])[2*i+1] = (short) ((sqrt(sigma2/2)*gaussdouble(0.0,1.0)));
1093 1094 1095 1096 1097
            }
          }

          // multipath channel

1098
          for (i=0; i<eNB->frame_parms.samples_per_tti; i++) {
1099
            for (aa=0; aa<1; aa++) {
1100 1101
              s_re[aa][i] = ((double)(((short *)&UE->common_vars.txdata[aa][eNB->frame_parms.samples_per_tti*subframe]))[(i<<1)]);
              s_im[aa][i] = ((double)(((short *)&UE->common_vars.txdata[aa][eNB->frame_parms.samples_per_tti*subframe]))[(i<<1)+1]);
1102 1103 1104 1105 1106 1107
            }
          }

          if (awgn_flag == 0) {
            if (UE2eNB->max_Doppler == 0) {
              multipath_channel(UE2eNB,s_re,s_im,r_re,r_im,
1108
                                eNB->frame_parms.samples_per_tti,hold_channel);
1109 1110
            } else {
              multipath_tv_channel(UE2eNB,s_re,s_im,r_re,r_im,
1111
                                   2*eNB->frame_parms.samples_per_tti,hold_channel);
1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129
            }
          }

          if(abstx) {
            if(saving_bler==0)
              if (trials==0 && round==0) {
                // calculate freq domain representation to compute SINR
                freq_channel(UE2eNB, N_RB_DL,12*N_RB_DL + 1);
                // snr=pow(10.0,.1*SNR);
                fprintf(csv_fdUL,"%f,%d,%d,%f,%f,%f,",SNR,tx_lev,tx_lev_dB,sigma2_dB,tx_gain,SNR2);

                //fprintf(csv_fdUL,"%f,",SNR);
                for (u=0; u<12*nb_rb; u++) {
                  for (aarx=0; aarx<UE2eNB->nb_rx; aarx++) {
                    for (aatx=0; aatx<UE2eNB->nb_tx; aatx++) {
                      // abs_channel = (eNB2UE->chF[aarx+(aatx*eNB2UE->nb_rx)][u].x*eNB2UE->chF[aarx+(aatx*eNB2UE->nb_rx)][u].x + eNB2UE->chF[aarx+(aatx*eNB2UE->nb_rx)][u].y*eNB2UE->chF[aarx+(aatx*eNB2UE->nb_rx)][u].y);
                      channelx = UE2eNB->chF[aarx+(aatx*UE2eNB->nb_rx)][u].x;
                      channely = UE2eNB->chF[aarx+(aatx*UE2eNB->nb_rx)][u].y;
Thomas Laurent's avatar
Thomas Laurent committed
1130
                      // if(transmission_m==5){
1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145
                      fprintf(csv_fdUL,"%e+i*(%e),",channelx,channely);
                      // }
                      // else{
                      //  pilot_sinr = 10*log10(snr*abs_channel);
                      //  fprintf(csv_fd,"%e,",pilot_sinr);
                      // }
                    }
                  }
                }
              }
          }

          if (n_frames==1)
            printf("Sigma2 %f (sigma2_dB %f), tx_gain %f (%f dB)\n",sigma2,sigma2_dB,tx_gain,20*log10(tx_gain));

1146 1147
          for (i=0; i<eNB->frame_parms.samples_per_tti; i++) {
            for (aa=0; aa<eNB->frame_parms.nb_antennas_rx; aa++) {
Thomas Laurent's avatar
Thomas Laurent committed
1148 1149 1150 1151 1152
              ((short *) &ru->common.rxdata[aa][eNB->frame_parms.samples_per_tti*subframe])[2*i] =
                (short) ((tx_gain*r_re[aa][i]) + sqrt(sigma2/2)*gaussdouble(0.0,1.0));
              ((short *) &ru->common.rxdata[aa][eNB->frame_parms.samples_per_tti*subframe])[2*i+1] =
                (short) ((tx_gain*r_im[aa][i]) + (iqim*tx_gain*r_re[aa][i]) +
                         sqrt(sigma2/2)*gaussdouble(0.0,1.0));
1153 1154 1155
            }
          }

1156
          if (n_frames<=10) {
Thomas Laurent's avatar
Thomas Laurent committed
1157 1158 1159 1160 1161 1162
            printf("rx_level Null symbol %f\n",10*log10((double)signal_energy((int *)
                   &ru->common.rxdata[0][(eNB->frame_parms.samples_per_tti<<1) -
                                         eNB->frame_parms.ofdm_symbol_size],OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2)));
            printf("rx_level data symbol %f\n",
                   10*log10(signal_energy((int *)&ru->common.rxdata[0][160+(eNB->frame_parms.samples_per_tti*subframe)],
                                          OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2)));
1163 1164
          }

Thomas Laurent's avatar
Thomas Laurent committed
1165 1166
          SNRmeas = 10*log10(((double)signal_energy((int *)&ru->common.rxdata[0][160+(eNB->frame_parms.samples_per_tti*subframe)],
                              OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2))/((double)signal_energy((int *)
1167
                                  &ru->common.rxdata[0][(eNB->frame_parms.samples_per_tti<<1) -eNB->frame_parms.ofdm_symbol_size],
1168
                                  OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2)) - 1)+10*log10(eNB->frame_parms.N_RB_UL/nb_rb);
1169

1170
          if (n_frames<=10) {
1171
            printf("SNRmeas %f\n",SNRmeas);
Thomas Laurent's avatar
Thomas Laurent committed
1172
            LOG_M("rxsig0UL.m","rxs0", &ru->common.rxdata[0][eNB->frame_parms.samples_per_tti*subframe],eNB->frame_parms.samples_per_tti,1,1);
1173

Thomas Laurent's avatar
Thomas Laurent committed
1174
            if (eNB->frame_parms.nb_antennas_rx>1) LOG_M("rxsig1UL.m","rxs1", &ru->common.rxdata[1][eNB->frame_parms.samples_per_tti*subframe],eNB->frame_parms.samples_per_tti,1,1);
1175 1176
          }

Thomas Laurent's avatar
Thomas Laurent committed
1177
	  start_meas(&eNB->phy_proc_rx);
Thomas Laurent's avatar
Thomas Laurent committed
1178 1179 1180 1181
          ru->feprx = (parallel_flag == 1) ? ru_fep_full_2thread        : fep_full;
          eNB->td  = (parallel_flag == 1) ? ulsch_decoding_data_2thread : ulsch_decoding_data;
          ru->feprx(ru);
          phy_procedures_eNB_uespec_RX(eNB,proc_rxtx);
Thomas Laurent's avatar
Thomas Laurent committed
1182 1183
	  stop_meas(&eNB->phy_proc_rx);
	  
1184 1185 1186
          if (cqi_flag > 0) {
            cqi_error = 0;

1187
            if (eNB->ulsch[0]->harq_processes[harq_pid]->Or1 < 32) {
1188
              for (i=2; i<4; i++) {
1189 1190
                //                printf("cqi %d : %d (%d)\n",i,eNB->ulsch[0]->o[i],UE->ulsch[0]->o[i]);
                if (eNB->ulsch[0]->harq_processes[harq_pid]->o[i] != UE->ulsch[0]->o[i])
1191 1192 1193 1194 1195 1196 1197 1198
                  cqi_error = 1;
              }
            } else {
            }

            if (cqi_error == 1) {
              cqi_errors++;

1199
              if (eNB->ulsch[0]->harq_processes[harq_pid]->cqi_crc_status == 1)
1200 1201
                cqi_crc_falsepositives++;
            } else {
1202
              if (eNB->ulsch[0]->harq_processes[harq_pid]->cqi_crc_status == 0)
1203 1204 1205 1206
                cqi_crc_falsenegatives++;
            }
          }

1207
          if (eNB->ulsch[0]->harq_processes[harq_pid]->o_ACK[0] != UE->ulsch[0]->o_ACK[0])
1208 1209
            ack_errors++;

1210
          //    printf("ulsch_coding: O[%d] %d\n",i,o_flip[i]);
Thomas Laurent's avatar
Thomas Laurent committed
1211 1212
          //          if (ret <= eNB->ulsch[0]->max_turbo_iterations) {
          iter_trials++;
1213

Thomas Laurent's avatar
Thomas Laurent committed
1214
          if (eNB->ulsch[0]->harq_processes[harq_pid]->status == SCH_IDLE) {
1215
            if (n_frames==1) {
1216
              printf("No ULSCH errors found, o_ACK[0]= %d, cqi_crc_status=%d\n",eNB->ulsch[0]->harq_processes[harq_pid]->o_ACK[0],eNB->ulsch[0]->harq_processes[harq_pid]->cqi_crc_status);
1217

1218 1219 1220
              if (eNB->ulsch[0]->harq_processes[harq_pid]->cqi_crc_status==1)
                print_CQI(eNB->ulsch[0]->harq_processes[harq_pid]->o,
                          eNB->ulsch[0]->harq_processes[harq_pid]->uci_format,0,eNB->frame_parms.N_RB_DL);
1221

1222
              dump_ulsch(eNB,eNB->proc.frame_rx,subframe,0,round);
1223 1224 1225 1226 1227 1228 1229 1230
              exit(-1);
            }

            round=5;
          } else {
            errs[round]++;

            if (n_frames==1) {
1231
              printf("ULSCH errors found o_ACK[0]= %d\n",eNB->ulsch[0]->harq_processes[harq_pid]->o_ACK[0]);
1232

1233 1234 1235
              for (s=0; s<eNB->ulsch[0]->harq_processes[harq_pid]->C; s++) {
                if (s<eNB->ulsch[0]->harq_processes[harq_pid]->Cminus)
                  Kr = eNB->ulsch[0]->harq_processes[harq_pid]->Kminus;
1236
                else
1237
                  Kr = eNB->ulsch[0]->harq_processes[harq_pid]->Kplus;
1238 1239 1240 1241 1242

                Kr_bytes = Kr>>3;
                printf("Decoded_output (Segment %d):\n",s);

                for (i=0; i<Kr_bytes; i++)
1243 1244
                  printf("%d : %x (%x)\n",i,eNB->ulsch[0]->harq_processes[harq_pid]->c[s][i],
                         eNB->ulsch[0]->harq_processes[harq_pid]->c[s][i]^UE->ulsch[0]->harq_processes[harq_pid]->c[s][i]);
1245 1246
              }

1247
              dump_ulsch(eNB,eNB->proc.frame_rx,subframe,0,round);
Thomas Laurent's avatar
Thomas Laurent committed
1248

1249
              if (round == 4) exit(-1);
1250 1251
            }

Thomas Laurent's avatar
Thomas Laurent committed
1252 1253
            if (n_frames==1) printf("round %d errors %d/%d\n",round,errs[round],trials);

1254 1255 1256 1257 1258
            round++;

            if (n_frames==1) {
              printf("ULSCH in error in round %d\n",round);
            }
Cedric Roux's avatar