dlsim.c 199 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;
62

Elena Lukashova's avatar
Elena Lukashova committed
63
64
65
extern uint8_t interf_unaw_shift0;
extern uint8_t interf_unaw_shift1;
extern uint8_t interf_unaw_shift;
66

67
68
#include "PHY/TOOLS/lte_phy_scope.h"

69
#define PRINT_BYTES
70
71
72
73
74

PHY_VARS_eNB *PHY_vars_eNB;
PHY_VARS_UE *PHY_vars_UE;

int otg_enabled=0;
75
76
77
78
79
80
81
/*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 */
82

83
84
void handler(int sig)
{
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
  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);
}



//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
108
109
void do_OFDM_mod_l(mod_sym_t **txdataF, int32_t **txdata, uint16_t next_slot, LTE_DL_FRAME_PARMS *frame_parms)
{
110
111
112
113
114
115
116
117
118
119
120

  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
121
                   &txdata[aa][slot_offset],         // output
122
                   frame_parms->ofdm_symbol_size,
123
124
125
                   6,                 // number of symbols
                   frame_parms->nb_prefix_samples,               // number of prefix samples
                   CYCLIC_PREFIX);
126
127
    else {
      normal_prefix_mod(&txdataF[aa][slot_offset_F],
128
129
130
                        &txdata[aa][slot_offset],
                        7,
                        frame_parms);
131
    }
132

133
134
135
136
137

  }

}

