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

22
23
24
25
26
27
28
29
30
31
32
/*! \file PHY/LTE_TRANSPORT/ulsch_coding.c
* \brief Top-level routines for coding the ULSCH transport channel as described in 36.212 V8.6 2009-03
* \author R. Knopp
* \date 2011
* \version 0.1
* \company Eurecom
* \email: knopp@eurecom.fr
* \note
* \warning
*/

33
34
#include "PHY/defs_UE.h"
#include "PHY/phy_extern_ue.h"
35
#include "PHY/LTE_UE_TRANSPORT/transport_proto_ue.h"
36
37
#include "PHY/CODING/coding_defs.h"
#include "PHY/CODING/coding_extern.h"
38
#include "PHY/CODING/lte_interleaver_inline.h"
39
#include "PHY/LTE_UE_TRANSPORT/transport_ue.h"
40
#include "common/utils/LOG/vcd_signal_dumper.h"
41

42
//#define DEBUG_ULSCH_CODING
43
44
45
//#define DEBUG_ULSCH_FREE 1

/*
46
47
48
#define is_not_pilot(pilots,first_pilot,re) (pilots==0) || \
  ((pilots==1)&&(first_pilot==1)&&(((re>2)&&(re<6))||((re>8)&&(re<12)))) || \
  ((pilots==1)&&(first_pilot==0)&&(((re<3))||((re>5)&&(re<9)))) \
49
50
51
52
53
54
*/
#define is_not_pilot(pilots,first_pilot,re) (1)




55
void free_ue_ulsch(LTE_UE_ULSCH_t *ulsch) {
56
57
58
59
60
  int i;
  int r;

  if (ulsch) {
#ifdef DEBUG_ULSCH_FREE
61
    printf("Freeing ulsch %p\n",ulsch);
62
#endif
63

64
    for (i=0; i<8; i++) {
65
      if (ulsch->harq_processes[i]) {
66
67
68
69
        if (ulsch->harq_processes[i]->b) {
          free16(ulsch->harq_processes[i]->b,MAX_ULSCH_PAYLOAD_BYTES);
          ulsch->harq_processes[i]->b = NULL;
        }
70

71
        for (r=0; r<MAX_NUM_ULSCH_SEGMENTS; r++) {
72
73
74
75
76
          if (ulsch->harq_processes[i]->c[r]) {
            free16(ulsch->harq_processes[i]->c[r],((r==0)?8:0) + 3+768);
            ulsch->harq_processes[i]->c[r] = NULL;
          }
        }
77
78
79

        free16(ulsch->harq_processes[i],sizeof(LTE_UL_UE_HARQ_t));
        ulsch->harq_processes[i] = NULL;
80
81
      }
    }
82

83
84
85
86
    free16(ulsch,sizeof(LTE_UE_ULSCH_t));
  }
}

87
LTE_UE_ULSCH_t *new_ue_ulsch(unsigned char N_RB_UL, uint8_t abstraction_flag) {
88
89
90
  LTE_UE_ULSCH_t *ulsch;
  unsigned char exit_flag = 0,i,j,r;
  unsigned char bw_scaling =1;
91
92

  switch (N_RB_UL) {
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
    case 6:
      bw_scaling =16;
      break;

    case 25:
      bw_scaling =4;
      break;

    case 50:
      bw_scaling =2;
      break;

    default:
      bw_scaling =1;
      break;
108
  }
109

110
  ulsch = (LTE_UE_ULSCH_t *)malloc16(sizeof(LTE_UE_ULSCH_t));
111

112
  if (ulsch) {
113
    memset(ulsch,0,sizeof(LTE_UE_ULSCH_t));
Bilel's avatar
Bilel committed
114
115
    ulsch->Mlimit = 4;

116
    for (i=0; i<8; i++) {
117
      ulsch->harq_processes[i] = (LTE_UL_UE_HARQ_t *)malloc16(sizeof(LTE_UL_UE_HARQ_t));
118

119
120
      //      printf("ulsch->harq_processes[%d] %p\n",i,ulsch->harq_processes[i]);
      if (ulsch->harq_processes[i]) {
121
        memset(ulsch->harq_processes[i], 0, sizeof(LTE_UL_UE_HARQ_t));
122
        ulsch->harq_processes[i]->b = (unsigned char *)malloc16(MAX_ULSCH_PAYLOAD_BYTES/bw_scaling);
123
124
125
126
127
128
129
130
131
132

        if (ulsch->harq_processes[i]->b)
          memset(ulsch->harq_processes[i]->b,0,MAX_ULSCH_PAYLOAD_BYTES/bw_scaling);
        else {
          LOG_E(PHY,"Can't get b\n");
          exit_flag=1;
        }

        if (abstraction_flag==0) {
          for (r=0; r<MAX_NUM_ULSCH_SEGMENTS; r++) {
133
            ulsch->harq_processes[i]->c[r] = (unsigned char *)malloc16(((r==0)?8:0) + 3+768); // account for filler in first segment and CRCs for multiple segment case
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148

            if (ulsch->harq_processes[i]->c[r])
              memset(ulsch->harq_processes[i]->c[r],0,((r==0)?8:0) + 3+768);
            else {
              LOG_E(PHY,"Can't get c\n");
              exit_flag=2;
            }
          }
        }

        ulsch->harq_processes[i]->subframe_scheduling_flag = 0;
        ulsch->harq_processes[i]->first_tx = 1;
      } else {
        LOG_E(PHY,"Can't get harq_p %d\n",i);
        exit_flag=3;
149
150
151
152
      }
    }

    if ((abstraction_flag == 0) && (exit_flag==0)) {
153
      for (i=0; i<8; i++)
154
155
156
157
        for (j=0; j<96; j++)
          for (r=0; r<MAX_NUM_ULSCH_SEGMENTS; r++)
            ulsch->harq_processes[i]->d[r][j] = LTE_NULL;

158
      return(ulsch);
159
    } else if (abstraction_flag==1)
160
161
      return(ulsch);
  }
162

Cedric Roux's avatar
Cedric Roux committed
163
  LOG_E(PHY,"new_ue_ulsch exit flag, size of  %d ,   %zu\n",exit_flag, sizeof(LTE_UE_ULSCH_t));
164
165
166
167
168
  free_ue_ulsch(ulsch);
  return(NULL);
}


