ulsim.c 67.6 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
THREAD_STRUCT thread_struct;
100

101
void
Cedric Roux's avatar
Cedric Roux committed
102
fill_nfapi_ulsch_config_request(nfapi_ul_config_request_pdu_t *ul_config_pdu,
Thomas Laurent's avatar
Thomas Laurent committed
103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124
                                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) {
125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144
  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
145 146
  //printf("Filling ul_config_pdu : Q %d, TBS %d, rv %d, ndi %d\n", modulation_type,size,redundancy_version,new_data_indication);

147 148 149 150 151 152 153 154 155 156
  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
157
        && (tmode == 3 || tmode == 4 || tmode == 8 || tmode == 9 || tmode == 10))
158 159 160 161
      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
162
    for (int ri = 0;
Thomas Laurent's avatar
Thomas Laurent committed
163 164 165
         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
166

167 168 169 170
    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;
  }
}
171

172
void fill_ulsch_dci(PHY_VARS_eNB *eNB,
Thomas Laurent's avatar
Thomas Laurent committed
173 174 175 176 177 178 179 180 181 182
                    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
183
                    int TBS,
Thomas Laurent's avatar
Thomas Laurent committed
184 185 186 187
                    int cqi_flag,
                    uint8_t beta_CQI,
                    uint8_t beta_RI,
                    uint8_t cqi_size) {
188 189 190
  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
191
  //printf("ulsch in frame %d, subframe %d => harq_pid %d, mcs %d, ndi %d\n",frame,subframe,harq_pid,mcs,ndi);
192 193

  switch (eNB->frame_parms.N_RB_UL) {
Thomas Laurent's avatar
Thomas Laurent committed
194 195
    case 6:
      break;
196

Thomas Laurent's avatar
Thomas Laurent committed
197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217
    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;
      }
218

Thomas Laurent's avatar
Thomas Laurent committed
219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265
      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;
      }
266

Thomas Laurent's avatar
Thomas Laurent committed
267
      break;
268

Thomas Laurent's avatar
Thomas Laurent committed
269 270
    default:
      break;
271 272
  }

Cedric Roux's avatar
Cedric Roux committed
273
  fill_nfapi_ulsch_config_request(&ul_req->ul_config_pdu_list[0],
Thomas Laurent's avatar
Thomas Laurent committed
274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297
                                  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);
298 299 300
  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
301
}
Cedric Roux's avatar
Cedric Roux committed
302

Thomas Laurent's avatar
Thomas Laurent committed
303
enum eTypes { eBool, eInt, eFloat, eText };
304
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;
305

Thomas Laurent's avatar
Thomas Laurent committed
306
int main(int argc, char **argv) {
307
  int i,j,aa,u;
308 309 310
  PHY_VARS_eNB *eNB;
  PHY_VARS_UE *UE;
  RU_t *ru;
311 312
  int aarx,aatx;
  double channelx,channely;
Thomas Laurent's avatar
Thomas Laurent committed
313 314
  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;
315
  double blerr;
Thomas Laurent's avatar
Thomas Laurent committed
316
  int rvidx[8]= {0,2,3,1,0,2,3,1};
317 318
  int **txdata;
  LTE_DL_FRAME_PARMS *frame_parms;
319 320
  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];
321 322
  double r_re2[30720],r_im2[30720];
  double r_re3[30720],r_im3[30720];
Thomas Laurent's avatar
Thomas Laurent committed
323 324 325 326
  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};
327 328
  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
329
  int cqi_error,cqi_errors,ack_errors,cqi_crc_falsepositives,cqi_crc_falsenegatives;
330 331 332 333
  int ch_realization;
  int eNB_id = 0;
  int chMod = 0 ;
  int UE_id = 0;
Thomas Laurent's avatar
Thomas Laurent committed
334
  static int nb_rb=25,first_rb=0,mcs=0,round=0;
335
  //unsigned char l;
Thomas Laurent's avatar
Thomas Laurent committed
336
  static int awgn_flag = 0 ;
337
  SCM_t channel_model=Rice1;
Thomas Laurent's avatar
Thomas Laurent committed
338
  unsigned char *input_buffer=0,harq_pid;