138
139
int main(int argc, char **argv)
{
140

141
  int c;
142
143
144
145
  int k,i,aa,aarx,aatx;

  int s,Kr,Kr_bytes;

146
  double sigma2, sigma2_dB=10,SNR,snr0=-2.0,snr1,rate;
147
148
149
150
151
152
153
  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;

154
  uint8_t extended_prefix_flag=0,transmission_mode=1,n_tx=1,n_rx=1;
gauthier's avatar
gauthier committed
155
  uint16_t Nid_cell=0;
156

157
158
  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;
159
160
161
  unsigned char i_mod = 2;
  unsigned short NB_RB;
  unsigned char Ns,l,m;
gauthier's avatar
gauthier committed
162
163
  uint16_t tdd_config=3;
  uint16_t n_rnti=0x1234;
164
  int n_users = 1;
165
166
  RX_type_t rx_type=rx_standard;

167
168
169
170

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

171
172
  unsigned char *input_buffer0[2],*input_buffer1[2];
  unsigned short input_buffer_length0,input_buffer_length1;
173
  unsigned int ret[2];
174
  unsigned int coded_bits_per_codeword=0,nsymb,dci_cnt,tbs=0;
175

176
  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;
177

178
  int re_allocated;
knopp's avatar
knopp committed
179
  char fname[32],vname[32];
180
181
  FILE *bler_fd;
  char bler_fname[256];
182
183
  FILE *time_meas_fd;
  char time_meas_fname[256];
184

185
186
  //FILE *tikz_fd;
  //char tikz_fname[256];
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206

  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
207
208
209
  uint8_t num_pdcch_symbols=1,num_pdcch_symbols_2=0;
  uint8_t pilot1,pilot2,pilot3;
  uint8_t rx_sample_offset = 0;
210
211
  //char stats_buffer[4096];
  //int len;
212

213
  uint8_t num_rounds = 4;
gauthier's avatar
gauthier committed
214
  uint8_t subframe=7;
215
216
217
218
219
220
221
222
223
  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;
224

225
226
227
  // void *data;
  // int ii;
  int bler;
228
  double blerr[2][4],uncoded_ber,avg_ber;
229
  short *uncoded_ber_bit=NULL;
gauthier's avatar
gauthier committed
230
231
  uint8_t N_RB_DL=25,osf=1;
  uint8_t fdd_flag = 0;
Elena Lukashova's avatar
Elena Lukashova committed
232
  int xforms = 0;
233
  frame_t frame_type = FDD;
234
235
  FD_lte_phy_scope_ue *form_ue;
  char title[255];
gauthier's avatar
gauthier committed
236
  uint32_t DLSCH_RB_ALLOC = 0x1fff;
237
238
239
240
241
  int numCCE=0;
  int dci_length_bytes=0,dci_length=0;
  int common_flag=0,TPC=0;

  double cpu_freq_GHz;
242
  //  time_stats_t ts;//,sts,usts;
243
  int avg_iter[2],iter_trials[2];
244
245
  int rballocset=0;
  int print_perf=0;
246
  int test_perf=0;
nikaeinn's avatar
nikaeinn committed
247
  int dump_table=0;
248
249
  int llr8_flag=0;

250
  double effective_rate=0.0;
251
  char channel_model_input[17]="I";
252
253

  int TB0_active = 1;
knopp's avatar
   
knopp committed
254
  uint32_t perfect_ce = 0;
255

Elena Lukashova's avatar
Elena Lukashova committed
256
257
 LTE_DL_UE_HARQ_t *dlsch0_ue_harq;
 LTE_DL_eNB_HARQ_t *dlsch0_eNB_harq;
258
  uint8_t Kmimo;
knopp's avatar
knopp committed
259

260
  mod_sym_t **sic_buffer;
Elena Lukashova's avatar
Elena Lukashova committed
261
262
263
264
  int8_t cw_to_decode_interf;
  int8_t cw_to_decode_interf_free;
  int8_t  cw_non_sic;
  int8_t  cw_sic;
265
266
  FILE    *proc_fd = NULL;
  char buf[64];
267
268
  uint8_t ue_category=4;
  uint32_t Nsoft;
269
270

  int CCE_table[800];
knopp's avatar
knopp committed
271

272
273
  int threequarter_fs=0;

274

275

knopp's avatar
   
knopp committed
276
277
  opp_enabled=1; // to enable the time meas

278
#if defined(__arm__)
279
280
281
  FILE    *proc_fd = NULL;
  char buf[64];

282
283
284
285
286
287
288
289
290
291
292
293
  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
294
295
  printf("Detected cpu_freq %f GHz\n",cpu_freq_GHz);

296
297
  //signal(SIGSEGV, handler);
  //signal(SIGABRT, handler);
298
299

  logInit();
300
  set_glog(LOG_INFO, LOG_MED);
301
302
303
304

  // default parameters
  n_frames = 1000;
  snr0 = 0;
305
  //  num_layers = 1;
knopp's avatar
   
knopp committed
306
  perfect_ce = 0;
307

308
  while ((c = getopt (argc, argv, "ahdpZDe:Em: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:XYv:W:J:")) != -1) {
309

310
311
312
313
314
315
    switch (c) {
    case 'a':
      awgn_flag = 1;
      channel_model = AWGN;
      break;

316
317
318
319
    case 'A':
      abstx = 1;
      break;

320
321
322
323
324
325
326
327
    case 'b':
      tdd_config=atoi(optarg);
      break;

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

328
329
330
331
332
333
334
335
    case 'c':
      num_pdcch_symbols=atoi(optarg);
      break;

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

336
337
338
339
    case 'd':
      dci_flag = 1;
      break;

340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
    case 'D':
      frame_type=TDD;
      break;

    case 'e':
      num_rounds=1;
      common_flag = 1;
      TPC = atoi(optarg);
      break;
      
    case 'E':
      threequarter_fs=1;
      break;

    case 'f':
      input_snr_step= atof(optarg);
      break;

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

    case 'i':
      input_fd = fopen(optarg,"r");
      input_file=1;
      dci_flag = 1;
      break;

    case 'I':
      input_trch_fd = fopen(optarg,"r");
      input_trch_file=1;
      break;

    case 'L':
      llr8_flag=1;
      break;

    case 'l':
      offset_mumimo_llr_drange_fix=atoi(optarg);
      break;

381
382
383
384
385
386
    case 'm':
      mcs1 = atoi(optarg);
      break;

    case 'M':
      mcs2 = atoi(optarg);
387
      //i_mod = get_Qm(mcs2); /// think here again!!!
388
389
      break;

390
391
392
393
394
    case 'O':
      test_perf=atoi(optarg);
      //print_perf =1;
      break;

395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
    case 't':
      mcs_i = atoi(optarg);
      i_mod = get_Qm(mcs_i);
      break;

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

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

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

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

    case 'w':
      snr_int = atof(optarg);
      break;
420
 
421
422
423
424
425
426
427
428
429
    case 'N':
      n_ch_rlz= atof(optarg);
      break;

    case 'p':
      extended_prefix_flag=1;
      break;

    case 'g':
Elena Lukashova's avatar
Elena Lukashova committed
430
      memcpy(channel_model_input,optarg,17);
431
432
433
434
435

      switch((char)*optarg) {
      case 'A':
        channel_model=SCM_A;
        break;
436
      case 'B':
437
438
439
        channel_model=SCM_B;
        break;

440
      case 'C':
441
442
443
        channel_model=SCM_C;
        break;

444
      case 'D':
445
446
447
448
449
450
451
        channel_model=SCM_D;
        break;

      case 'E':
        channel_model=EPA;
        break;

452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
	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;
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
	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;
497
498
499
500
501
502
503
504
505
	default:
	  msg("Unsupported channel model!\n");
	  exit(-1);
	}
	break;
      case 'x':
	transmission_mode=atoi(optarg);
	if ((transmission_mode!=1) &&
	    (transmission_mode!=2) &&
506
	    (transmission_mode!=3) &&
507
	    (transmission_mode!=4) &&
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
	    (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;
531
    
532
533
534
535
536
537
538
539
540
541
      case 'R':
	num_rounds=atoi(optarg);
	break;
      case 'S':
	subframe=atoi(optarg);
	break;
      case 'T':
	n_rnti=atoi(optarg);
	break;	
      case 'u':
542
	rx_type = (RX_type_t) atoi(optarg);
Elena Lukashova's avatar
Elena Lukashova committed
543
	if (rx_type<rx_standard || rx_type>rx_SIC_dual_stream) {
544
	  printf("Unsupported rx type %d\n",rx_type);
545
546
547
548
549
550
551
552
553
554
555
556
557
	  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;
558
559
560
      case 'X':
	xforms = 1;
	break;
561
      case 'Z':
nikaeinn's avatar
nikaeinn committed
562
563
	dump_table=1;
	break;
knopp's avatar
   
knopp committed
564
565
566
      case 'Y':
	perfect_ce=1;
	break;
567
568
569
570
571
572
      case 'V':
	interf_unaw_shift0=atoi(optarg);
	break;	 
      case 'W':
	interf_unaw_shift1=atoi(optarg);
	break;
573
574
575
      case 'J':
	interf_unaw_shift=atoi(optarg);
	break;
576
577
      case 'h':
      default:
578
	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]);
579
580
581
	printf("-h This message\n");
	printf("-a Use AWGN channel and not multipath\n");
	printf("-c Number of PDCCH symbols\n");
582
583
	printf("-m MCS1 for TB 1\n");
	printf("-M MCS2 for TB 2\n");
584
585
586
587
588
589
590
	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");
591
	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");
592
593
594
595
596
597
	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");
598
	printf("-A Turns on calibration mode for abstraction.\n");
Florian Kaltenberger's avatar
   
Florian Kaltenberger committed
599
	printf("-N Determines the number of Channel Realizations in Abstraction mode. Default value is 1. \n");
600
	printf("-O Set the percenatge of effective rate to testbench the modem performance (typically 30 and 70, range 1-100) \n");
601
	printf("-I Input filename for TrCH data (binary)\n");
602
	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");
603
604
	exit(1);
	break;
Elena Lukashova's avatar
Elena Lukashova committed
605
    
606
   }
607
  }
Elena Lukashova's avatar
Elena Lukashova committed
608
    
609
  if (common_flag == 0) {
610
611
612
613
614
    switch (N_RB_DL) {
    case 6:
      if (rballocset==0) DLSCH_RB_ALLOC = 0x3f;
      num_pdcch_symbols = 3;
      break;
615

616
617
618
    case 25:
      if (rballocset==0) DLSCH_RB_ALLOC = 0x1fff;
      break;
619

620
621
622
    case 50:
      if (rballocset==0) DLSCH_RB_ALLOC = 0x1ffff;
      break;
623

624
625
626
627
    case 100:
      if (rballocset==0) DLSCH_RB_ALLOC = 0x1ffffff;
      break;
    }
628

629
    NB_RB=conv_nprb(0,DLSCH_RB_ALLOC,N_RB_DL);
630
  } else
631
632
    NB_RB = 4;

633
  NB_RB=conv_nprb(0,DLSCH_RB_ALLOC,N_RB_DL); 
634

635
  if ((transmission_mode > 1) && (n_tx != 2)) {
636

637
    printf("n_tx must be >1 for transmission_mode %d\n",transmission_mode);
638
639
640
641
642
643
644
    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);
  }
Elena Lukashova's avatar
Elena Lukashova committed
645
  if (((transmission_mode==5) || (transmission_mode==6)) && (rx_type > rx_IC_single_stream)) {
646
647
648
649
    printf("only standard rx or single stream IC available for TM5 and TM6\n");
    exit(-1);
  }

650

651
  if (xforms==1) {
652
653
654
  fl_initialize (&argc, argv, NULL, 0, 0);
  form_ue = create_lte_phy_scope_ue();
  sprintf (title, "LTE PHY SCOPE eNB");
655
  fl_show_form (form_ue->lte_phy_scope_ue, FL_PLACE_HOTSPOT, FL_FULLBORDER, title);
656

657
658
  /*
  if (rx_type==rx_IC_single_stream) {
659
660
661
662
663
    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);
  }
664

665
666
  */
  }
Elena Lukashova's avatar
Elena Lukashova committed
667

668
  if (transmission_mode==5) { 
669
    n_users = 2;
670
671
    eNB_id_i = PHY_vars_UE->n_connected_eNB;
    //eNB_id_i=1;
Elena Lukashova's avatar
Elena Lukashova committed
672
      
673
  }
674
675
676
  else 
    eNB_id_i = eNB_id;

677
678
679
680
681
682
683
684
685
686
687
  lte_param_init(n_tx,
		 n_rx,
		 transmission_mode,
		 extended_prefix_flag,
		 frame_type,
		 Nid_cell,
		 tdd_config,
		 N_RB_DL,
		 threequarter_fs,
		 osf,
		 perfect_ce);
688
689


690
691
  printf("Setting mcs1 = %d\n",mcs1);
  printf("Setting mcs2 = %d\n",mcs2);
692
693
694
  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);