gauthier's avatar
gauthier committed
169
uint32_t ulsch_encoding(uint8_t *a,
170
                        PHY_VARS_UE *ue,
171
172
                        uint8_t harq_pid,
                        uint8_t eNB_id,
Bilel's avatar
Bilel committed
173
                        uint8_t subframe_rx,
174
175
                        uint8_t tmode,
                        uint8_t control_only_flag,
176
                        uint8_t Nbundled) {
177
178
179
180
181
  time_stats_t *seg_stats=&ue->ulsch_segmentation_stats;
  time_stats_t *rm_stats=&ue->ulsch_rate_matching_stats;
  time_stats_t *te_stats=&ue->ulsch_turbo_encoding_stats;
  time_stats_t *i_stats=&ue->ulsch_interleaving_stats;
  time_stats_t *m_stats=&ue->ulsch_multiplexing_stats;
gauthier's avatar
gauthier committed
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
  //  uint16_t offset;
  uint32_t crc=1;
  uint32_t A;
  uint8_t Q_m=0;
  uint32_t Kr=0,Kr_bytes,r,r_offset=0;
  uint8_t y[6*14*1200],*yptr;;
  uint8_t *columnset;
  uint32_t sumKr=0;
  uint32_t Qprime,L,G,Q_CQI=0,Q_RI=0,Q_ACK=0,H=0,Hprime=0,Hpp=0,Cmux=0,Rmux=0,Rmux_prime=0;
  uint32_t Qprime_ACK=0,Qprime_CQI=0,Qprime_RI=0,len_ACK=0,len_RI=0;
  //  uint32_t E;
  uint8_t ack_parity;
  uint32_t i,q,j,iprime,j2;
  uint16_t o_RCC;
  uint8_t o_flip[8];
  uint32_t wACK_idx;
198
199
200
  LTE_DL_FRAME_PARMS *frame_parms=&ue->frame_parms;
  PHY_MEASUREMENTS *meas = &ue->measurements;
  LTE_UE_ULSCH_t *ulsch=ue->ulsch[eNB_id];
Bilel's avatar
Bilel committed
201
  LTE_UE_DLSCH_t **dlsch = ue->dlsch[0][eNB_id];
202
  uint16_t rnti = 0xffff;
203

204
205
206
207
208
  if (!ulsch) {
    LOG_E(PHY,"Null ulsch ptr %p\n",ulsch);
    return(-1);
  }

209
  if (harq_pid >= 8) {
210
211
212
213
    LOG_E(PHY,"Illegal harq_pid %d\n",harq_pid);
    return(-1);
  }

214
215
216
217
218
219
220
  if (ulsch->harq_processes[harq_pid]->O_ACK > 2) {
    LOG_E(PHY,"Illegal O_ACK %d\n",ulsch->harq_processes[harq_pid]->O_ACK);
    return(-1);
  }

  if (ulsch->O_RI > 1) {
    LOG_E(PHY,"Illegal O_RI %d\n",ulsch->O_RI);
221
222
223
    return(-1);
  }

gauthier's avatar
gauthier committed
224
  VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_IN);
225

226
227
  // fill CQI/PMI information
  if (ulsch->O>0) {
228
    VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING_FILL_CQI, VCD_FUNCTION_IN);
229
    rnti = ue->pdcch_vars[ue->current_thread_id[subframe_rx]][eNB_id]->crnti;
230
    fill_CQI(ulsch,meas,0,harq_pid,ue->frame_parms.N_RB_DL,rnti, tmode,ue->sinr_eff);
231
    LOG_D(PHY,"ULSCH Encoding rnti %x \n", rnti);
232
    print_CQI(ulsch->o,ulsch->uci_format,0,ue->frame_parms.N_RB_DL);
233
234

    // save PUSCH pmi for later (transmission modes 4,5,6)