339 340 341
  unsigned short input_buffer_length;
  unsigned int ret;
  unsigned int coded_bits_per_codeword,nsymb;
342
  unsigned int tx_lev=0,tx_lev_dB,trials,errs[4]= {0,0,0,0},round_trials[4]= {0,0,0,0};
343
  FILE *bler_fd=NULL;
344
  char bler_fname[512];
345 346
  FILE *time_meas_fd=NULL;
  char time_meas_fname[256];
347 348 349 350
  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;
351
  FILE *csv_fdUL=NULL;
352
  /*
Thomas Laurent's avatar
Thomas Laurent committed
353 354
    FILE *fperen=NULL;
    char fperen_name[512];
355

Thomas Laurent's avatar
Thomas Laurent committed
356 357
    FILE *fmageren=NULL;
    char fmageren_name[512];
358

Thomas Laurent's avatar
Thomas Laurent committed
359 360
    FILE *flogeren=NULL;
    char flogeren_name[512];
361
  */
362 363 364 365
  /* FILE *ftxlev;
     char ftxlev_name[512];
  */
  char csv_fname[512];
Thomas Laurent's avatar
Thomas Laurent committed
366 367 368
  static int n_frames=5000;
  static int n_ch_rlz = 1;
  static int abstx = 0;
369
  int hold_channel=0;
370
  channel_desc_t *UE2eNB;
371
  //uint8_t control_only_flag = 0;
Thomas Laurent's avatar
Thomas Laurent committed
372 373 374 375
  static int delay = 0;
  static double maxDoppler = 0.0;
  static int srs_flag = 0;
  static int N_RB_DL=25,osf=1;
376
  //uint8_t cyclic_shift = 0;
Thomas Laurent's avatar
Thomas Laurent committed
377 378 379 380
  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;
381
  double cpu_freq_GHz;
Thomas Laurent's avatar
Thomas Laurent committed
382
  int iter_trials;
383 384 385
  uint32_t UL_alloc_pdu;
  int s,Kr,Kr_bytes;
  int dump_perf=0;
Thomas Laurent's avatar
Thomas Laurent committed
386
  static int dump_table =0;
387
  double effective_rate=0.0;
Thomas Laurent's avatar
Thomas Laurent committed
388 389
  char channel_model_input[10]= {0};
  static int max_turbo_iterations=4;
390
  int nb_rb_set = 0;
391
  int sf;
Thomas Laurent's avatar
Thomas Laurent committed
392
  static int threequarter_fs=0;
393
  int ndi;
knopp's avatar
knopp committed
394
  opp_enabled=1; // to enable the time meas
395 396 397 398
  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
399 400 401 402
  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));
403 404
  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;
405
  cpu_freq_GHz = (double)get_cpu_freq_GHz();
406
  cpuf = cpu_freq_GHz;
407 408
  set_parallel_conf("PARALLEL_SINGLE_THREAD");

409
  printf("Detected cpu_freq %f GHz\n",cpu_freq_GHz);
410
  AssertFatal(load_configmodule(argc,argv,CONFIG_ENABLECMDLINEONLY) != NULL,
Thomas Laurent's avatar
Thomas Laurent committed
411
              "cannot load configuration module, exiting\n");
412
  logInit();
Calisson's avatar
Calisson committed
413 414 415
  // 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
416 417 418
  //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
419 420
  //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[] = {