695
  printf("Using receiver type %d\n", rx_type);
Elena Lukashova's avatar
Elena Lukashova committed
696
697
698
  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);
699
700
701
702
703
704
705
  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);
706

707
708
709
710
711
712
713
    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;

714
  s_re = malloc(2*sizeof(double*)); //transmitted signal (Sent)
715
  s_im = malloc(2*sizeof(double*));
716
  r_re = malloc(2*sizeof(double*)); //received signal
717
718
719
720
721
722
  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;

723
  printf("Channel Model= (%s,%d)\n",channel_model_input, channel_model);
724
725
  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);
726
  if(transmission_mode==5)
727
728
729
    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)
730
	 sprintf(bler_fname,"bler_tx%d_r%d_ch%d_%d_nrx%d_mcs%d_mcsi%d_ab_pce_sh%d_pmnew.csv",transmission_mode,rx_type,channel_model,n_frames, n_rx,mcs1, mcs2,interf_unaw_shift );
731
  	else
732
    	sprintf(bler_fname,"bler_tx%d_r%d_ch%d_%d_nrx%d_mcs%d_mcsi%d_ab_sh%d_pmnew.csv",transmission_mode,rx_type,channel_model, n_frames, n_rx,mcs1, mcs2,interf_unaw_shift );
Elena Lukashova's avatar
Elena Lukashova committed
733
  else //abstx=0