235
236
237
238
    if (dlsch[0]) {
      //LOG_I(PHY,"XXX saving pmi for DL %x\n",pmi2hex_2Ar1(((wideband_cqi_rank1_2A_5MHz *)ulsch->o)->pmi));
      dlsch[0]->pmi_alloc = ((wideband_cqi_rank1_2A_5MHz *)ulsch->o)->pmi;
    }
239

240
    VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING_FILL_CQI, VCD_FUNCTION_OUT);
241
242
243
244
245
246
  }

  if (ulsch->O<=32) {
    o_flip[0] = ulsch->o[3];
    o_flip[1] = ulsch->o[2];
    o_flip[2] = ulsch->o[1];
247
248
    o_flip[3] = ulsch->o[0];
  } else {
249
250
251
252
253
254
255
256
257
    o_flip[0] = ulsch->o[7];
    o_flip[1] = ulsch->o[6];
    o_flip[2] = ulsch->o[5];
    o_flip[3] = ulsch->o[4];
    o_flip[4] = ulsch->o[3];
    o_flip[5] = ulsch->o[2];
    o_flip[6] = ulsch->o[1];
    o_flip[7] = ulsch->o[0];
  }
258

259
260
261
262
263
  if (control_only_flag == 0) {
    A=ulsch->harq_processes[harq_pid]->TBS;
    Q_m = get_Qm_ul(ulsch->harq_processes[harq_pid]->mcs);
    ulsch->harq_processes[harq_pid]->control_only = 0;
#ifdef DEBUG_ULSCH_CODING
264
    printf("[PHY][UE] ULSCH coding : A %d, Qm %d, mcs %d, harq_pid %d, round %d, RV %d\n",
265
266
267
268
269
270
           ulsch->harq_processes[harq_pid]->TBS,
           Q_m,
           ulsch->harq_processes[harq_pid]->mcs,
           harq_pid,
           ulsch->harq_processes[harq_pid]->round,
           ulsch->harq_processes[harq_pid]->rvidx);
271
272

    for (i=0; i<ulsch->harq_processes[harq_pid]->O_ACK; i++)
273
      printf("ulsch_coding: o_ACK[%d] %d\n",i,ulsch->o_ACK[i]);
274
275

    for (i=0; i<ulsch->O_RI; i++)
276
      printf("ulsch_coding: o_RI[%d] %d\n",i,ulsch->o_RI[i]);
277

278
    printf("ulsch_coding: O=%d\n",ulsch->O);
279
280
281

    for (i=0; i<1+((8+ulsch->O)/8); i++) {
      //    ulsch->o[i] = i;
282
      printf("ulsch_coding: O[%d] %d\n",i,ulsch->o[i]);
283
284
285
    }

    if ((tmode != 4))
286
      print_CQI(ulsch->o,wideband_cqi_rank1_2A,0,ue->frame_parms.N_RB_DL);
287
    else
288
      print_CQI(ulsch->o,HLC_subband_cqi_rank1_2A,0,ue->frame_parms.N_RB_DL);
289

290
#endif
291

292
    if (ulsch->harq_processes[harq_pid]->round == 0) {  // this is a new packet
293
294
295
      start_meas(seg_stats);
      // Add 24-bit crc (polynomial A) to payload
      crc = crc24a(a,
296
                   A)>>8;
297
298
299
      a[A>>3] = ((uint8_t *)&crc)[2];
      a[1+(A>>3)] = ((uint8_t *)&crc)[1];
      a[2+(A>>3)] = ((uint8_t *)&crc)[0];
300
301
302
      ulsch->harq_processes[harq_pid]->B = A+24;
      ulsch->harq_processes[harq_pid]->b = a;
      lte_segmentation(ulsch->harq_processes[harq_pid]->b,
303
304
305
306
307
308
309
310
                       ulsch->harq_processes[harq_pid]->c,
                       ulsch->harq_processes[harq_pid]->B,
                       &ulsch->harq_processes[harq_pid]->C,
                       &ulsch->harq_processes[harq_pid]->Cplus,
                       &ulsch->harq_processes[harq_pid]->Cminus,
                       &ulsch->harq_processes[harq_pid]->Kplus,
                       &ulsch->harq_processes[harq_pid]->Kminus,
                       &ulsch->harq_processes[harq_pid]->F);
311
      stop_meas(seg_stats);
312
313
314
315
316
317
318
319

      for (r=0; r<ulsch->harq_processes[harq_pid]->C; r++) {
        if (r<ulsch->harq_processes[harq_pid]->Cminus)
          Kr = ulsch->harq_processes[harq_pid]->Kminus;
        else
          Kr = ulsch->harq_processes[harq_pid]->Kplus;

        Kr_bytes = Kr>>3;
320
#ifdef DEBUG_ULSCH_CODING
321
        printf("Generating Code Segment %d (%d bits)\n",r,Kr);
322
        // generate codewords
323
324
325
326
        printf("bits_per_codeword (Kr)= %d\n",Kr);
        printf("N_RB = %d\n",ulsch->harq_processes[harq_pid]->nb_rb);
        printf("Ncp %d\n",frame_parms->Ncp);
        printf("Qm %d\n",Q_m);
327
#endif
328
329
        //  offset=0;
        start_meas(te_stats);
330
        encoder(ulsch->harq_processes[harq_pid]->c[r],
331
332
333
                Kr>>3,
                &ulsch->harq_processes[harq_pid]->d[r][96],
                (r==0) ? ulsch->harq_processes[harq_pid]->F : 0
334
               );
335
        stop_meas(te_stats);
336
#ifdef DEBUG_ULSCH_CODING
337
338

        if (r==0)
Mongazon's avatar
Mongazon committed
339
          LOG_M("enc_output0.m","enc0",&ulsch->harq_processes[harq_pid]->d[r][96],(3*8*Kr_bytes)+12,1,4);
340

341
#endif
342
343
344
345
346
347
        start_meas(i_stats);
        ulsch->harq_processes[harq_pid]->RTC[r] =
          sub_block_interleaving_turbo(4+(Kr_bytes*8),
                                       &ulsch->harq_processes[harq_pid]->d[r][96],
                                       ulsch->harq_processes[harq_pid]->w[r]);
        stop_meas(i_stats);
348
349
      }
    }
350

351
352
    if (ulsch->harq_processes[harq_pid]->C == 0) {
      LOG_E(PHY,"null segment\n");
gauthier's avatar
gauthier committed
353
      VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
354
355
      return(-1);
    }
356

357
    sumKr = 0;
358
359

    for (r=0; r<ulsch->harq_processes[harq_pid]->C; r++) {
360
      if (r<ulsch->harq_processes[harq_pid]->Cminus)
361
        Kr = ulsch->harq_processes[harq_pid]->Kminus;
362
      else
363
364
        Kr = ulsch->harq_processes[harq_pid]->Kplus;

365
366
      sumKr += Kr;
    }
367
  } else { // This is a control-only PUSCH, set sumKr to O_CQI-MIN
368
369
370
    ulsch->harq_processes[harq_pid]->control_only = 1;
    sumKr = ulsch->O_CQI_MIN;
  }