421
    { "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
422 423 424 425
    { "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 },
426
    { "wsnrInterrupt", "snr int ?", 0, dblptr:&snr_int,  defdblval:30, TYPE_DOUBLE, 0 },
427
    { "e_snr_step", "step increasing snr",0, dblptr:&input_snr_step,  defdblval:0.2, TYPE_DOUBLE, 0 },
Thomas Laurent's avatar
Thomas Laurent committed
428 429 430 431 432 433 434 435
    { "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 },
436
    { "Lparallel", "Enable parallel execution",0, strptr:NULL,  defstrval:NULL, TYPE_STRING,  0 },
Thomas Laurent's avatar
Thomas Laurent committed
437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453
    { "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 },
454
    { "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
455
    { "verbose", "display debug text", PARAMFLAG_BOOL,  iptr:&verbose, defintval:0, TYPE_INT, 0 },
456
    { "help", "display help and exit", PARAMFLAG_BOOL,  iptr:&help, defintval:0, TYPE_INT, 0 },
Thomas Laurent's avatar
Thomas Laurent committed
457 458
    { "", "",0,  iptr:NULL, defintval:0, TYPE_INT, 0 },
  };
459
  struct option * long_options = parse_oai_options(options); 
460

Thomas Laurent's avatar
Thomas Laurent committed
461
  int option_index;
462

Thomas Laurent's avatar
Thomas Laurent committed
463
  int res;
464

Thomas Laurent's avatar
Thomas Laurent committed
465 466 467 468 469
  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) {
470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488
	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
489
        }
490

Thomas Laurent's avatar
Thomas Laurent committed
491 492
      continue;
    }
493

Thomas Laurent's avatar
Thomas Laurent committed
494
    switch (long_options[option_index].name[0]) {
495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 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
    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;
550

Thomas Laurent's avatar
Thomas Laurent committed
551 552 553
      //case 'c':
      //  cyclic_shift = atoi(optarg);
      //  break;
554

555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580
    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;

581 582 583 584
    case 'L':
      set_parallel_conf(optarg);
      break;
      
585 586 587 588
    default:
      printf("Wrong option: %s\n",long_options[option_index].name);
      exit(1);
      break;
Thomas Laurent's avatar
Thomas Laurent committed
589 590
    }
  }
591

Thomas Laurent's avatar
Thomas Laurent committed
592 593 594 595
  if ( res != -1 ) {
    printf("A wrong option has been found\n");
    exit(1);
  }
596 597


598 599 600 601 602
  if (help || verbose )
     display_options_values(options, true);
  if (help)
    exit(0);
  
603 604
  if (thread_struct.parallel_conf != PARALLEL_SINGLE_THREAD)
    set_worker_conf("WORKER_ENABLE");
605 606 607
  RC.nb_L1_inst = 1;
  RC.nb_RU = 1;
  lte_param_init(&eNB,&UE,&ru,
608
                 1,
Thomas Laurent's avatar
Thomas Laurent committed
609 610 611 612 613 614 615 616 617 618 619 620 621
                 1,
                 n_rx,
                 1,
                 1,
                 extended_prefix_flag,
                 frame_type,
                 0,
                 tdd_config,
                 N_RB_DL,
                 4,
                 threequarter_fs,
                 osf,
                 0);
622 623 624 625 626 627
  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
628
  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
629

Thomas Laurent's avatar
Thomas Laurent committed
630
  memset((void *)&eNB->UL_INFO,0,sizeof(eNB->UL_INFO));
631 632 633 634 635 636 637 638
  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
639
  // 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
640 641
  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
642 643
  PHY_vars_UE_g[0][0] = UE;

644
  if (nb_rb_set == 0)
645
    nb_rb = eNB->frame_parms.N_RB_UL;
646

647
  printf("1 . rxdataF_comp[0] %p\n",eNB->pusch_vars[0]->rxdataF_comp[0]);
648
  printf("Setting mcs = %d\n",mcs);
649
  printf("n_frames = %d\n", n_frames);
knopp's avatar
knopp committed
650
  snr1 = snr0+snr_int;
651
  printf("SNR0 %f, SNR1 %f\n",snr0,snr1);
652 653 654
  frame_parms = &eNB->frame_parms;
  txdata = UE->common_vars.txdata;
  nsymb = (eNB->frame_parms.Ncp == NORMAL) ? 14 : 12;
655
  sprintf(bler_fname,"ULbler_mcs%d_nrb%d_ChannelModel%d_nsim%d.csv",mcs,nb_rb,chMod,n_frames);
656
  bler_fd = fopen(bler_fname,"w");
Thomas Laurent's avatar
Thomas Laurent committed
657

658 659 660 661
  if (bler_fd==NULL) {
    fprintf(stderr,"Problem creating file %s\n",bler_fname);
    exit(-1);
  }
662

663
  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");
664 665

  if (test_perf != 0) {
666 667 668 669
    char hostname[1024];
    hostname[1023] = '\0';
    gethostname(hostname, 1023);
    printf("Hostname: %s\n", hostname);
670 671 672 673
    //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
674
            N_RB_DL,mcs,n_rx,channel_model_input,transmission_m);
675
    time_meas_fd = fopen(time_meas_fname,"w");
Thomas Laurent's avatar
Thomas Laurent committed
676

677 678 679 680
    if (time_meas_fd==NULL) {
      fprintf(stderr,"Cannot create file %s!\n",time_meas_fname);
      exit(-1);
    }
681
  }
