ulsch_decoding.c 40.1 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
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
 * except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *      http://www.openairinterface.org/?page_id=698
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 *-------------------------------------------------------------------------------
 * For more information about the OpenAirInterface (OAI) Software Alliance:
 *      contact@openairinterface.org
 */
21 22 23 24 25 26 27 28 29 30 31 32

/*! \file PHY/LTE_TRANSPORT/ulsch_decoding.c
* \brief Top-level routines for decoding  the ULSCH transport channel from 36.212 V8.6 2009-03
* \author R. Knopp
* \date 2011
* \version 0.1
* \company Eurecom
* \email: knopp@eurecom.fr
* \note
* \warning
*/

33

Wang Tsu-Han's avatar
Wang Tsu-Han committed
34
#include <syscall.h>
35 36 37 38 39 40 41
#include "PHY/defs_eNB.h"
#include "PHY/phy_extern.h"
#include "PHY/CODING/coding_extern.h"
#include "SCHED/sched_eNB.h"
#include "LAYER2/MAC/mac.h"
#include "RRC/LTE/rrc_extern.h"
#include "PHY_INTERFACE/phy_interface.h"
42 43 44

#include "UTIL/LOG/vcd_signal_dumper.h"
//#define DEBUG_ULSCH_DECODING
45
#include "targets/RT/USER/rt_wrapper.h"
46
#include "transport_proto.h"
47

Wang Tsu-Han's avatar
Wang Tsu-Han committed
48
extern int codingw;
49

50 51
void free_eNB_ulsch(LTE_eNB_ULSCH_t *ulsch)
{
52 53 54 55

  int i,r;

  if (ulsch) {
56
    for (i=0; i<8; i++) {
57
      if (ulsch->harq_processes[i]) {
58 59 60 61 62 63
        if (ulsch->harq_processes[i]->b) {
          free16(ulsch->harq_processes[i]->b,MAX_ULSCH_PAYLOAD_BYTES);
          ulsch->harq_processes[i]->b = NULL;
        }

        for (r=0; r<MAX_NUM_ULSCH_SEGMENTS; r++) {
64 65 66
          free16(ulsch->harq_processes[i]->c[r],((r==0)?8:0) + 768);
          ulsch->harq_processes[i]->c[r] = NULL;
        }
67 68 69 70 71 72 73 74 75

        for (r=0; r<MAX_NUM_ULSCH_SEGMENTS; r++)
          if (ulsch->harq_processes[i]->d[r]) {
            free16(ulsch->harq_processes[i]->d[r],((3*8*6144)+12+96)*sizeof(short));
            ulsch->harq_processes[i]->d[r] = NULL;
          }

        free16(ulsch->harq_processes[i],sizeof(LTE_UL_eNB_HARQ_t));
        ulsch->harq_processes[i] = NULL;
76 77
      }
    }
78 79 80

    free16(ulsch,sizeof(LTE_eNB_ULSCH_t));
    ulsch = NULL;
81 82 83
  }
}

84
LTE_eNB_ULSCH_t *new_eNB_ulsch(uint8_t max_turbo_iterations,uint8_t N_RB_UL, uint8_t abstraction_flag)
85
{
86 87 88 89

  LTE_eNB_ULSCH_t *ulsch;
  uint8_t exit_flag = 0,i,r;
  unsigned char bw_scaling =1;
90 91 92

  switch (N_RB_UL) {
  case 6:
93 94
    bw_scaling =16;
    break;
95

96 97 98
  case 25:
    bw_scaling =4;
    break;
99 100

  case 50:
101 102
    bw_scaling =2;
    break;
103

104 105 106 107
  default:
    bw_scaling =1;
    break;
  }
108

109
  ulsch = (LTE_eNB_ULSCH_t *)malloc16(sizeof(LTE_eNB_ULSCH_t));
110

111
  if (ulsch) {
112
    memset(ulsch,0,sizeof(LTE_eNB_ULSCH_t));
113
    ulsch->max_turbo_iterations = max_turbo_iterations;
114
    ulsch->Mlimit = 4;
115

116
    for (i=0; i<8; i++) {
Raymond Knopp's avatar
Raymond Knopp committed
117
      //      printf("new_ue_ulsch: Harq process %d\n",i);
118
      ulsch->harq_processes[i] = (LTE_UL_eNB_HARQ_t *)malloc16(sizeof(LTE_UL_eNB_HARQ_t));
119

120
      if (ulsch->harq_processes[i]) {
121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146
        memset(ulsch->harq_processes[i],0,sizeof(LTE_UL_eNB_HARQ_t));
        ulsch->harq_processes[i]->b = (uint8_t*)malloc16(MAX_ULSCH_PAYLOAD_BYTES/bw_scaling);

        if (ulsch->harq_processes[i]->b)
          memset(ulsch->harq_processes[i]->b,0,MAX_ULSCH_PAYLOAD_BYTES/bw_scaling);
        else
          exit_flag=3;

        if (abstraction_flag==0) {
          for (r=0; r<MAX_NUM_ULSCH_SEGMENTS/bw_scaling; r++) {
            ulsch->harq_processes[i]->c[r] = (uint8_t*)malloc16(((r==0)?8:0) + 3+768);
            if (ulsch->harq_processes[i]->c[r])
              memset(ulsch->harq_processes[i]->c[r],0,((r==0)?8:0) + 3+768);
            else
              exit_flag=2;

            ulsch->harq_processes[i]->d[r] = (short*)malloc16(((3*8*6144)+12+96)*sizeof(short));

            if (ulsch->harq_processes[i]->d[r])
              memset(ulsch->harq_processes[i]->d[r],0,((3*8*6144)+12+96)*sizeof(short));
            else
              exit_flag=2;
          }
        }
      } else {
        exit_flag=1;
147 148 149 150 151 152
      }
    }

    if (exit_flag==0)
      return(ulsch);
  }
153

154 155 156 157 158 159
  LOG_E(PHY,"new_ue_ulsch: exit_flag = %d\n",exit_flag);
  free_eNB_ulsch(ulsch);

  return(NULL);
}

160
void clean_eNb_ulsch(LTE_eNB_ULSCH_t *ulsch)
161
{
162 163 164 165 166 167

  unsigned char i;

  //ulsch = (LTE_eNB_ULSCH_t *)malloc16(sizeof(LTE_eNB_ULSCH_t));
  if (ulsch) {
    ulsch->rnti = 0;
168
    ulsch->harq_mask = 0;
169

170
    for (i=0; i<8; i++) {
171
      if (ulsch->harq_processes[i]) {
172 173 174 175 176
        //    ulsch->harq_processes[i]->Ndi = 0;
        ulsch->harq_processes[i]->status = 0;
        //ulsch->harq_processes[i]->phich_active = 0; //this will be done later after transmission of PHICH
        ulsch->harq_processes[i]->phich_ACK = 0;
        ulsch->harq_processes[i]->round = 0;
177 178 179 180 181 182 183 184 185 186 187 188 189 190
        ulsch->harq_processes[i]->rar_alloc = 0;
        ulsch->harq_processes[i]->first_rb = 0;
        ulsch->harq_processes[i]->nb_rb = 0;
        ulsch->harq_processes[i]->TBS = 0;
        ulsch->harq_processes[i]->Or1 = 0;
        ulsch->harq_processes[i]->Or2 = 0;
        for ( int j = 0; j < 2; j++ ) {
          ulsch->harq_processes[i]->o_RI[j] = 0;
        }
        ulsch->harq_processes[i]->O_ACK = 0;
        ulsch->harq_processes[i]->srs_active = 0;
        ulsch->harq_processes[i]->rvidx = 0;
        ulsch->harq_processes[i]->Msc_initial = 0;
        ulsch->harq_processes[i]->Nsymb_initial = 0;
191 192
      }
    }
193 194 195 196
    ulsch->beta_offset_cqi_times8 = 0;
    ulsch->beta_offset_ri_times8 = 0;
    ulsch->beta_offset_harqack_times8 = 0;
    ulsch->Msg3_active = 0;
197 198 199 200
  }
}


