dlsim.c 190 KB
Newer Older
1
/*******************************************************************************
2
    OpenAirInterface
ghaddab's avatar
ghaddab committed
3
    Copyright(c) 1999 - 2014 Eurecom
4

ghaddab's avatar
ghaddab committed
5
6
7
8
    OpenAirInterface is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
9
10


ghaddab's avatar
ghaddab committed
11
12
13
14
    OpenAirInterface is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.
15

ghaddab's avatar
ghaddab committed
16
    You should have received a copy of the GNU General Public License
17
18
    along with OpenAirInterface.The full GNU General Public License is
   included in this distribution in the file called "COPYING". If not,
ghaddab's avatar
ghaddab committed
19
   see <http://www.gnu.org/licenses/>.
20
21

  Contact Information
ghaddab's avatar
ghaddab committed
22
23
  OpenAirInterface Admin: openair_admin@eurecom.fr
  OpenAirInterface Tech : openair_tech@eurecom.fr
24
  OpenAirInterface Dev  : openair4g-devel@lists.eurecom.fr
25

ghaddab's avatar
ghaddab committed
26
  Address      : Eurecom, Campus SophiaTech, 450 Route des Chappes, CS 50193 - 06904 Biot Sophia Antipolis cedex, FRANCE
27

ghaddab's avatar
ghaddab committed
28
 *******************************************************************************/
29
30
31
32

/*! \file dlsim.c
 \brief Top-level DL simulator
 \author R. Knopp
33
 \date 2011 - 2014
34
35
36
37
38
39
40
 \version 0.1
 \company Eurecom
 \email: knopp@eurecom.fr
 \note
 \warning
*/

41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
#include <string.h>
#include <math.h>
#include <unistd.h>
#include <execinfo.h>
#include <signal.h>

#include "SIMULATION/TOOLS/defs.h"
#include "PHY/types.h"
#include "PHY/defs.h"
#include "PHY/vars.h"
#include "MAC_INTERFACE/vars.h"

#include "SCHED/defs.h"
#include "SCHED/vars.h"
#include "LAYER2/MAC/vars.h"
#include "OCG_vars.h"
#include "UTIL/LOG/log.h"
58
#include "UTIL/LISTS/list.h"
59
60

extern unsigned int dlsch_tbs25[27][25],TBStable[27][110];
61
extern unsigned char offset_mumimo_llr_drange_fix;
Elena Lukashova's avatar
Elena Lukashova committed
62
63
64
extern uint8_t interf_unaw_shift0;
extern uint8_t interf_unaw_shift1;
extern uint8_t interf_unaw_shift;
65
66
67
68
#ifdef XFORMS
#include "PHY/TOOLS/lte_phy_scope.h"
#endif

69
#define PRINT_BYTES
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88

//#define AWGN
//#define NO_DCI



//#define ABSTRACTION

/*
  #define RBmask0 0x00fc00fc
  #define RBmask1 0x0
  #define RBmask2 0x0
  #define RBmask3 0x0
*/

PHY_VARS_eNB *PHY_vars_eNB;
PHY_VARS_UE *PHY_vars_UE;

int otg_enabled=0;
89
90
91
92
93
94
95
/*the following parameters are used to control the processing times calculations*/
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 rx */
int n_tx_dropped = 0; /*!< \brief initial max process time for tx */
int n_rx_dropped = 0; /*!< \brief initial max process time for rx */
96

97
98
void handler(int sig)
{
99
100
101
102
103
104
105
106
107
108
109
110
  void *array[10];
  size_t size;

  // get void*'s for all entries on the stack
  size = backtrace(array, 10);

  // print out all the frames to stderr
  fprintf(stderr, "Error: signal %d:\n", sig);
  backtrace_symbols_fd(array, size, 2);
  exit(1);
}

111
112
113
void lte_param_init(unsigned char N_tx, unsigned char N_rx,unsigned char transmission_mode,uint8_t extended_prefix_flag,uint8_t fdd_flag, uint16_t Nid_cell,uint8_t tdd_config,uint8_t N_RB_DL,
                    uint8_t osf,uint32_t perfect_ce)
{
114
115
116
117
118
119
120
121
122
123
124
125
126
127

  LTE_DL_FRAME_PARMS *lte_frame_parms;
  int i;


  printf("Start lte_param_init\n");
  PHY_vars_eNB = malloc(sizeof(PHY_VARS_eNB));
  PHY_vars_UE = malloc(sizeof(PHY_VARS_UE));
  //PHY_config = malloc(sizeof(PHY_CONFIG));
  mac_xface = malloc(sizeof(MAC_xface));

  srand(0);
  randominit(0);
  set_taus_seed(0);
128

129
130
131
  lte_frame_parms = &(PHY_vars_eNB->lte_frame_parms);

  lte_frame_parms->N_RB_DL            = N_RB_DL;   //50 for 10MHz and 25 for 5 MHz
132
  lte_frame_parms->N_RB_UL            = N_RB_DL;
133
134
135
136
137
138
  lte_frame_parms->Ncp                = extended_prefix_flag;
  lte_frame_parms->Nid_cell           = Nid_cell;
  lte_frame_parms->nushift            = Nid_cell%6;
  lte_frame_parms->nb_antennas_tx     = N_tx;
  lte_frame_parms->nb_antennas_rx     = N_rx;
  lte_frame_parms->nb_antennas_tx_eNB = N_tx;
139
  lte_frame_parms->phich_config_common.phich_resource         = one;
140
141
142
143
144
145
146
147
148
  lte_frame_parms->tdd_config         = tdd_config;
  lte_frame_parms->frame_type         = (fdd_flag==1)?0 : 1;
  //  lte_frame_parms->Csrs = 2;
  //  lte_frame_parms->Bsrs = 0;
  //  lte_frame_parms->kTC = 0;44
  //  lte_frame_parms->n_RRC = 0;
  lte_frame_parms->mode1_flag = (transmission_mode == 1)? 1 : 0;

  init_frame_parms(lte_frame_parms,osf);
149

150
  //copy_lte_parms_to_phy_framing(lte_frame_parms, &(PHY_config->PHY_framing));
151

152
  //  phy_init_top(lte_frame_parms); //allocation
153

Elena Lukashova's avatar
Elena Lukashova committed
154
  PHY_vars_UE->is_secondary_ue = 0;
155
156
157
158
159
160
161
162
163
164
  PHY_vars_UE->lte_frame_parms = *lte_frame_parms;
  PHY_vars_eNB->lte_frame_parms = *lte_frame_parms;

  phy_init_lte_top(lte_frame_parms);
  dump_frame_parms(lte_frame_parms);

  PHY_vars_UE->PHY_measurements.n_adj_cells=0;
  PHY_vars_UE->PHY_measurements.adj_cell_id[0] = Nid_cell+1;
  PHY_vars_UE->PHY_measurements.adj_cell_id[1] = Nid_cell+2;

165
166
  for (i=0; i<3; i++)
    lte_gold(lte_frame_parms,PHY_vars_UE->lte_gold_table[i],Nid_cell+i);
167
168
169
170
171
172

  phy_init_lte_ue(PHY_vars_UE,1,0);
  phy_init_lte_eNB(PHY_vars_eNB,0,0,0);

  generate_pcfich_reg_mapping(&PHY_vars_UE->lte_frame_parms);
  generate_phich_reg_mapping(&PHY_vars_UE->lte_frame_parms);
173

174
  // DL power control init
175
  if (transmission_mode == 1) {
176
177
178
    PHY_vars_eNB->pdsch_config_dedicated->p_a  = dB0; // 4 = 0dB
    ((PHY_vars_eNB->lte_frame_parms).pdsch_config_common).p_b = 0;
    PHY_vars_UE->pdsch_config_dedicated->p_a  = dB0; // 4 = 0dB
179
180
    ((PHY_vars_UE->lte_frame_parms).pdsch_config_common).p_b = 0;
  } else { // rho_a = rhob
181
182
183
184
185
    PHY_vars_eNB->pdsch_config_dedicated->p_a  = dB0; // 4 = 0dB
    ((PHY_vars_eNB->lte_frame_parms).pdsch_config_common).p_b = 1;
    PHY_vars_UE->pdsch_config_dedicated->p_a  = dB0; // 4 = 0dB
    ((PHY_vars_UE->lte_frame_parms).pdsch_config_common).p_b = 1;
  }
186

knopp's avatar
   
knopp committed
187
  PHY_vars_UE->perfect_ce = perfect_ce;
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203

  printf("Done lte_param_init\n");


}


