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

22 23 24 25 26 27 28
/*! \file PHY/NR_TRANSPORT/nr_prach.c
 * \brief Routines for UE PRACH physical channel
 * \author R. Knopp, G. Casati
 * \date 2019
 * \version 0.2
 * \company Eurecom, Fraunhofer IIS
 * \email: knopp@eurecom.fr, guido.casati@iis.fraunhofer.de
Agustin's avatar
Agustin committed
29 30 31 32 33 34 35 36
 * \note
 * \warning
 */
#include "PHY/sse_intrin.h"
#include "common/utils/LOG/vcd_signal_dumper.h"

#include "PHY/impl_defs_nr.h"
#include "PHY/defs_nr_UE.h"
37
#include "PHY/NR_UE_TRANSPORT/nr_transport_proto_ue.h"
38
#include "PHY/NR_TRANSPORT/nr_transport_common_proto.h"
Agustin's avatar
Agustin committed
39 40 41 42 43 44

#include "common/utils/LOG/log.h"
#include "common/utils/LOG/vcd_signal_dumper.h"

#include "T.h"

45
//#define NR_PRACH_DEBUG 1
46

47 48 49 50
extern uint16_t prach_root_sequence_map_0_3[838];
extern uint16_t prach_root_sequence_map_abc[138];
extern uint16_t nr_du[838];
extern int16_t nr_ru[2*839];
51
extern const char *prachfmt[];
Agustin's avatar
Agustin committed
52