734
    if (perfect_ce==1)
735
	sprintf(bler_fname,"bler_tx%d_r%d_ch%d_%d_nrx%d_mcs%d_mcsi%d_pce_sh%d_pmnew.csv",transmission_mode,rx_type,channel_model,n_frames, n_rx,mcs1, mcs2, interf_unaw_shift);
Elena Lukashova's avatar
Elena Lukashova committed
736
   else
737
        sprintf(bler_fname,"bler_tx%d_r%d_ch%d_%d_nrx%d_mcs%d_mcsi%d_sh%d_pmnew.csv",transmission_mode,rx_type,channel_model,n_frames,n_rx,mcs1, mcs2, interf_unaw_shift);
Elena Lukashova's avatar
Elena Lukashova committed
738
  
739
  bler_fd = fopen(bler_fname,"w");
740
741
742
743
  if (bler_fd==NULL) {
    fprintf(stderr,"Cannot create file %s!\n",bler_fname);
    exit(-1);
  }
744

745
  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
746
    if (test_perf != 0) {
747
    
748
    char hostname[1024];
Elena Lukashova's avatar
Elena Lukashova committed
749

750
751
752
    hostname[1023] = '\0';
    gethostname(hostname, 1023);
    printf("Hostname: %s\n", hostname);
753
754
755
756
757
    //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);
758
    time_meas_fd = fopen(time_meas_fname,"w");
759
760
761
762
    if (time_meas_fd==NULL) {
      fprintf(stderr,"Cannot create file %s!\n",time_meas_fname);
      exit(-1);
    }
763
  }