371

372
373
374
  ulsch->harq_processes[harq_pid]->sumKr = sumKr;
  // Compute Q_ri (p. 23 36-212)
  Qprime = ulsch->O_RI*ulsch->harq_processes[harq_pid]->Msc_initial*ulsch->harq_processes[harq_pid]->Nsymb_initial * ulsch->beta_offset_ri_times8;
375

376
377
378
379
380
  if (Qprime > 0) {
    if ((Qprime % (8*sumKr)) > 0)
      Qprime = 1+(Qprime/(8*sumKr));
    else
      Qprime = Qprime/(8*sumKr);
381

382
383
384
385
386
387
388
389
    if (Qprime > 4*ulsch->harq_processes[harq_pid]->nb_rb * 12)
      Qprime = 4*ulsch->harq_processes[harq_pid]->nb_rb * 12;
  }

  Q_RI = Q_m*Qprime;
  Qprime_RI = Qprime;
  // Compute Q_ack (p. 23 36-212)
  Qprime = ulsch->harq_processes[harq_pid]->O_ACK*ulsch->harq_processes[harq_pid]->Msc_initial*ulsch->harq_processes[harq_pid]->Nsymb_initial * ulsch->beta_offset_harqack_times8;
390

391
392
393
394
395
  if (Qprime > 0) {
    if ((Qprime % (8*sumKr)) > 0)
      Qprime = 1+(Qprime/(8*sumKr));
    else
      Qprime = Qprime/(8*sumKr);
396

397
398
399
400
401
402
    if (Qprime > 4*ulsch->harq_processes[harq_pid]->nb_rb * 12)
      Qprime = 4*ulsch->harq_processes[harq_pid]->nb_rb * 12;
  }

  Q_ACK = Qprime * Q_m;
  Qprime_ACK = Qprime;
403
  LOG_D(PHY,"UE (%x/%d) O_ACK %d, Mcs_initial %d, Nsymb_initial %d, beta_offset_harqack*8 %d, sum Kr %d, Qprime_ACK %d, Q_ACK %d\n",
404
405
406
407
408
409
410
411
        rnti, harq_pid,
        ulsch->harq_processes[harq_pid]->O_ACK,
        ulsch->harq_processes[harq_pid]->Msc_initial,
        ulsch->harq_processes[harq_pid]->Nsymb_initial,
        ulsch->beta_offset_harqack_times8,
        sumKr,
        Qprime_ACK,
        Q_ACK);
412