53 54 55 56 57 58 59 60
// Note:
// - prach_fmt_id is an ID used to map to the corresponding PRACH format value in prachfmt
// WIP todo:
// - take prach start symbol into account
// - idft for short sequence assumes we are transmitting starting in symbol 0 of a PRACH slot
// - Assumes that PRACH SCS is same as PUSCH SCS @ 30 kHz, take values for formats 0-2 and adjust for others below
// - Preamble index different from 0 is not detected by gNB
int32_t generate_nr_prach(PHY_VARS_NR_UE *ue, uint8_t gNB_id, uint8_t slot){
Agustin's avatar
Agustin committed
61

62
  NR_DL_FRAME_PARMS *fp=&ue->frame_parms;
63 64
  fapi_nr_config_request_t *nrUE_config = &ue->nrUE_config;
  NR_PRACH_RESOURCES_t *prach_resources = ue->prach_resources[gNB_id];
65
  fapi_nr_ul_config_prach_pdu *prach_pdu = &ue->prach_vars[gNB_id]->prach_pdu;
Agustin's avatar
Agustin committed
66

67
  uint8_t Mod_id, fd_occasion, preamble_index, restricted_set, not_found;
68 69
  uint16_t rootSequenceIndex, prach_fmt_id, NCS, *prach_root_sequence_map, preamble_offset = 0;
  uint16_t preamble_shift = 0, preamble_index0, n_shift_ra, n_shift_ra_bar, d_start, numshift, N_ZC, u, offset, offset2, first_nonzero_root_idx;
70
  int16_t prach_tmp[98304*2*4] __attribute__((aligned(32)));
Agustin's avatar
Agustin committed
71

72
  int16_t Ncp = 0, amp, *prach, *prach2, *prachF, *Xu;
73
  int32_t Xu_re, Xu_im;
74
  int prach_start, prach_sequence_length, i, prach_len, dftlen, mu, kbar, K, n_ra_prb, k, prachStartSymbol, sample_offset_slot;
75 76 77 78 79 80 81 82 83
  //int restricted_Type;

  prach                   = prach_tmp;
  prachF                  = ue->prach_vars[gNB_id]->prachF;
  amp                     = ue->prach_vars[gNB_id]->amp;
  Mod_id                  = ue->Mod_id;
  prach_sequence_length   = nrUE_config->prach_config.prach_sequence_length;
  N_ZC                    = (prach_sequence_length == 0) ? 839:139;
  mu                      = nrUE_config->prach_config.prach_sub_c_spacing;
84 85 86 87 88
  restricted_set          = prach_pdu->restricted_set;
  rootSequenceIndex       = prach_pdu->root_seq_id;
  n_ra_prb                = prach_pdu->freq_msg1;
  NCS                     = prach_pdu->num_cs;
  prach_fmt_id            = prach_pdu->prach_format;
89
  preamble_index          = prach_resources->ra_PreambleIndex;
90 91 92 93 94 95 96
  fd_occasion             = 0;
  prach_len               = 0;
  dftlen                  = 0;
  first_nonzero_root_idx  = 0;
  kbar                    = 1;
  K                       = 24;
  k                       = 12*n_ra_prb - 6*fp->N_RB_UL;
97
  prachStartSymbol        = prach_pdu->prach_start_symbol;
98 99 100 101 102 103 104
  //restricted_Type         = 0;

  compute_nr_prach_seq(nrUE_config->prach_config.prach_sequence_length,
                       nrUE_config->prach_config.num_prach_fd_occasions_list[fd_occasion].num_root_sequences,
                       nrUE_config->prach_config.num_prach_fd_occasions_list[fd_occasion].prach_root_sequence_index,
                       ue->X_u);

105 106
  sample_offset_slot = (prachStartSymbol==0?0:fp->ofdm_symbol_size*prachStartSymbol+fp->nb_prefix_samples0+fp->nb_prefix_samples*(prachStartSymbol-1));
  prach_start = fp->get_samples_slot_timestamp(slot, fp, 0) + sample_offset_slot;
107

Agustin's avatar
Agustin committed
108 109 110 111 112 113 114 115
  // First compute physical root sequence
  /************************************************************************
  * 4G and NR NCS tables are slightly different and depend on prach format
  * Table 6.3.3.1-5:  for preamble formats with delta_f_RA = 1.25 Khz (formats 0,1,2)
  * Table 6.3.3.1-6:  for preamble formats with delta_f_RA = 5 Khz (formats 3)
  * NOTE: Restricted set type B is not implemented
  *************************************************************************/

116
  prach_root_sequence_map = (prach_sequence_length == 0) ? prach_root_sequence_map_0_3 : prach_root_sequence_map_abc;
Agustin's avatar
Agustin committed
117 118

  if (restricted_set == 0) {
119 120
    // This is the relative offset (for unrestricted case) in the root sequence table (5.7.2-4 from 36.211) for the given preamble index
    preamble_offset = ((NCS==0)? preamble_index : (preamble_index/(N_ZC/NCS)));
Agustin's avatar
Agustin committed
121 122 123 124 125
    // This is the \nu corresponding to the preamble index
    preamble_shift  = (NCS==0)? 0 : (preamble_index % (N_ZC/NCS));
    preamble_shift *= NCS;
  } else { // This is the high-speed case

126 127 128
    #ifdef NR_PRACH_DEBUG
      LOG_I(PHY, "PRACH [UE %d] High-speed mode, NCS %d\n", Mod_id, NCS);
    #endif
Agustin's avatar
Agustin committed
129 130 131 132 133 134 135 136 137 138

    not_found = 1;
    preamble_index0 = preamble_index;
    // set preamble_offset to initial rootSequenceIndex and look if we need more root sequences for this
    // preamble index and find the corresponding cyclic shift
    preamble_offset = 0; // relative rootSequenceIndex;

    while (not_found == 1) {
      // current root depending on rootSequenceIndex and preamble_offset
      int index = (rootSequenceIndex + preamble_offset) % N_ZC;
139
      uint16_t n_group_ra = 0;
Agustin's avatar
Agustin committed
140

141
      if (prach_fmt_id<4) {
Agustin's avatar
Agustin committed
142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182
        // prach_root_sequence_map points to prach_root_sequence_map0_3
        DevAssert( index < sizeof(prach_root_sequence_map_0_3) / sizeof(prach_root_sequence_map_0_3[0]) );
      } else {
        // prach_root_sequence_map points to prach_root_sequence_map4
        DevAssert( index < sizeof(prach_root_sequence_map_abc) / sizeof(prach_root_sequence_map_abc[0]) );
      }

      u = prach_root_sequence_map[index];

      if ( (nr_du[u]<(N_ZC/3)) && (nr_du[u]>=NCS) ) {
        n_shift_ra     = nr_du[u]/NCS;
        d_start        = (nr_du[u]<<1) + (n_shift_ra * NCS);
        n_group_ra     = N_ZC/d_start;
        n_shift_ra_bar = max(0,(N_ZC-(nr_du[u]<<1)-(n_group_ra*d_start))/N_ZC);
      } else if  ( (nr_du[u]>=(N_ZC/3)) && (nr_du[u]<=((N_ZC - NCS)>>1)) ) {
        n_shift_ra     = (N_ZC - (nr_du[u]<<1))/NCS;
        d_start        = N_ZC - (nr_du[u]<<1) + (n_shift_ra * NCS);
        n_group_ra     = nr_du[u]/d_start;
        n_shift_ra_bar = min(n_shift_ra,max(0,(nr_du[u]- (n_group_ra*d_start))/NCS));
      } else {
        n_shift_ra     = 0;
        n_shift_ra_bar = 0;
      }

      // This is the number of cyclic shifts for the current root u
      numshift = (n_shift_ra*n_group_ra) + n_shift_ra_bar;

      if (numshift>0 && preamble_index0==preamble_index)
        first_nonzero_root_idx = preamble_offset;

      if (preamble_index0 < numshift) {
        not_found      = 0;
        preamble_shift = (d_start * (preamble_index0/n_shift_ra)) + ((preamble_index0%n_shift_ra)*NCS);

      } else { // skip to next rootSequenceIndex and recompute parameters
        preamble_offset++;
        preamble_index0 -= numshift;
      }
    }
  }

183
  // now generate PRACH signal
184
#ifdef NR_PRACH_DEBUG
185
    if (NCS>0)
cig's avatar
cig committed
186 187 188 189 190
      LOG_I(PHY, "PRACH [UE %d] generate PRACH in slot %d for RootSeqIndex %d, Preamble Index %d, PRACH Format %s, NCS %d (N_ZC %d): Preamble_offset %d, Preamble_shift %d msg1 frequency start %d\n",
        Mod_id,
        slot,
        rootSequenceIndex,
        preamble_index,
191
        prachfmt[prach_fmt_id],
192 193 194
        NCS,
        N_ZC,
        preamble_offset,
cig's avatar
cig committed
195 196
        preamble_shift,
        n_ra_prb);
197
  #endif
Agustin's avatar
Agustin committed
198 199

  //  nsymb = (frame_parms->Ncp==0) ? 14:12;
200 201 202 203 204 205 206 207
  //  subframe_offset = (unsigned int)frame_parms->ofdm_symbol_size*slot*nsymb;

  if (prach_sequence_length == 0 && prach_fmt_id == 3) {
    K = 4;
    kbar = 10;
  } else if (prach_sequence_length == 1) {
    K = 1;
    kbar = 2;
208
  }
Agustin's avatar
Agustin committed
209

210 211
  if (k<0)
    k += fp->ofdm_symbol_size;
Agustin's avatar
Agustin committed
212

213 214 215
  k *= K;
  k += kbar;
  k *= 2;
Agustin's avatar
Agustin committed
216

217 218 219 220 221 222
  LOG_I(PHY, "PRACH [UE %d] in slot %d, placing PRACH in position %d, msg1 frequency start %d, preamble_offset %d, first_nonzero_root_idx %d\n", Mod_id,
    slot,
    k,
    n_ra_prb,
    preamble_offset,
    first_nonzero_root_idx);
Agustin's avatar
Agustin committed
223 224 225

  Xu = (int16_t*)ue->X_u[preamble_offset-first_nonzero_root_idx];

226 227 228 229
  #if defined (PRACH_WRITE_OUTPUT_DEBUG)
    LOG_M("X_u.m", "X_u", (int16_t*)ue->X_u[preamble_offset-first_nonzero_root_idx], N_ZC, 1, 1);
  #endif

230 231 232 233 234 235 236 237 238 239 240 241 242 243
  /********************************************************
   *
   * In function init_prach_tables:
   * to compute quantized roots of unity ru(n) = 32767 * exp j*[ (2 * PI * n) / N_ZC ]
   *
   * In compute_prach_seq:
   * to calculate Xu = DFT xu = xu (inv_u*k) * Xu[0] (This is a Zadoff-Chou sequence property: DFT ZC sequence is another ZC sequence)
   *
   * In generate_prach:
   * to do the cyclic-shifted DFT by multiplying Xu[k] * ru[k*preamble_shift] as:
   * If X[k] = DFT x(n) -> X_shifted[k] = DFT x(n+preamble_shift) = X[k] * exp -j*[ (2*PI*k*preamble_shift) / N_ZC ]
   *
   *********************************************************/

244
  if (mu==1) {
245 246 247 248 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 274 275 276 277 278 279 280 281 282 283 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 321 322 323 324
  if (fp->N_RB_UL <= 100)
    AssertFatal(1 == 0, "N_RB_UL %d not support for NR PRACH yet\n", fp->N_RB_UL);
  else if (fp->N_RB_UL < 137) {
    if (fp->threequarter_fs == 0) {
      //40 MHz @ 61.44 Ms/s
      //50 MHz @ 61.44 Ms/s
      if (prach_sequence_length == 0) {
        if (prach_fmt_id == 0 || prach_fmt_id == 1 || prach_fmt_id == 2)
            dftlen = 49152;
        if (prach_fmt_id == 3)
            dftlen = 12288;
      } // 839 sequence
      else {
        switch (mu){
          case 1:
            dftlen = 2048;
            break;
          default:
            AssertFatal(1 == 0, "Shouldn't get here\n");
            break;
        }
      }
    } else { // threequarter sampling
      //  40 MHz @ 46.08 Ms/s
      if (prach_sequence_length == 0) {
        AssertFatal(fp->N_RB_UL <= 107, "cannot do 108..136 PRBs with 3/4 sampling\n");
        if (prach_fmt_id == 0 || prach_fmt_id == 1 || prach_fmt_id == 2)
          dftlen = 36864;
        if (prach_fmt_id == 3)
          dftlen = 9216;
      } else {
        switch (mu){
          case 1:
            dftlen = 1536;
          break;
          default:
            AssertFatal(1 == 0, "Shouldn't get here\n");
            break;
        }
      } // short format
    } // 3/4 sampling
  } // <=50 MHz BW
  else if (fp->N_RB_UL <= 273) {
    if (fp->threequarter_fs == 0) {
    //80,90,100 MHz @ 122.88 Ms/s
      if (prach_sequence_length == 0) {
        if (prach_fmt_id == 0 || prach_fmt_id == 1 || prach_fmt_id == 2)
          dftlen = 98304;
        if (prach_fmt_id == 3)
          dftlen = 24576;
      }
    } else { // threequarter sampling
      switch (mu){
        case 1:
          dftlen = 4096;
          break;
        default:
          AssertFatal(1 == 0, "Shouldn't get here\n");
          break;
      }
    }
  } else {
    AssertFatal(fp->N_RB_UL <= 217, "cannot do more than 217 PRBs with 3/4 sampling\n");
    //  80 MHz @ 92.16 Ms/s
    if (prach_sequence_length == 0) {
      if (prach_fmt_id == 0 || prach_fmt_id == 1 || prach_fmt_id == 2)
        dftlen = 73728;
      if (prach_fmt_id == 3)
        dftlen = 18432;
    } else {
      switch (mu){
        case 1:
          dftlen = 3072;
          break;
        default:
          AssertFatal(1 == 0, "Shouldn't get here\n");
          break;
      }
    }
  }
325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340
  }
  else if (mu==3) {
    if (fp->threequarter_fs) 
      AssertFatal(1==0,"3/4 sampling not supported for numerology %d\n",mu);
    
    if (prach_sequence_length == 0) 
	AssertFatal(1==0,"long prach not supported for numerology %d\n",mu);

    if (fp->N_RB_UL == 32) 
      dftlen=512;
    else if (fp->N_RB_UL == 66) 
      dftlen=1024;
    else 
      AssertFatal(1==0,"N_RB_UL %d not support for numerology %d\n",fp->N_RB_UL,mu);
  }

341 342 343 344 345 346 347 348 349 350

  for (offset=0,offset2=0; offset<N_ZC; offset++,offset2+=preamble_shift) {

    if (offset2 >= N_ZC)
      offset2 -= N_ZC;

    Xu_re = (((int32_t)Xu[offset<<1]*amp)>>15);
    Xu_im = (((int32_t)Xu[1+(offset<<1)]*amp)>>15);
    prachF[k++]= ((Xu_re*nr_ru[offset2<<1]) - (Xu_im*nr_ru[1+(offset2<<1)]))>>15;
    prachF[k++]= ((Xu_im*nr_ru[offset2<<1]) + (Xu_re*nr_ru[1+(offset2<<1)]))>>15;
351

352 353
    if (k==dftlen) k=0;
  }
Agustin's avatar
Agustin committed
354

355 356 357 358 359 360 361 362 363
  #if defined (PRACH_WRITE_OUTPUT_DEBUG)
    LOG_M("prachF.m", "prachF", &prachF[1804], 1024, 1, 1);
    LOG_M("Xu.m", "Xu", Xu, N_ZC, 1, 1);
  #endif

  if (prach_sequence_length == 0) {

    AssertFatal(prach_fmt_id < 4, "Illegal PRACH format %d for sequence length 839\n", prach_fmt_id);

364
    // Ncp here is given in terms of T_s wich is 30.72MHz sampling
365 366 367 368 369 370 371 372 373 374 375 376 377 378
    switch (prach_fmt_id) {
    case 0:
      Ncp = 3168;
      break;
    case 1:
      Ncp = 21024;
      break;
    case 2:
      Ncp = 4688;
      break;
    case 3:
      Ncp = 3168;
      break;
    }
Agustin's avatar
Agustin committed
379

380 381 382
  } else {

    switch (prach_fmt_id) {
383
    case 4: //A1
384 385
      Ncp = 288/(1<<mu);
      break;
386
    case 5: //A2
387 388
      Ncp = 576/(1<<mu);
      break;
389
    case 6: //A3
390 391
      Ncp = 864/(1<<mu);
      break;
392
    case 7: //B1
393
      Ncp = 216/(1<<mu);
Agustin's avatar
Agustin committed
394
    break;
395
    /*
396 397 398 399 400 401
    case 4: //B2
      Ncp = 360/(1<<mu);
      break;
    case 5: //B3
      Ncp = 504/(1<<mu);
      break;
402 403
    */
    case 8: //B4
404 405
      Ncp = 936/(1<<mu);
      break;
406
    case 9: //C0
407 408
      Ncp = 1240/(1<<mu);
      break;
409
    case 10: //C2
410 411 412 413 414 415
      Ncp = 2048/(1<<mu);
      break;
    default:
      AssertFatal(1==0,"Unknown PRACH format ID %d\n", prach_fmt_id);
      break;
    }
Agustin's avatar
Agustin committed
416 417
  }

cig's avatar
cig committed
418 419 420
  #ifdef NR_PRACH_DEBUG
    LOG_D(PHY, "PRACH [UE %d] Ncp %d, dftlen %d \n", Mod_id, Ncp, dftlen);
  #endif
Agustin's avatar
Agustin committed
421

422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 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 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 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559
    //actually what we should be checking here is how often the current prach crosses a 0.5ms boundary. I am not quite sure for which paramter set this would be the case, so I will ignore it for now and just check if the prach starts on a 0.5ms boundary
    uint8_t  use_extended_prach_prefix = 0;
    if(fp->numerology_index == 0) {
      if (prachStartSymbol == 0 || prachStartSymbol == 7)
	use_extended_prach_prefix = 1;
    }
    else {
      if (slot%(fp->slots_per_subframe/2)==0 && prachStartSymbol == 0)
	use_extended_prach_prefix = 1;
    }
    
  if (fp->N_RB_UL <= 34) { //32 PRB case 61.44Msps
    if (fp->threequarter_fs == 0) {
      Ncp<<=1; //to account for 61.44Mbps
      // This is after cyclic prefix 
      prach2 = prach+(Ncp<<1); //times 2 for complex samples
      if (prach_sequence_length == 0)
	AssertFatal(1==0,"no long PRACH for this PRACH size %d\n",fp->N_RB_UL);
      else {
	if (use_extended_prach_prefix) 
          Ncp+=32;  // 16*kappa, kappa=2 for 61.44Msps
	prach2 = prach+(Ncp<<1); //times 2 for complex samples
        if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
          idft(IDFT_512,prachF,prach2,1);
          // here we have |empty | Prach512 |
          if (prach_fmt_id != 9) {
            memmove(prach2+(512<<1),prach2,(512<<2));
            prach_len = (512*2)+Ncp;
          }
          else prach_len = (512*1)+Ncp;
          memmove(prach,prach+(512<<1),(Ncp<<2));
          // here we have |Prefix | Prach512 | Prach512 (if ! 0xc0)  |
        } else if (prach_fmt_id == 5) { // 6x512
          idft(IDFT_512,prachF,prach2,1);
          // here we have |empty | Prach512 |
          memmove(prach2+(512<<1),prach2,(512<<2));
          // here we have |empty | Prach512 | Prach512| empty512 | empty512 |
          memmove(prach2+(512<<2),prach2,(512<<3));
          // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 |
          memmove(prach,prach+(512<<1),(Ncp<<2));
          // here we have |Prefix | Prach512 |
          prach_len = (512*4)+Ncp;
        } else if (prach_fmt_id == 6) { // 6x512
          idft(IDFT_512,prachF,prach2,1);
          // here we have |empty | Prach512 |
          memmove(prach2+(512<<1),prach2,(512<<2));
          // here we have |empty | Prach512 | Prach512| empty512 | empty512 | empty512 | empty512
          memmove(prach2+(512<<2),prach2,(512<<3));
          // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 | empty512 | empty512
          memmove(prach2+(512<<3),prach2,(512<<3));
          // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512
          memmove(prach,prach+(512<<1),(Ncp<<2));
          // here we have |Prefix | Prach512 |
          prach_len = (512*6)+Ncp;
        } else if (prach_fmt_id == 8) { // 12x512
          idft(IDFT_512,prachF,prach2,1);
          // here we have |empty | Prach512 |
          memmove(prach2+(512<<1),prach2,(512<<2));
          // here we have |empty | Prach512 | Prach512| empty512 | empty512 | empty512 | empty512
          memmove(prach2+(512<<2),prach2,(512<<3));
          // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 | empty512 | empty512
          memmove(prach2+(512<<3),prach2,(512<<3));
          // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512
          memmove(prach2+(512<<1)*6,prach2,(512<<2)*6);
          // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512|
          memmove(prach,prach+(512<<1),(Ncp<<2));
          // here we have |Prefix | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512|
          prach_len = (512*12)+Ncp;
	}		
      }
    }
    else
      AssertFatal(1==0,"3/4 sampling not supported for this PRACH size %d\n",fp->N_RB_UL);
   
  }
  else if (fp->N_RB_UL <= 68) {//66 PRB case, 122.88 Msps 
    if (fp->threequarter_fs == 0) {
      Ncp<<=2; //to account for 122.88Mbps
      // This is after cyclic prefix 
      prach2 = prach+(Ncp<<1); //times 2 for complex samples
      if (prach_sequence_length == 0)
	AssertFatal(1==0,"no long PRACH for this PRACH size %d\n",fp->N_RB_UL);
      else {
	if (use_extended_prach_prefix) 
          Ncp+=64;  // 16*kappa, kappa=4 for 122.88Msps
	prach2 = prach+(Ncp<<1); //times 2 for complex samples
        if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
          idft(IDFT_1024,prachF,prach2,1);
          // here we have |empty | Prach1024 |
          if (prach_fmt_id != 9) {
            memmove(prach2+(1024<<1),prach2,(1024<<2));
            prach_len = (1024*2)+Ncp;
          }
          else prach_len = (1024*1)+Ncp;
          memmove(prach,prach+(1024<<1),(Ncp<<2));
          // here we have |Prefix | Prach1024 | Prach1024 (if ! 0xc0)  |
        } else if (prach_fmt_id == 5) { // 6x1024
          idft(IDFT_1024,prachF,prach2,1);
          // here we have |empty | Prach1024 |
          memmove(prach2+(1024<<1),prach2,(1024<<2));
          // here we have |empty | Prach1024 | Prach1024| empty1024 | empty1024 |
          memmove(prach2+(1024<<2),prach2,(1024<<3));
          // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 |
          memmove(prach,prach+(1024<<1),(Ncp<<2));
          // here we have |Prefix | Prach1024 |
          prach_len = (1024*4)+Ncp;
        } else if (prach_fmt_id == 6) { // 6x1024
          idft(IDFT_1024,prachF,prach2,1);
          // here we have |empty | Prach1024 |
          memmove(prach2+(1024<<1),prach2,(1024<<2));
          // here we have |empty | Prach1024 | Prach1024| empty1024 | empty1024 | empty1024 | empty1024
          memmove(prach2+(1024<<2),prach2,(1024<<3));
          // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 | empty1024 | empty1024
          memmove(prach2+(1024<<3),prach2,(1024<<3));
          // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024
          memmove(prach,prach+(1024<<1),(Ncp<<2));
          // here we have |Prefix | Prach1024 |
          prach_len = (1024*6)+Ncp;
        } else if (prach_fmt_id == 8) { // 12x1024
          idft(IDFT_1024,prachF,prach2,1);
          // here we have |empty | Prach1024 |
          memmove(prach2+(1024<<1),prach2,(1024<<2));
          // here we have |empty | Prach1024 | Prach1024| empty1024 | empty1024 | empty1024 | empty1024
          memmove(prach2+(1024<<2),prach2,(1024<<3));
          // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 | empty1024 | empty1024
          memmove(prach2+(1024<<3),prach2,(1024<<3));
          // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024
          memmove(prach2+(1024<<1)*6,prach2,(1024<<2)*6);
          // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024|
          memmove(prach,prach+(1024<<1),(Ncp<<2));
          // here we have |Prefix | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024|
          prach_len = (1024*12)+Ncp;
	}	
      }
    }
    else
      AssertFatal(1==0,"3/4 sampling not supported for this PRACH size %d\n",fp->N_RB_UL);
  }