//DCI2_5MHz_2A_M10PRB_TDD_t DLSCH_alloc_pdu2_2A[2];

DCI1E_5MHz_2A_M10PRB_TDD_t  DLSCH_alloc_pdu2_1E[2];
uint64_t DLSCH_alloc_pdu_1[2];

#define UL_RB_ALLOC 0x1ff;
#define CCCH_RB_ALLOC computeRIV(PHY_vars_eNB->lte_frame_parms.N_RB_UL,0,2)
//#define DLSCH_RB_ALLOC 0x1fbf // igore DC component,RB13
//#define DLSCH_RB_ALLOC 0x0001
204
205
void do_OFDM_mod_l(mod_sym_t **txdataF, int32_t **txdata, uint16_t next_slot, LTE_DL_FRAME_PARMS *frame_parms)
{
206
207
208
209
210
211
212
213
214
215
216

  int aa, slot_offset, slot_offset_F;

  slot_offset_F = (next_slot)*(frame_parms->ofdm_symbol_size)*((frame_parms->Ncp==1) ? 6 : 7);
  slot_offset = (next_slot)*(frame_parms->samples_per_tti>>1);

  for (aa=0; aa<frame_parms->nb_antennas_tx; aa++) {
    //    printf("Thread %d starting ... aa %d (%llu)\n",omp_get_thread_num(),aa,rdtsc());

    if (frame_parms->Ncp == 1)
      PHY_ofdm_mod(&txdataF[aa][slot_offset_F],        // input
217
218
219
220
221
                   &txdata[aa][slot_offset],         // output
                   frame_parms->log2_symbol_size,                // log2_fft_size
                   6,                 // number of symbols
                   frame_parms->nb_prefix_samples,               // number of prefix samples
                   CYCLIC_PREFIX);
222
223
    else {
      normal_prefix_mod(&txdataF[aa][slot_offset_F],
224
225
226
                        &txdata[aa][slot_offset],
                        7,
                        frame_parms);
227
    }
228

229
230
231
232
233

  }

}