201 202
uint8_t extract_cqi_crc(uint8_t *cqi,uint8_t CQI_LENGTH)
{
203 204 205 206

  uint8_t crc;

  crc = cqi[CQI_LENGTH>>3];
Raymond Knopp's avatar
Raymond Knopp committed
207
  //  printf("crc1: %x, shift %d\n",crc,CQI_LENGTH&0x7);
208 209 210
  crc = (crc<<(CQI_LENGTH&0x7));
  // clear crc bits
  //  ((char *)cqi)[CQI_LENGTH>>3] &= 0xff>>(8-(CQI_LENGTH&0x7));
Raymond Knopp's avatar
Raymond Knopp committed
211
  //  printf("crc2: %x, cqi0 %x\n",crc,cqi[1+(CQI_LENGTH>>3)]);
212 213 214 215 216 217 218 219 220 221 222
  crc |= (cqi[1+(CQI_LENGTH>>3)])>>(8-(CQI_LENGTH&0x7));
  // clear crc bits
  //(((char *)cqi)[1+(CQI_LENGTH>>3)]) = 0;

  //  printf("crc : %x\n",crc);
  return(crc);

}



223 224 225 226 227 228 229 230

int ulsch_decoding_data_2thread0(td_params* tdp) {

  PHY_VARS_eNB *eNB = tdp->eNB;
  int UE_id         = tdp->UE_id;
  int harq_pid      = tdp->harq_pid;
  int llr8_flag     = tdp->llr8_flag;

231
  unsigned int r,r_offset=0,Kr,Kr_bytes;
232 233 234 235 236 237
  uint8_t crc_type;
  int offset = 0;
  int ret = 1;
  int16_t dummy_w[MAX_NUM_ULSCH_SEGMENTS][3*(6144+64)];
  LTE_eNB_ULSCH_t *ulsch = eNB->ulsch[UE_id];
  LTE_UL_eNB_HARQ_t *ulsch_harq = ulsch->harq_processes[harq_pid];
238
  int Q_m = ulsch_harq->Qm;
239
  int G = ulsch_harq->G;
Raymond Knopp's avatar
Raymond Knopp committed
240
  uint32_t E=0;
241 242
  uint32_t Gp,GpmodC,Nl=1;
  uint32_t C = ulsch_harq->C;
243
  decoder_if_t *tc;
244 245

  if (llr8_flag == 0)
246
    tc = decoder16;
247
  else
248
    tc = decoder8;
249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273



  // go through first half of segments to get r_offset
  for (r=0; r<(ulsch_harq->C/2); r++) {

    // Get Turbo interleaver parameters
    if (r<ulsch_harq->Cminus)
      Kr = ulsch_harq->Kminus;
    else
      Kr = ulsch_harq->Kplus;

    Kr_bytes = Kr>>3;

    // This is stolen from rate-matching algorithm to get the value of E
    
    Gp = G/Nl/Q_m;
    GpmodC = Gp%C;
    
    if (r < (C-(GpmodC)))
      E = Nl*Q_m * (Gp/C);
    else
      E = Nl*Q_m * ((GpmodC==0?0:1) + (Gp/C));
    
    r_offset += E;
274 275 276 277 278 279

    if (r==0) {
      offset = Kr_bytes - (ulsch_harq->F>>3) - ((ulsch_harq->C>1)?3:0);
    } else {
      offset += (Kr_bytes- ((ulsch_harq->C>1)?3:0));
    }
280 281 282
  }

  // go through second half of segments
283
  for (; r<(ulsch_harq->C); r++) {
284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320


    //    printf("before subblock deinterleaving c[%d] = %p\n",r,ulsch_harq->c[r]);
    // Get Turbo interleaver parameters
    if (r<ulsch_harq->Cminus)
      Kr = ulsch_harq->Kminus;
    else
      Kr = ulsch_harq->Kplus;

    Kr_bytes = Kr>>3;

    memset(&dummy_w[r][0],0,3*(6144+64)*sizeof(short));
    ulsch_harq->RTC[r] = generate_dummy_w(4+(Kr_bytes*8),
                                          (uint8_t*)&dummy_w[r][0],
                                          (r==0) ? ulsch_harq->F : 0);

#ifdef DEBUG_ULSCH_DECODING
    printf("Rate Matching Segment %d (coded bits (G) %d,unpunctured/repeated bits %d, Q_m %d, nb_rb %d, Nl %d)...\n",
        r, G,
        Kr*3,
        Q_m,
        nb_rb,
        ulsch_harq->Nl);
#endif


    if (lte_rate_matching_turbo_rx(ulsch_harq->RTC[r],
                                   G,
                                   ulsch_harq->w[r],
                                   (uint8_t*) &dummy_w[r][0],
                                   ulsch_harq->e+r_offset,
                                   ulsch_harq->C,
                                   NSOFT,
                                   0,   //Uplink
                                   1,
                                   ulsch_harq->rvidx,
                                   (ulsch_harq->round==0)?1:0,  // clear
321
                                   ulsch_harq->Qm,
322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339
                                   1,
                                   r,
                                   &E)==-1) {
      LOG_E(PHY,"ulsch_decoding.c: Problem in rate matching\n");
      return(-1);
    }

    r_offset += E;

    sub_block_deinterleaving_turbo(4+Kr,
                                   &ulsch_harq->d[r][96],
                                   ulsch_harq->w[r]);

    if (ulsch_harq->C == 1)
      crc_type = CRC24_A;
    else
      crc_type = CRC24_B;
    
Raymond Knopp's avatar
Raymond Knopp committed
340
   
341
    ret = tc(&ulsch_harq->d[r][96],
342
             NULL,
343
	     ulsch_harq->c[r],
344
             NULL,
345 346 347 348 349 350 351 352 353 354 355 356
	     Kr,
	     ulsch->max_turbo_iterations,//MAX_TURBO_ITERATIONS,
	     crc_type,
	     (r==0) ? ulsch_harq->F : 0,
	     &eNB->ulsch_tc_init_stats,
	     &eNB->ulsch_tc_alpha_stats,
	     &eNB->ulsch_tc_beta_stats,
	     &eNB->ulsch_tc_gamma_stats,
	     &eNB->ulsch_tc_ext_stats,
	     &eNB->ulsch_tc_intl1_stats,
	     &eNB->ulsch_tc_intl2_stats);
    
357
    // Reassembly of Transport block here
358 359 360 361 362 363 364 365 366

    if (ret != (1+ulsch->max_turbo_iterations)) {
      if (r<ulsch_harq->Cminus)
	Kr = ulsch_harq->Kminus;
      else
	Kr = ulsch_harq->Kplus;
      
      Kr_bytes = Kr>>3;
      
367 368 369 370 371
      memcpy(ulsch_harq->b+offset,
	     ulsch_harq->c[r],
	     Kr_bytes - ((ulsch_harq->C>1)?3:0));
      offset += (Kr_bytes- ((ulsch_harq->C>1)?3:0));
      
372 373 374 375 376 377 378 379 380 381 382 383 384 385
      
    } else {
      break;
    }
    
  }

  return(ret);
}

extern int oai_exit;
void *td_thread(void *param) {
  PHY_VARS_eNB *eNB = ((td_params*)param)->eNB;
  eNB_proc_t *proc  = &eNB->proc;
Wang Tsu-Han's avatar
Wang Tsu-Han committed
386 387
  cpu_set_t cpuset;
  CPU_ZERO(&cpuset);
388
  
Wang Tsu-Han's avatar
Wang Tsu-Han committed
389
  thread_top_init("td_thread",1,200000,250000,500000);
390 391
  pthread_setname_np( pthread_self(),"td processing");
  LOG_I(PHY,"thread td created id=%ld\n", syscall(__NR_gettid));
392
  //wait_sync("td_thread");
393 394 395 396

  while (!oai_exit) {

    if (wait_on_condition(&proc->mutex_td,&proc->cond_td,&proc->instance_cnt_td,"td thread")<0) break;  
Wang Tsu-Han's avatar
Wang Tsu-Han committed
397
    if(oai_exit) break;
398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415

    ((td_params*)param)->ret = ulsch_decoding_data_2thread0((td_params*)param);

    if (release_thread(&proc->mutex_td,&proc->instance_cnt_td,"td thread")<0) break;

    if (pthread_cond_signal(&proc->cond_td) != 0) {
      printf("[eNB] ERROR pthread_cond_signal for td thread exit\n");
      exit_fun( "ERROR pthread_cond_signal" );
      return(NULL);
    }
  }

  return(NULL);
}