682 683 684

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

688 689 690 691
    if (csv_fdUL == NULL) {
      fprintf(stderr,"Problem opening file %s\n",csv_fname);
      exit(-1);
    }
Thomas Laurent's avatar
Thomas Laurent committed
692

693 694
    fprintf(csv_fdUL,"data_all%d=[",mcs);
  }
695

696 697 698 699 700 701
  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);
  }
702

Bilel's avatar
Bilel committed
703
  UE->pdcch_vars[0][0]->crnti = 14;
704
  UE->frame_parms.soundingrs_ul_config_common.enabled_flag = srs_flag;
705
  UE->frame_parms.soundingrs_ul_config_common.srs_BandwidthConfig = 2;
706 707 708 709
  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;
710 711 712
  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;
713
  UE->soundingrs_ul_config_dedicated[eNB_id].cyclicShift = 0;
714
  eNB->frame_parms.soundingrs_ul_config_common.enabled_flag = srs_flag;
715
  eNB->frame_parms.soundingrs_ul_config_common.srs_BandwidthConfig = 2;
716 717 718 719
  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;
720 721 722
  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;
723
  eNB->soundingrs_ul_config_dedicated[UE_id].cyclicShift = 0;
724 725 726 727 728 729 730
  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;
731 732 733
  // 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;
734
  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);
735 736
  UE2eNB = new_channel_desc_scm(1,
                                n_rx,
737
                                channel_model,
Thomas Laurent's avatar
Thomas Laurent committed
738 739
                                N_RB2sampling_rate(eNB->frame_parms.N_RB_UL),
                                N_RB2channel_bandwidth(eNB->frame_parms.N_RB_UL),
740 741 742 743 744 745 746
                                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
747 748
  for (int k=0; k<NUMBER_OF_UE_MAX; k++) eNB->ulsch[k] = new_eNB_ulsch(max_turbo_iterations,N_RB_DL,0);

749
  UE->ulsch[0]   = new_ue_ulsch(N_RB_DL,0);
750
  printf("ULSCH %p\n",UE->ulsch[0]);
751

752
  if(get_thread_worker_conf() == WORKER_ENABLE) {
753
    extern void init_fep_thread(PHY_VARS_eNB *, pthread_attr_t *);
754
    extern void init_td_thread(PHY_VARS_eNB *);
755
    init_fep_thread(eNB,NULL);
756
    init_td_thread(eNB);
757
  }
Thomas Laurent's avatar
Thomas Laurent committed
758

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

763
    if (!eNB->dlsch[0][i]) {
764 765 766
      printf("Can't get eNB dlsch structures\n");
      exit(-1);
    }
Thomas Laurent's avatar
Thomas Laurent committed
767

768
    eNB->dlsch[0][i]->rnti = 14;
769
  }
Thomas Laurent's avatar
Thomas Laurent committed
770

771
  /* allocate memory for both subframes (only one is really used
Thomas Laurent's avatar
Thomas Laurent committed
772 773 774 775
     but there is now "copy_harq_proc_struct" which needs both
     to be valid)
     TODO: refine this somehow (necessary?)
  */
776 777 778
  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
779

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

785 786 787
      UE->dlsch[sf][0][i]->rnti   = 14;
    }
  }
788

789 790
  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);
791 792
  UE->measurements.rank[0] = 0;
  UE->transmission_mode[0] = 2;
Thomas Laurent's avatar
Thomas Laurent committed
793
  UE->pucch_config_dedicated[0].tdd_AckNackFeedbackMode = disable_bundling == 0 ? bundling : multiplexing;