764
765
  
  if(abstx){
766
767
768
769
770
771
772
    // 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)
773
		sprintf(csv_fname,"dout_tx%d_r%d_mcs%d_mcsi%d_ch%d_ns%d_R%d_pce_sh%d_%d_pmnew.m",transmission_mode,rx_type,mcs1,mcs2,channel_model,n_frames,num_rounds, interf_unaw_shift, n_ch_rlz);
774
 	else
775
		sprintf(csv_fname,"dout_tx%d_r%d_mcs%d_mcsi%d_ch%d_ns%d_R%d_sh%d_%d_pmnew.m",transmission_mode,rx_type,mcs1,mcs2,channel_model,n_frames,num_rounds, interf_unaw_shift, n_ch_rlz);
776

777
    csv_fd = fopen(csv_fname,"w");
778
    fprintf(csv_fd,"data_all%d=[",mcs1);
779
780
781
782
    if (csv_fd==NULL) {
      fprintf(stderr,"Cannot create file %s!\n",csv_fname);
      exit(-1);
    }
783
784
  }

785
  /*
786
787
  //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);
788
789
  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");
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
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
  //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;
    }
881
882
  */

883
  for (i=0; i<2; i++) {
884
885
886
887
888
889
890
891
892
893
894
895
    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
896
897
898
899
900
  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;
901
902
903
904
905
906
907
908
  // 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;
909
  
910
911
912
  CCCH_alloc_pdu.type               = 0;
  CCCH_alloc_pdu.vrb_type           = 0;
  CCCH_alloc_pdu.rballoc            = CCCH_RB_ALLOC;
913
914
915
  CCCH_alloc_pdu.ndi                = 1;
  CCCH_alloc_pdu.mcs                = 1;
  CCCH_alloc_pdu.harq_pid           = 0;
916
917
918
919
920
921

  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;
922
923
924
925
  //DLSCH_alloc_pdu2_1E[0].tb_swap        = 0;
  DLSCH_alloc_pdu2_1E[0].mcs              = mcs1;
  DLSCH_alloc_pdu2_1E[0].ndi              = 1;
  DLSCH_alloc_pdu2_1E[0].rv               = 0;
926
927
928
929
930
931
932
933
934
935
  // 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;
936
937
938
  DLSCH_alloc_pdu2_1E[1].mcs              = mcs_i;
  DLSCH_alloc_pdu2_1E[1].ndi              = 1;
  DLSCH_alloc_pdu2_1E[1].rv               = 0;
939
940
941
942
943
  // 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,
944
945
                                   PHY_vars_UE->lte_frame_parms.nb_antennas_rx,
                                   channel_model,
Florian Kaltenberger's avatar
Florian Kaltenberger committed
946
947
                                   N_RB2sampling_rate(PHY_vars_eNB->lte_frame_parms.N_RB_DL),
				   N_RB2channel_bandwidth(PHY_vars_eNB->lte_frame_parms.N_RB_DL),
948
949
950
951
952
953
                                   forgetting_factor,
                                   rx_sample_offset,
                                   0);

  if(num_rounds>1) {
    for(n=1; n<4; n++)
954
      eNB2UE[n] = new_channel_desc_scm(PHY_vars_eNB->lte_frame_parms.nb_antennas_tx,
955
956
                                       PHY_vars_UE->lte_frame_parms.nb_antennas_rx,
                                       channel_model,
Florian Kaltenberger's avatar
Florian Kaltenberger committed
957
958
959
				       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,
960
961
                                       rx_sample_offset,
                                       0);
962
  }