560
  else if (fp->N_RB_UL < 137) { // 46.08 or 61.44 Ms/s
561
    if (fp->threequarter_fs == 0) { // full sampling @ 61.44 Ms/s
562 563 564
      Ncp<<=1; //to account for 61.44Mbps 
      // This is after cyclic prefix 
      prach2 = prach+(Ncp<<1); //times 2 for complex samples
565 566
      if (prach_sequence_length == 0){
        if (prach_fmt_id == 0) { //24576 samples @ 30.72 Ms/s, 49152 samples @ 61.44 Ms/s
567
          idft(IDFT_49152,prachF,prach2,1);
568 569 570 571 572
          // here we have |empty | Prach49152|
          memmove(prach,prach+(49152<<1),(Ncp<<2));
          // here we have |Prefix | Prach49152|
          prach_len = 49152+Ncp;
        } else if (prach_fmt_id == 1) { //24576 samples @ 30.72 Ms/s, 49152 samples @ 61.44 Ms/s
573
          idft(IDFT_49152,prachF,prach2,1);
574 575 576 577 578 579
          memmove(prach2+(49152<<1),prach2,(49152<<2));
          // here we have |empty | Prach49152 | Prach49152|
          memmove(prach,prach+(49152<<2),(Ncp<<2));
          // here we have |Prefix | Prach49152 | Prach49152|
          prach_len = (49152*2)+Ncp;
        } else if (prach_fmt_id == 2) { //24576 samples @ 30.72 Ms/s, 49152 samples @ 61.44 Ms/s
580
          idft(IDFT_49152,prachF,prach2,1);
581 582 583 584 585 586 587 588
          memmove(prach2+(49152<<1),prach2,(49152<<2));
          // here we have |empty | Prach49152 | Prach49152| empty49152 | empty49152
          memmove(prach2+(49152<<2),prach2,(49152<<3));
          // here we have |empty | Prach49152 | Prach49152| Prach49152 | Prach49152
          memmove(prach,prach+(49152<<3),(Ncp<<2));
          // here we have |Prefix | Prach49152 | Prach49152| Prach49152 | Prach49152
          prach_len = (49152*4)+Ncp;
        } else if (prach_fmt_id == 3) { // //6144 samples @ 30.72 Ms/s, 12288 samples @ 61.44 Ms/s
589
          idft(IDFT_12288,prachF,prach2,1);
590 591 592 593 594 595 596 597 598
          memmove(prach2+(12288<<1),prach2,(12288<<2));
          // here we have |empty | Prach12288 | Prach12288| empty12288 | empty12288
          memmove(prach2+(12288<<2),prach2,(12288<<3));
          // here we have |empty | Prach12288 | Prach12288| Prach12288 | Prach12288
          memmove(prach,prach+(12288<<3),(Ncp<<2));
          // here we have |Prefix | Prach12288 | Prach12288| Prach12288 | Prach12288
          prach_len = (12288*4)+Ncp;
        }
      } else { // short PRACH sequence
599 600 601
	if (use_extended_prach_prefix) 
	  Ncp+=32; // 16*kappa, kappa=2 for 61.44Msps 
	prach2 = prach+(Ncp<<1); //times 2 for complex samples
602
        if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
603
          idft(IDFT_2048,prachF,prach2,1);
604
          // here we have |empty | Prach2048 |
605
          if (prach_fmt_id != 9) {
606 607 608 609 610 611
            memmove(prach2+(2048<<1),prach2,(2048<<2));
            prach_len = (2048*2)+Ncp;
          }
          else prach_len = (2048*1)+Ncp;
          memmove(prach,prach+(2048<<1),(Ncp<<2));
          // here we have |Prefix | Prach2048 | Prach2048 (if ! 0xc0)  |
612
        } else if (prach_fmt_id == 5) { // 6x2048
613
          idft(IDFT_2048,prachF,prach2,1);
614 615 616 617 618 619 620 621
          // here we have |empty | Prach2048 |
          memmove(prach2+(2048<<1),prach2,(2048<<2));
          // here we have |empty | Prach2048 | Prach2048| empty2048 | empty2048 |
          memmove(prach2+(2048<<2),prach2,(2048<<3));
          // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 |
          memmove(prach,prach+(2048<<1),(Ncp<<2));
          // here we have |Prefix | Prach2048 |
          prach_len = (2048*4)+Ncp;
622
        } else if (prach_fmt_id == 6) { // 6x2048
623
          idft(IDFT_2048,prachF,prach2,1);
624 625 626 627 628 629 630 631 632 633
          // here we have |empty | Prach2048 |
          memmove(prach2+(2048<<1),prach2,(2048<<2));
          // here we have |empty | Prach2048 | Prach2048| empty2048 | empty2048 | empty2048 | empty2048
          memmove(prach2+(2048<<2),prach2,(2048<<3));
          // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 | empty2048 | empty2048
          memmove(prach2+(2048<<3),prach2,(2048<<3));
          // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048
          memmove(prach,prach+(2048<<1),(Ncp<<2));
          // here we have |Prefix | Prach2048 |
          prach_len = (2048*6)+Ncp;
634
        } else if (prach_fmt_id == 8) { // 12x2048
635
          idft(IDFT_2048,prachF,prach2,1);
636 637 638 639 640 641 642 643 644 645 646 647 648 649 650
          // here we have |empty | Prach2048 |
          memmove(prach2+(2048<<1),prach2,(2048<<2));
          // here we have |empty | Prach2048 | Prach2048| empty2048 | empty2048 | empty2048 | empty2048
          memmove(prach2+(2048<<2),prach2,(2048<<3));
          // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 | empty2048 | empty2048
          memmove(prach2+(2048<<3),prach2,(2048<<3));
          // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048
          memmove(prach2+(2048<<1)*6,prach2,(2048<<2)*6);
          // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048 | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048|
          memmove(prach,prach+(2048<<1),(Ncp<<2));
          // here we have |Prefix | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048 | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048|
          prach_len = (2048*12)+Ncp;
        }
      }
    } else {  // threequarter sampling @ 46.08 Ms/s
651 652
      Ncp = (Ncp*3)/2;
      prach2 = prach+(Ncp<<1);
653 654
      if (prach_sequence_length == 0){
        if (prach_fmt_id == 0) {
655
          idft(IDFT_36864,prachF,prach2,1);
656 657 658 659 660
          // here we have |empty | Prach73728|
          memmove(prach,prach+(36864<<1),(Ncp<<2));
          // here we have |Prefix | Prach73728|
          prach_len = (36864*1)+Ncp;
        } else if (prach_fmt_id == 1) {
661
          idft(IDFT_36864,prachF,prach2,1);
662 663 664 665 666 667
          memmove(prach2+(36864<<1),prach2,(36864<<2));
          // here we have |empty | Prach73728 | Prach73728|
          memmove(prach,prach+(36864<<2),(Ncp<<2));
          // here we have |Prefix | Prach73728 | Prach73728|
          prach_len = (36864*2)+Ncp;
        } else if (prach_fmt_id == 2) {
668
          idft(IDFT_36864,prachF,prach2,1);
669 670 671 672 673 674 675 676
          memmove(prach2+(36864<<1),prach2,(36864<<2));
          // here we have |empty | Prach73728 | Prach73728| empty73728 | empty73728
          memmove(prach2+(36864<<2),prach2,(36864<<3));
          // here we have |empty | Prach73728 | Prach73728| Prach73728 | Prach73728
          memmove(prach,prach+(36864<<3),(Ncp<<2));
          // here we have |Prefix | Prach73728 | Prach73728| Prach73728 | Prach73728
          prach_len = (36864*4)+Ncp;
        } else if (prach_fmt_id == 3) {
677
          idft(IDFT_9216,prachF,prach2,1);
678 679 680 681 682 683 684 685 686
          memmove(prach2+(9216<<1),prach2,(9216<<2));
          // here we have |empty | Prach9216 | Prach9216| empty9216 | empty9216
          memmove(prach2+(9216<<2),prach2,(9216<<3));
          // here we have |empty | Prach9216 | Prach9216| Prach9216 | Prach9216
          memmove(prach,prach+(9216<<3),(Ncp<<2));
          // here we have |Prefix | Prach9216 | Prach9216| Prach9216 | Prach9216
          prach_len = (9216*4)+Ncp;
        }
      } else { // short sequence
687 688 689
	if (use_extended_prach_prefix) 
	  Ncp+=24; // 16*kappa, kappa=1.5 for 46.08Msps 
	prach2 = prach+(Ncp<<1); //times 2 for complex samples
690
        if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
691
          idft(IDFT_1536,prachF,prach2,1);
692
          // here we have |empty | Prach1536 |
693
          if (prach_fmt_id != 9) {
694 695 696 697 698 699 700
            memmove(prach2+(1536<<1),prach2,(1536<<2));
            prach_len = (1536*2)+Ncp;
          }	else prach_len = (1536*1)+Ncp;

          memmove(prach,prach+(1536<<1),(Ncp<<2));
          // here we have |Prefix | Prach1536 | Prach1536 (if ! 0xc0) |

701
        } else if (prach_fmt_id == 5) { // 6x1536
702
          idft(IDFT_1536,prachF,prach2,1);
703 704 705 706 707 708 709 710
          // here we have |empty | Prach1536 |
          memmove(prach2+(1536<<1),prach2,(1536<<2));
          // here we have |empty | Prach1536 | Prach1536| empty1536 | empty1536 |
          memmove(prach2+(1536<<2),prach2,(1536<<3));
          // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 |
          memmove(prach,prach+(1536<<1),(Ncp<<2));
          // here we have |Prefix | Prach1536 |
          prach_len = (1536*4)+Ncp;
711
        } else if (prach_fmt_id == 6) { // 6x1536
712
          idft(IDFT_1536,prachF,prach2,1);
713 714 715 716 717 718 719 720 721 722
          // here we have |empty | Prach1536 |
          memmove(prach2+(1536<<1),prach2,(1536<<2));
          // here we have |empty | Prach1536 | Prach1536| empty1536 | empty1536 | empty1536 | empty1536
          memmove(prach2+(1536<<2),prach2,(1536<<3));
          // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 | empty1536 | empty1536
          memmove(prach2+(1536<<3),prach2,(1536<<3));
          // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536
          memmove(prach,prach+(1536<<1),(Ncp<<2));
          // here we have |Prefix | Prach1536 |
          prach_len = (1536*6)+Ncp;
723
        } else if (prach_fmt_id == 8) { // 12x1536
724
          idft(IDFT_1536,prachF,prach2,1);
725 726 727 728 729 730 731 732 733 734 735 736 737
          // here we have |empty | Prach1536 |
          memmove(prach2+(1536<<1),prach2,(1536<<2));
          // here we have |empty | Prach1536 | Prach1536| empty1536 | empty1536 | empty1536 | empty1536
          memmove(prach2+(1536<<2),prach2,(1536<<3));
          // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 | empty1536 | empty1536
          memmove(prach2+(1536<<3),prach2,(1536<<3));
          // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536
          memmove(prach2+(1536<<1)*6,prach2,(1536<<2)*6);
          // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536 | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536|
          memmove(prach,prach+(1536<<1),(Ncp<<2));
          // here we have |Prefix | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536 | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536|
          prach_len = (1536*12)+Ncp;
        }
738
      }