234
235
int main(int argc, char **argv)
{
236

237
  int c;
238
239
240
241
  int k,i,aa,aarx,aatx;

  int s,Kr,Kr_bytes;

242
  double sigma2, sigma2_dB=10,SNR,snr0=-2.0,snr1,rate;
243
244
245
246
247
248
249
  double snr_step=1,input_snr_step=1, snr_int=30;

  LTE_DL_FRAME_PARMS *frame_parms;
  double **s_re,**s_im,**r_re,**r_im;
  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;

250
  uint8_t extended_prefix_flag=0,transmission_mode=1,n_tx=1,n_rx=1;
gauthier's avatar
gauthier committed
251
  uint16_t Nid_cell=0;
252

253
254
  int8_t eNB_id = 0, eNB_id_i = 1;
  unsigned char mcs1=0,mcs2=0,mcs_i=0,awgn_flag=0,round,dci_flag=0;
255
256
257
  unsigned char i_mod = 2;
  unsigned short NB_RB;
  unsigned char Ns,l,m;
gauthier's avatar
gauthier committed
258
259
  uint16_t tdd_config=3;
  uint16_t n_rnti=0x1234;
260
  int n_users = 1;
261
262
  RX_type_t rx_type=rx_standard;

263
264
265
266

  SCM_t channel_model=Rayleigh1;
  //  unsigned char *input_data,*decoded_output;

267
268
  unsigned char *input_buffer0[2],*input_buffer1[2];
  unsigned short input_buffer_length0,input_buffer_length1;
269
270
  unsigned int ret;
  unsigned int coded_bits_per_codeword=0,nsymb,dci_cnt,tbs=0;
271

272
 
273
  unsigned int tx_lev=0,tx_lev_dB=0,trials,errs[2][4],round_trials[4]={0,0,0,0},dci_errors=0,dlsch_active=0,num_layers;
274
  int re_allocated;
knopp's avatar
knopp committed
275
  char fname[32],vname[32];
276
277
  FILE *bler_fd;
  char bler_fname[256];
278
279
  FILE *time_meas_fd;
  char time_meas_fname[256];
280
281
  //FILE *tikz_fd;
  //char tikz_fname[256];
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301

  FILE *input_trch_fd=NULL;
  unsigned char input_trch_file=0;
  FILE *input_fd=NULL;
  unsigned char input_file=0;
  //  char input_val_str[50],input_val_str2[50];

  char input_trch_val[16];
  double channelx,channely;

  //  unsigned char pbch_pdu[6];

  DCI_ALLOC_t dci_alloc[8],dci_alloc_rx[8];
  int num_common_dci=0,num_ue_spec_dci=0,num_dci=0;

  //  FILE *rx_frame_file;

  int n_frames;
  int n_ch_rlz = 1;
  channel_desc_t *eNB2UE[4];
gauthier's avatar
gauthier committed
302
303
304
  uint8_t num_pdcch_symbols=1,num_pdcch_symbols_2=0;
  uint8_t pilot1,pilot2,pilot3;
  uint8_t rx_sample_offset = 0;
305
306
  //char stats_buffer[4096];
  //int len;
307
  uint8_t num_rounds = 4;
gauthier's avatar
gauthier committed
308
  uint8_t subframe=7;
309
310
311
312
313
314
315
316
317
  int u;
  int n=0;
  int abstx=0;
  int iii;
  FILE *csv_fd=NULL;
  char csv_fname[512];
  int ch_realization;
  int pmi_feedback=0;
  int hold_channel=0;
318

319
320
321
  // void *data;
  // int ii;
  int bler;
322
  double blerr[2][4],uncoded_ber,avg_ber;
323
  short *uncoded_ber_bit=NULL;
gauthier's avatar
gauthier committed
324
325
  uint8_t N_RB_DL=25,osf=1;
  uint8_t fdd_flag = 0;
Elena Lukashova's avatar
Elena Lukashova committed
326
  int xforms = 0;
327
328
329
330
#ifdef XFORMS
  FD_lte_phy_scope_ue *form_ue;
  char title[255];
#endif
gauthier's avatar
gauthier committed
331
  uint32_t DLSCH_RB_ALLOC = 0x1fff;
332
333
334
335
336
  int numCCE=0;
  int dci_length_bytes=0,dci_length=0;
  int common_flag=0,TPC=0;

  double cpu_freq_GHz;
337
  time_stats_t ts;//,sts,usts;
338
339
340
  int avg_iter,iter_trials;
  int rballocset=0;
  int print_perf=0;
341
  int test_perf=0;
nikaeinn's avatar
nikaeinn committed
342
  int dump_table=0;
343
344
  int llr8_flag=0;

345
  double effective_rate=0.0;
346
  char channel_model_input[17]="I";
347
348

  int TB0_active = 1;
knopp's avatar
   
knopp committed
349
  uint32_t perfect_ce = 0;
350

knopp's avatar
knopp committed
351
352
  LTE_DL_UE_HARQ_t *dlsch0_ue_harq;
  LTE_DL_eNB_HARQ_t *dlsch0_eNB_harq;
353
  uint8_t Kmimo;
354
  mod_sym_t **sic_buffer;
knopp's avatar
knopp committed
355

356
357
  FILE    *proc_fd = NULL;
  char buf[64];
knopp's avatar
   
knopp committed
358

359

knopp's avatar
   
knopp committed
360
361
  opp_enabled=1; // to enable the time meas

362
363
364
365
366
367
368
369
370
371
372
373
374
#if defined(__arm__)
  proc_fd = fopen("/sys/devices/system/cpu/cpu4/cpufreq/cpuinfo_cur_freq", "r");
  if(!proc_fd)
     printf("cannot open /sys/devices/system/cpu/cpu4/cpufreq/cpuinfo_cur_freq");
  else {
     while(fgets(buf, 63, proc_fd))
        printf("%s", buf);
  }
  fclose(proc_fd);
  cpu_freq_GHz = ((double)atof(buf))/1e6;
#else
  cpu_freq_GHz = get_cpu_freq_GHz();
#endif
375
376
  printf("Detected cpu_freq %f GHz\n",cpu_freq_GHz);

377
378
  //signal(SIGSEGV, handler);
  //signal(SIGABRT, handler);
379
380

  logInit();
381
  set_glog(LOG_INFO, LOG_MED);
382
383
384
385
386

  // default parameters
  n_frames = 1000;
  snr0 = 0;
  num_layers = 1;
knopp's avatar
   
knopp committed
387
  perfect_ce = 0;
388

389
  while ((c = getopt (argc, argv, "ahdpZDe:m:n:o:s:f:t:c:g:r:F:x:y:z:AM:N:I:i:O:R:S:C:T:b:u:v:w:B:PLl:YXV:W:J:")) != -1) {
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
    switch (c) {
    case 'a':
      awgn_flag = 1;
      channel_model = AWGN;
      break;

    case 'b':
      tdd_config=atoi(optarg);
      break;

    case 'B':
      N_RB_DL=atoi(optarg);
      break;

    case 'd':
      dci_flag = 1;
      break;

    case 'm':
      mcs1 = atoi(optarg);
      break;

    case 'M':
      mcs2 = atoi(optarg);
414
      //i_mod = get_Qm(mcs2); /// think here again!!!
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
      break;

    case 't':
      mcs_i = atoi(optarg);
      i_mod = get_Qm(mcs_i);
      break;

    case 'n':
      n_frames = atoi(optarg);
      break;

    case 'C':
      Nid_cell = atoi(optarg);
      break;

    case 'o':
      rx_sample_offset = atoi(optarg);
      break;

    case 'D':
      fdd_flag = 1;
      break;

    case 'r':
      DLSCH_RB_ALLOC = atoi(optarg);
      rballocset = 1;
      break;

    case 'F':
      forgetting_factor = atof(optarg);
      break;

    case 's':
      snr0 = atof(optarg);
      break;

    case 'w':
      snr_int = atof(optarg);
      break;
Elena Lukashova's avatar
Elena Lukashova committed
454
      
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
    case 'f':
      input_snr_step= atof(optarg);
      break;

    case 'A':
      abstx = 1;
      break;

    case 'N':
      n_ch_rlz= atof(optarg);
      break;

    case 'p':
      extended_prefix_flag=1;
      break;

    case 'c':
      num_pdcch_symbols=atoi(optarg);
      break;

    case 'g':
Elena Lukashova's avatar
Elena Lukashova committed
476
      memcpy(channel_model_input,optarg,17);
477
478
479
480
481

      switch((char)*optarg) {
      case 'A':
        channel_model=SCM_A;
        break;
482
      case 'B':
483
484
485
        channel_model=SCM_B;
        break;

486
      case 'C':
487
488
489
        channel_model=SCM_C;
        break;

490
      case 'D':
491
492
493
494
495
496
497
        channel_model=SCM_D;
        break;

      case 'E':
        channel_model=EPA;
        break;

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
	case 'F': 
	  channel_model=EVA;
	  break;
	case 'G': 
	  channel_model=ETU;
	  break;
	case 'H':
	  channel_model=Rayleigh8;
	  break;
	case 'I':
	  channel_model=Rayleigh1;
	  break;
	case 'J':
	  channel_model=Rayleigh1_corr;
	  break;
	case 'K':
	  channel_model=Rayleigh1_anticorr;
	  break;
	case 'L':
	  channel_model=Rice8;
	  break;
	case 'M':
	  channel_model=Rice1;
	  break;
	case 'N':
	  channel_model=AWGN;
	  break;
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
	case 'P':
	  channel_model=Rayleigh1_orthogonal;
	  break;
	case 'Q':
	  channel_model=Rayleigh1_orth_eff_ch_TM4_prec_real; // for DUALSTREAM_UNIFORM_PRECODING1 when interf is precancelled
	break;
	case 'R':
	  channel_model=Rayleigh1_orth_eff_ch_TM4_prec_imag; // for DUALSTREAM_UNIFORM_PRECODINGj when interf is precancelled
	  break;
	case 'S':
	  channel_model=Rayleigh8_orth_eff_ch_TM4_prec_real;//when interf is precancelled
	  break;
	case 'T':
	  channel_model=Rayleigh8_orth_eff_ch_TM4_prec_imag;//when interf is precancelled
	  break;
	case 'U':
	  channel_model = TS_SHIFT;
	  break;
543
544
545
546
547
548
549
550
551
	default:
	  msg("Unsupported channel model!\n");
	  exit(-1);
	}
	break;
      case 'x':
	transmission_mode=atoi(optarg);
	if ((transmission_mode!=1) &&
	    (transmission_mode!=2) &&
552
	    (transmission_mode!=3) &&
553
	    (transmission_mode!=4) &&
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
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
	    (transmission_mode!=5) &&
	    (transmission_mode!=6)) {
	  msg("Unsupported transmission mode %d\n",transmission_mode);
	  exit(-1);
	}
	if (transmission_mode>1) {
	  n_tx = 2;
	}
	break;
      case 'y':
	n_tx=atoi(optarg);
	if ((n_tx==0) || (n_tx>2)) {
	  msg("Unsupported number of tx antennas %d\n",n_tx);
	  exit(-1);
	}
	break;
      case 'z':
	n_rx=atoi(optarg);
	if ((n_rx==0) || (n_rx>2)) {
	  msg("Unsupported number of rx antennas %d\n",n_rx);
	  exit(-1);
	}
	break;
      case 'I':
	input_trch_fd = fopen(optarg,"r");
	input_trch_file=1;
	break;
      case 'i':
	input_fd = fopen(optarg,"r");
	input_file=1;
	dci_flag = 1;	
	break;
      case 'e':
	num_rounds=1;
	common_flag = 1;
	TPC = atoi(optarg);
	break;
      case 'R':
	num_rounds=atoi(optarg);
	break;
      case 'S':
	subframe=atoi(optarg);
	break;
      case 'T':
	n_rnti=atoi(optarg);
	break;	
      case 'u':
601
602
603
	rx_type = (RX_type_t) atoi(optarg);
	if (rx_type<rx_standard || rx_type>rx_IC_dual_stream) {
	  printf("Unsupported rx type %d\n",rx_type);
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
	  exit(-1);
	}
	break;
      case 'v':
	i_mod = atoi(optarg);
	if (i_mod!=2 && i_mod!=4 && i_mod!=6) {
	  msg("Wrong i_mod %d, should be 2,4 or 6\n",i_mod);
	  exit(-1);
	}
	break;
      case 'P':
	print_perf=1;
	break;
      case 'L':
	llr8_flag=1;
	break;
620
621
622
      case 'l':
	offset_mumimo_llr_drange_fix=atoi(optarg);
	break;
623
624
625
626
      case 'O':
	test_perf=atoi(optarg);
	//print_perf =1;
	break;
627
628
629
      case 'X':
	xforms = 1;
	break;
630
      case 'Z':
nikaeinn's avatar
nikaeinn committed
631
632
	dump_table=1;
	break;
knopp's avatar
   
knopp committed
633
634
635
      case 'Y':
	perfect_ce=1;
	break;
636
637
638
639
640
641
      case 'V':
	interf_unaw_shift0=atoi(optarg);
	break;	 
      case 'W':
	interf_unaw_shift1=atoi(optarg);
	break;
642
643
644
      case 'J':
	interf_unaw_shift=atoi(optarg);
	break;
645
646
      case 'h':
      default:
647
	printf("%s -h(elp) -a(wgn on) -d(ci decoding on) -p(extended prefix on) -m mcs1 -M mcs2 -n n_frames -s snr0 -x transmission mode (1,2,3,5,6) -y TXant -z RXant -I trch_file\n",argv[0]);
648
649
650
	printf("-h This message\n");
	printf("-a Use AWGN channel and not multipath\n");
	printf("-c Number of PDCCH symbols\n");
651
652
	printf("-m MCS1 for TB 1\n");
	printf("-M MCS2 for TB 2\n");
653
654
655
656
657
658
659
	printf("-d Transmit the DCI and compute its error statistics and the overall throughput\n");
	printf("-p Use extended prefix mode\n");
	printf("-n Number of frames to simulate\n");
	printf("-o Sample offset for receiver\n");
	printf("-s Starting SNR, runs from SNR to SNR+%.1fdB in steps of %.1fdB. If n_frames is 1 then just SNR is simulated and MATLAB/OCTAVE output is generated\n", snr_int, snr_step);
	printf("-f step size of SNR, default value is 1.\n");
	printf("-r ressource block allocation (see  section 7.1.6.3 in 36.213\n");
660
	printf("-g Channel model, possible values are 3GPP 25.814 SCM-A/B/C/D('A','B','C','D'), 36-101 EPA('E'), EVA ('F'),ETU('G'), Rayghleigh8 ('H'), Rayleigh1('I'), Rayleigh1_corr('J'), Rayleigh1_anticorr('K'), Rice8('L'), Rice1('M'), AWGN('N'), Rayleigh1_orthogonal('P'), Rayleigh1_orth_eff_ch_TM4_prec_real ('Q'), Rayleigh1_orth_eff_ch_TM4_prec_imag ('R'), Rayleigh8_orth_eff_ch_TM4_prec_real ('S'),Rayleigh8_orth_eff_ch_TM4_prec_imag ('T')   \n");
661
662
663
664
665
666
	printf("-F forgetting factor (0 new channel every trial, 1 channel constant\n");
	printf("-x Transmission mode (1,2,6 for the moment)\n");
	printf("-y Number of TX antennas used in eNB\n");
	printf("-z Number of RX antennas used in UE\n");
	printf("-t MCS of interfering UE\n");
	printf("-R Number of HARQ rounds (fixed)\n");
667
	printf("-A Turns on calibration mode for abstraction.\n");
Florian Kaltenberger's avatar
   
Florian Kaltenberger committed
668
	printf("-N Determines the number of Channel Realizations in Abstraction mode. Default value is 1. \n");
669
	printf("-O Set the percenatge of effective rate to testbench the modem performance (typically 30 and 70, range 1-100) \n");
670
	printf("-I Input filename for TrCH data (binary)\n");
671
	printf("-u Receiver type: 0=standard, 1 = single stream IC (for TM3,4,5,6), 2 = dual stream IC (for TM3,4), 3 = SIC (for TM3,4) \n");
672
673
	exit(1);
	break;
Elena Lukashova's avatar
Elena Lukashova committed
674
    
675
    }
Elena Lukashova's avatar
Elena Lukashova committed
676
    
677
  }
Elena Lukashova's avatar
Elena Lukashova committed
678
    
679
  if (common_flag == 0) {
680
681
682
683
684
    switch (N_RB_DL) {
    case 6:
      if (rballocset==0) DLSCH_RB_ALLOC = 0x3f;
      num_pdcch_symbols = 3;
      break;
685

686
687
688
    case 25:
      if (rballocset==0) DLSCH_RB_ALLOC = 0x1fff;
      break;
689

690
691
692
    case 50:
      if (rballocset==0) DLSCH_RB_ALLOC = 0x1ffff;
      break;
693

694
695
696
697
    case 100:
      if (rballocset==0) DLSCH_RB_ALLOC = 0x1ffffff;
      break;
    }
698

699
    NB_RB=conv_nprb(0,DLSCH_RB_ALLOC,N_RB_DL);
700
  } else
701
702
    NB_RB = 4;

703
  NB_RB=conv_nprb(0,DLSCH_RB_ALLOC,N_RB_DL); 
704

705
  if ((transmission_mode > 1) && (n_tx != 2)) {
706

707
    printf("n_tx must be >1 for transmission_mode %d\n",transmission_mode);
708
709
710
711
712
713
714
715
716
717
718
719
    exit(-1);
  }

  if (((transmission_mode==1) || (transmission_mode==2)) && (rx_type != rx_standard)) {
    printf("only standard rx available for TM1 and TM2\n");
    exit(-1);
  }
  if (((transmission_mode==5) || (transmission_mode==6)) && (rx_type == rx_IC_dual_stream)) {
    printf("only standard rx or single stream IC available for TM5 and TM6\n");
    exit(-1);
  }

720

721
#ifdef XFORMS
722
  if (xforms==1) {
723
724
725
  fl_initialize (&argc, argv, NULL, 0, 0);
  form_ue = create_lte_phy_scope_ue();
  sprintf (title, "LTE PHY SCOPE eNB");
726
  fl_show_form (form_ue->lte_phy_scope_ue, FL_PLACE_HOTSPOT, FL_FULLBORDER, title);
727

728
729
  /*
  if (rx_type==rx_IC_single_stream) {
730
731
732
733
734
    openair_daq_vars.use_ia_receiver = 1;
    fl_set_button(form_ue->button_0,1);
    fl_set_object_label(form_ue->button_0, "IA Receiver ON");
    fl_set_object_color(form_ue->button_0, FL_GREEN, FL_GREEN);
  }
735

736
737
  */
  }
Elena Lukashova's avatar
Elena Lukashova committed
738

739
740
#endif

741
  if (transmission_mode==5) { 
742
    n_users = 2;
Elena Lukashova's avatar
Elena Lukashova committed
743
    //eNB_id_i = PHY_vars_UE->n_connected_eNB;
Florian Kaltenberger's avatar
Florian Kaltenberger committed
744
    eNB_id_i=1;
Elena Lukashova's avatar
Elena Lukashova committed
745
      
746
  }
747
748
749
  else 
    eNB_id_i = eNB_id;

750

Elena Lukashova's avatar
Elena Lukashova committed
751
  lte_param_init(n_tx,n_rx,transmission_mode,extended_prefix_flag,fdd_flag,Nid_cell,tdd_config,N_RB_DL,osf,perfect_ce);
752

753
754
  printf("Setting mcs1 = %d\n",mcs1);
  printf("Setting mcs2 = %d\n",mcs2);
755
756
757
  printf("NPRB = %d\n",NB_RB);
  printf("n_frames = %d\n",n_frames);
  printf("Transmission mode %d with %dx%d antenna configuration, Extended Prefix %d\n",transmission_mode,n_tx,n_rx,extended_prefix_flag);
758
  printf("Using receiver type %d\n", rx_type);
Elena Lukashova's avatar
Elena Lukashova committed
759
760
761
  printf("TM1 shift %d\n", interf_unaw_shift);
  //printf("Using I_UA rec shift layer 1  %d\n", interf_unaw_shift0);
  //printf("Using I_UA rec shift layer 2  %d\n", interf_unaw_shift1);
762
763
764
765
766
767
768
  snr1 = snr0+snr_int;
  printf("SNR0 %f, SNR1 %f\n",snr0,snr1);

  /*
    txdataF    = (int **)malloc16(2*sizeof(int*));
    txdataF[0] = (int *)malloc16(FRAME_LENGTH_BYTES);
    txdataF[1] = (int *)malloc16(FRAME_LENGTH_BYTES);
769

770
771
772
773
774
775
776
    txdata    = (int **)malloc16(2*sizeof(int*));
    txdata[0] = (int *)malloc16(FRAME_LENGTH_BYTES);
    txdata[1] = (int *)malloc16(FRAME_LENGTH_BYTES);
  */

  frame_parms = &PHY_vars_eNB->lte_frame_parms;

777
  s_re = malloc(2*sizeof(double*)); //transmitted signal (Sent)
778
  s_im = malloc(2*sizeof(double*));
779
  r_re = malloc(2*sizeof(double*)); //received signal
780
781
782
783
784
785
  r_im = malloc(2*sizeof(double*));
  //  r_re0 = malloc(2*sizeof(double*));
  //  r_im0 = malloc(2*sizeof(double*));

  nsymb = (PHY_vars_eNB->lte_frame_parms.Ncp == 0) ? 14 : 12;

786
  printf("Channel Model= (%s,%d)\n",channel_model_input, channel_model);
787
788
  printf("SCM-A=%d, SCM-B=%d, SCM-C=%d, SCM-D=%d, EPA=%d, EVA=%d, ETU=%d, Rayleigh8=%d, Rayleigh1=%d, Rayleigh1_corr=%d, Rayleigh1_anticorr=%d, Rice1=%d, Rice8=%d, Rayleigh1_orthogonal=%d, Rayleigh1_orth_eff_ch_TM4_prec_real22=%d, Rayleigh1_orth_eff_ch_TM4_prec_imag=%d, Rayleigh8_orth_eff_ch_TM4_prec_real=%d,  Rayleigh8_orth_eff_ch_TM4_prec_imag=%d , TS_SHIFT=%d\n",
	 SCM_A, SCM_B, SCM_C, SCM_D, EPA, EVA, ETU, Rayleigh8, Rayleigh1, Rayleigh1_corr, Rayleigh1_anticorr, Rice1, Rice8, Rayleigh1_orthogonal, Rayleigh1_orth_eff_ch_TM4_prec_real, Rayleigh1_orth_eff_ch_TM4_prec_imag, Rayleigh8_orth_eff_ch_TM4_prec_real, Rayleigh8_orth_eff_ch_TM4_prec_imag, TS_SHIFT);
789
  if(transmission_mode==5)
790
791
792
    sprintf(bler_fname,"bler_tx%d_rec%d_chan%d_nrx%d_mcs%d_mcsi%d_u%d_imod%d.csv",transmission_mode,rx_type,channel_model,n_rx,mcs1,mcs_i,rx_type,i_mod);
  else if (abstx == 1) 
    if (perfect_ce==1)
793
	 sprintf(bler_fname,"bler_tx%d_r%d_ch%d_%d_nrx%d_mcs%d_mcsi%d_ab_pce_sh%d_maxamant.csv",transmission_mode,rx_type,channel_model,n_frames, n_rx,mcs1, mcs2,interf_unaw_shift );
794
  	else
795
    	sprintf(bler_fname,"bler_tx%d_r%d_ch%d_%d_nrx%d_mcs%d_mcsi%d_ab_sh%d_d2_maxamant.csv",transmission_mode,rx_type,channel_model, n_frames, n_rx,mcs1, mcs2,interf_unaw_shift );
Elena Lukashova's avatar
Elena Lukashova committed
796
  else //abstx=0
797
    if (perfect_ce==1)
798
	sprintf(bler_fname,"bler_tx%d_r%d_ch%d_%d_nrx%d_mcs%d_mcsi%d_pce_sh%d_d2_maxamant.csv",transmission_mode,rx_type,channel_model,n_frames, n_rx,mcs1, mcs2, interf_unaw_shift);
Elena Lukashova's avatar
Elena Lukashova committed
799
   else
800
        sprintf(bler_fname,"bler_tx%d_r%d_ch%d_%d_nrx%d_mcs%d_mcsi%d_sh%d_d2_maxamant.csv",transmission_mode,rx_type,channel_model,n_frames,n_rx,mcs1, mcs2, interf_unaw_shift);
Elena Lukashova's avatar
Elena Lukashova committed
801
  
802
  bler_fd = fopen(bler_fname,"w");
803
804
805
806
  if (bler_fd==NULL) {
    fprintf(stderr,"Cannot create file %s!\n",bler_fname);
    exit(-1);
  }
807

808
  fprintf(bler_fd,"SNR; MCS1; MCS2; TBS1; TBS2; rate; err0_st1; err0_st2 trials0; err1_st1; err1_st2; trials1; err2_st1; err2_st2; trials2; err3_st1; err3_st2; trials3; dci_err\n");
Elena Lukashova's avatar
Elena Lukashova committed
809
    if (test_perf != 0) {
810
    
811
    char hostname[1024];
Elena Lukashova's avatar
Elena Lukashova committed
812

813
814
815
    hostname[1023] = '\0';
    gethostname(hostname, 1023);
    printf("Hostname: %s\n", hostname);
816
817
818
819
820
    //char dirname[FILENAME_MAX];
    //sprintf(dirname, "%s/SIMU/USER/pre-ci-logs-%s", getenv("OPENAIR_TARGETS"),hostname );
    sprintf(time_meas_fname,"time_meas_prb%d_mcs%d_anttx%d_antrx%d_pdcch%d_channel%s_tx%d.csv",
            N_RB_DL,mcs1,n_tx,n_rx,num_pdcch_symbols,channel_model_input,transmission_mode);
    //mkdir(dirname,0777);
821
    time_meas_fd = fopen(time_meas_fname,"w");
822
823
824
825
    if (time_meas_fd==NULL) {
      fprintf(stderr,"Cannot create file %s!\n",time_meas_fname);
      exit(-1);
    }
826
  }
827
828
  
  if(abstx){
829
830
831
832
833
834
835
    // CSV file  // add here second stream mcs
    
    if (transmission_mode == 5)
      sprintf(csv_fname,"dataout_tx%d_u2%d_mcs%d_chan%d_nsimus%d_R%d_abstr.m",transmission_mode,rx_type,mcs1,channel_model,n_frames,num_rounds);
    
    else 
      if (perfect_ce==1)
836
		sprintf(csv_fname,"dout_tx%d_r%d_mcs%d_mcsi%d_ch%d_ns%d_R%d_ab_pce_sh%d_d2_%d_maxamant.m",transmission_mode,rx_type,mcs1,mcs2,channel_model,n_frames,num_rounds, interf_unaw_shift, n_ch_rlz);
837
 	else
838
		sprintf(csv_fname,"dout_tx%d_r%d_mcs%d_mcsi%d_ch%d_ns%d_R%d_ab_sh%d_d2_%d_maxamant.m",transmission_mode,rx_type,mcs1,mcs2,channel_model,n_frames,num_rounds, interf_unaw_shift, n_ch_rlz);
839

840
    csv_fd = fopen(csv_fname,"w");
841
    fprintf(csv_fd,"data_all%d=[",mcs1);
842
843
844
845
    if (csv_fd==NULL) {
      fprintf(stderr,"Cannot create file %s!\n",csv_fname);
      exit(-1);
    }
846
847
  }

848
  /*
849
850
  //sprintf(tikz_fname, "second_bler_tx%d_u2=%d_mcs%d_chan%d_nsimus%d.tex",transmission_mode,dual_stream_UE,mcs,channel_model,n_frames);
  sprintf(tikz_fname, "second_bler_tx%d_u2%d_mcs%d_chan%d_nsimus%d",transmission_mode,dual_stream_UE,mcs,channel_model,n_frames);
851
852
  tikz_fd = fopesprintf(csv_fname,"dataout_tx%d_mcs%d_mcs_interf%d_chan%d_nsimus%d_R%d_abstr_old_perf_ce_llr.m",transmission_mode,mcs1,mcs2,channel_model,n_frames,num_rounds);
n(tikz_fname,"w");
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
  //fprintf(tikz_fd,"\\addplot[color=red, mark=o] plot coordinates {");
  switch (mcs)
    {
    case 0:
      fprintf(tikz_fd,"\\addplot[color=blue, mark=star] plot coordinates {");
      break;
    case 1:
      fprintf(tikz_fd,"\\addplot[color=red, mark=star] plot coordinates {");
      break;
    case 2:
      fprintf(tikz_fd,"\\addplot[color=green, mark=star] plot coordinates {");
      break;
    case 3:
      fprintf(tikz_fd,"\\addplot[color=yellow, mark=star] plot coordinates {");
      break;
    case 4:
      fprintf(tikz_fd,"\\addplot[color=black, mark=star] plot coordinates {");
      break;
    case 5:
      fprintf(tikz_fd,"\\addplot[color=blue, mark=o] plot coordinates {");
      break;
    case 6:
      fprintf(tikz_fd,"\\addplot[color=red, mark=o] plot coordinates {");
      break;
    case 7:
      fprintf(tikz_fd,"\\addplot[color=green, mark=o] plot coordinates {");
      break;
    case 8:
      fprintf(tikz_fd,"\\addplot[color=yellow, mark=o] plot coordinates {");
      break;
    case 9:
      fprintf(tikz_fd,"\\addplot[color=black, mark=o] plot coordinates {");
      break;
    case 10:
      fprintf(tikz_fd,"\\addplot[color=blue, mark=square] plot coordinates {");
      break;
    case 11:
      fprintf(tikz_fd,"\\addplot[color=red, mark=square] plot coordinates {");
      break;
    case 12:
      fprintf(tikz_fd,"\\addplot[color=green, mark=square] plot coordinates {");
      break;
    case 13:
      fprintf(tikz_fd,"\\addplot[color=yellow, mark=square] plot coordinates {");
      break;
    case 14:
      fprintf(tikz_fd,"\\addplot[color=black, mark=square] plot coordinates {");
      break;
    case 15:
      fprintf(tikz_fd,"\\addplot[color=blue, mark=diamond] plot coordinates {");
      break;
    case 16:
      fprintf(tikz_fd,"\\addplot[color=red, mark=diamond] plot coordinates {");
      break;
    case 17:
      fprintf(tikz_fd,"\\addplot[color=green, mark=diamond] plot coordinates {");
      break;
    case 18:
      fprintf(tikz_fd,"\\addplot[color=yellow, mark=diamond] plot coordinates {");
      break;
    case 19:
      fprintf(tikz_fd,"\\addplot[color=black, mark=diamond] plot coordinates {");
      break;
    case 20:
      fprintf(tikz_fd,"\\addplot[color=blue, mark=x] plot coordinates {");
      break;
    case 21:
      fprintf(tikz_fd,"\\addplot[color=red, mark=x] plot coordinates {");
      break;
    case 22:
      fprintf(tikz_fd,"\\addplot[color=green, mark=x] plot coordinates {");
      break;
    case 23:
      fprintf(tikz_fd,"\\addplot[color=yellow, mark=x] plot coordinates {");
      break;
    case 24:
      fprintf(tikz_fd,"\\addplot[color=black, mark=x] plot coordinates {");
      break;
    case 25:
      fprintf(tikz_fd,"\\addplot[color=blue, mark=x] plot coordinates {");
      break;
    case 26:
      fprintf(tikz_fd,"\\addplot[color=red, mark=+] plot coordinates {");
      break;
    case 27:
      fprintf(tikz_fd,"\\addplot[color=green, mark=+] plot coordinates {");
      break;
    case 28:
      fprintf(tikz_fd,"\\addplot[color=yellow, mark=+] plot coordinates {");
      break;
    }
944
945
  */

946
  for (i=0; i<2; i++) {
947
948
949
950
951
952
953
954
955
956
957
958
    s_re[i] = malloc(FRAME_LENGTH_COMPLEX_SAMPLES*sizeof(double));
    s_im[i] = malloc(FRAME_LENGTH_COMPLEX_SAMPLES*sizeof(double));
    r_re[i] = malloc(FRAME_LENGTH_COMPLEX_SAMPLES*sizeof(double));
    r_im[i] = malloc(FRAME_LENGTH_COMPLEX_SAMPLES*sizeof(double));
    //    r_re0[i] = malloc(FRAME_LENGTH_COMPLEX_SAMPLES*sizeof(double));
    //    bzero(r_re0[i],FRAME_LENGTH_COMPLEX_SAMPLES*sizeof(double));
    //    r_im0[i] = malloc(FRAME_LENGTH_COMPLEX_SAMPLES*sizeof(double));
    //    bzero(r_im0[i],FRAME_LENGTH_COMPLEX_SAMPLES*sizeof(double));
  }


  PHY_vars_UE->lte_ue_pdcch_vars[0]->crnti = n_rnti;
Elena Lukashova's avatar
Elena Lukashova committed
959
960
961
962
963
  PHY_vars_UE->transmission_mode[eNB_id]=transmission_mode;
  if (PHY_vars_UE->transmission_mode[eNB_id] !=4)
      PHY_vars_UE->PHY_measurements.rank[eNB_id]=0;
  else 
      PHY_vars_UE->PHY_measurements.rank[eNB_id]=1;
964
965
966
967
968
969
970
971
  // Fill in UL_alloc
  UL_alloc_pdu.type    = 0;
  UL_alloc_pdu.hopping = 0;
  UL_alloc_pdu.rballoc = UL_RB_ALLOC;
  UL_alloc_pdu.mcs     = 1;
  UL_alloc_pdu.ndi     = 1;
  UL_alloc_pdu.TPC     = 0;
  UL_alloc_pdu.cqi_req = 1;
972
  
973
974
975
976
977
978
979
980
981
982
983
984
985
  CCCH_alloc_pdu.type               = 0;
  CCCH_alloc_pdu.vrb_type           = 0;
  CCCH_alloc_pdu.rballoc            = CCCH_RB_ALLOC;
  CCCH_alloc_pdu.ndi      = 1;
  CCCH_alloc_pdu.mcs      = 1;
  CCCH_alloc_pdu.harq_pid = 0;

  DLSCH_alloc_pdu2_1E[0].rah              = 0;
  DLSCH_alloc_pdu2_1E[0].rballoc          = DLSCH_RB_ALLOC;
  DLSCH_alloc_pdu2_1E[0].TPC              = 0;
  DLSCH_alloc_pdu2_1E[0].dai              = 0;
  DLSCH_alloc_pdu2_1E[0].harq_pid         = 0;
  //DLSCH_alloc_pdu2_1E[0].tb_swap          = 0;
986
  DLSCH_alloc_pdu2_1E[0].mcs             = mcs1;
987
988
989
990
991
992
993
994
995
996
997
998
  DLSCH_alloc_pdu2_1E[0].ndi             = 1;
  DLSCH_alloc_pdu2_1E[0].rv              = 0;
  // Forget second codeword
  DLSCH_alloc_pdu2_1E[0].tpmi             = (transmission_mode>=5 ? 5 : 0);  // precoding
  DLSCH_alloc_pdu2_1E[0].dl_power_off     = (transmission_mode==5 ? 0 : 1);

  DLSCH_alloc_pdu2_1E[1].rah              = 0;
  DLSCH_alloc_pdu2_1E[1].rballoc          = DLSCH_RB_ALLOC;
  DLSCH_alloc_pdu2_1E[1].TPC              = 0;
  DLSCH_alloc_pdu2_1E[1].dai              = 0;
  DLSCH_alloc_pdu2_1E[1].harq_pid         = 0;
  //DLSCH_alloc_pdu2_1E[1].tb_swap          = 0;
999
  DLSCH_alloc_pdu2_1E[1].mcs             = mcs_i;
1000
1001
1002
1003
1004
1005
1006
  DLSCH_alloc_pdu2_1E[1].ndi             = 1;
  DLSCH_alloc_pdu2_1E[1].rv              = 0;
  // Forget second codeword
  DLSCH_alloc_pdu2_1E[1].tpmi             = (transmission_mode>=5 ? 5 : 0) ;  // precoding
  DLSCH_alloc_pdu2_1E[1].dl_power_off     = (transmission_mode==5 ? 0 : 1);

  eNB2UE[0] = new_channel_desc_scm(PHY_vars_eNB->lte_frame_parms.nb_antennas_tx,
1007
1008
                                   PHY_vars_UE->lte_frame_parms.nb_antennas_rx,
                                   channel_model,
Florian Kaltenberger's avatar
Florian Kaltenberger committed
1009
1010
                                   N_RB2sampling_rate(PHY_vars_eNB->lte_frame_parms.N_RB_DL),
				   N_RB2channel_bandwidth(PHY_vars_eNB->lte_frame_parms.N_RB_DL),
1011
1012
1013
1014
1015
1016
                                   forgetting_factor,
                                   rx_sample_offset,
                                   0);

  if(num_rounds>1) {
    for(n=1; n<4; n++)
1017
      eNB2UE[n] = new_channel_desc_scm(PHY_vars_eNB->lte_frame_parms.nb_antennas_tx,
1018
1019
                                       PHY_vars_UE->lte_frame_parms.nb_antennas_rx,
                                       channel_model,
Florian Kaltenberger's avatar
Florian Kaltenberger committed
1020
1021
1022
				       N_RB2sampling_rate(PHY_vars_eNB->lte_frame_parms.N_RB_DL),
				       N_RB2channel_bandwidth(PHY_vars_eNB->lte_frame_parms.N_RB_DL),
				       forgetting_factor,
1023
1024
                                       rx_sample_offset,
                                       0);
1025
  }
1026

1027
1028
1029
1030
1031
  if (eNB2UE[0]==NULL) {
    msg("Problem generating channel model. Exiting.\n");
    exit(-1);
  }

1032
  if ((transmission_mode == 3) || (transmission_mode==4) || (transmission_mode==8))  
1033
1034
1035
1036
    Kmimo=2;
  else
    Kmimo=1;

1037
  for (k=0; k<n_users; k++) {
1038
    // Create transport channel structures for 2 transport blocks (MIMO)
Elena Lukashova's avatar
Elena Lukashova committed
1039
1040

    for (i=0;i<2;i++) { 
1041
      PHY_vars_eNB->dlsch_eNB[k][i] = new_eNB_dlsch(Kmimo,8,N_RB_DL,0);
1042

1043
      if (!PHY_vars_eNB->dlsch_eNB[k][i]) {
1044
1045
        printf("Can't get eNB dlsch structures\n");
        exit(-1);
1046
      }
1047

1048
1049
1050
      PHY_vars_eNB->dlsch_eNB[k][i]->rnti = n_rnti+k;
    }
  }
1051
1052

  for (i=0; i<2; i++) {
1053
    PHY_vars_UE->dlsch_ue[0][i]  = new_ue_dlsch(Kmimo,8,MAX_TURBO_ITERATIONS,N_RB_DL,0);
1054

1055
1056
1057
    if (!PHY_vars_UE->dlsch_ue[0][i]) {
      printf("Can't get ue dlsch structures\n");
      exit(-1);
1058
1059
    }

1060
1061
    PHY_vars_UE->dlsch_ue[0][i]->rnti   = n_rnti;
  }
1062

1063
1064
1065
  // structure for SIC at UE
  PHY_vars_UE->dlsch_eNB[0] = new_eNB_dlsch(Kmimo,8,N_RB_DL,0);

1066
1067
1068
  if (DLSCH_alloc_pdu2_1E[0].tpmi == 5) {

    PHY_vars_eNB->eNB_UE_stats[0].DL_pmi_single = (unsigned short)(taus()&0xffff);
1069

1070
    if (n_users>1)
1071
1072
      PHY_vars_eNB->eNB_UE_stats[1].DL_pmi_single = (PHY_vars_eNB->eNB_UE_stats[0].DL_pmi_single ^ 0x1555); //opposite PMI
  } else {
1073
    PHY_vars_eNB->eNB_UE_stats[0].DL_pmi_single = 0;
1074

1075
1076
1077
1078
1079
    if (n_users>1)
      PHY_vars_eNB->eNB_UE_stats[1].DL_pmi_single = 0;
  }


1080
1081
1082
1083
1084
  sic_buffer = (mod_sym_t **) malloc16(frame_parms->nb_antennas_tx*sizeof(mod_sym_t*) );
  for (i=0; i<frame_parms->nb_antennas_tx; i++) {
    sic_buffer[i] = malloc16_clear(FRAME_LENGTH_COMPLEX_SAMPLES_NO_PREFIX*sizeof(mod_sym_t));
  }

1085
1086
1087
1088
  if (input_fd==NULL) {


    /*
1089
    // common DCI
1090
1091
1092
1093
1094
1095
1096
1097
1098
    memcpy(&dci_alloc[num_dci].dci_pdu[0],&CCCH_alloc_pdu,sizeof(DCI1A_5MHz_TDD_1_6_t));
    dci_alloc[num_dci].dci_length = sizeof_DCI1A_5MHz_TDD_1_6_t;
    dci_alloc[num_dci].L          = 2;
    dci_alloc[num_dci].rnti       = SI_RNTI;
    num_dci++;
    num_common_dci++;
    */

    // UE specific DCI
1099
    for(k=0; k<n_users; k++) {
1100
1101
1102
1103
1104
1105
1106
1107
1108
1109
1110
1111
1112
1113
1114
1115
      switch(transmission_mode) {
      case 1:
      case 2:
	if (common_flag == 0) {
	  
	  if (PHY_vars_eNB->lte_frame_parms.frame_type == TDD) {
	    
	    switch (PHY_vars_eNB->lte_frame_parms.N_RB_DL) {
	    case 6:
	      dci_length = sizeof_DCI1_1_5MHz_TDD_t;
	      dci_length_bytes = sizeof(DCI1_1_5MHz_TDD_t);
	      ((DCI1_1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rah              = 0;
	      ((DCI1_1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = DLSCH_RB_ALLOC;
	      ((DCI1_1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->TPC              = 0;
	      ((DCI1_1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->dai              = 0;
	      ((DCI1_1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1116
	      ((DCI1_1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
1117
1118
1119
1120
1121
1122
1123
1124
1125
1126
1127
	      ((DCI1_1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->ndi             = 1;
	      ((DCI1_1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rv              = 0;
	      break;
	    case 25:
	      dci_length = sizeof_DCI1_5MHz_TDD_t;
	      dci_length_bytes = sizeof(DCI1_5MHz_TDD_t);
	      ((DCI1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rah              = 0;
	      ((DCI1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = DLSCH_RB_ALLOC;
	      ((DCI1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->TPC              = 0;
	      ((DCI1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->dai              = 0;
	      ((DCI1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1128
	      ((DCI1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
1139
	      ((DCI1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->ndi             = 1;
	      ((DCI1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rv              = 0;
	      break;
	    case 50:
	      dci_length = sizeof_DCI1_10MHz_TDD_t;
	      dci_length_bytes = sizeof(DCI1_10MHz_TDD_t);
	      ((DCI1_10MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rah              = 0;
	      ((DCI1_10MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = DLSCH_RB_ALLOC;
	      ((DCI1_10MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->TPC              = 0;
	      ((DCI1_10MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->dai              = 0;
	      ((DCI1_10MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1140
	      ((DCI1_10MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
1141
1142
1143
1144
1145
1146
1147
1148
1149
	      ((DCI1_10MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->ndi             = 1;
	      ((DCI1_10MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rv              = 0;
	      break;
	    case 100:
	      ((DCI1_20MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rah              = 0;
	      ((DCI1_20MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = DLSCH_RB_ALLOC;
	      ((DCI1_20MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->TPC              = 0;
	      ((DCI1_20MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->dai              = 0;
	      ((DCI1_20MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1150
	      ((DCI1_20MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
1151
1152
1153
1154
1155
1156
1157
	      ((DCI1_20MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->ndi             = 1;
	      ((DCI1_20MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rv              = 0;
	      dci_length = sizeof_DCI1_20MHz_TDD_t;
	      dci_length_bytes = sizeof(DCI1_20MHz_TDD_t);
	      break;
	    }
	  }
1158
	  else { //fdd
1159
1160
1161
1162
1163
1164
1165
1166
	    switch (PHY_vars_eNB->lte_frame_parms.N_RB_DL) {
	    case 6:
	      dci_length = sizeof_DCI1_1_5MHz_FDD_t;
	      dci_length_bytes = sizeof(DCI1_1_5MHz_FDD_t);
	      ((DCI1_1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rah              = 0;
	      ((DCI1_1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = DLSCH_RB_ALLOC;
	      ((DCI1_1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->TPC              = 0;
	      ((DCI1_1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1167
	      ((DCI1_1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
1168
1169
1170
1171
1172
1173
1174
1175
1176
1177
	      ((DCI1_1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->ndi             = 1;
	      ((DCI1_1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rv              = 0;
	      break;
	    case 25:
	      dci_length = sizeof_DCI1_5MHz_FDD_t;
	      dci_length_bytes = sizeof(DCI1_5MHz_FDD_t);
	      ((DCI1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rah              = 0;
	      ((DCI1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = DLSCH_RB_ALLOC;
	      ((DCI1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->TPC              = 0;
	      ((DCI1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1178
	      ((DCI1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
1179
1180
1181
1182
1183
1184
1185
1186
1187
1188
	      ((DCI1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->ndi             = 1;
	      ((DCI1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rv              = 0;
	      break;
	    case 50:
	      dci_length = sizeof_DCI1_10MHz_FDD_t;
	      dci_length_bytes = sizeof(DCI1_10MHz_FDD_t);
	      ((DCI1_10MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rah              = 0;
	      ((DCI1_10MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = DLSCH_RB_ALLOC;
	      ((DCI1_10MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->TPC              = 0;
	      ((DCI1_10MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1189
	      ((DCI1_10MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
1190
1191
1192
1193
1194
1195
1196
1197
1198
1199
	      ((DCI1_10MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->ndi             = 1;
	      ((DCI1_10MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rv              = 0;
	      break;
	    case 100:
	      dci_length = sizeof_DCI1_20MHz_FDD_t;
	      dci_length_bytes = sizeof(DCI1_20MHz_FDD_t);
	      ((DCI1_20MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rah              = 0;
	      ((DCI1_20MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = DLSCH_RB_ALLOC;
	      ((DCI1_20MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->TPC              = 0;
	      ((DCI1_20MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1200
	      ((DCI1_20MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
1201
1202
1203
1204
1205
1206
1207
	      ((DCI1_20MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->ndi             = 1;
	      ((DCI1_20MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rv              = 0;
	      break;
	    }	  
	  }
	  memcpy(&dci_alloc[num_dci].dci_pdu[0],&DLSCH_alloc_pdu_1[k],dci_length_bytes);
	  dci_alloc[num_dci].dci_length = dci_length;
1208
	  dci_alloc[num_dci].L          = 1;
1209
1210
	  dci_alloc[num_dci].rnti       = n_rnti+k;
	  dci_alloc[num_dci].format     = format1;
knopp's avatar
   
knopp committed
1211
1212
1213
1214
	  dump_dci(&PHY_vars_eNB->lte_frame_parms,&dci_alloc[num_dci]);	

	  printf("Generating dlsch params for user %d\n",k);
	  generate_eNB_dlsch_params_from_dci(0,
kaltenbe's avatar
kaltenbe committed
1215
					     subframe,
knopp's avatar
   
knopp committed
1216
1217
1218
1219
1220
1221
1222
1223
1224
1225
1226
1227
1228
					     &DLSCH_alloc_pdu_1[0],
					     n_rnti+k,
					     format1,
					     PHY_vars_eNB->dlsch_eNB[0],
					     &PHY_vars_eNB->lte_frame_parms,
					     PHY_vars_eNB->pdsch_config_dedicated,
					     SI_RNTI,
					     0,
					     P_RNTI,
					     PHY_vars_eNB->eNB_UE_stats[0].DL_pmi_single);
	  num_dci++;
	  num_ue_spec_dci++;
	}
1229
	else { //common flag =1
knopp's avatar
   
knopp committed
1230
1231
1232
1233
1234
1235
1236
1237
	  if (PHY_vars_eNB->lte_frame_parms.frame_type == TDD) {
	    
	    switch (PHY_vars_eNB->lte_frame_parms.N_RB_DL) {
	    case 6:
	      dci_length = sizeof_DCI1A_1_5MHz_TDD_1_6_t;
	      dci_length_bytes = sizeof(DCI1A_1_5MHz_TDD_1_6_t);
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->type             = 1;
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->vrb_type         = 0;
1238
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = computeRIV(PHY_vars_eNB->lte_frame_parms.N_RB_DL,0,9);
knopp's avatar
   
knopp committed
1239
1240
1241
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->TPC              = TPC;
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->dai              = 0;
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1242
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
knopp's avatar
   
knopp committed
1243
1244
1245
1246
1247
1248
1249
1250
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->ndi             = 0;
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->rv              = 0;
	      break;
	    case 25:
	      dci_length = sizeof_DCI1A_5MHz_TDD_1_6_t;
	      dci_length_bytes = sizeof(DCI1A_5MHz_TDD_1_6_t);
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->type             = 1;
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->vrb_type         = 0;
1251
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = computeRIV(PHY_vars_eNB->lte_frame_parms.N_RB_DL,0,9);
knopp's avatar
   
knopp committed
1252
1253
1254
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->TPC              = TPC;
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->dai              = 0;
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1255
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
knopp's avatar
   
knopp committed
1256
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->ndi             = 0;
1257
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->rv              = 0;
knopp's avatar
   
knopp committed
1258
1259
1260
1261
1262
1263
	      break;
	    case 50:
	      dci_length = sizeof_DCI1A_10MHz_TDD_1_6_t;
	      dci_length_bytes = sizeof(DCI1A_10MHz_TDD_1_6_t);
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->type             = 1;
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->vrb_type         = 1;
1264
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = computeRIV(PHY_vars_eNB->lte_frame_parms.N_RB_DL,0,9);
knopp's avatar
   
knopp committed
1265
1266
1267
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->TPC              = TPC;
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->dai              = 0;
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1268
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
knopp's avatar
   
knopp committed
1269
1270
1271
1272
1273
1274
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->ndi             = 0;
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->rv              = 0;
	      break;
	    case 100:
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->type             = 1;
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->vrb_type         = 1;
1275
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->rballoc          = computeRIV(PHY_vars_eNB->lte_frame_parms.N_RB_DL,0,9);
knopp's avatar
   
knopp committed
1276
1277
1278
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->TPC              = TPC;
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->dai              = 0;
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->harq_pid         = 0;
1279
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->mcs             = mcs1;  
knopp's avatar
   
knopp committed