963

964
965
966
967
968
  if (eNB2UE[0]==NULL) {
    msg("Problem generating channel model. Exiting.\n");
    exit(-1);
  }

969
  if ((transmission_mode == 3) || (transmission_mode==4) || (transmission_mode==8))  
970
971
972
973
    Kmimo=2;
  else
    Kmimo=1;

974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
  switch (ue_category) {
  case 1:
    Nsoft = 250368;
    break;
  case 2:
  case 3:
    Nsoft = 1237248;
    break;
  case 4:
    Nsoft = 1827072;
    break;
  default:
    printf("Unsupported UE category %d\n",ue_category);
    exit(-1);
    break;
  }

991
  for (k=0; k<n_users; k++) {
992
    // Create transport channel structures for 2 transport blocks (MIMO)
993
    for (i=0; i<2; i++) {
994
      PHY_vars_eNB->dlsch_eNB[k][i] = new_eNB_dlsch(Kmimo,8,Nsoft,N_RB_DL,0);
995

996
      if (!PHY_vars_eNB->dlsch_eNB[k][i]) {
997
998
        printf("Can't get eNB dlsch structures\n");
        exit(-1);
999
      }
1000

1001
1002
1003
      PHY_vars_eNB->dlsch_eNB[k][i]->rnti = n_rnti+k;
    }
  }
1004
1005

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

1008
1009
1010
    if (!PHY_vars_UE->dlsch_ue[0][i]) {
      printf("Can't get ue dlsch structures\n");
      exit(-1);
1011
1012
    }

1013
1014
    PHY_vars_UE->dlsch_ue[0][i]->rnti   = n_rnti;
  }
1015

1016
  // structure for SIC at UE
1017
  PHY_vars_UE->dlsch_eNB[0] = new_eNB_dlsch(Kmimo,8,Nsoft,N_RB_DL,0);
1018

1019
1020
1021
  if (DLSCH_alloc_pdu2_1E[0].tpmi == 5) {

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

1023
    if (n_users>1)
1024
1025
      PHY_vars_eNB->eNB_UE_stats[1].DL_pmi_single = (PHY_vars_eNB->eNB_UE_stats[0].DL_pmi_single ^ 0x1555); //opposite PMI
  } else {
1026
    PHY_vars_eNB->eNB_UE_stats[0].DL_pmi_single = 0;
1027

1028
1029
1030
1031
1032
    if (n_users>1)
      PHY_vars_eNB->eNB_UE_stats[1].DL_pmi_single = 0;
  }


1033
1034
1035
1036
1037
  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));
  }