int ulsch_decoding_data_2thread(PHY_VARS_eNB *eNB,int UE_id,int harq_pid,int llr8_flag) {

  eNB_proc_t *proc = &eNB->proc;
416
  unsigned int r,r_offset=0,Kr,Kr_bytes;
417 418 419 420 421 422
  uint8_t crc_type;
  int offset = 0;
  int ret = 1;
  int16_t dummy_w[MAX_NUM_ULSCH_SEGMENTS][3*(6144+64)];
  LTE_eNB_ULSCH_t *ulsch = eNB->ulsch[UE_id];
  LTE_UL_eNB_HARQ_t *ulsch_harq = ulsch->harq_processes[harq_pid];
423

424
  int G = ulsch_harq->G;
425
  unsigned int E;
426
  int Cby2;
427
  decoder_if_t *tc;
428 429 430 431 432 433 434 435

  struct timespec wait;

  wait.tv_sec=0;
  wait.tv_nsec=5000000L;


  if (llr8_flag == 0)
436
    tc = decoder16;
437
  else
438
    tc = decoder8;
439

440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467
  if (ulsch_harq->C>1) { // wakeup worker if more than 1 segment
    if (pthread_mutex_timedlock(&proc->mutex_td,&wait) != 0) {
      printf("[eNB] ERROR pthread_mutex_lock for TD thread (IC %d)\n", proc->instance_cnt_td);
      exit_fun( "error locking mutex_fep" );
      return -1;
    }
    
    if (proc->instance_cnt_td==0) {
      printf("[eNB] TD thread busy\n");
      exit_fun("TD thread busy");
      pthread_mutex_unlock( &proc->mutex_td );
      return -1;
    }
    
    ++proc->instance_cnt_td;
    
    proc->tdp.eNB       = eNB;
    proc->tdp.UE_id     = UE_id;
    proc->tdp.harq_pid  = harq_pid;
    proc->tdp.llr8_flag = llr8_flag;
    
    
    // wakeup worker to do second half segments 
    if (pthread_cond_signal(&proc->cond_td) != 0) {
      printf("[eNB] ERROR pthread_cond_signal for td thread exit\n");
      exit_fun( "ERROR pthread_cond_signal" );
      return (1+ulsch->max_turbo_iterations);
    }
468 469

    pthread_mutex_unlock( &proc->mutex_td );
470
    Cby2 = ulsch_harq->C/2;
471
  }
472 473
  else {
    Cby2 = 1;
474 475 476
  }

  // go through first half of segments in main thread
477
  for (r=0; r<Cby2; r++) {
478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514

    //    printf("before subblock deinterleaving c[%d] = %p\n",r,ulsch_harq->c[r]);
    // Get Turbo interleaver parameters
    if (r<ulsch_harq->Cminus)
      Kr = ulsch_harq->Kminus;
    else
      Kr = ulsch_harq->Kplus;

    Kr_bytes = Kr>>3;

    memset(&dummy_w[r][0],0,3*(6144+64)*sizeof(short));
    ulsch_harq->RTC[r] = generate_dummy_w(4+(Kr_bytes*8),
                                          (uint8_t*)&dummy_w[r][0],
                                          (r==0) ? ulsch_harq->F : 0);

#ifdef DEBUG_ULSCH_DECODING
    printf("Rate Matching Segment %d (coded bits (G) %d,unpunctured/repeated bits %d, Q_m %d, nb_rb %d, Nl %d)...\n",
        r, G,
        Kr*3,
        Q_m,
        nb_rb,
        ulsch_harq->Nl);
#endif

    start_meas(&eNB->ulsch_rate_unmatching_stats);

    if (lte_rate_matching_turbo_rx(ulsch_harq->RTC[r],
                                   G,
                                   ulsch_harq->w[r],
                                   (uint8_t*) &dummy_w[r][0],
                                   ulsch_harq->e+r_offset,
                                   ulsch_harq->C,
                                   NSOFT,
                                   0,   //Uplink
                                   1,
                                   ulsch_harq->rvidx,
                                   (ulsch_harq->round==0)?1:0,  // clear
515
                                   ulsch_harq->Qm,
516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539
                                   1,
                                   r,
                                   &E)==-1) {
      LOG_E(PHY,"ulsch_decoding.c: Problem in rate matching\n");
      return(-1);
    }

    stop_meas(&eNB->ulsch_rate_unmatching_stats);
    r_offset += E;

    start_meas(&eNB->ulsch_deinterleaving_stats);
    sub_block_deinterleaving_turbo(4+Kr,
                                   &ulsch_harq->d[r][96],
                                   ulsch_harq->w[r]);
    stop_meas(&eNB->ulsch_deinterleaving_stats);

    if (ulsch_harq->C == 1)
      crc_type = CRC24_A;
    else
      crc_type = CRC24_B;

    start_meas(&eNB->ulsch_turbo_decoding_stats);
    
    ret = tc(&ulsch_harq->d[r][96],
540
             NULL,
541
	     ulsch_harq->c[r],
542
             NULL,
543 544 545 546 547 548 549 550 551 552 553
	     Kr,
	     ulsch->max_turbo_iterations,//MAX_TURBO_ITERATIONS,
	     crc_type,
	     (r==0) ? ulsch_harq->F : 0,
	     &eNB->ulsch_tc_init_stats,
	     &eNB->ulsch_tc_alpha_stats,
	     &eNB->ulsch_tc_beta_stats,
	     &eNB->ulsch_tc_gamma_stats,
	     &eNB->ulsch_tc_ext_stats,
	     &eNB->ulsch_tc_intl1_stats,
	     &eNB->ulsch_tc_intl2_stats);
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
  // Reassembly of Transport block here

    if (ret != (1+ulsch->max_turbo_iterations)) {
      if (r<ulsch_harq->Cminus)
	Kr = ulsch_harq->Kminus;
      else
	Kr = ulsch_harq->Kplus;
      
      Kr_bytes = Kr>>3;
      
      if (r==0) {
	memcpy(ulsch_harq->b,
	       &ulsch_harq->c[0][(ulsch_harq->F>>3)],
	       Kr_bytes - (ulsch_harq->F>>3) - ((ulsch_harq->C>1)?3:0));
	offset = Kr_bytes - (ulsch_harq->F>>3) - ((ulsch_harq->C>1)?3:0);
      } else {
	memcpy(ulsch_harq->b+offset,
	       ulsch_harq->c[r],
	       Kr_bytes - ((ulsch_harq->C>1)?3:0));
	offset += (Kr_bytes- ((ulsch_harq->C>1)?3:0));
      }
      
    } else {
      break;
    }
    stop_meas(&eNB->ulsch_turbo_decoding_stats);    
581
  //printf("/////////////////////////////////////////**************************loop for %d time in ulsch_decoding main\n",r);
582 583 584 585 586 587 588 589 590 591 592
  }

   // wait for worker to finish

  wait_on_busy_condition(&proc->mutex_td,&proc->cond_td,&proc->instance_cnt_td,"td thread");  

  return( (ret>proc->tdp.ret) ? ret : proc->tdp.ret );
}