Agustin's avatar
Agustin committed
739
    }
740 741
  } else if (fp->N_RB_UL <= 273) {// 92.16 or 122.88 Ms/s
    if (fp->threequarter_fs == 0) { // full sampling @ 122.88 Ms/s
742 743 744
      Ncp<<=2; //to account for 122.88Mbps
      // This is after cyclic prefix
      prach2 = prach+(Ncp<<1); //times 2 for complex samples
745 746
      if (prach_sequence_length == 0){
        if (prach_fmt_id == 0) { //24576 samples @ 30.72 Ms/s, 98304 samples @ 122.88 Ms/s
747
          idft(IDFT_98304,prachF,prach2,1);
748 749 750 751 752
          // here we have |empty | Prach98304|
          memmove(prach,prach+(98304<<1),(Ncp<<2));
          // here we have |Prefix | Prach98304|
          prach_len = (98304*1)+Ncp;
        } else if (prach_fmt_id == 1) {
753
          idft(IDFT_98304,prachF,prach2,1);
754 755 756 757 758 759
          memmove(prach2+(98304<<1),prach2,(98304<<2));
          // here we have |empty | Prach98304 | Prach98304|
          memmove(prach,prach+(98304<<2),(Ncp<<2));
          // here we have |Prefix | Prach98304 | Prach98304|
          prach_len = (98304*2)+Ncp;
        } else if (prach_fmt_id == 2) {
760
          idft(IDFT_98304,prachF,prach2,1);
761 762 763 764 765 766 767 768
          memmove(prach2+(98304<<1),prach2,(98304<<2));
          // here we have |empty | Prach98304 | Prach98304| empty98304 | empty98304
          memmove(prach2+(98304<<2),prach2,(98304<<3));
          // here we have |empty | Prach98304 | Prach98304| Prach98304 | Prach98304
          memmove(prach,prach+(98304<<3),(Ncp<<2));
          // here we have |Prefix | Prach98304 | Prach98304| Prach98304 | Prach98304
          prach_len = (98304*4)+Ncp;
        } else if (prach_fmt_id == 3) { // 4x6144, Ncp 3168
769
          idft(IDFT_24576,prachF,prach2,1);
770 771 772 773 774 775 776 777 778
          memmove(prach2+(24576<<1),prach2,(24576<<2));
          // here we have |empty | Prach24576 | Prach24576| empty24576 | empty24576
          memmove(prach2+(24576<<2),prach2,(24576<<3));
          // here we have |empty | Prach24576 | Prach24576| Prach24576 | Prach24576
          memmove(prach,prach+(24576<<3),(Ncp<<2));
          // here we have |Prefix | Prach24576 | Prach24576| Prach24576 | Prach24576
          prach_len = (24576*4)+Ncp;
        }
      } else { // short sequence
779 780 781
	if (use_extended_prach_prefix) 
          Ncp+=64; // 16*kappa, kappa=4 for 122.88Msps
	prach2 = prach+(Ncp<<1); //times 2 for complex samples
782
        if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
783
          idft(IDFT_4096,prachF,prach2,1);
784
          // here we have |empty | Prach4096 |
785
          if (prach_fmt_id != 9) {
786 787 788 789 790
            memmove(prach2+(4096<<1),prach2,(4096<<2));
            prach_len = (4096*2)+Ncp; 
          }	else 	prach_len = (4096*1)+Ncp;
          memmove(prach,prach+(4096<<1),(Ncp<<2));
          // here we have |Prefix | Prach4096 | Prach4096 (if ! 0xc0) |
791
        } else if (prach_fmt_id == 5) { // 4x4096
792
          idft(IDFT_4096,prachF,prach2,1);
793 794 795 796 797 798 799 800
          // here we have |empty | Prach4096 |
          memmove(prach2+(4096<<1),prach2,(4096<<2));
          // here we have |empty | Prach4096 | Prach4096| empty4096 | empty4096 |
          memmove(prach2+(4096<<2),prach2,(4096<<3));
          // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 |
          memmove(prach,prach+(4096<<1),(Ncp<<2));
          // here we have |Prefix | Prach4096 |
          prach_len = (4096*4)+Ncp;
801
        } else if (prach_fmt_id == 6) { // 6x4096
802
          idft(IDFT_4096,prachF,prach2,1);
803 804 805 806 807 808 809 810 811 812
          // here we have |empty | Prach4096 |
          memmove(prach2+(4096<<1),prach2,(4096<<2));
          // here we have |empty | Prach4096 | Prach4096| empty4096 | empty4096 | empty4096 | empty4096
          memmove(prach2+(4096<<2),prach2,(4096<<3));
          // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 | empty4096 | empty4096
          memmove(prach2+(4096<<3),prach2,(4096<<3));
          // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096
          memmove(prach,prach+(4096<<1),(Ncp<<2));
          // here we have |Prefix | Prach4096 |
          prach_len = (4096*6)+Ncp;
813
        } else if (prach_fmt_id == 8) { // 12x4096
814
          idft(IDFT_4096,prachF,prach2,1);
815 816 817 818 819 820 821 822 823 824 825 826 827 828 829
          // here we have |empty | Prach4096 |
          memmove(prach2+(4096<<1),prach2,(4096<<2));
          // here we have |empty | Prach4096 | Prach4096| empty4096 | empty4096 | empty4096 | empty4096
          memmove(prach2+(4096<<2),prach2,(4096<<3));
          // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 | empty4096 | empty4096
          memmove(prach2+(4096<<3),prach2,(4096<<3));
          // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096
          memmove(prach2+(4096<<1)*6,prach2,(4096<<2)*6);
          // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096 | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096|
          memmove(prach,prach+(4096<<1),(Ncp<<2));
          // here we have |Prefix | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096 | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096|
          prach_len = (4096*12)+Ncp;
        }
      }
    } else { // three quarter sampling @ 92.16 Ms/s
830 831
      Ncp = (Ncp*3); //to account for 92.16 Msps
      prach2 = prach+(Ncp<<1); //times 2 for complex samples
832 833
      if (prach_sequence_length == 0){
        if (prach_fmt_id == 0) {
834
          idft(IDFT_73728,prachF,prach2,1);
835 836 837 838 839
          // here we have |empty | Prach73728|
          memmove(prach,prach+(73728<<1),(Ncp<<2));
          // here we have |Prefix | Prach73728|
          prach_len = (73728*1)+Ncp;
        } else if (prach_fmt_id == 1) {
840
          idft(IDFT_73728,prachF,prach2,1);
841 842 843 844 845 846
          memmove(prach2+(73728<<1),prach2,(73728<<2));
          // here we have |empty | Prach73728 | Prach73728|
          memmove(prach,prach+(73728<<2),(Ncp<<2));
          // here we have |Prefix | Prach73728 | Prach73728|
          prach_len = (73728*2)+Ncp;
        } if (prach_fmt_id == 2) {
847
          idft(IDFT_73728,prachF,prach2,1);
848 849 850 851 852 853 854 855
          memmove(prach2+(73728<<1),prach2,(73728<<2));
          // here we have |empty | Prach73728 | Prach73728| empty73728 | empty73728
          memmove(prach2+(73728<<2),prach2,(73728<<3));
          // here we have |empty | Prach73728 | Prach73728| Prach73728 | Prach73728
          memmove(prach,prach+(73728<<3),(Ncp<<2));
          // here we have |Prefix | Prach73728 | Prach73728| Prach73728 | Prach73728
          prach_len = (73728*4)+Ncp;
        } else if (prach_fmt_id == 3) {
856
          idft(IDFT_18432,prachF,prach2,1);
857 858 859 860 861 862 863 864 865
          memmove(prach2+(18432<<1),prach2,(18432<<2));
          // here we have |empty | Prach18432 | Prach18432| empty18432 | empty18432
          memmove(prach2+(18432<<2),prach2,(18432<<3));
          // here we have |empty | Prach18432 | Prach18432| Prach18432 | Prach18432
          memmove(prach,prach+(18432<<3),(Ncp<<2));
          // here we have |Prefix | Prach18432 | Prach18432| Prach18432 | Prach18432
          prach_len = (18432*4)+Ncp;
        }
      } else { // short sequence
866 867 868 869
	if (use_extended_prach_prefix) 
          Ncp+=48; // 16*kappa, kappa=3 for 92.16Msps 
	prach2 = prach+(Ncp<<1); //times 2 for complex samples
	if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
870
          idft(IDFT_3072,prachF,prach2,1);
871
          // here we have |empty | Prach3072 |
872
          if (prach_fmt_id != 9) {
873 874 875
            memmove(prach2+(3072<<1),prach2,(3072<<2));
            prach_len = (3072*2)+Ncp;
          } else 	  prach_len = (3072*1)+Ncp;
876 877
	  memmove(prach,prach+(3072<<1),(Ncp<<2));
	  // here we have |Prefix | Prach3072 | Prach3072 (if ! 0xc0) |
878
        } else if (prach_fmt_id == 6) { // 6x3072
879
          idft(IDFT_3072,prachF,prach2,1);
880 881 882 883 884 885 886 887 888 889
          // here we have |empty | Prach3072 |
          memmove(prach2+(3072<<1),prach2,(3072<<2));
          // here we have |empty | Prach3072 | Prach3072| empty3072 | empty3072 | empty3072 | empty3072
          memmove(prach2+(3072<<2),prach2,(3072<<3));
          // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 | empty3072 | empty3072
          memmove(prach2+(3072<<3),prach2,(3072<<3));
          // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072
          memmove(prach,prach+(3072<<1),(Ncp<<2));
          // here we have |Prefix | Prach3072 |
          prach_len = (3072*6)+Ncp;
890
        } else if (prach_fmt_id == 5) { // 4x3072
891
          idft(IDFT_3072,prachF,prach2,1);
892 893 894 895 896 897 898 899 900
          // here we have |empty | Prach3072 |
          memmove(prach2+(3072<<1),prach2,(3072<<2));
          // here we have |empty | Prach3072 | Prach3072| empty3072 | empty3072 |
          memmove(prach2+(3072<<2),prach2,(3072<<3));
          // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 |
          memmove(prach,prach+(3072<<1),(Ncp<<2));
          // here we have |Prefix | Prach3072 |
          prach_len = (3072*4)+Ncp;
        } else if (prach_fmt_id == 6) { // 12x3072
901
          idft(IDFT_3072,prachF,prach2,1);
902 903 904 905 906 907 908 909 910 911 912 913 914
          // here we have |empty | Prach3072 |
          memmove(prach2+(3072<<1),prach2,(3072<<2));
          // here we have |empty | Prach3072 | Prach3072| empty3072 | empty3072 | empty3072 | empty3072
          memmove(prach2+(3072<<2),prach2,(3072<<3));
          // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 | empty3072 | empty3072
          memmove(prach2+(3072<<3),prach2,(3072<<3));
          // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072
          memmove(prach2+(3072<<1)*6,prach2,(3072<<2)*6);
          // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072 | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072|
          memmove(prach,prach+(3072<<1),(Ncp<<2));
          // here we have |Prefix | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072 | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072|
          prach_len = (3072*12)+Ncp;
        }
915
      }