794
  eNB->transmission_mode[0] = 2;
Thomas Laurent's avatar
Thomas Laurent committed
795
  eNB->pucch_config_dedicated[0].tdd_AckNackFeedbackMode = disable_bundling == 0 ? bundling : multiplexing;
796 797 798 799 800 801
  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;
802 803 804
  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];
805 806 807 808
  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);
809
  proc_rxtx_ue->frame_tx = proc_rxtx->frame_rx;
810
  proc_rxtx_ue->frame_rx = (subframe<4)?(proc_rxtx->frame_tx-1):(proc_rxtx->frame_tx);
811
  proc_rxtx_ue->subframe_tx = proc_rxtx->subframe_rx;
812
  proc_rxtx_ue->subframe_rx = (proc_rxtx->subframe_tx+6)%10;
813
  printf("Init UL hopping UE\n");
814
  init_ul_hopping(&UE->frame_parms);
815
  printf("Init UL hopping eNB\n");
816
  init_ul_hopping(&eNB->frame_parms);
817
  UE->dlsch[subframe&1][0][0]->harq_ack[ul_subframe2pdcch_alloc_subframe(&eNB->frame_parms,subframe)].send_harq_status = 1;
818 819
  UE->ulsch_Msg3_active[eNB_id] = 0;
  UE->ul_power_control_dedicated[eNB_id].accumulationEnabled=1;
820 821
  coded_bits_per_codeword = nb_rb * (12 * get_Qm_ul(mcs)) * nsymb;

822
  if (cqi_flag == 1) coded_bits_per_codeword-=UE->ulsch[0]->O;
823 824 825

  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);
826

827
  for (ch_realization=0; ch_realization<n_ch_rlz; ch_realization++) {
828 829 830 831 832 833
    /*
      if(abstx){
      int ulchestim_f[300*12];
      int ulchestim_t[2*(frame_parms->ofdm_symbol_size)];
      }
    */
834
    if(abstx) {
835 836 837
      printf("**********************Channel Realization Index = %d **************************\n", ch_realization);
      saving_bler=1;
    }
838

knopp's avatar
knopp committed
839
    //    if ((subframe>5) || (subframe < 4))
840
    //      UE->frame++;
841 842

    for (SNR=snr0; SNR<snr1; SNR+=input_snr_step) {
843 844 845 846 847 848 849 850 851 852 853 854 855 856
      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);
857
      harq_pid = subframe2harq_pid(&UE->frame_parms,proc_rxtx_ue->frame_tx,subframe);
858
      input_buffer_length = UE->ulsch[0]->harq_processes[harq_pid]->TBS/8;
Thomas Laurent's avatar
Thomas Laurent committed
859 860 861 862

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

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

865
      //      printf("UL frame %d/subframe %d, harq_pid %d\n",UE->frame,subframe,harq_pid);
866
      if (input_fdUL == NULL) {
867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886
        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)) {
887
          ret=fscanf(input_fdUL,"%s %s",input_val_str,input_val_str2);//&input_val1,&input_val2);
Thomas Laurent's avatar
Thomas Laurent committed
888

889
          if (ret != 2) printf("ERROR: error reading file\n");
890 891

          if ((i%4)==0) {
Thomas Laurent's avatar
Thomas Laurent committed
892 893
            ((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));
894 895

            if ((i/4)<100)
Thomas Laurent's avatar
Thomas Laurent committed
896
              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,);
897 898 899 900 901 902 903 904 905
          }

          i++;

          if (i>(FRAME_LENGTH_SAMPLES))
            break;
        }

        printf("Read in %d samples\n",i/4);
Mongazon's avatar
Mongazon committed
906 907
        //      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);
908 909 910
        tx_lev = signal_energy(&txdata[0][0],
                               OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES);
        tx_lev_dB = (unsigned int) dB_fixed(tx_lev);
911 912
      }

913
      iter_trials=0;
914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931
      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
932
      reset_meas(&eNB->ulsch_demodulation_stats);
933 934 935 936 937 938 939
      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
940
      // initialization
Thomas Laurent's avatar
Thomas Laurent committed
941 942 943 944 945 946 947 948
      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));