int ulsch_decoding_data(PHY_VARS_eNB *eNB,int UE_id,int harq_pid,int llr8_flag) {

593
  unsigned int r,r_offset=0,Kr,Kr_bytes;
594 595 596 597 598 599
  uint8_t crc_type;
  int offset = 0;
  int ret = 1;
  int16_t dummy_w[MAX_NUM_ULSCH_SEGMENTS][3*(6144+64)];
  LTE_eNB_ULSCH_t *ulsch = eNB->ulsch[UE_id];
  LTE_UL_eNB_HARQ_t *ulsch_harq = ulsch->harq_processes[harq_pid];
600

601
  int G = ulsch_harq->G;
602
  unsigned int E;
603
  decoder_if_t *tc;
604 605

  if (llr8_flag == 0)
606
    tc = *decoder16;
607
  else
608
    tc = *decoder8;
609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636


  for (r=0; r<ulsch_harq->C; r++) {

    //    printf("before subblock deinterleaving c[%d] = %p\n",r,ulsch_harq->c[r]);
    // Get Turbo interleaver parameters
    if (r<ulsch_harq->Cminus)
      Kr = ulsch_harq->Kminus;
    else
      Kr = ulsch_harq->Kplus;

    Kr_bytes = Kr>>3;

    memset(&dummy_w[r][0],0,3*(6144+64)*sizeof(short));
    ulsch_harq->RTC[r] = generate_dummy_w(4+(Kr_bytes*8),
                                          (uint8_t*)&dummy_w[r][0],
                                          (r==0) ? ulsch_harq->F : 0);

#ifdef DEBUG_ULSCH_DECODING
    printf("Rate Matching Segment %d (coded bits (G) %d,unpunctured/repeated bits %d, Q_m %d, nb_rb %d, Nl %d)...\n",
        r, G,
        Kr*3,
        Q_m,
        nb_rb,
        ulsch_harq->Nl);
#endif

    start_meas(&eNB->ulsch_rate_unmatching_stats);
Raymond Knopp's avatar
Raymond Knopp committed
637

638 639 640 641 642 643 644 645 646 647 648
    if (lte_rate_matching_turbo_rx(ulsch_harq->RTC[r],
                                   G,
                                   ulsch_harq->w[r],
                                   (uint8_t*) &dummy_w[r][0],
                                   ulsch_harq->e+r_offset,
                                   ulsch_harq->C,
                                   NSOFT,
                                   0,   //Uplink
                                   1,
                                   ulsch_harq->rvidx,
                                   (ulsch_harq->round==0)?1:0,  // clear
649
                                   ulsch_harq->Qm,
650 651 652 653 654 655
                                   1,
                                   r,
                                   &E)==-1) {
      LOG_E(PHY,"ulsch_decoding.c: Problem in rate matching\n");
      return(-1);
    }
Raymond Knopp's avatar
Raymond Knopp committed
656

657 658
    stop_meas(&eNB->ulsch_rate_unmatching_stats);
    r_offset += E;
Raymond Knopp's avatar
Raymond Knopp committed
659

660 661 662 663 664
    start_meas(&eNB->ulsch_deinterleaving_stats);
    sub_block_deinterleaving_turbo(4+Kr,
                                   &ulsch_harq->d[r][96],
                                   ulsch_harq->w[r]);
    stop_meas(&eNB->ulsch_deinterleaving_stats);
Raymond Knopp's avatar
Raymond Knopp committed
665

666 667 668 669
    if (ulsch_harq->C == 1)
      crc_type = CRC24_A;
    else
      crc_type = CRC24_B;
Raymond Knopp's avatar
Raymond Knopp committed
670 671
  
 
672 673
    start_meas(&eNB->ulsch_turbo_decoding_stats);
    
Raymond Knopp's avatar
Raymond Knopp committed
674
    ret = tc(&ulsch_harq->d[r][96],
675
             NULL,
Raymond Knopp's avatar
Raymond Knopp committed
676
	     ulsch_harq->c[r],
677
             NULL,
Raymond Knopp's avatar
Raymond Knopp committed
678 679 680 681 682 683 684 685 686 687 688 689 690 691
	     Kr,
	     ulsch->max_turbo_iterations,//MAX_TURBO_ITERATIONS,
	     crc_type,
	     (r==0) ? ulsch_harq->F : 0,
	     &eNB->ulsch_tc_init_stats,
	     &eNB->ulsch_tc_alpha_stats,
	     &eNB->ulsch_tc_beta_stats,
	     &eNB->ulsch_tc_gamma_stats,
	     &eNB->ulsch_tc_ext_stats,
	     &eNB->ulsch_tc_intl1_stats,
	     &eNB->ulsch_tc_intl2_stats);
    
    stop_meas(&eNB->ulsch_turbo_decoding_stats);
    
692
  // Reassembly of Transport block here
Raymond Knopp's avatar
Raymond Knopp committed
693

694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715
    if (ret != (1+ulsch->max_turbo_iterations)) {
      if (r<ulsch_harq->Cminus)
	Kr = ulsch_harq->Kminus;
      else
	Kr = ulsch_harq->Kplus;
      
      Kr_bytes = Kr>>3;
      
      if (r==0) {
	memcpy(ulsch_harq->b,
	       &ulsch_harq->c[0][(ulsch_harq->F>>3)],
	       Kr_bytes - (ulsch_harq->F>>3) - ((ulsch_harq->C>1)?3:0));
	offset = Kr_bytes - (ulsch_harq->F>>3) - ((ulsch_harq->C>1)?3:0);
      } else {
	memcpy(ulsch_harq->b+offset,
	       ulsch_harq->c[r],
	       Kr_bytes - ((ulsch_harq->C>1)?3:0));
	offset += (Kr_bytes- ((ulsch_harq->C>1)?3:0));
      }
      
    } else {
      break;
Raymond Knopp's avatar
Raymond Knopp committed
716
    }
717
    
Raymond Knopp's avatar
Raymond Knopp committed
718 719
  }

720 721 722
  return(ret);
}

723 724 725
int ulsch_decoding_data_all(PHY_VARS_eNB *eNB,int UE_id,int harq_pid,int llr8_flag) 
{
  int ret = 0;
726
  /*if(codingw)
727 728 729
  {
    ret = ulsch_decoding_data_2thread(eNB,UE_id,harq_pid,llr8_flag);
  }
730
  else*/
731 732 733 734 735 736
  {
    ret = ulsch_decoding_data(eNB,UE_id,harq_pid,llr8_flag);
  }
  return ret;
}

737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754
static inline unsigned int lte_gold_unscram(unsigned int *x1, unsigned int *x2, unsigned char reset) __attribute__((always_inline));
static inline unsigned int lte_gold_unscram(unsigned int *x1, unsigned int *x2, unsigned char reset)
{
  int n;

  if (reset) {
    *x1 = 1+ (1<<31);
    *x2=*x2 ^ ((*x2 ^ (*x2>>1) ^ (*x2>>2) ^ (*x2>>3))<<31);

    // skip first 50 double words (1600 bits)
    //      printf("n=0 : x1 %x, x2 %x\n",x1,x2);
    for (n=1; n<50; n++) {
      *x1 = (*x1>>1) ^ (*x1>>4);
      *x1 = *x1 ^ (*x1<<31) ^ (*x1<<28);
      *x2 = (*x2>>1) ^ (*x2>>2) ^ (*x2>>3) ^ (*x2>>4);
      *x2 = *x2 ^ (*x2<<31) ^ (*x2<<30) ^ (*x2<<29) ^ (*x2<<28);
    }
  }
755

756 757 758 759 760 761
  *x1 = (*x1>>1) ^ (*x1>>4);
  *x1 = *x1 ^ (*x1<<31) ^ (*x1<<28);
  *x2 = (*x2>>1) ^ (*x2>>2) ^ (*x2>>3) ^ (*x2>>4);
  *x2 = *x2 ^ (*x2<<31) ^ (*x2<<30) ^ (*x2<<29) ^ (*x2<<28);
  return(*x1^*x2);
  //  printf("n=%d : c %x\n",n,x1^x2);
Raymond Knopp's avatar
Raymond Knopp committed
762 763 764

}
  
765
unsigned int  ulsch_decoding(PHY_VARS_eNB *eNB,eNB_rxtx_proc_t *proc,
766 767 768 769 770
                             uint8_t UE_id,
                             uint8_t control_only_flag,
                             uint8_t Nbundled,
                             uint8_t llr8_flag)
{
771 772


773 774 775
  int16_t *ulsch_llr = eNB->pusch_vars[UE_id]->llr;
  LTE_DL_FRAME_PARMS *frame_parms = &eNB->frame_parms;
  LTE_eNB_ULSCH_t *ulsch = eNB->ulsch[UE_id];
776 777
  uint8_t harq_pid;
  unsigned short nb_rb;
778
  unsigned int A;
779 780 781
  uint8_t Q_m;
  unsigned int i,i2,q,j,j2;
  int iprime;
782 783
  unsigned int ret=0;

784
  //  uint8_t dummy_channel_output[(3*8*block_length)+12];
785
  int r,Kr;
786 787 788 789

  uint8_t *columnset;
  unsigned int sumKr=0;
  unsigned int Qprime,L,G,Q_CQI,Q_RI,H,Hprime,Hpp,Cmux,Rmux_prime,O_RCC;
790
  unsigned int Qprime_ACK,Qprime_RI,len_ACK=0,len_RI=0;
Raymond Knopp's avatar
 
Raymond Knopp committed
791
  //  uint8_t q_ACK[MAX_ACK_PAYLOAD],q_RI[MAX_RI_PAYLOAD];
792 793 794 795 796
  int metric,metric_new;
  uint32_t x1, x2, s=0;
  int16_t ys,c;
  uint32_t wACK_idx;
  uint8_t dummy_w_cc[3*(MAX_CQI_BITS+8+32)];
797
  int16_t y[6*14*1200] __attribute__((aligned(32)));
798 799
  uint8_t ytag[14*1200];
  //  uint8_t ytag2[6*14*1200],*ytag2_ptr;
Cedric Roux's avatar
Cedric Roux committed
800
  int16_t cseq[6*14*1200] __attribute__((aligned(32)));
801
  int off;
802

803
  int frame = proc->frame_rx;
804
  int subframe = proc->subframe_rx;
Raymond Knopp's avatar
Raymond Knopp committed
805
  LTE_UL_eNB_HARQ_t *ulsch_harq;
806

807

808

809
  harq_pid = subframe2harq_pid(frame_parms,frame,subframe);
810 811

  VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_PHY_ENB_ULSCH_DECODING0+harq_pid,1);