413
414
415
416
  // Compute Q_cqi, assume O>11, p. 26 36-212
  if (control_only_flag == 0) {
    if (ulsch->O < 12)
      L=0;
417
    else
418
      L=8;
419

420
421
422
423
    if (ulsch->O > 0)
      Qprime = (ulsch->O + L) * ulsch->harq_processes[harq_pid]->Msc_initial*ulsch->harq_processes[harq_pid]->Nsymb_initial * ulsch->beta_offset_cqi_times8;
    else
      Qprime = 0;
424

425
426
427
428
429
430
    LOG_D(PHY,"Qprime %d, O_RI %d + %d, Msc %d, Nym %d beta %d\n",
          Qprime,
          ulsch->O, L,
          ulsch->harq_processes[harq_pid]->Msc_initial,
          ulsch->harq_processes[harq_pid]->Nsymb_initial,
          ulsch->beta_offset_cqi_times8);
431

432
433
    if (Qprime > 0) {
      if ((Qprime % (8*sumKr)) > 0)
434
        Qprime = 1+(Qprime/(8*sumKr));
435
      else
436
        Qprime = Qprime/(8*sumKr);
437
    }
438

439
    G = ulsch->harq_processes[harq_pid]->nb_rb * (12 * Q_m) * (ulsch->Nsymb_pusch);
440
    LOG_D(PHY,"G: rb %d * ( 12 * Qm %d ) * nsymb %d, Qprime %d, O_RI %d\n", ulsch->harq_processes[harq_pid]->nb_rb, Q_m, ulsch->Nsymb_pusch, Qprime, ulsch->O_RI);
441

442
443
    if (Qprime > (G - ulsch->O_RI))
      Qprime = G - ulsch->O_RI;
444

445
446
447
    Q_CQI = Q_m * Qprime;
    Qprime_CQI = Qprime;
    G = G - Q_RI - Q_CQI;
448
    LOG_D(PHY,"new G: %d, Q_RI %d Q_CQI %d\n",  G, Q_RI, Q_CQI);
Gabriel's avatar
Gabriel committed
449
450
    ulsch->harq_processes[harq_pid]->G = G;

451
452
453
454
455
456
457
458
459
460
461
462
463
    /*
        LOG_I(PHY,"ULSCH Encoding G %d, Q_RI %d (O_RI%d, Msc_initial %d, Nsymb_initial%d, beta_offset_ri_times8 %d), Q_CQI %d, Q_ACK %d \n",G,Q_RI,ulsch->O_RI,ulsch->harq_processes[harq_pid]->Msc_initial,ulsch->harq_processes[harq_pid]->Nsymb_initial,ulsch->beta_offset_ri_times8,Q_CQI,Q_ACK);

        LOG_I(PHY,"ULSCH Encoding (Nid_cell %d, rnti %x): harq_pid %d round %d, RV %d, mcs %d, O_RI %d, O_ACK %d, G %d\n",
              frame_parms->Nid_cell,ulsch->rnti,
              harq_pid,
              ulsch->harq_processes[harq_pid]->round,
              ulsch->harq_processes[harq_pid]->rvidx,
              ulsch->harq_processes[harq_pid]->mcs,
              ulsch->O_RI,
              ulsch->harq_processes[harq_pid]->O_ACK,
              G);
    */
464

465
466
    if ((int)G < 0) {
      LOG_E(PHY,"FATAL: ulsch_coding.c G < 0 (%d) : Q_RI %d, Q_CQI %d, O %d, betaCQI_times8 %d)\n",G,Q_RI,Q_CQI,ulsch->O,ulsch->beta_offset_cqi_times8);
gauthier's avatar
gauthier committed
467
      VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
468
469
470
471
472
473
474
475
476
      return(-1);
    }

    // Data and control multiplexing (5.2.2.7 36-212)
    H = G + Q_CQI;
    Hprime = H/Q_m;

    // Fill in the "e"-sequence from 36-212, V8.6 2009-03, p. 16-17 (for each "e") and concatenate the
    // outputs for each code segment, see Section 5.1.5 p.20
477
478

    for (r=0; r<ulsch->harq_processes[harq_pid]->C; r++) {
479
#ifdef DEBUG_ULSCH_CODING
480
      printf("Rate Matching, Code segment %d (coded bits (G) %d,unpunctured/repeated bits per code segment %d,mod_order %d, nb_rb %d)...\n",
481
482
483
484
             r,
             G,
             Kr*3,
             Q_m,ulsch->harq_processes[harq_pid]->nb_rb);
485
#endif
486
      start_meas(rm_stats);
487
      r_offset += lte_rate_matching_turbo(ulsch->harq_processes[harq_pid]->RTC[r],
488
489
490
491
492
                                          G,
                                          ulsch->harq_processes[harq_pid]->w[r],
                                          ulsch->e+r_offset,
                                          ulsch->harq_processes[harq_pid]->C, // C
                                          NSOFT,                    // Nsoft,
493
                                          0,  // this means UL
494
495
496
497
498
                                          1,
                                          ulsch->harq_processes[harq_pid]->rvidx,
                                          get_Qm_ul(ulsch->harq_processes[harq_pid]->mcs),
                                          1,
                                          r,
499
                                          ulsch->harq_processes[harq_pid]->nb_rb);
500
      //ulsch->harq_processes[harq_pid]->mcs);                       // r
501
502
      stop_meas(rm_stats);
#ifdef DEBUG_ULSCH_CODING
503

504
      if (r==ulsch->harq_processes[harq_pid]->C-1)
Mongazon's avatar
Mongazon committed
505
        LOG_M("enc_output.m","enc",ulsch->e,r_offset,1,4);
506

507
508
#endif
    }
509
  } else { //control-only PUSCH
510
511
512
513
514
515
516
517
    Q_CQI = (ulsch->harq_processes[harq_pid]->nb_rb * (12 * Q_m) * (ulsch->Nsymb_pusch)) - Q_RI;
    H = Q_CQI;
    Hprime = H/Q_m;
  }

  //  Do CQI coding
  if ((ulsch->O>1) && (ulsch->O < 12)) {
    LOG_E(PHY,"short CQI sizes not supported yet\n");
gauthier's avatar
gauthier committed
518
    VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
519
    return(-1);
520
  } else {
521
522
    // add 8-bit CRC
    crc = crc8(o_flip,
523
               ulsch->O)>>24;
524
525
526
527
528
#ifdef DEBUG_ULSCH_CODING
    printf("crc(cqi) tx : %x\n",crc);
#endif
    memset((void *)&ulsch->o_d[0],LTE_NULL,96);
    ccodelte_encode(ulsch->O,
529
530
531
532
533
534
535
                    1,
                    o_flip,
                    &ulsch->o_d[96],
                    0);
    o_RCC = sub_block_interleaving_cc(8+ulsch->O,
                                      &ulsch->o_d[96],
                                      ulsch->o_w);
536
    lte_rate_matching_cc(o_RCC,
537
538
539
                         Q_CQI,
                         ulsch->o_w,
                         ulsch->q);
540
541
542
  }

  i=0;