949
      ndi=0;
950
      phy_reset_ue(0,0,0);
951
      UE->UE_mode[eNB_id]=PUSCH;
Thomas Laurent's avatar
Thomas Laurent committed
952
      SET_LOG_DEBUG(UE_TIMING);
953

954 955
      for (trials = 0; trials<n_frames; trials++) {
        //      printf("*");
956 957
        //        UE->frame++;
        //        eNB->frame++;
Thomas Laurent's avatar
Thomas Laurent committed
958
        ndi = (1-ndi);
959 960 961 962
        fflush(stdout);
        round=0;

        while (round < 4) {
Thomas Laurent's avatar
Thomas Laurent committed
963 964 965 966 967 968 969 970
          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;
971 972
          eNB->ulsch[0]->harq_processes[harq_pid]->round=round;
          UE->ulsch[0]->harq_processes[harq_pid]->round=round;
973

Thomas Laurent's avatar
Thomas Laurent committed
974
          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);
975

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

Thomas Laurent's avatar
Thomas Laurent committed
979
          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);
980

Thomas Laurent's avatar
Thomas Laurent committed
981
          int modulation_type;
982

Thomas Laurent's avatar
Thomas Laurent committed
983 984 985 986 987 988 989
          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);
          }
990

Thomas Laurent's avatar
Thomas Laurent committed
991 992 993 994 995 996
          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;
997

Thomas Laurent's avatar
Thomas Laurent committed
998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015
          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);
1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033

          /////////////////////
          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
1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044
            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
1045

1046
            if (n_frames==1) {
Mongazon's avatar
Mongazon committed
1047
              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
1048
                    1);
Mongazon's avatar
Mongazon committed
1049
              //LOG_M("txsigF1.m","txsF1", UE->common_vars.txdataF[0],FRAME_LENGTH_COMPLEX_SAMPLES_NO_PREFIX,1,1);
1050
            }
Thomas Laurent's avatar
Thomas Laurent committed
1051
          }  // input_fd == NULL
Cedric Roux's avatar
Cedric Roux committed
1052

1053 1054 1055
          tx_lev_dB = (unsigned int) dB_fixed_times10(tx_lev);

          if (n_frames==1) {
Mongazon's avatar
Mongazon committed
1056 1057
            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);
1058 1059 1060 1061
          }

          //AWGN
          //Set target wideband RX noise level to N0
1062
          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;
1063 1064
          sigma2 = pow(10,sigma2_dB/10);
          // compute tx_gain to achieve target SNR (per resource element!)
1065
          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);
1066

Thomas Laurent's avatar
Thomas Laurent committed
1067
          if (n_frames==1)
1068 1069 1070 1071
            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++) {
1072
            for (aa=0; aa<eNB->frame_parms.nb_antennas_rx; aa++) {
Thomas Laurent's avatar
Thomas Laurent committed
1073 1074
              ((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)));
1075 1076 1077 1078 1079
            }
          }

          // multipath channel

1080
          for (i=0; i<eNB->frame_parms.samples_per_tti; i++) {
1081
            for (aa=0; aa<1; aa++) {
1082 1083
              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]);
1084 1085 1086 1087 1088 1089
            }
          }

          if (awgn_flag == 0) {
            if (UE2eNB->max_Doppler == 0) {
              multipath_channel(UE2eNB,s_re,s_im,r_re,r_im,
1090
                                eNB->frame_parms.samples_per_tti,hold_channel);
1091 1092
            } else {
              multipath_tv_channel(UE2eNB,s_re,s_im,r_re,r_im,
1093
                                   2*eNB->frame_parms.samples_per_tti,hold_channel);
1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111
            }
          }

          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
1112
                      // if(transmission_m==5){
1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127
                      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));

1128 1129
          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
1130 1131 1132 1133 1134
              ((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));
1135 1136 1137
            }
          }

1138
          if (n_frames<=10) {
Thomas Laurent's avatar
Thomas Laurent committed
1139 1140 1141 1142 1143 1144
            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)));
1145 1146
          }

Thomas Laurent's avatar
Thomas Laurent committed
1147 1148
          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 *)