812 813 814

  // x1 is set in lte_gold_generic
  x2 = ((uint32_t)ulsch->rnti<<14) + ((uint32_t)subframe<<9) + frame_parms->Nid_cell; //this is c_init in 36.211 Sec 6.3.1
Raymond Knopp's avatar
Raymond Knopp committed
815
  ulsch_harq = ulsch->harq_processes[harq_pid];
816

817 818
  AssertFatal(harq_pid!=255,
              "FATAL ERROR: illegal harq_pid, returning\n");
819

820 821
  AssertFatal(ulsch_harq->Nsymb_pusch != 0,
              "FATAL ERROR: harq_pid %d, Nsymb 0!\n",harq_pid);
822

823

Raymond Knopp's avatar
Raymond Knopp committed
824
  nb_rb = ulsch_harq->nb_rb;
825

Raymond Knopp's avatar
Raymond Knopp committed
826
  A = ulsch_harq->TBS;
827

828

829
  Q_m = ulsch_harq->Qm;
Raymond Knopp's avatar
Raymond Knopp committed
830
  G = nb_rb * (12 * Q_m) * ulsch_harq->Nsymb_pusch;
831 832


Raymond Knopp's avatar
Raymond Knopp committed
833
  //#ifdef DEBUG_ULSCH_DECODING
834 835
  LOG_D(PHY,"[PUSCH %d] Frame %d, Subframe %d: ulsch_decoding (Nid_cell %d, rnti %x, x2 %x): A %d, round %d, RV %d, O_r1 %d, O_RI %d, O_ACK %d, G %d, Q_m %d Nsymb_pusch %d nb_rb %d\n",
      harq_pid,
Raymond Knopp's avatar
Raymond Knopp committed
836 837 838 839 840 841 842 843
	proc->frame_rx,subframe,
	frame_parms->Nid_cell,ulsch->rnti,x2,
	A,
	ulsch_harq->round,
	ulsch_harq->rvidx,
	ulsch_harq->Or1,
	ulsch_harq->O_RI,
	ulsch_harq->O_ACK,
844 845 846 847
	G,
        ulsch_harq->Qm,
        ulsch_harq->Nsymb_pusch,
        nb_rb);
Raymond Knopp's avatar
Raymond Knopp committed
848 849
	
  //#endif
850
  //if (ulsch_harq->round == 0) { // delete for RB shortage pattern
851
    // This is a new packet, so compute quantities regarding segmentation
Raymond Knopp's avatar
Raymond Knopp committed
852
    ulsch_harq->B = A+24;
853
    lte_segmentation(NULL,
854 855 856 857 858 859 860 861
                     NULL,
                     ulsch_harq->B,
                     &ulsch_harq->C,
                     &ulsch_harq->Cplus,
                     &ulsch_harq->Cminus,
                     &ulsch_harq->Kplus,
                     &ulsch_harq->Kminus,
                     &ulsch_harq->F);
862
    //  CLEAR LLR's HERE for first packet in process
863
  //}
864
  //  printf("after segmentation c[%d] = %p\n",0,ulsch_harq->c[0]);
865 866

  sumKr = 0;
867 868

  for (r=0; r<ulsch_harq->C; r++) {
Raymond Knopp's avatar
Raymond Knopp committed
869 870
    if (r<ulsch_harq->Cminus)
      Kr = ulsch_harq->Kminus;
871
    else
Raymond Knopp's avatar
Raymond Knopp committed
872
      Kr = ulsch_harq->Kplus;
873

874 875
    sumKr += Kr;
  }
876

877
  AssertFatal(sumKr>0,
Raymond Knopp's avatar
Raymond Knopp committed
878
	      "[eNB] ulsch_decoding.c: FATAL sumKr is 0! (Nid_cell %d, rnti %x, x2 %x): harq_pid %d round %d, RV %d, O_RI %d, O_ACK %d, G %d, subframe %d\n",
879 880 881 882 883 884 885 886 887
	      frame_parms->Nid_cell,ulsch->rnti,x2,
	      harq_pid,
	      ulsch_harq->round,
	      ulsch_harq->rvidx,
	      ulsch_harq->O_RI,
	      ulsch_harq->O_ACK,
	      G,
	      subframe);

888

889
  // Compute Q_ri
Raymond Knopp's avatar
Raymond Knopp committed
890
  Qprime = ulsch_harq->O_RI*ulsch_harq->Msc_initial*ulsch_harq->Nsymb_initial * ulsch->beta_offset_ri_times8;
891 892 893 894 895 896

  if (Qprime > 0 ) {
    if ((Qprime % (8*sumKr)) > 0)
      Qprime = 1+(Qprime/(8*sumKr));
    else
      Qprime = Qprime/(8*sumKr);
897

898 899 900 901 902 903 904 905 906 907
    if (Qprime > 4*nb_rb * 12)
      Qprime = 4*nb_rb * 12;
  }

  Q_RI = Q_m*Qprime;
  Qprime_RI = Qprime;


  // Compute Q_ack

Raymond Knopp's avatar
Raymond Knopp committed
908
  Qprime = ulsch_harq->O_ACK*ulsch_harq->Msc_initial*ulsch_harq->Nsymb_initial * ulsch->beta_offset_harqack_times8;
909

910 911 912 913 914 915 916 917 918 919 920 921 922
  if (Qprime > 0) {
    if ((Qprime % (8*sumKr)) > 0)
      Qprime = 1+(Qprime/(8*sumKr));
    else
      Qprime = Qprime/(8*sumKr);

    if (Qprime > (4*nb_rb * 12))
      Qprime = 4*nb_rb * 12;
  }

  //  Q_ACK = Qprime * Q_m;
  Qprime_ACK = Qprime;
#ifdef DEBUG_ULSCH_DECODING
Raymond Knopp's avatar
Raymond Knopp committed
923
  printf("ulsch_decoding.c: Qprime_ACK %d, Msc_initial %d, Nsymb_initial %d, sumKr %d\n",
Raymond Knopp's avatar
Raymond Knopp committed
924
      Qprime_ACK,ulsch_harq->Msc_initial,ulsch_harq->Nsymb_initial,sumKr);
925
#endif
926

927
  // Compute Q_cqi
Raymond Knopp's avatar
Raymond Knopp committed
928
  if (ulsch_harq->Or1 < 12)
929
    L=0;
930
  else
931
    L=8;
932

933
  // NOTE: we have to handle the case where we have a very small number of bits (condition on pg. 26 36.212)
Raymond Knopp's avatar
Raymond Knopp committed
934 935
  if (ulsch_harq->Or1 > 0)
    Qprime = (ulsch_harq->Or1 + L) * ulsch_harq->Msc_initial*ulsch_harq->Nsymb_initial * ulsch->beta_offset_cqi_times8;
936 937 938
  else
    Qprime=0;

939
  if (Qprime > 0) {  // check if ceiling is larger than floor in Q' expression
940 941 942 943 944 945
    if ((Qprime % (8*sumKr)) > 0)
      Qprime = 1+(Qprime/(8*sumKr));
    else
      Qprime = Qprime/(8*sumKr);
  }

Raymond Knopp's avatar
Raymond Knopp committed
946
  G = nb_rb * (12 * Q_m) * (ulsch_harq->Nsymb_pusch);
947

948

949 950

  Q_CQI = Q_m * Qprime;
Raymond Knopp's avatar
Raymond Knopp committed
951
#ifdef DEBUG_ULSCH_DECODING
Raymond Knopp's avatar
Raymond Knopp committed
952
  printf("ulsch_decoding: G %d, Q_RI %d, Q_CQI %d (L %d, Or1 %d) O_ACK %d\n",G,Q_RI,Q_CQI,L,ulsch_harq->Or1,ulsch_harq->O_ACK);
Raymond Knopp's avatar
Raymond Knopp committed
953
#endif
954 955

  G = G - Q_RI - Q_CQI;
956
  ulsch_harq->G = G;
957

958 959
  AssertFatal((int)G > 0,
              "FATAL: ulsch_decoding.c G < 0 (%d) : Q_RI %d, Q_CQI %d\n",G,Q_RI,Q_CQI);
960 961 962 963

  H = G + Q_CQI;
  Hprime = H/Q_m;

964

965
  // Demultiplexing/Deinterleaving of PUSCH/ACK/RI/CQI
966
  start_meas(&eNB->ulsch_demultiplexing_stats);
967
  Hpp = Hprime + Qprime_RI;
968

Raymond Knopp's avatar
Raymond Knopp committed
969
  Cmux       = ulsch_harq->Nsymb_pusch;
970
  Rmux_prime = Hpp/Cmux;
971

972 973 974
  // Clear "tag" interleaving matrix to allow for CQI/DATA identification
  memset(ytag,0,Cmux*Rmux_prime);