543

544
545
546
  //  Do RI coding
  if (ulsch->O_RI == 1) {
    switch (Q_m) {
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
      case 2:
        ulsch->q_RI[0] = ulsch->o_RI[0];
        ulsch->q_RI[1] = PUSCH_y;//ulsch->o_RI[0];
        len_RI=2;
        break;

      case 4:
        ulsch->q_RI[0] = ulsch->o_RI[0];
        ulsch->q_RI[1] = PUSCH_y;//1;
        ulsch->q_RI[2] = PUSCH_x;//ulsch->o_RI[0];
        ulsch->q_RI[3] = PUSCH_x;//1;
        len_RI=4;
        break;

      case 6:
        ulsch->q_RI[0] = ulsch->o_RI[0];
        ulsch->q_RI[1] = PUSCH_y;//1;
        ulsch->q_RI[2] = PUSCH_x;//1;
        ulsch->q_RI[3] = PUSCH_x;//ulsch->o_RI[0];
        ulsch->q_RI[4] = PUSCH_x;//1;
        ulsch->q_RI[5] = PUSCH_x;//1;
        len_RI=6;
        break;
570
    }
571
  } else if (ulsch->O_RI>1) {
572
    LOG_E(PHY,"RI cannot be more than 1 bit yet\n");
gauthier's avatar
gauthier committed
573
    VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
574
575
    return(-1);
  }
576

577
578
579
  //  Do ACK coding, Section 5.2.2.6 36.213 (p.23-24 in v8.6)
  wACK_idx = (ulsch->bundling==0) ? 4 : ((Nbundled-1)&3);
#ifdef DEBUG_ULSCH_CODING
580
  printf("ulsch_coding.c: Bundling %d, Nbundled %d, wACK_idx %d\n",
581
         ulsch->bundling,Nbundled,wACK_idx);
582
#endif
583

584
585
586
  // 1-bit ACK/NAK
  if (ulsch->harq_processes[harq_pid]->O_ACK == 1) {
    switch (Q_m) {
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
      case 2:
        ulsch->q_ACK[0] = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[1] = (ulsch->bundling==0)? PUSCH_y : ((ulsch->o_ACK[0]+wACK[wACK_idx][1])&1);//ulsch->o_ACK[0];
        len_ACK = 2;
        break;

      case 4:
        ulsch->q_ACK[0] = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[1] = (ulsch->bundling==0)? PUSCH_y : ((ulsch->o_ACK[0]+wACK[wACK_idx][1])&1);
        ulsch->q_ACK[2] = PUSCH_x;
        ulsch->q_ACK[3] = PUSCH_x;
        len_ACK = 4;
        break;

      case 6:
        ulsch->q_ACK[0] = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[1] = (ulsch->bundling==0)? PUSCH_y : ((ulsch->o_ACK[0]+wACK[wACK_idx][1])&1);
        ulsch->q_ACK[2] = PUSCH_x;
        ulsch->q_ACK[3] = PUSCH_x;
        ulsch->q_ACK[4] = PUSCH_x;
        ulsch->q_ACK[6] = PUSCH_x;
        len_ACK = 6;
        break;
610
611
    }
  }
612