Agustin's avatar
Agustin committed
916 917 918
    }
  }

919
  #ifdef NR_PRACH_DEBUG
920
    LOG_I(PHY, "PRACH [UE %d] N_RB_UL %d prach_start %d, prach_len %d\n", Mod_id,
921 922
      fp->N_RB_UL,
      prach_start,
923
      prach_len);
924 925
  #endif

926 927 928 929
  for (i=0; i<prach_len; i++) {
    ((int16_t*)(&ue->common_vars.txdata[0][prach_start]))[2*i] = prach[2*i];
    ((int16_t*)(&ue->common_vars.txdata[0][prach_start]))[2*i+1] = prach[2*i+1];
  }
Agustin's avatar
Agustin committed
930

931 932 933 934 935
  //printf("----------------------\n");
  //for(int ii = prach_start; ii<2*(prach_start + prach_len); ii++){
  //  printf("PRACH rx data[%d] = %d\n", ii, ue->common_vars.txdata[0][ii]);
  //}
  //printf(" \n");
Agustin's avatar
Agustin committed
936

cig's avatar
cig committed
937
  #ifdef PRACH_WRITE_OUTPUT_DEBUG
938 939 940
    LOG_M("prach_tx0.m", "prachtx0", prach+(Ncp<<1), prach_len-Ncp, 1, 1);
    LOG_M("Prach_txsig.m","txs",(int16_t*)(&ue->common_vars.txdata[0][prach_start]), 2*(prach_start+prach_len), 1, 1)
  #endif
Agustin's avatar
Agustin committed
941

942
  return signal_energy((int*)prach, 256);
Agustin's avatar
Agustin committed
943 944
}