975

976 977

  i=0;
Raymond Knopp's avatar
 
Raymond Knopp committed
978
  memset(y,LTE_NULL,Q_m*Hpp);
979 980 981 982 983

  // read in buffer and unscramble llrs for everything but placeholder bits
  // llrs stored per symbol correspond to columns of interleaving matrix


984
  s = lte_gold_unscram(&x1, &x2, 1);
985
  i2=0;
986 987

  for (i=0; i<((Hpp*Q_m)>>5); i++) {
988
    /*
989
    for (j=0; j<32; j++) {
990 991
      cseq[i2++] = (int16_t)((((s>>j)&1)<<1)-1);
    }
992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011
    */
#if defined(__x86_64__) || defined(__i386__)
#ifndef __AVX2__
    ((__m128i*)cseq)[i2++] = ((__m128i*)unscrambling_lut)[(s&65535)<<1];
    ((__m128i*)cseq)[i2++] = ((__m128i*)unscrambling_lut)[1+((s&65535)<<1)];
    s>>=16;
    ((__m128i*)cseq)[i2++] = ((__m128i*)unscrambling_lut)[(s&65535)<<1];
    ((__m128i*)cseq)[i2++] = ((__m128i*)unscrambling_lut)[1+((s&65535)<<1)];
#else
    ((__m256i*)cseq)[i2++] = ((__m256i*)unscrambling_lut)[s&65535];
    ((__m256i*)cseq)[i2++] = ((__m256i*)unscrambling_lut)[(s>>16)&65535];
#endif
#elif defined(__arm__)
    ((int16x8_t*)cseq)[i2++] = ((int16x8_t*)unscrambling_lut)[(s&65535)<<1];
    ((int16x8_t*)cseq)[i2++] = ((int16x8_t*)unscrambling_lut)[1+((s&65535)<<1)];
    s>>=16;
    ((int16x8_t*)cseq)[i2++] = ((int16x8_t*)unscrambling_lut)[(s&65535)<<1];
    ((int16x8_t*)cseq)[i2++] = ((int16x8_t*)unscrambling_lut)[1+((s&65535)<<1)];
#endif
    s = lte_gold_unscram(&x1, &x2, 0);
1012 1013
  }

1014

1015 1016
  //  printf("after unscrambling c[%d] = %p\n",0,ulsch_harq->c[0]);

1017 1018 1019 1020
  if (frame_parms->Ncp == 0)
    columnset = cs_ri_normal;
  else
    columnset = cs_ri_extended;
1021 1022 1023 1024

  j=0;

  for (i=0; i<Qprime_RI; i++) {
1025
    r = Rmux_prime - 1 - (i>>2);
1026 1027
    /*
    for (q=0;q<Q_m;q++)
1028 1029
      ytag2[q+(Q_m*((r*Cmux) + columnset[j]))]  = q_RI[(q+(Q_m*i))%len_RI];
    */
1030
    off =((Rmux_prime*Q_m*columnset[j])+(r*Q_m));
1031
    cseq[off+1] = cseq[off];  // PUSCH_y
1032 1033

    for (q=2; q<Q_m; q++)
1034 1035 1036 1037 1038 1039
      cseq[off+q] = -1;    // PUSCH_x

    j=(j+3)&3;

  }

1040 1041
  //  printf("after RI c[%d] = %p\n",0,ulsch_harq->c[0]);

1042 1043 1044 1045 1046 1047 1048
  // HARQ-ACK Bits (Note these overwrite some bits)
  if (frame_parms->Ncp == 0)
    columnset = cs_ack_normal;
  else
    columnset = cs_ack_extended;

  j=0;
1049 1050

  for (i=0; i<Qprime_ACK; i++) {
1051
    r = Rmux_prime - 1 - (i>>2);
1052 1053
    off =((Rmux_prime*Q_m*columnset[j])+(r*Q_m));

Raymond Knopp's avatar
Raymond Knopp committed
1054
    if (ulsch_harq->O_ACK == 1) {
1055
      if (ulsch->bundling==0)
1056 1057 1058 1059 1060 1061 1062
        cseq[off+1] = cseq[off];  // PUSCH_y

      for (q=2; q<Q_m; q++)
        cseq[off+q] = -1;    // PUSCH_x
    } else if (ulsch_harq->O_ACK == 2) {
      for (q=2; q<Q_m; q++)
        cseq[off+q] = -1;    // PUSCH_x
1063
    }
1064

1065
#ifdef DEBUG_ULSCH_DECODING
Raymond Knopp's avatar
Raymond Knopp committed
1066
    printf("ulsch_decoding.c: ACK i %d, r %d, j %d, ColumnSet[j] %d\n",i,r,j,columnset[j]);
1067 1068 1069 1070 1071 1072 1073
#endif
    j=(j+3)&3;
  }



  i=0;
1074

1075 1076
  switch (Q_m) {
  case 2:
1077
    for (j=0; j<Cmux; j++) {
1078
      i2=j<<1;
1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089

      for (r=0; r<Rmux_prime; r++) {
        c = cseq[i];
        //  printf("ulsch %d: %d * ",i,c);
        y[i2++] = c*ulsch_llr[i++];
        //  printf("%d\n",ulsch_llr[i-1]);
        c = cseq[i];
        //  printf("ulsch %d: %d * ",i,c);
        y[i2] = c*ulsch_llr[i++];
        //  printf("%d\n",ulsch_llr[i-1]);
        i2=(i2+(Cmux<<1)-1);
1090 1091
      }
    }
1092

1093
    break;
1094

1095
  case 4:
1096
    for (j=0; j<Cmux; j++) {
1097
      i2=j<<2;
1098 1099

      for (r=0; r<Rmux_prime; r++) {
1100
	/*
1101 1102 1103 1104 1105 1106 1107 1108 1109
        c = cseq[i];
        y[i2++] = c*ulsch_llr[i++];
        c = cseq[i];
        y[i2++] = c*ulsch_llr[i++];
        c = cseq[i];
        y[i2++] = c*ulsch_llr[i++];
        c = cseq[i];
        y[i2] = c*ulsch_llr[i++];
        i2=(i2+(Cmux<<2)-3);
1110
	*/
1111
	// slightly more optimized version (equivalent to above) for 16QAM to improve computational performance
1112 1113 1114
	*(__m64 *)&y[i2] = _mm_sign_pi16(*(__m64*)&ulsch_llr[i],*(__m64*)&cseq[i]);i+=4;i2+=(Cmux<<2);


1115 1116
      }
    }
1117 1118 1119

    break;

1120
  case 6:
1121
    for (j=0; j<Cmux; j++) {
1122
      i2=j*6;
1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137

      for (r=0; r<Rmux_prime; r++) {
        c = cseq[i];
        y[i2++] = c*ulsch_llr[i++];
        c = cseq[i];
        y[i2++] = c*ulsch_llr[i++];
        c = cseq[i];
        y[i2++] = c*ulsch_llr[i++];
        c = cseq[i];
        y[i2++] = c*ulsch_llr[i++];
        c = cseq[i];
        y[i2++] = c*ulsch_llr[i++];
        c = cseq[i];
        y[i2] = c*ulsch_llr[i++];
        i2=(i2+(Cmux*6)-5);
1138 1139
      }
    }
1140

1141 1142 1143 1144
    break;
  }


1145

1146 1147

  if (i!=(H+Q_RI))
1148
    LOG_D(PHY,"ulsch_decoding.c: Error in input buffer length (j %d, H+Q_RI %d)\n",i,H+Q_RI);
1149 1150 1151 1152 1153 1154 1155 1156 1157 1158

  // HARQ-ACK Bits (LLRs are nulled in overwritten bits after copying HARQ-ACK LLR)

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

  j=0;

Raymond Knopp's avatar
Raymond Knopp committed
1159
  if (ulsch_harq->O_ACK == 1) {
1160 1161 1162 1163
    switch (Q_m) {
    case 2:
      len_ACK = 2;
      break;
1164

1165 1166 1167
    case 4:
      len_ACK = 4;
      break;
1168

1169 1170 1171 1172 1173
    case 6:
      len_ACK = 6;
      break;
    }
  }
1174

Raymond Knopp's avatar
Raymond Knopp committed
1175
  if (ulsch_harq->O_ACK == 2) {
1176 1177 1178 1179
    switch (Q_m) {
    case 2:
      len_ACK = 6;
      break;
1180

1181 1182 1183
    case 4:
      len_ACK = 12;
      break;
1184

1185 1186 1187 1188 1189
    case 6:
      len_ACK = 18;
      break;
    }
  }