1038
1039
1040
1041
  if (input_fd==NULL) {


    /*
1042
    // common DCI
1043
1044
1045
1046
1047
1048
1049
1050
1051
    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
1052
    for(k=0; k<n_users; k++) {
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065
1066
1067
1068
      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;
1069
1070
1071
	      ((DCI1_1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1_1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1_1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
1072
1073
1074
1075
1076
1077
1078
1079
1080
	      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;
1081
1082
1083
	      ((DCI1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1_5MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
1084
1085
1086
1087
1088
1089
1090
1091
1092
	      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;
1093
1094
1095
	      ((DCI1_10MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1_10MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1_10MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
1096
1097
1098
1099
1100
1101
1102
	      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;
1103
1104
1105
	      ((DCI1_20MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1_20MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1_20MHz_TDD_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
1106
1107
1108
1109
1110
	      dci_length = sizeof_DCI1_20MHz_TDD_t;
	      dci_length_bytes = sizeof(DCI1_20MHz_TDD_t);
	      break;
	    }
	  }
1111
	  else { //fdd
1112
1113
1114
1115
1116
1117
1118
1119
	    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;
1120
1121
1122
	      ((DCI1_1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1_1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1_1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
1123
1124
1125
1126
1127
1128
1129
1130
	      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;
1131
1132
1133
	      ((DCI1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
1134
1135
1136
1137
1138
1139
1140
1141
	      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;
1142
1143
1144
	      ((DCI1_10MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1_10MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1_10MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
1145
1146
1147
1148
1149
1150
1151
1152
	      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;
1153
1154
1155
	      ((DCI1_20MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1_20MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1_20MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
1156
1157
1158
1159
1160
	      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;
1161
	  dci_alloc[num_dci].L          = 1;
1162
1163
	  dci_alloc[num_dci].rnti       = n_rnti+k;
	  dci_alloc[num_dci].format     = format1;
knopp's avatar
   
knopp committed
1164
1165
1166
1167
	  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
1168
					     subframe,
knopp's avatar
   
knopp committed
1169
1170
1171
1172
1173
1174
1175
1176
1177
1178
1179
1180
1181
					     &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++;
	}
1182
	else { //common flag =1
knopp's avatar
   
knopp committed
1183
1184
1185
1186
1187
1188
1189
1190
	  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;
1191
	      ((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
1192
1193
1194
	      ((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;
1195
1196
1197
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1A_1_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
knopp's avatar
   
knopp committed
1198
1199
1200
1201
1202
1203
	      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;
1204
	      ((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
1205
1206
1207
	      ((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;
1208
1209
1210
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1A_5MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
knopp's avatar
   
knopp committed
1211
1212
1213
1214
1215
	      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;
1216
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->vrb_type         = 0;
1217
	      ((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
1218
1219
1220
	      ((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;
1221
1222
1223
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1A_10MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
knopp's avatar
   
knopp committed
1224
1225
1226
	      break;
	    case 100:
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->type             = 1;
1227
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->vrb_type         = 0;
1228
	      ((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
1229
1230
1231
	      ((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;
1232
1233
1234
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->mcs              = mcs1;  
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->ndi              = 1;
	      ((DCI1A_20MHz_TDD_1_6_t *)&DLSCH_alloc_pdu_1[k])->rv               = 0;
knopp's avatar
   
knopp committed
1235
1236
1237
1238
1239
	      dci_length = sizeof_DCI1A_20MHz_TDD_1_6_t;
	      dci_length_bytes = sizeof(DCI1A_20MHz_TDD_1_6_t);
	      break;
	    }
	  }
1240
	  else { // fdd
knopp's avatar
   
knopp committed
1241
1242
1243
1244
1245
	    switch (PHY_vars_eNB->lte_frame_parms.N_RB_DL) {
	    case 6:
	      dci_length = sizeof_DCI1A_1_5MHz_FDD_t;
	      dci_length_bytes = sizeof(DCI1A_1_5MHz_FDD_t);
	      ((DCI1A_1_5MHz_FDD_t *)&DLSCH_alloc_pdu_1[k])->type             = 1;