ulsch_coding.c 30.5 KB
Newer Older
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.0  (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 29 30 31 32 33 34 35 36 37 38 39 40 41 42
/*! \file PHY/LTE_TRANSPORT/ulsch_coding.c
* \brief Top-level routines for coding the ULSCH transport channel as described in 36.212 V8.6 2009-03
* \author R. Knopp
* \date 2011
* \version 0.1
* \company Eurecom
* \email: knopp@eurecom.fr
* \note
* \warning
*/

#include "PHY/defs.h"
#include "PHY/extern.h"

#include "PHY/CODING/defs.h"
#include "PHY/CODING/extern.h"
#include "PHY/CODING/lte_interleaver_inline.h"
#include "PHY/LTE_TRANSPORT/defs.h"
#include "defs.h"
#include "extern.h"
#include "SIMULATION/ETH_TRANSPORT/extern.h"
43
#include "UTIL/LOG/vcd_signal_dumper.h"
44

45
//#define DEBUG_ULSCH_CODING
46 47 48
//#define DEBUG_ULSCH_FREE 1

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




58 59
void free_ue_ulsch(LTE_UE_ULSCH_t *ulsch)
{
60 61 62 63 64
  int i;
  int r;

  if (ulsch) {
#ifdef DEBUG_ULSCH_FREE
65
    printf("Freeing ulsch %p\n",ulsch);
66
#endif
67

68
    for (i=0; i<8; i++) {
69
#ifdef DEBUG_ULSCH_FREE
70
      printf("Freeing ulsch process %d\n",i);
71
#endif
72

73 74
      if (ulsch->harq_processes[i]) {
#ifdef DEBUG_ULSCH_FREE
75
        printf("Freeing ulsch process %d (%p)\n",i,ulsch->harq_processes[i]);
76
#endif
77 78 79 80

        if (ulsch->harq_processes[i]->b) {
          free16(ulsch->harq_processes[i]->b,MAX_ULSCH_PAYLOAD_BYTES);
          ulsch->harq_processes[i]->b = NULL;
81
#ifdef DEBUG_ULSCH_FREE
82
          printf("Freeing ulsch process %d b (%p)\n",i,ulsch->harq_processes[i]->b);
83
#endif
84 85
        }

86
#ifdef DEBUG_ULSCH_FREE
87
        printf("Freeing ulsch process %d c (%p)\n",i,ulsch->harq_processes[i]->c);
88
#endif
89 90

        for (r=0; r<MAX_NUM_ULSCH_SEGMENTS; r++) {
91 92

#ifdef DEBUG_ULSCH_FREE
93
          printf("Freeing ulsch process %d c[%d] (%p)\n",i,r,ulsch->harq_processes[i]->c[r]);
94
#endif
95

96 97 98 99 100
          if (ulsch->harq_processes[i]->c[r]) {
            free16(ulsch->harq_processes[i]->c[r],((r==0)?8:0) + 3+768);
            ulsch->harq_processes[i]->c[r] = NULL;
          }
        }
101 102 103

        free16(ulsch->harq_processes[i],sizeof(LTE_UL_UE_HARQ_t));
        ulsch->harq_processes[i] = NULL;
104 105
      }
    }
106

107
    free16(ulsch,sizeof(LTE_UE_ULSCH_t));
108
    ulsch = NULL;
109
  }
110

111 112
}

113
LTE_UE_ULSCH_t *new_ue_ulsch(unsigned char N_RB_UL, uint8_t abstraction_flag)
114
{
115 116 117 118

  LTE_UE_ULSCH_t *ulsch;
  unsigned char exit_flag = 0,i,j,r;
  unsigned char bw_scaling =1;
119 120 121

  switch (N_RB_UL) {
  case 6:
122 123
    bw_scaling =16;
    break;
124

125 126 127
  case 25:
    bw_scaling =4;
    break;
128 129

  case 50:
130 131
    bw_scaling =2;
    break;
132

133 134 135 136
  default:
    bw_scaling =1;
    break;
  }
137

138
  ulsch = (LTE_UE_ULSCH_t *)malloc16(sizeof(LTE_UE_ULSCH_t));
139

140
  if (ulsch) {
141
    memset(ulsch,0,sizeof(LTE_UE_ULSCH_t));
142

143 144
    ulsch->Mlimit = 4;

145
    for (i=0; i<8; i++) {
knopp's avatar
knopp committed
146

147
      ulsch->harq_processes[i] = (LTE_UL_UE_HARQ_t *)malloc16(sizeof(LTE_UL_UE_HARQ_t));
148

149 150
      //      printf("ulsch->harq_processes[%d] %p\n",i,ulsch->harq_processes[i]);
      if (ulsch->harq_processes[i]) {
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
        memset(ulsch->harq_processes[i], 0, sizeof(LTE_UL_UE_HARQ_t));
        ulsch->harq_processes[i]->b = (unsigned char*)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 {
          LOG_E(PHY,"Can't get b\n");
          exit_flag=1;
        }

        if (abstraction_flag==0) {
          for (r=0; r<MAX_NUM_ULSCH_SEGMENTS; r++) {
            ulsch->harq_processes[i]->c[r] = (unsigned char*)malloc16(((r==0)?8:0) + 3+768);  // account for filler in first segment and CRCs for multiple segment case

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

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

    if ((abstraction_flag == 0) && (exit_flag==0)) {
183
      for (i=0; i<8; i++)
184 185 186 187
        for (j=0; j<96; j++)
          for (r=0; r<MAX_NUM_ULSCH_SEGMENTS; r++)
            ulsch->harq_processes[i]->d[r][j] = LTE_NULL;

188
      return(ulsch);
189
    } else if (abstraction_flag==1)
190 191
      return(ulsch);
  }
192

193 194 195
  LOG_E(PHY,"new_ue_ulsch exit flag, size of  %d ,   %d\n",exit_flag, sizeof(LTE_UE_ULSCH_t));
  free_ue_ulsch(ulsch);
  return(NULL);
196 197


198 199 200
}


201
uint32_t ulsch_encoding(uint8_t *a,
202
                        PHY_VARS_UE *ue,
203 204 205 206 207 208
                        uint8_t harq_pid,
                        uint8_t eNB_id,
                        uint8_t tmode,
                        uint8_t control_only_flag,
                        uint8_t Nbundled)
{
209

210 211 212 213 214
  time_stats_t *seg_stats=&ue->ulsch_segmentation_stats;
  time_stats_t *rm_stats=&ue->ulsch_rate_matching_stats;
  time_stats_t *te_stats=&ue->ulsch_turbo_encoding_stats;
  time_stats_t *i_stats=&ue->ulsch_interleaving_stats;
  time_stats_t *m_stats=&ue->ulsch_multiplexing_stats;
215

216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232
  //  uint16_t offset;
  uint32_t crc=1;
  uint16_t iind;
  uint32_t A;
  uint8_t Q_m=0;
  uint32_t Kr=0,Kr_bytes,r,r_offset=0;
  uint8_t y[6*14*1200],*yptr;;
  uint8_t *columnset;
  uint32_t sumKr=0;
  uint32_t Qprime,L,G,Q_CQI=0,Q_RI=0,Q_ACK=0,H=0,Hprime=0,Hpp=0,Cmux=0,Rmux=0,Rmux_prime=0;
  uint32_t Qprime_ACK=0,Qprime_CQI=0,Qprime_RI=0,len_ACK=0,len_RI=0;
  //  uint32_t E;
  uint8_t ack_parity;
  uint32_t i,q,j,iprime,j2;
  uint16_t o_RCC;
  uint8_t o_flip[8];
  uint32_t wACK_idx;
233 234 235 236
  LTE_DL_FRAME_PARMS *frame_parms=&ue->frame_parms;
  PHY_MEASUREMENTS *meas = &ue->measurements;
  LTE_UE_ULSCH_t *ulsch=ue->ulsch[eNB_id];
  LTE_UE_DLSCH_t **dlsch = ue->dlsch[eNB_id];
237
  uint16_t rnti;
238

239 240 241 242 243 244 245 246 247 248
  if (!ulsch) {
    LOG_E(PHY,"Null ulsch ptr %p\n",ulsch);
    return(-1);
  }

  if (harq_pid > 7) {
    LOG_E(PHY,"Illegal harq_pid %d\n",harq_pid);
    return(-1);
  }

249 250 251 252 253 254 255
  if (ulsch->harq_processes[harq_pid]->O_ACK > 2) {
    LOG_E(PHY,"Illegal O_ACK %d\n",ulsch->harq_processes[harq_pid]->O_ACK);
    return(-1);
  }

  if (ulsch->O_RI > 1) {
    LOG_E(PHY,"Illegal O_RI %d\n",ulsch->O_RI);
256 257 258
    return(-1);
  }

gauthier's avatar
gauthier committed
259
  VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_IN);
260

261 262
  // fill CQI/PMI information
  if (ulsch->O>0) {
263 264
    rnti = ue->pdcch_vars[eNB_id]->crnti;
    fill_CQI(ulsch,meas,0,harq_pid,ue->frame_parms.N_RB_DL,rnti, tmode,ue->sinr_eff);
265

266
    LOG_D(PHY,"UE CQI\n");
267
    print_CQI(ulsch->o,ulsch->uci_format,0,ue->frame_parms.N_RB_DL);
268 269

    // save PUSCH pmi for later (transmission modes 4,5,6)
270 271 272 273
    if (dlsch[0]) {
      //LOG_I(PHY,"XXX saving pmi for DL %x\n",pmi2hex_2Ar1(((wideband_cqi_rank1_2A_5MHz *)ulsch->o)->pmi));
      dlsch[0]->pmi_alloc = ((wideband_cqi_rank1_2A_5MHz *)ulsch->o)->pmi;
    }
274 275 276 277 278 279
  }

  if (ulsch->O<=32) {
    o_flip[0] = ulsch->o[3];
    o_flip[1] = ulsch->o[2];
    o_flip[2] = ulsch->o[1];
280 281
    o_flip[3] = ulsch->o[0];
  } else {
282 283 284 285 286 287 288 289 290
    o_flip[0] = ulsch->o[7];
    o_flip[1] = ulsch->o[6];
    o_flip[2] = ulsch->o[5];
    o_flip[3] = ulsch->o[4];
    o_flip[4] = ulsch->o[3];
    o_flip[5] = ulsch->o[2];
    o_flip[6] = ulsch->o[1];
    o_flip[7] = ulsch->o[0];
  }
291

292 293 294 295 296
  if (control_only_flag == 0) {
    A=ulsch->harq_processes[harq_pid]->TBS;
    Q_m = get_Qm_ul(ulsch->harq_processes[harq_pid]->mcs);

    ulsch->harq_processes[harq_pid]->control_only = 0;
297

298
#ifdef DEBUG_ULSCH_CODING
299
    printf("[PHY][UE] ULSCH coding : A %d, Qm %d, mcs %d, harq_pid %d, round %d, RV %d\n",
300 301 302 303 304 305 306 307
        ulsch->harq_processes[harq_pid]->TBS,
        Q_m,
        ulsch->harq_processes[harq_pid]->mcs,
        harq_pid,
        ulsch->harq_processes[harq_pid]->round,
        ulsch->harq_processes[harq_pid]->rvidx);

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

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

313
    printf("ulsch_coding: O=%d\n",ulsch->O);
314 315 316

    for (i=0; i<1+((8+ulsch->O)/8); i++) {
      //    ulsch->o[i] = i;
317
      printf("ulsch_coding: O[%d] %d\n",i,ulsch->o[i]);
318 319 320
    }

    if ((tmode != 4))
321
      print_CQI(ulsch->o,wideband_cqi_rank1_2A,0,ue->frame_parms.N_RB_DL);
322
    else
323
      print_CQI(ulsch->o,HLC_subband_cqi_rank1_2A,0,ue->frame_parms.N_RB_DL);
324

325
#endif
326

327
    if (ulsch->harq_processes[harq_pid]->round == 0) {  // this is a new packet
328

329 330 331
      start_meas(seg_stats);
      // Add 24-bit crc (polynomial A) to payload
      crc = crc24a(a,
332 333
                   A)>>8;

334 335 336
      a[A>>3] = ((uint8_t*)&crc)[2];
      a[1+(A>>3)] = ((uint8_t*)&crc)[1];
      a[2+(A>>3)] = ((uint8_t*)&crc)[0];
337

338 339 340
      ulsch->harq_processes[harq_pid]->B = A+24;
      ulsch->harq_processes[harq_pid]->b = a;
      lte_segmentation(ulsch->harq_processes[harq_pid]->b,
341 342 343 344 345 346 347 348
                       ulsch->harq_processes[harq_pid]->c,
                       ulsch->harq_processes[harq_pid]->B,
                       &ulsch->harq_processes[harq_pid]->C,
                       &ulsch->harq_processes[harq_pid]->Cplus,
                       &ulsch->harq_processes[harq_pid]->Cminus,
                       &ulsch->harq_processes[harq_pid]->Kplus,
                       &ulsch->harq_processes[harq_pid]->Kminus,
                       &ulsch->harq_processes[harq_pid]->F);
349 350

      stop_meas(seg_stats);
351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370

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

        Kr_bytes = Kr>>3;

        // get interleaver index for Turbo code (lookup in Table 5.1.3-3 36-212, V8.6 2009-03, p. 13-14)
        if (Kr_bytes<=64)
          iind = (Kr_bytes-5);
        else if (Kr_bytes <=128)
          iind = 59 + ((Kr_bytes-64)>>1);
        else if (Kr_bytes <= 256)
          iind = 91 + ((Kr_bytes-128)>>2);
        else if (Kr_bytes <= 768)
          iind = 123 + ((Kr_bytes-256)>>3);
        else {
          LOG_E(PHY,"ulsch_coding: Illegal codeword size %d!!!\n",Kr_bytes);
gauthier's avatar
gauthier committed
371
          VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
372 373 374 375
          return(-1);
        }


376
#ifdef DEBUG_ULSCH_CODING
377
        printf("Generating Code Segment %d (%d bits)\n",r,Kr);
378 379
        // generate codewords

380 381 382 383
        printf("bits_per_codeword (Kr)= %d\n",Kr);
        printf("N_RB = %d\n",ulsch->harq_processes[harq_pid]->nb_rb);
        printf("Ncp %d\n",frame_parms->Ncp);
        printf("Qm %d\n",Q_m);
384
#endif
385 386 387 388 389

        //  offset=0;


#ifdef DEBUG_ULSCH_CODING
390
        printf("Encoding ... iind %d f1 %d, f2 %d\n",iind,f1f2mat_old[iind*2],f1f2mat_old[(iind*2)+1]);
391
#endif
392 393 394 395 396 397 398 399 400
        start_meas(te_stats);
        threegpplte_turbo_encoder(ulsch->harq_processes[harq_pid]->c[r],
                                  Kr>>3,
                                  &ulsch->harq_processes[harq_pid]->d[r][96],
                                  (r==0) ? ulsch->harq_processes[harq_pid]->F : 0,
                                  f1f2mat_old[iind*2],   // f1 (see 36212-820, page 14)
                                  f1f2mat_old[(iind*2)+1]  // f2 (see 36212-820, page 14)
                                 );
        stop_meas(te_stats);
401
#ifdef DEBUG_ULSCH_CODING
402 403 404 405

        if (r==0)
          write_output("enc_output0.m","enc0",&ulsch->harq_processes[harq_pid]->d[r][96],(3*8*Kr_bytes)+12,1,4);

406
#endif
407 408 409 410 411 412
        start_meas(i_stats);
        ulsch->harq_processes[harq_pid]->RTC[r] =
          sub_block_interleaving_turbo(4+(Kr_bytes*8),
                                       &ulsch->harq_processes[harq_pid]->d[r][96],
                                       ulsch->harq_processes[harq_pid]->w[r]);
        stop_meas(i_stats);
413
      }
414

415
    }
416

417 418
    if (ulsch->harq_processes[harq_pid]->C == 0) {
      LOG_E(PHY,"null segment\n");
gauthier's avatar
gauthier committed
419
      VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
420 421
      return(-1);
    }
422

423
    sumKr = 0;
424 425

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

431 432
      sumKr += Kr;
    }
433
  } else { // This is a control-only PUSCH, set sumKr to O_CQI-MIN
434 435 436
    ulsch->harq_processes[harq_pid]->control_only = 1;
    sumKr = ulsch->O_CQI_MIN;
  }
437

438 439 440 441
  ulsch->harq_processes[harq_pid]->sumKr = sumKr;
  // Compute Q_ri (p. 23 36-212)

  Qprime = ulsch->O_RI*ulsch->harq_processes[harq_pid]->Msc_initial*ulsch->harq_processes[harq_pid]->Nsymb_initial * ulsch->beta_offset_ri_times8;
442

443 444 445 446 447
  if (Qprime > 0) {
    if ((Qprime % (8*sumKr)) > 0)
      Qprime = 1+(Qprime/(8*sumKr));
    else
      Qprime = Qprime/(8*sumKr);
448

449 450 451 452 453 454 455 456 457
    if (Qprime > 4*ulsch->harq_processes[harq_pid]->nb_rb * 12)
      Qprime = 4*ulsch->harq_processes[harq_pid]->nb_rb * 12;
  }

  Q_RI = Q_m*Qprime;
  Qprime_RI = Qprime;

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

459 460 461 462 463
  if (Qprime > 0) {
    if ((Qprime % (8*sumKr)) > 0)
      Qprime = 1+(Qprime/(8*sumKr));
    else
      Qprime = Qprime/(8*sumKr);
464

465 466 467 468 469 470 471 472 473
    if (Qprime > 4*ulsch->harq_processes[harq_pid]->nb_rb * 12)
      Qprime = 4*ulsch->harq_processes[harq_pid]->nb_rb * 12;
  }

  Q_ACK = Qprime * Q_m;
  Qprime_ACK = Qprime;

  // Compute Q_cqi, assume O>11, p. 26 36-212
  if (control_only_flag == 0) {
474

475 476
    if (ulsch->O < 12)
      L=0;
477
    else
478
      L=8;
479

480 481 482 483
    if (ulsch->O > 0)
      Qprime = (ulsch->O + L) * ulsch->harq_processes[harq_pid]->Msc_initial*ulsch->harq_processes[harq_pid]->Nsymb_initial * ulsch->beta_offset_cqi_times8;
    else
      Qprime = 0;
484

485 486
    if (Qprime > 0) {
      if ((Qprime % (8*sumKr)) > 0)
487
        Qprime = 1+(Qprime/(8*sumKr));
488
      else
489
        Qprime = Qprime/(8*sumKr);
490
    }
491

492 493 494 495
    G = ulsch->harq_processes[harq_pid]->nb_rb * (12 * Q_m) * (ulsch->Nsymb_pusch);

    if (Qprime > (G - ulsch->O_RI))
      Qprime = G - ulsch->O_RI;
496

497 498 499 500
    Q_CQI = Q_m * Qprime;
    Qprime_CQI = Qprime;


501

502
    G = G - Q_RI - Q_CQI;
503 504 505
    ulsch->harq_processes[harq_pid]->G = G;

/*
506
    LOG_I(PHY,"ULSCH Encoding G %d, Q_RI %d (O_RI%d, Msc_initial %d, Nsymb_initial%d, beta_offset_ri_times8 %d), Q_CQI %d, Q_ACK %d \n",G,Q_RI,ulsch->O_RI,ulsch->harq_processes[harq_pid]->Msc_initial,ulsch->harq_processes[harq_pid]->Nsymb_initial,ulsch->beta_offset_ri_times8,Q_CQI,Q_ACK);
507

Gabriel's avatar
Gabriel committed
508
    LOG_I(PHY,"ULSCH Encoding (Nid_cell %d, rnti %x): harq_pid %d round %d, RV %d, mcs %d, O_RI %d, O_ACK %d, G %d\n",
509 510 511 512 513 514 515 516 517
          frame_parms->Nid_cell,ulsch->rnti,
          harq_pid,
          ulsch->harq_processes[harq_pid]->round,
          ulsch->harq_processes[harq_pid]->rvidx,
          ulsch->harq_processes[harq_pid]->mcs,
          ulsch->O_RI,
          ulsch->harq_processes[harq_pid]->O_ACK,
          G);
*/
518

519 520
    if ((int)G < 0) {
      LOG_E(PHY,"FATAL: ulsch_coding.c G < 0 (%d) : Q_RI %d, Q_CQI %d, O %d, betaCQI_times8 %d)\n",G,Q_RI,Q_CQI,ulsch->O,ulsch->beta_offset_cqi_times8);
gauthier's avatar
gauthier committed
521
      VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
522 523 524 525 526 527 528 529 530 531 532 533 534
      return(-1);
    }


    // Data and control multiplexing (5.2.2.7 36-212)

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



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

    for (r=0; r<ulsch->harq_processes[harq_pid]->C; r++) {
537
#ifdef DEBUG_ULSCH_CODING
538
      printf("Rate Matching, Code segment %d (coded bits (G) %d,unpunctured/repeated bits per code segment %d,mod_order %d, nb_rb %d)...\n",
539 540 541 542
          r,
          G,
          Kr*3,
          Q_m,ulsch->harq_processes[harq_pid]->nb_rb);
543
#endif
544 545

      start_meas(rm_stats);
546
      r_offset += lte_rate_matching_turbo(ulsch->harq_processes[harq_pid]->RTC[r],
547 548 549 550 551
                                          G,
                                          ulsch->harq_processes[harq_pid]->w[r],
                                          ulsch->e+r_offset,
                                          ulsch->harq_processes[harq_pid]->C, // C
                                          NSOFT,                    // Nsoft,
552
                                          0,  // this means UL
553 554 555 556 557 558 559
                                          1,
                                          ulsch->harq_processes[harq_pid]->rvidx,
                                          get_Qm_ul(ulsch->harq_processes[harq_pid]->mcs),
                                          1,
                                          r,
                                          ulsch->harq_processes[harq_pid]->nb_rb,
                                          ulsch->harq_processes[harq_pid]->mcs);                       // r
560 561
      stop_meas(rm_stats);
#ifdef DEBUG_ULSCH_CODING
562

563
      if (r==ulsch->harq_processes[harq_pid]->C-1)
564 565
        write_output("enc_output.m","enc",ulsch->e,r_offset,1,4);

566 567
#endif
    }
568
  } else { //control-only PUSCH
569 570 571 572 573 574 575 576 577
    Q_CQI = (ulsch->harq_processes[harq_pid]->nb_rb * (12 * Q_m) * (ulsch->Nsymb_pusch)) - Q_RI;
    H = Q_CQI;
    Hprime = H/Q_m;
  }


  //  Do CQI coding
  if ((ulsch->O>1) && (ulsch->O < 12)) {
    LOG_E(PHY,"short CQI sizes not supported yet\n");
gauthier's avatar
gauthier committed
578
    VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
579
    return(-1);
580
  } else {
581 582
    // add 8-bit CRC
    crc = crc8(o_flip,
583
               ulsch->O)>>24;
584 585 586 587
#ifdef DEBUG_ULSCH_CODING
    printf("crc(cqi) tx : %x\n",crc);
#endif
    memset((void *)&ulsch->o_d[0],LTE_NULL,96);
588

589
    ccodelte_encode(ulsch->O,
590 591 592 593
                    1,
                    o_flip,
                    &ulsch->o_d[96],
                    0);
594 595


596 597 598
    o_RCC = sub_block_interleaving_cc(8+ulsch->O,
                                      &ulsch->o_d[96],
                                      ulsch->o_w);
599 600

    lte_rate_matching_cc(o_RCC,
601 602 603 604
                         Q_CQI,
                         ulsch->o_w,
                         ulsch->q);

605 606 607
  }

  i=0;
608

609 610 611 612 613 614 615 616
  //  Do RI coding
  if (ulsch->O_RI == 1) {
    switch (Q_m) {
    case 2:
      ulsch->q_RI[0] = ulsch->o_RI[0];
      ulsch->q_RI[1] = PUSCH_y;//ulsch->o_RI[0];
      len_RI=2;
      break;
617

618 619 620 621 622 623 624
    case 4:
      ulsch->q_RI[0] = ulsch->o_RI[0];
      ulsch->q_RI[1] = PUSCH_y;//1;
      ulsch->q_RI[2] = PUSCH_x;//ulsch->o_RI[0];
      ulsch->q_RI[3] = PUSCH_x;//1;
      len_RI=4;
      break;
625

626 627 628 629 630 631 632 633 634 635
    case 6:
      ulsch->q_RI[0] = ulsch->o_RI[0];
      ulsch->q_RI[1] = PUSCH_y;//1;
      ulsch->q_RI[2] = PUSCH_x;//1;
      ulsch->q_RI[3] = PUSCH_x;//ulsch->o_RI[0];
      ulsch->q_RI[4] = PUSCH_x;//1;
      ulsch->q_RI[5] = PUSCH_x;//1;
      len_RI=6;
      break;
    }
636
  } else if (ulsch->O_RI>1) {
637
    LOG_E(PHY,"RI cannot be more than 1 bit yet\n");
gauthier's avatar
gauthier committed
638
    VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
639 640
    return(-1);
  }
641

642 643 644
  //  Do ACK coding, Section 5.2.2.6 36.213 (p.23-24 in v8.6)
  wACK_idx = (ulsch->bundling==0) ? 4 : ((Nbundled-1)&3);
#ifdef DEBUG_ULSCH_CODING
645
  printf("ulsch_coding.c: Bundling %d, Nbundled %d, wACK_idx %d\n",
646 647
      ulsch->bundling,Nbundled,wACK_idx);
#endif
648

649 650 651 652 653 654 655 656
  // 1-bit ACK/NAK
  if (ulsch->harq_processes[harq_pid]->O_ACK == 1) {
    switch (Q_m) {
    case 2:
      ulsch->q_ACK[0] = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[1] = (ulsch->bundling==0)? PUSCH_y : ((ulsch->o_ACK[0]+wACK[wACK_idx][1])&1);//ulsch->o_ACK[0];
      len_ACK = 2;
      break;
657

658 659 660 661 662 663 664
    case 4:
      ulsch->q_ACK[0] = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[1] = (ulsch->bundling==0)? PUSCH_y : ((ulsch->o_ACK[0]+wACK[wACK_idx][1])&1);
      ulsch->q_ACK[2] = PUSCH_x;
      ulsch->q_ACK[3] = PUSCH_x;
      len_ACK = 4;
      break;
665

666 667 668 669 670 671 672 673 674 675 676
    case 6:
      ulsch->q_ACK[0] = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[1] = (ulsch->bundling==0)? PUSCH_y : ((ulsch->o_ACK[0]+wACK[wACK_idx][1])&1);
      ulsch->q_ACK[2] = PUSCH_x;
      ulsch->q_ACK[3] = PUSCH_x;
      ulsch->q_ACK[4] = PUSCH_x;
      ulsch->q_ACK[6] = PUSCH_x;
      len_ACK = 6;
      break;
    }
  }
677

678 679 680
  // two-bit ACK/NAK
  if (ulsch->harq_processes[harq_pid]->O_ACK == 2) {
    ack_parity = (ulsch->o_ACK[0]+ulsch->o_ACK[1])&1;
681

682 683 684 685 686 687 688 689 690 691
    switch (Q_m) {
    case 2:
      ulsch->q_ACK[0] = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[1] = (ulsch->o_ACK[1]+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[2] = (ack_parity+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[3] = (ulsch->o_ACK[0]+wACK[wACK_idx][1])&1;
      ulsch->q_ACK[4] = (ulsch->o_ACK[1]+wACK[wACK_idx][1])&1;
      ulsch->q_ACK[5] = (ack_parity+wACK[wACK_idx][1])&1;
      len_ACK = 6;
      break;
692

693 694 695 696 697 698 699 700 701 702 703 704 705 706 707
    case 4:
      ulsch->q_ACK[0]  = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[1]  = (ulsch->o_ACK[1]+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[2]  = PUSCH_x;
      ulsch->q_ACK[3]  = PUSCH_x;//1;
      ulsch->q_ACK[4]  = (ack_parity+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[5]  = (ulsch->o_ACK[0]+wACK[wACK_idx][1])&1;
      ulsch->q_ACK[6]  = PUSCH_x;
      ulsch->q_ACK[7]  = PUSCH_x;//1;
      ulsch->q_ACK[8]  = (ulsch->o_ACK[1]+wACK[wACK_idx][1])&1;
      ulsch->q_ACK[9]  = (ack_parity+wACK[wACK_idx][1])&1;
      ulsch->q_ACK[10] = PUSCH_x;
      ulsch->q_ACK[11] = PUSCH_x;//1;
      len_ACK = 12;
      break;
708

709 710 711 712 713 714 715 716 717 718 719
    case 6:
      ulsch->q_ACK[0] = (ulsch->o_ACK[0]+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[1] = (ulsch->o_ACK[1]+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[2] = PUSCH_x;
      ulsch->q_ACK[3] = PUSCH_x;
      ulsch->q_ACK[4] = PUSCH_x;
      ulsch->q_ACK[5] = PUSCH_x;

      ulsch->q_ACK[6] = (ack_parity+wACK[wACK_idx][0])&1;
      ulsch->q_ACK[7] = (ulsch->o_ACK[0]+wACK[wACK_idx][1])&1;
      ulsch->q_ACK[8] = PUSCH_x;
720
      ulsch->q_ACK[9] = PUSCH_x;
721 722 723 724 725 726 727 728 729 730 731 732 733 734
      ulsch->q_ACK[10] = PUSCH_x;
      ulsch->q_ACK[11] = PUSCH_x;

      ulsch->q_ACK[12] = (ulsch->o_ACK[1]+wACK[wACK_idx][1])&1;
      ulsch->q_ACK[13] = (ack_parity+wACK[wACK_idx][1])&1;
      ulsch->q_ACK[14] = PUSCH_x;
      ulsch->q_ACK[15] = PUSCH_x;
      ulsch->q_ACK[16] = PUSCH_x;
      ulsch->q_ACK[17] = PUSCH_x;
      len_ACK = 18;

      break;
    }
  }
735

736 737
  if (ulsch->harq_processes[harq_pid]->O_ACK > 2) {
    LOG_E(PHY,"ACK cannot be more than 2 bits yet\n");
gauthier's avatar
gauthier committed
738
    VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
739 740 741 742 743 744 745 746
    return(-1);
  }


  // channel multiplexing/interleaving

  start_meas(m_stats);
  Hpp = Hprime + Q_RI;
747

748 749 750 751 752 753 754 755 756
  Cmux       = ulsch->Nsymb_pusch;
  Rmux       = Hpp*Q_m/Cmux;
  Rmux_prime = Rmux/Q_m;

  Qprime_RI  = Q_RI / Q_m;
  Qprime_ACK = Q_ACK / Q_m;
  Qprime_CQI = Q_CQI / Q_m;

  //  printf("Qprime_CQI = %d\n",Qprime_CQI);
757
  // RI BITS
758 759 760 761 762 763 764

  memset(y,LTE_NULL,Q_m*Hpp);

  if (frame_parms->Ncp == 0)
    columnset = cs_ri_normal;
  else
    columnset = cs_ri_extended;
765 766 767 768

  j=0;

  for (i=0; i<Qprime_RI; i++) {
769
    r = Rmux_prime - 1 - (i>>2);
770 771

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

776 777 778 779 780 781 782 783 784 785 786
    j=(j+3)&3;

  }


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

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

788 789
    if (i<Qprime_CQI) {
      for (q=0;q<Q_m;q++) {
790 791
  y[q+(Q_m*j)] = ulsch->q[q+(Q_m*i)];
  //printf("cqi[%d] %d => y[%d]\n",q+(Q_m*i),ulsch->q[q+(Q_m*i)],q+(Q_m*j));
792 793 794 795
      }
    }
    else {
      for (q=0;q<Q_m;q++) {
796 797
  y[q+(Q_m*j)] = ulsch->e[q+(Q_m*iprime)];
  //  printf("e[%d] %d => y[%d]\n",q+(Q_m*iprime),ulsch->e[q+(Q_m*iprime)],q+(Q_m*j));
798 799 800 801 802 803
      }
    }
    j++;
  }
  */

804
  for (i=0; i<Qprime_CQI; i++) {
805 806 807

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

808
    for (q=0; q<Q_m; q++) {
809 810 811
      y[q+(Q_m*j)] = ulsch->q[q+(Q_m*i)];
      //        printf("cqi[%d] %d => y[%d] (j %d)\n",q+(Q_m*i),ulsch->q[q+(Q_m*i)],q+(Q_m*j),j);
    }
812

813 814 815 816
    j++;
  }

  j2 = j*Q_m;
817

818 819 820 821
  switch (Q_m) {

  case 2:

822
    for (iprime=0; iprime<(Hprime-Qprime_CQI)<<1; iprime+=2) {
823
      while (y[j2] != LTE_NULL) j2+=2;
824

825 826 827 828
      y[j2]   = ulsch->e[iprime];
      y[1+j2] = ulsch->e[1+iprime];
      j2+=2;
    }
829

830 831 832
    break;

  case 4:
833
    for (iprime=0; iprime<(Hprime-Qprime_CQI)<<2; iprime+=4) {
834
      while (y[j2] != LTE_NULL) j2+=4;
835

836 837 838 839 840
      y[j2]   = ulsch->e[iprime];
      y[1+j2] = ulsch->e[1+iprime];
      y[2+j2] = ulsch->e[2+iprime];
      y[3+j2] = ulsch->e[3+iprime];
      j2+=4;
841 842
    }

843 844 845
    break;

  case 6:
846
    for (iprime=0; iprime<(Hprime-Qprime_CQI)*6; iprime+=6) {
847
      while (y[j2] != LTE_NULL) j2+=6;
848

849 850 851 852 853 854 855
      y[j2]   = ulsch->e[iprime];
      y[1+j2] = ulsch->e[1+iprime];
      y[2+j2] = ulsch->e[2+iprime];
      y[3+j2] = ulsch->e[3+iprime];
      y[4+j2] = ulsch->e[4+iprime];
      y[5+j2] = ulsch->e[5+iprime];
      j2+=6;
856 857
    }

858 859 860
    break;

  }
861

862 863 864 865 866 867 868 869
  // HARQ-ACK Bits (Note these overwrite some bits)

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

  j=0;
870 871

  for (i=0; i<Qprime_ACK; i++) {
872
    r = Rmux_prime - 1 - (i>>2);
873 874

    for (q=0; q<Q_m; q++) {
875 876
      y[q+(Q_m*((r*Cmux) + columnset[j]))]  = ulsch->q_ACK[(q+(Q_m*i))%len_ACK];
#ifdef DEBUG_ULSCH_CODING
877
      printf("ulsch_coding.c: ACK %d => y[%d]=%d (i %d, r*Cmux %d, columnset %d)\n",q+(Q_m*i),
878 879
          q+(Q_m*((r*Cmux) + columnset[j])),ulsch->q_ACK[(q+(Q_m*i))%len_ACK],
          i,r*Cmux,columnset[j]);
880 881
#endif
    }
882

883 884 885 886 887 888
    j=(j+3)&3;

  }

  // write out buffer
  j=0;
889

890 891
  switch (Q_m) {
  case 2:
892 893 894 895 896
    for (i=0; i<Cmux; i++)
      for (r=0; r<Rmux_prime; r++) {
        yptr=&y[((r*Cmux)+i)<<1];
        ulsch->h[j++] = *yptr++;
        ulsch->h[j++] = *yptr++;
897
      }
898

899
    break;
900

901
  case 4:
902 903 904 905 906 907 908
    for (i=0; i<Cmux; i++)
      for (r=0; r<Rmux_prime; r++) {
        yptr = &y[((r*Cmux)+i)<<2];
        ulsch->h[j++] = *yptr++;
        ulsch->h[j++] = *yptr++;
        ulsch->h[j++] = *yptr++;
        ulsch->h[j++] = *yptr++;
909
      }
910

911
    break;
912

913
  case 6:
914 915 916 917 918 919 920 921 922
    for (i=0; i<Cmux; i++)
      for (r=0; r<Rmux_prime; r++) {
        yptr = &y[((r*Cmux)+i)*6];
        ulsch->h[j++] = *yptr++;
        ulsch->h[j++] = *yptr++;
        ulsch->h[j++] = *yptr++;
        ulsch->h[j++] = *yptr++;
        ulsch->h[j++] = *yptr++;
        ulsch->h[j++] = *yptr++;
923
      }
924

925
    break;
926

927 928 929
  default:
    break;
  }
930

931 932 933
  stop_meas(m_stats);

  if (j!=(H+Q_RI)) {
934
    LOG_E(PHY,"Error in output buffer length (j %d, H+Q_RI %d)\n",j,H+Q_RI);
gauthier's avatar
gauthier committed
935
    VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
936 937 938
    return(-1);
  }

gauthier's avatar
gauthier committed
939
  VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_UE_ULSCH_ENCODING, VCD_FUNCTION_OUT);
940 941 942 943 944 945 946 947 948
  return(0);
}


#ifdef PHY_ABSTRACTION
#ifdef OPENAIR2
#include "LAYER2/MAC/extern.h"
#include "LAYER2/MAC/defs.h"
#endif
949
int ulsch_encoding_emul(uint8_t *ulsch_buffer,
950
                        PHY_VARS_UE *ue,
951 952 953 954
                        uint8_t eNB_id,
                        uint8_t harq_pid,
                        uint8_t control_only_flag)
{
955

956 957 958 959 960
  LTE_UE_ULSCH_t *ulsch = ue->ulsch[eNB_id];
  LTE_UE_DLSCH_t **dlsch = ue->dlsch[eNB_id];
  PHY_MEASUREMENTS *meas = &ue->measurements;
  uint8_t tmode = ue->transmission_mode[eNB_id];
  uint16_t rnti=ue->pdcch_vars[eNB_id]->crnti;
961
  LOG_D(PHY,"EMUL UE ulsch_encoding for eNB %d,mod_id %d, harq_pid %d rnti %x, ACK(%d,%d) \n",
962
        eNB_id,ue->Mod_id, harq_pid, rnti,ulsch->o_ACK[0],ulsch->o_ACK[1]);
963 964

  if (ulsch->O>0) {
965 966
    /*
    if(flag_LA==1)
967
      sinr_eff = sinr_eff_cqi_calc(ue, eNB_id);
968 969 970
    else
      sinr_eff = meas->wideband_cqi_avg[eNB_id];
    */
971

972
    fill_CQI(ulsch,meas,eNB_id,harq_pid,ue->frame_parms.N_RB_DL,rnti,tmode,ue->sinr_eff);
973
    //LOG_D(PHY,"UE CQI\n");
974 975 976
    //    print_CQI(ulsch->o,ulsch->uci_format,eNB_id);

    // save PUSCH pmi for later (transmission modes 4,5,6)
977
    //    printf("ulsch: saving pmi for DL %x\n",pmi2hex_2Ar1(((wideband_cqi_rank1_2A_5MHz *)ulsch->o)->pmi));
978 979 980 981
    // if (ulsch->uci_format != HLC_subband_cqi_mcs_CBA)
    dlsch[0]->harq_processes[harq_pid]->pmi_alloc = ((wideband_cqi_rank1_2A_5MHz *)ulsch->o)->pmi;
  }

982
  memcpy(ue->ulsch[eNB_id]->harq_processes[harq_pid]->b,
983
         ulsch_buffer,
984
         ue->ulsch[eNB_id]->harq_processes[harq_pid]->TBS>>3);
985

986

987 988
  //memcpy(&UE_transport_info[ue->Mod_id].transport_blocks[UE_transport_info_TB_index[ue->Mod_id]],
  memcpy(&UE_transport_info[ue->Mod_id][ue->CC_id].transport_blocks,
989
         ulsch_buffer,
990 991
         ue->ulsch[eNB_id]->harq_processes[harq_pid]->TBS>>3);
  //UE_transport_info_TB_index[ue->Mod_id]+=ue->ulsch[eNB_id]->harq_processes[harq_pid]->TBS>>3;
992
  // navid: currently more than one eNB is not supported in the code
993 994 995 996 997
  UE_transport_info[ue->Mod_id][ue->CC_id].num_eNB = 1;
  UE_transport_info[ue->Mod_id][ue->CC_id].rnti[0] = ue->pdcch_vars[0]->crnti;
  UE_transport_info[ue->Mod_id][ue->CC_id].eNB_id[0]  = eNB_id;
  UE_transport_info[ue->Mod_id][ue->CC_id].harq_pid[0] = harq_pid;
  UE_transport_info[ue->Mod_id][ue->CC_id].tbs[0]     = ue->ulsch[eNB_id]->harq_processes[harq_pid]->TBS>>3 ;
998
  // printf("\nue->Mod_id%d\n",ue->Mod_id);
999 1000 1001 1002

  UE_transport_info[ue->Mod_id][ue->CC_id].cntl.pusch_flag = 1;
  //UE_transport_info[ue->Mod_id].cntl.pusch_uci = *(uint32_t *)ulsch->o;
  memcpy(UE_transport_info[ue->Mod_id][ue->CC_id].cntl.pusch_uci,
1003 1004
         ulsch->o,
         MAX_CQI_BYTES);
1005
  // printf("[UE]cqi is %d \n", ((HLC_subband_cqi_rank1_2A_5MHz *)ulsch->o)->cqi1);
1006

1007 1008 1009 1010
  UE_transport_info[ue->Mod_id][ue->CC_id].cntl.length_uci = ulsch->O;
  UE_transport_info[ue->Mod_id][ue->CC_id].cntl.uci_format = ulsch->uci_format;
  UE_transport_info[ue->Mod_id][ue->CC_id].cntl.pusch_ri = (ulsch->o_RI[0]&1)+((ulsch->o_RI[1]&1)<<1);
  UE_transport_info[ue->Mod_id][ue->CC_id].cntl.pusch_ack =   (ulsch->o_ACK[0]&1) + ((ulsch->o_ACK[1]&1)<<1);
1011
  //printf("ack is %d %d %d\n",UE_transport_info[ue->Mod_id].cntl.pusch_ack, (ulsch->o_ACK[1]&1)<<1, ulsch->o_ACK[0]&1);
1012
  return(0);
1013

1014 1015
}
#endif