1190

Raymond Knopp's avatar
Raymond Knopp committed
1191
  if (ulsch_harq->O_ACK > 2) {
1192
    LOG_E(PHY,"ulsch_decoding: FATAL, ACK cannot be more than 2 bits yet O_ACK:%d SFN/SF:%04d%d UE_id:%d rnti:%x\n",ulsch_harq->O_ACK,proc->frame_rx,proc->subframe_rx,UE_id,ulsch->rnti);
1193 1194 1195
    return(-1);
  }

1196
  for (i=0; i<len_ACK; i++)
Raymond Knopp's avatar
Raymond Knopp committed
1197
    ulsch_harq->q_ACK[i] = 0;
1198 1199


1200
  for (i=0; i<Qprime_ACK; i++) {
1201
    r = Rmux_prime -1 - (i>>2);
1202 1203

    for (q=0; q<Q_m; q++) {
1204
      if (y[q+(Q_m*((r*Cmux) + columnset[j]))]!=0)
1205
        ulsch_harq->q_ACK[(q+(Q_m*i))%len_ACK] += y[q+(Q_m*((r*Cmux) + columnset[j]))];
1206 1207
      y[q+(Q_m*((r*Cmux) + columnset[j]))]=0;  // NULL LLRs in ACK positions
    }
1208

1209 1210 1211
    j=(j+3)&3;
  }

1212
  //  printf("after ACKNAK c[%d] = %p\n",0,ulsch_harq->c[0]);
1213

1214
  // RI BITS
1215

Raymond Knopp's avatar
Raymond Knopp committed
1216
  if (ulsch_harq->O_RI == 1) {
1217 1218
    switch (Q_m) {
    case 2:
1219
      len_RI=2;
1220
      break;
1221

1222 1223 1224
    case 4:
      len_RI=4;
      break;
1225

1226 1227 1228 1229 1230 1231
    case 6:
      len_RI=6;
      break;
    }
  }

Raymond Knopp's avatar
Raymond Knopp committed
1232
  if (ulsch_harq->O_RI > 1) {
1233 1234 1235 1236
    LOG_E(PHY,"ulsch_decoding: FATAL, RI cannot be more than 1 bit yet\n");
    return(-1);
  }

1237
  for (i=0; i<len_RI; i++)
Raymond Knopp's avatar
Raymond Knopp committed
1238
    ulsch_harq->q_RI[i] = 0;
1239

1240 1241 1242 1243
  if (frame_parms->Ncp == 0)
    columnset = cs_ri_normal;
  else
    columnset = cs_ri_extended;
1244 1245 1246 1247

  j=0;

  for (i=0; i<Qprime_RI; i++) {
1248
    r = Rmux_prime -1 - (i>>2);
1249 1250

    for (q=0; q<Q_m; q++)
Raymond Knopp's avatar
Raymond Knopp committed
1251
      ulsch_harq->q_RI[(q+(Q_m*i))%len_RI] += y[q+(Q_m*((r*Cmux) + columnset[j]))];
1252

1253 1254 1255 1256
    ytag[(r*Cmux) + columnset[j]] = LTE_NULL;
    j=(j+3)&3;
  }

1257 1258
  //  printf("after RI2 c[%d] = %p\n",0,ulsch_harq->c[0]);

1259
  // CQI and Data bits
1260 1261 1262
  j=0;
  j2=0;

1263
  //  r=0;
1264 1265 1266 1267 1268 1269
  if (Q_RI>0) {
    for (i=0; i<(Q_CQI/Q_m); i++) {
      
      while (ytag[j]==LTE_NULL) {
	j++;
	j2+=Q_m;
1270
      }
1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283
      
      for (q=0; q<Q_m; q++) {
	//      ys = y[q+(Q_m*((r*Cmux)+j))];
	ys = y[q+j2];
	
	if (ys>127)
	  ulsch_harq->q[q+(Q_m*i)] = 127;
	else if (ys<-128)
	  ulsch_harq->q[q+(Q_m*i)] = -128;
	else
	  ulsch_harq->q[q+(Q_m*i)] = ys;
      }
      
1284 1285
      j2+=Q_m;
    }
1286
    
1287
        
1288 1289
    switch (Q_m) {
    case 2:
1290
      for (iprime=0; iprime<G;) {
1291 1292 1293 1294 1295 1296 1297 1298
	while (ytag[j]==LTE_NULL) {
	  j++;
	  j2+=2;
	}
	
	ulsch_harq->e[iprime++] = y[j2++];
	ulsch_harq->e[iprime++] = y[j2++];
	
1299
      }
1300 1301 1302 1303 1304
      
      
      break;
      
    case 4:
1305
      for (iprime=0; iprime<G;) {
1306 1307 1308 1309 1310 1311 1312 1313 1314 1315
	while (ytag[j]==LTE_NULL) {
	  j++;
	  j2+=4;
	}
	
	ulsch_harq->e[iprime++] = y[j2++];
	ulsch_harq->e[iprime++] = y[j2++];
	ulsch_harq->e[iprime++] = y[j2++];
	ulsch_harq->e[iprime++] = y[j2++];
	
1316
      }
1317 1318 1319 1320
      
      break;
      
    case 6:
1321
      for (iprime=0; iprime<G;) {
1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333
	while (ytag[j]==LTE_NULL) {
	  j++;
	  j2+=6;
	}
	
	ulsch_harq->e[iprime++] = y[j2++];
	ulsch_harq->e[iprime++] = y[j2++];
	ulsch_harq->e[iprime++] = y[j2++];
	ulsch_harq->e[iprime++] = y[j2++];
	ulsch_harq->e[iprime++] = y[j2++];
	ulsch_harq->e[iprime++] = y[j2++];
	
1334
      }
1335 1336
      
      break;
1337

1338
    }
1339 1340 1341
    

  } // Q_RI>0
1342
  else {
1343

1344 1345 1346 1347 1348 1349 1350 1351 1352 1353
    for (i=0; i<(Q_CQI/Q_m); i++) {
      
      for (q=0; q<Q_m; q++) {
	ys = y[q+j2];
	if (ys>127)
	  ulsch_harq->q[q+(Q_m*i)] = 127;
	else if (ys<-128)
	  ulsch_harq->q[q+(Q_m*i)] = -128;
	else
	  ulsch_harq->q[q+(Q_m*i)] = ys;
1354
      }
1355 1356
      
      j2+=Q_m;
1357
    }
1358
    /* To be improved according to alignment of j2
1359
#if defined(__x86_64__)||defined(__i386__)
1360
#ifndef __AVX2__
1361 1362 1363 1364 1365 1366 1367 1368 1369
    for (iprime=0; iprime<G;iprime+=8,j2+=8)
      *((__m128i *)&ulsch_harq->e[iprime]) = *((__m128i *)&y[j2]);
#else
    for (iprime=0; iprime<G;iprime+=16,j2+=16)
      *((__m256i *)&ulsch_harq->e[iprime]) = *((__m256i *)&y[j2]);
#endif
#elif defined(__arm__)
    for (iprime=0; iprime<G;iprime+=8,j2+=8)
      *((int16x8_t *)&ulsch_harq->e[iprime]) = *((int16x8_t *)&y[j2]);
1370 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380 1381 1382 1383 1384
#endif
    */
    int16_t *yp,*ep;
    for (iprime=0,yp=&y[j2],ep=&ulsch_harq->e[0]; 
	 iprime<G;
	 iprime+=8,j2+=8,ep+=8,yp+=8) {
      ep[0] = yp[0];
      ep[1] = yp[1];
      ep[2] = yp[2];
      ep[3] = yp[3];
      ep[4] = yp[4];
      ep[5] = yp[5];
      ep[6] = yp[6];
      ep[7] = yp[7];
    }
1385
  }
1386 1387
    
   
1388 1389
  stop_meas(&eNB->ulsch_demultiplexing_stats);

1390 1391
  //  printf("after ACKNAK2 c[%d] = %p (iprime %d, G %d)\n",0,ulsch_harq->c[0],iprime,G);

1392 1393
  // Do CQI/RI/HARQ-ACK Decoding first and pass to MAC

1394
  // HARQ-ACK
1395
  wACK_idx = (ulsch->bundling==0) ? 4 : ((Nbundled-1)&3);
1396