1149
                                  &ru->common.rxdata[0][(eNB->frame_parms.samples_per_tti<<1) -eNB->frame_parms.ofdm_symbol_size],
1150
                                  OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2)) - 1)+10*log10(eNB->frame_parms.N_RB_UL/nb_rb);
1151

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

Thomas Laurent's avatar
Thomas Laurent committed
1156
            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);
1157 1158
          }

Thomas Laurent's avatar
Thomas Laurent committed
1159
	  start_meas(&eNB->phy_proc_rx);
1160
          ru->feprx = (get_thread_worker_conf() == WORKER_ENABLE) ? ru_fep_full_2thread        : fep_full;
Thomas Laurent's avatar
Thomas Laurent committed
1161 1162
          ru->feprx(ru);
          phy_procedures_eNB_uespec_RX(eNB,proc_rxtx);
Thomas Laurent's avatar
Thomas Laurent committed
1163 1164
	  stop_meas(&eNB->phy_proc_rx);
	  
1165 1166 1167
          if (cqi_flag > 0) {
            cqi_error = 0;

1168
            if (eNB->ulsch[0]->harq_processes[harq_pid]->Or1 < 32) {
1169
              for (i=2; i<4; i++) {
1170 1171
                //                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])
1172 1173 1174 1175 1176 1177 1178 1179
                  cqi_error = 1;
              }
            } else {
            }

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

1180
              if (eNB->ulsch[0]->harq_processes[harq_pid]->cqi_crc_status == 1)
1181 1182
                cqi_crc_falsepositives++;
            } else {
1183
              if (eNB->ulsch[0]->harq_processes[harq_pid]->cqi_crc_status == 0)
1184 1185 1186 1187
                cqi_crc_falsenegatives++;
            }
          }

1188
          if (eNB->ulsch[0]->harq_processes[harq_pid]->o_ACK[0] != UE->ulsch[0]->o_ACK[0])
1189 1190
            ack_errors++;

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

Thomas Laurent's avatar
Thomas Laurent committed
1195
          if (eNB->ulsch[0]->harq_processes[harq_pid]->status == SCH_IDLE) {
1196
            if (n_frames==1) {
1197
              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);
1198

1199 1200 1201
              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);
1202

1203
              dump_ulsch(eNB,eNB->proc.frame_rx,subframe,0,round);
1204 1205 1206 1207 1208 1209 1210 1211
              exit(-1);
            }

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

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

1214 1215 1216
              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;
1217
                else
1218
                  Kr = eNB->ulsch[0]->harq_processes[harq_pid]->Kplus;
1219 1220 1221 1222 1223

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

                for (i=0; i<Kr_bytes; i++)
1224 1225
                  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]);
1226 1227
              }

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

1230
              if (round == 4) exit(-1);
1231 1232
            }

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

1235 1236 1237 1238 1239
            round++;

            if (n_frames==1) {
              printf("ULSCH in error in round %d\n",round);
            }
Cedric Roux's avatar
Cedric Roux committed
1240
          }  // ulsch error
1241 1242 1243 1244 1245 1246
        } // round

        //      printf("\n");
        if ((errs[0]>=100) && (trials>(n_frames/2)))
          break;

Thomas Laurent's avatar
Thomas Laurent committed
1247 1248
        if (xforms==1)
          phy_scope_eNB(form_enb,eNB,0);
1249

Thomas Laurent's avatar
Thomas Laurent committed
1250 1251 1252 1253 1254 1255 1256 1257
        double t_tx = inMicroS(UE->phy_proc_tx.p_time);
        double t_tx_ifft = inMicroS(UE->ofdm_mod_stats.p_time);
        double t_tx_mod = inMicroS(UE->ulsch_modulation_stats.p_time);
        double t_tx_enc = inMicroS(UE->ulsch_encoding_stats.p_time);
        double t_rx = inMicroS(eNB->phy_proc_rx.p_time);
        double t_rx_fft = inMicroS(ru->ofdm_demod_stats.p_time);
        double t_rx_demod = inMicroS(eNB->ulsch_demodulation_stats.p_time);
        double t_rx_dec = inMicroS(e