613
614
615
  // two-bit ACK/NAK
  if (ulsch->harq_processes[harq_pid]->O_ACK == 2) {
    ack_parity = (ulsch->o_ACK[0]+ulsch->o_ACK[1])&1;
616

617
    switch (Q_m) {
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
      case 2:
        ulsch->q_ACK[0] = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[1] = (ulsch->o_ACK[1]+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[2] = (ack_parity+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[3] = (ulsch->o_ACK[0]+wACK[wACK_idx][1])&1;
        ulsch->q_ACK[4] = (ulsch->o_ACK[1]+wACK[wACK_idx][1])&1;
        ulsch->q_ACK[5] = (ack_parity+wACK[wACK_idx][1])&1;
        len_ACK = 6;
        break;

      case 4:
        ulsch->q_ACK[0]  = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[1]  = (ulsch->o_ACK[1]+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[2]  = PUSCH_x;
        ulsch->q_ACK[3]  = PUSCH_x;//1;
        ulsch->q_ACK[4]  = (ack_parity+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[5]  = (ulsch->o_ACK[0]+wACK[wACK_idx][1])&1;
        ulsch->q_ACK[6]  = PUSCH_x;
        ulsch->q_ACK[7]  = PUSCH_x;//1;
        ulsch->q_ACK[8]  = (ulsch->o_ACK[1]+wACK[wACK_idx][1])&1;
        ulsch->q_ACK[9]  = (ack_parity+wACK[wACK_idx][1])&1;
        ulsch->q_ACK[10] = PUSCH_x;
        ulsch->q_ACK[11] = PUSCH_x;//1;
        len_ACK = 12;
        break;

      case 6:
        ulsch->q_ACK[0] = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[1] = (ulsch->o_ACK[1]+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[2] = PUSCH_x;
        ulsch->q_ACK[3] = PUSCH_x;
        ulsch->q_ACK[4] = PUSCH_x;
        ulsch->q_ACK[5] = PUSCH_x;
        ulsch->q_ACK[6] = (ack_parity+wACK[wACK_idx][0])&1;
        ulsch->q_ACK[7] = (ulsch->o_ACK[0]+wACK[wACK_idx][1])&1;
        ulsch->q_ACK[8] = PUSCH_x;
        ulsch->q_ACK[9] = PUSCH_x;
        ulsch->q_ACK[10] = PUSCH_x;
        ulsch->q_ACK[11] = PUSCH_x;
        ulsch->q_ACK[12] = (ulsch->o_ACK[1]+wACK[wACK_idx][1])&1;
        ulsch->q_ACK[13] = (ack_parity+wACK[wACK_idx][1])&1;
        ulsch->q_ACK[14] = PUSCH_x;
        ulsch->q_ACK[15] = PUSCH_x;
        ulsch->q_ACK[16] = PUSCH_x;
        ulsch->q_ACK[17] = PUSCH_x;
        len_ACK = 18;
        break;
665
666
    }
  }
667

668
669
  if (ulsch->harq_processes[harq_pid]->O_ACK > 2) {
    LOG_E(PHY,"ACK cannot be more than 2 bits yet\n");
gauthier's avatar
gauthier committed
670
    VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
671
672
673
674
675
676
677
678
679
680
681
682
683
    return(-1);
  }

  // channel multiplexing/interleaving
  start_meas(m_stats);
  Hpp = Hprime + Q_RI;
  Cmux       = ulsch->Nsymb_pusch;
  Rmux       = Hpp*Q_m/Cmux;
  Rmux_prime = Rmux/Q_m;
  Qprime_RI  = Q_RI / Q_m;
  Qprime_ACK = Q_ACK / Q_m;
  Qprime_CQI = Q_CQI / Q_m;
  //  printf("Qprime_CQI = %d\n",Qprime_CQI);
684
  // RI BITS
685
686
687
688
689
690
  memset(y,LTE_NULL,Q_m*Hpp);

  if (frame_parms->Ncp == 0)
    columnset = cs_ri_normal;
  else
    columnset = cs_ri_extended;
691
692
693
694

  j=0;

  for (i=0; i<Qprime_RI; i++) {
695
    r = Rmux_prime - 1 - (i>>2);
696
697

    for (q=0; q<Q_m; q++)  {
698
699
700
      y[q+(Q_m*((r*Cmux) + columnset[j]))]  = ulsch->q_RI[(q+(Q_m*i))%len_RI];
      //      printf("ri[%d] %d => y[%d]\n",q+(Q_m*i)%len_RI,ulsch->q_RI[(q+(Q_m*i))%len_RI],q+(Q_m*((r*Cmux) + columnset[j])),y[q+(Q_m*((r*Cmux) + columnset[j]))]);
    }
701

702
703
704
705
706
707
708
709
710
    j=(j+3)&3;
  }

  // CQI and Data bits
  j=0;
  /*
  for (i=0,iprime=-Qprime_CQI;i<Hprime;i++,iprime++) {

    while (y[Q_m*j] != LTE_NULL) j++;
711

712
713
    if (i<Qprime_CQI) {
      for (q=0;q<Q_m;q++) {
714
715
  y[q+(Q_m*j)] = ulsch->q[q+(Q_m*i)];
  //printf("cqi[%d] %d => y[%d]\n",q+(Q_m*i),ulsch->q[q+(Q_m*i)],q+(Q_m*j));
716
717
718
719
      }
    }
    else {
      for (q=0;q<Q_m;q++) {
720
721
  y[q+(Q_m*j)] = ulsch->e[q+(Q_m*iprime)];
  //  printf("e[%d] %d => y[%d]\n",q+(Q_m*iprime),ulsch->e[q+(Q_m*iprime)],q+(Q_m*j));
722
723
724
725
726
727
      }
    }
    j++;
  }
  */

728
  for (i=0; i<Qprime_CQI; i++) {
729
730
    while (y[Q_m*j] != LTE_NULL) j++;

731
    for (q=0; q<Q_m; q++) {
732
733
734
      y[q+(Q_m*j)] = ulsch->q[q+(Q_m*i)];
      //        printf("cqi[%d] %d => y[%d] (j %d)\n",q+(Q_m*i),ulsch->q[q+(Q_m*i)],q+(Q_m*j),j);
    }
735

736
737
738
739
    j++;
  }

  j2 = j*Q_m;
740

741
  switch (Q_m) {
742
743
744
    case 2:
      for (iprime=0; iprime<(Hprime-Qprime_CQI)<<1; iprime+=2) {
        while (y[j2] != LTE_NULL) j2+=2;
745

746
747
748
749
        y[j2]   = ulsch->e[iprime];
        y[1+j2] = ulsch->e[1+iprime];
        j2+=2;
      }
750

751
      break;
752

753
754
755
756
757
758
759
760
761
762
    case 4:
      for (iprime=0; iprime<(Hprime-Qprime_CQI)<<2; iprime+=4) {
        while (y[j2] != LTE_NULL) j2+=4;

        y[j2]   = ulsch->e[iprime];
        y[1+j2] = ulsch->e[1+iprime];
        y[2+j2] = ulsch->e[2+iprime];
        y[3+j2] = ulsch->e[3+iprime];
        j2+=4;
      }
763

764
      break;
765

766
767
768
769
770
771
772
773
774
775
776
777
    case 6:
      for (iprime=0; iprime<(Hprime-Qprime_CQI)*6; iprime+=6) {
        while (y[j2] != LTE_NULL) j2+=6;

        y[j2]   = ulsch->e[iprime];
        y[1+j2] = ulsch->e[1+iprime];
        y[2+j2] = ulsch->e[2+iprime];
        y[3+j2] = ulsch->e[3+iprime];
        y[4+j2] = ulsch->e[4+iprime];
        y[5+j2] = ulsch->e[5+iprime];
        j2+=6;
      }
778

779
      break;
780
  }
781

782
783
784
785
786
787
788
789
  // HARQ-ACK Bits (Note these overwrite some bits)

  if (frame_parms->Ncp == 0)
    columnset = cs_ack_normal;
  else
    columnset = cs_ack_extended;

  j=0;
790
791

  for (i=0; i<Qprime_ACK; i++) {
792
    r = Rmux_prime - 1 - (i>>2);
793
794

    for (q=0; q<Q_m; q++) {
795
796
      y[q+(Q_m*((r*Cmux) + columnset[j]))]  = ulsch->q_ACK[(q+(Q_m*i))%len_ACK];
#ifdef DEBUG_ULSCH_CODING
797
      printf("ulsch_coding.c: ACK %d => y[%d]=%d (i %d, r*Cmux %d, columnset %d)\n",q+(Q_m*i),
798
799
             q+(Q_m*((r*Cmux) + columnset[j])),ulsch->q_ACK[(q+(Q_m*i))%len_ACK],
             i,r*Cmux,columnset[j]);
800
801
#endif
    }
802

803
804
805
806
807
    j=(j+3)&3;
  }

  // write out buffer
  j=0;
808

809
  switch (Q_m) {
810
811
812
813
814
815
816
    case 2:
      for (i=0; i<Cmux; i++)
        for (r=0; r<Rmux_prime; r++) {
          yptr=&y[((r*Cmux)+i)<<1];
          ulsch->h[j++] = *yptr++;
          ulsch->h[j++] = *yptr++;
        }
817

818
      break;
819

820
821
822
823
824
825
826
827
828
    case 4:
      for (i=0; i<Cmux; i++)
        for (r=0; r<Rmux_prime; r++) {
          yptr = &y[((r*Cmux)+i)<<2];
          ulsch->h[j++] = *yptr++;
          ulsch->h[j++] = *yptr++;
          ulsch->h[j++] = *yptr++;
          ulsch->h[j++] = *yptr++;
        }
829

830
      break;
831

832
833
834
835
836
837
838
839
840
841
842
    case 6:
      for (i=0; i<Cmux; i++)
        for (r=0; r<Rmux_prime; r++) {
          yptr = &y[((r*Cmux)+i)*6];
          ulsch->h[j++] = *yptr++;
          ulsch->h[j++] = *yptr++;
          ulsch->h[j++] = *yptr++;
          ulsch->h[j++] = *yptr++;
          ulsch->h[j++] = *yptr++;
          ulsch->h[j++] = *yptr++;
        }
843

844
845
846
847
      break;

    default:
      break;
848
  }
849

850
851
852
  stop_meas(m_stats);

  if (j!=(H+Q_RI)) {
853
    LOG_E(PHY,"Error in output buffer length (j %d, H+Q_RI %d)\n",j,H+Q_RI);
gauthier's avatar
gauthier committed
854
    VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
855
856
857
    return(-1);
  }

gauthier's avatar
gauthier committed
858
  VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
859
860
861
  return(0);
}