Raymond Knopp's avatar
Raymond Knopp committed
1397
  if (ulsch_harq->O_ACK == 1) {
1398 1399 1400 1401 1402 1403 1404
    ulsch_harq->q_ACK[0] *= wACK_RX[wACK_idx][0];
    ulsch_harq->q_ACK[0] += (ulsch->bundling==0) ? ulsch_harq->q_ACK[1]*wACK_RX[wACK_idx][0] : ulsch_harq->q_ACK[1]*wACK_RX[wACK_idx][1];

    if (ulsch_harq->q_ACK[0] < 0)
      ulsch_harq->o_ACK[0] = 0;
    else
      ulsch_harq->o_ACK[0] = 1;
1405
  }
1406

Raymond Knopp's avatar
Raymond Knopp committed
1407
  if (ulsch_harq->O_ACK == 2) {
1408 1409 1410
    switch (Q_m) {

    case 2:
Raymond Knopp's avatar
Raymond Knopp committed
1411 1412 1413
      ulsch_harq->q_ACK[0] = ulsch_harq->q_ACK[0]*wACK_RX[wACK_idx][0] + ulsch_harq->q_ACK[3]*wACK_RX[wACK_idx][1];
      ulsch_harq->q_ACK[1] = ulsch_harq->q_ACK[1]*wACK_RX[wACK_idx][0] + ulsch_harq->q_ACK[4]*wACK_RX[wACK_idx][1];
      ulsch_harq->q_ACK[2] = ulsch_harq->q_ACK[2]*wACK_RX[wACK_idx][0] + ulsch_harq->q_ACK[5]*wACK_RX[wACK_idx][1];
1414 1415
      break;
    case 4:
Raymond Knopp's avatar
Raymond Knopp committed
1416 1417 1418
      ulsch_harq->q_ACK[0] = ulsch_harq->q_ACK[0]*wACK_RX[wACK_idx][0] + ulsch_harq->q_ACK[5]*wACK_RX[wACK_idx][1];
      ulsch_harq->q_ACK[1] = ulsch_harq->q_ACK[1]*wACK_RX[wACK_idx][0] + ulsch_harq->q_ACK[8]*wACK_RX[wACK_idx][1];
      ulsch_harq->q_ACK[2] = ulsch_harq->q_ACK[4]*wACK_RX[wACK_idx][0] + ulsch_harq->q_ACK[9]*wACK_RX[wACK_idx][1];
1419 1420
      break;
    case 6:
Raymond Knopp's avatar
Raymond Knopp committed
1421
      ulsch_harq->q_ACK[0] =  ulsch_harq->q_ACK[0]*wACK_RX[wACK_idx][0] + ulsch_harq->q_ACK[7]*wACK_RX[wACK_idx][1];
1422 1423
      ulsch_harq->q_ACK[1] =  ulsch_harq->q_ACK[1]*wACK_RX[wACK_idx][0] + ulsch_harq->q_ACK[12]*wACK_RX[wACK_idx][1];
      ulsch_harq->q_ACK[2] =  ulsch_harq->q_ACK[6]*wACK_RX[wACK_idx][0] + ulsch_harq->q_ACK[13]*wACK_RX[wACK_idx][1];
1424 1425
      break;
    }
1426

Raymond Knopp's avatar
Raymond Knopp committed
1427 1428 1429 1430
    ulsch_harq->o_ACK[0] = 1;
    ulsch_harq->o_ACK[1] = 1;
    metric     = ulsch_harq->q_ACK[0]+ulsch_harq->q_ACK[1]-ulsch_harq->q_ACK[2];
    metric_new = -ulsch_harq->q_ACK[0]+ulsch_harq->q_ACK[1]+ulsch_harq->q_ACK[2];
1431 1432

    if (metric_new > metric) {
Raymond Knopp's avatar
Raymond Knopp committed
1433 1434
      ulsch_harq->o_ACK[0]=0;
      ulsch_harq->o_ACK[1]=1;
1435 1436
      metric = metric_new;
    }
1437

Raymond Knopp's avatar
Raymond Knopp committed
1438
    metric_new = ulsch_harq->q_ACK[0]-ulsch_harq->q_ACK[1]+ulsch_harq->q_ACK[2];
1439 1440 1441


    if (metric_new > metric) {
Raymond Knopp's avatar
Raymond Knopp committed
1442 1443
      ulsch_harq->o_ACK[0] = 1;
      ulsch_harq->o_ACK[1] = 0;
1444 1445
      metric = metric_new;
    }
1446

Raymond Knopp's avatar
Raymond Knopp committed
1447
    metric_new = -ulsch_harq->q_ACK[0]-ulsch_harq->q_ACK[1]-ulsch_harq->q_ACK[2];
1448 1449

    if (metric_new > metric) {
Raymond Knopp's avatar
Raymond Knopp committed
1450 1451
      ulsch_harq->o_ACK[0] = 0;
      ulsch_harq->o_ACK[1] = 0;
1452 1453 1454 1455 1456 1457
      metric = metric_new;
    }
  }

  // RI

1458
  // rank 1
Raymond Knopp's avatar
Raymond Knopp committed
1459
  if ((ulsch_harq->O_RI == 1) && (Qprime_RI > 0)) {
1460
    ulsch_harq->o_RI[0] = ((ulsch_harq->q_RI[0] + ulsch_harq->q_RI[Q_m/2]) > 0) ? 0 : 1;
1461
  }
1462

1463 1464
  // CQI

1465
  //  printf("before cqi c[%d] = %p\n",0,ulsch_harq->c[0]);
1466
  ulsch_harq->cqi_crc_status = 0;
1467
  if (Q_CQI>0) {
Raymond Knopp's avatar
Raymond Knopp committed
1468
    memset((void *)&dummy_w_cc[0],0,3*(ulsch_harq->Or1+8+32));
1469

Raymond Knopp's avatar
Raymond Knopp committed
1470
    O_RCC = generate_dummy_w_cc(ulsch_harq->Or1+8,
1471 1472 1473
                                &dummy_w_cc[0]);


1474
    lte_rate_matching_cc_rx(O_RCC,
1475 1476 1477 1478 1479
                            Q_CQI,
                            ulsch_harq->o_w,
                            dummy_w_cc,
                            ulsch_harq->q);

Raymond Knopp's avatar
Raymond Knopp committed
1480
    sub_block_deinterleaving_cc((unsigned int)(ulsch_harq->Or1+8),
1481 1482 1483
                                &ulsch_harq->o_d[96],
                                &ulsch_harq->o_w[0]);

Cedric Roux's avatar
Cedric Roux committed
1484 1485
    memset(ulsch_harq->o,0,(7+8+ulsch_harq->Or1) / 8);
    phy_viterbi_lte_sse2(ulsch_harq->o_d+96,ulsch_harq->o,8+ulsch_harq->Or1);
1486

Cedric Roux's avatar
Cedric Roux committed
1487
    if (extract_cqi_crc(ulsch_harq->o,ulsch_harq->Or1) == (crc8(ulsch_harq->o,ulsch_harq->Or1)>>24))
Raymond Knopp's avatar
Raymond Knopp committed
1488
      ulsch_harq->cqi_crc_status = 1;
1489 1490

#ifdef DEBUG_ULSCH_DECODING
Raymond Knopp's avatar
Raymond Knopp committed
1491
    printf("ulsch_decoding: Or1=%d\n",ulsch_harq->Or1);
1492 1493

    for (i=0; i<1+((8+ulsch_harq->Or1)/8); i++)
Raymond Knopp's avatar
Raymond Knopp committed
1494
      printf("ulsch_decoding: O[%d] %d\n",i,ulsch_harq->o[i]);
1495

Raymond Knopp's avatar
Raymond Knopp committed
1496
    if (ulsch_harq->cqi_crc_status == 1)
Cedric Roux's avatar
Cedric Roux committed
1497
      printf("RX CQI CRC OK (%x)\n",extract_cqi_crc(ulsch_harq->o,ulsch_harq->Or1));
1498
    else
Cedric Roux's avatar
Cedric Roux committed
1499
      printf("RX CQI CRC NOT OK (%x)\n",extract_cqi_crc(ulsch_harq->o,ulsch_harq->Or1));
1500

1501 1502 1503
#endif
  }

Cedric Roux's avatar
Cedric Roux committed
1504
  LOG_D(PHY,"frame %d subframe %d O_ACK:%d o_ACK[]=%d:%d:%d:%d\n",frame,subframe,ulsch_harq->O_ACK,ulsch_harq->o_ACK[0],ulsch_harq->o_ACK[1],ulsch_harq->o_ACK[2],ulsch_harq->o_ACK[3]);
1505

1506
  // Do ULSCH Decoding for data portion
1507

1508
  ret = eNB->td(eNB,UE_id,harq_pid,llr8_flag);
1509

1510
  VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_PHY_ENB_ULSCH_DECODING0+harq_pid,0);
1511

1512 1513
  return(ret);
}