1600 lines
85 KiB
C
1600 lines
85 KiB
C
/*
|
|
* 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
|
|
*/
|
|
|
|
/*! \file PHY/NR_UE_TRANSPORT/pucch_nr.c
|
|
* \brief Top-level routines for generating and decoding the PUCCH physical channel
|
|
* \author A. Mico Pereperez
|
|
* \date 2018
|
|
* \version 0.1
|
|
* \company Eurecom
|
|
* \email:
|
|
* \note
|
|
* \warning
|
|
*/
|
|
//#include "PHY/defs.h"
|
|
#include "PHY/impl_defs_nr.h"
|
|
#include "PHY/defs_nr_common.h"
|
|
#include "PHY/defs_nr_UE.h"
|
|
//#include "PHY/extern.h"
|
|
//#include "LAYER2/MAC/extern.h"
|
|
#include "PHY/NR_UE_TRANSPORT/pucch_nr.h"
|
|
#include "PHY/NR_UE_TRANSPORT/nr_transport_proto_ue.h"
|
|
#include "PHY/NR_TRANSPORT/nr_transport_common_proto.h"
|
|
#include <openair1/PHY/CODING/nrSmallBlock/nr_small_block_defs.h>
|
|
#include "common/utils/LOG/log.h"
|
|
#include "common/utils/LOG/vcd_signal_dumper.h"
|
|
|
|
#include "T.h"
|
|
//#define NR_UNIT_TEST 1
|
|
#ifdef NR_UNIT_TEST
|
|
#define DEBUG_PUCCH_TX
|
|
#define DEBUG_NR_PUCCH_TX
|
|
#endif
|
|
|
|
//#define ONE_OVER_SQRT2 23170 // 32767/sqrt(2) = 23170 (ONE_OVER_SQRT2)
|
|
|
|
void nr_generate_pucch0(PHY_VARS_NR_UE *ue,
|
|
int32_t **txdataF,
|
|
NR_DL_FRAME_PARMS *frame_parms,
|
|
int16_t amp,
|
|
int nr_slot_tx,
|
|
fapi_nr_ul_config_pucch_pdu *pucch_pdu) {
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch0] start function at slot(nr_slot_tx)=%d\n",nr_slot_tx);
|
|
#endif
|
|
/*
|
|
* Implement TS 38.211 Subclause 6.3.2.3.1 Sequence generation
|
|
*
|
|
*/
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch0] sequence generation\n");
|
|
#endif
|
|
/*
|
|
* Defining cyclic shift hopping TS 38.211 Subclause 6.3.2.2.2
|
|
*/
|
|
// alpha is cyclic shift
|
|
double alpha;
|
|
// lnormal is the OFDM symbol number in the PUCCH transmission where l=0 corresponds to the first OFDM symbol of the PUCCH transmission
|
|
//uint8_t lnormal;
|
|
// lprime is the index of the OFDM symbol in the slot that corresponds to the first OFDM symbol of the PUCCH transmission in the slot given by [5, TS 38.213]
|
|
//uint8_t lprime;
|
|
// mcs is provided by TC 38.213 subclauses 9.2.3, 9.2.4, 9.2.5 FIXME!
|
|
//uint8_t mcs;
|
|
/*
|
|
* in TS 38.213 Subclause 9.2.1 it is said that:
|
|
* for PUCCH format 0 or PUCCH format 1, the index of the cyclic shift
|
|
* is indicated by higher layer parameter PUCCH-F0-F1-initial-cyclic-shift
|
|
*/
|
|
/*
|
|
* Implementing TS 38.211 Subclause 6.3.2.3.1, the sequence x(n) shall be generated according to:
|
|
* x(l*12+n) = r_u_v_alpha_delta(n)
|
|
*/
|
|
// the value of u,v (delta always 0 for PUCCH) has to be calculated according to TS 38.211 Subclause 6.3.2.2.1
|
|
uint8_t u[2]={0,0},v[2]={0,0};
|
|
|
|
LOG_D(PHY,"pucch0: nr_symbols %d, start_symbol %d, prb_start %d, second_hop_prb %d, group_hop_flag %d, sequence_hop_flag %d, mcs %d\n",pucch_pdu->nr_of_symbols,pucch_pdu->start_symbol_index,pucch_pdu->prb_start,pucch_pdu->second_hop_prb,pucch_pdu->group_hop_flag,pucch_pdu->sequence_hop_flag,pucch_pdu->mcs);
|
|
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch0] sequence generation: variable initialization for test\n");
|
|
#endif
|
|
// x_n contains the sequence r_u_v_alpha_delta(n)
|
|
int16_t x_n_re[2][24],x_n_im[2][24];
|
|
|
|
uint16_t startingPRB = pucch_pdu->prb_start + pucch_pdu->bwp_start;
|
|
pucch_GroupHopping_t pucch_GroupHopping = pucch_pdu->group_hop_flag + (pucch_pdu->sequence_hop_flag<<1);
|
|
|
|
// we proceed to calculate alpha according to TS 38.211 Subclause 6.3.2.2.2
|
|
int prb_offset[2]={startingPRB,startingPRB};
|
|
nr_group_sequence_hopping(pucch_GroupHopping,pucch_pdu->hopping_id,0,nr_slot_tx,&u[0],&v[0]); // calculating u and v value
|
|
if (pucch_pdu->freq_hop_flag == 1) {
|
|
nr_group_sequence_hopping(pucch_GroupHopping,pucch_pdu->hopping_id,1,nr_slot_tx,&u[1],&v[1]); // calculating u and v value
|
|
prb_offset[1] = pucch_pdu->second_hop_prb;
|
|
}
|
|
for (int l=0; l<pucch_pdu->nr_of_symbols; l++) {
|
|
alpha = nr_cyclic_shift_hopping(pucch_pdu->hopping_id,
|
|
pucch_pdu->initial_cyclic_shift,
|
|
pucch_pdu->mcs,l,
|
|
pucch_pdu->start_symbol_index,
|
|
nr_slot_tx);
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch0] sequence generation \tu=%d \tv=%d \talpha=%lf \t(for symbol l=%d)\n",u[l],v[l],alpha,l);
|
|
#endif
|
|
|
|
for (int n=0; n<12; n++) {
|
|
x_n_re[l][n] = (int16_t)(((((int32_t)(round(32767*cos(alpha*n))) * table_5_2_2_2_2_Re[u[l]][n])>>15)
|
|
- (((int32_t)(round(32767*sin(alpha*n))) * table_5_2_2_2_2_Im[u[l]][n])>>15))); // Re part of base sequence shifted by alpha
|
|
x_n_im[l][n] = (int16_t)(((((int32_t)(round(32767*cos(alpha*n))) * table_5_2_2_2_2_Im[u[l]][n])>>15)
|
|
+ (((int32_t)(round(32767*sin(alpha*n))) * table_5_2_2_2_2_Re[u[l]][n])>>15))); // Im part of base sequence shifted by alpha
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch0] sequence generation \tu=%d \tv=%d \talpha=%lf \tx_n(l=%d,n=%d)=(%d,%d)\n",
|
|
u[l],v[l],alpha,l,n,x_n_re[l][n],x_n_im[l][n]);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
/*
|
|
* Implementing TS 38.211 Subclause 6.3.2.3.2 Mapping to physical resources FIXME!
|
|
*/
|
|
//int32_t *txptr;
|
|
uint32_t re_offset=0;
|
|
uint8_t l2;
|
|
|
|
for (int l=0; l<pucch_pdu->nr_of_symbols; l++) {
|
|
l2=l+pucch_pdu->start_symbol_index;
|
|
re_offset = (12*prb_offset[l]) + frame_parms->first_carrier_offset;
|
|
if (re_offset>= frame_parms->ofdm_symbol_size)
|
|
re_offset-=frame_parms->ofdm_symbol_size;
|
|
|
|
//txptr = &txdataF[0][re_offset];
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch0] symbol %d PRB %d (%d)\n",l,prb_offset[l],re_offset);
|
|
#endif
|
|
for (int n=0; n<12; n++) {
|
|
|
|
((int16_t *)&txdataF[0][(l2*frame_parms->ofdm_symbol_size) + re_offset])[0] = (int16_t)(((int32_t)(amp) * x_n_re[l][n])>>15);
|
|
((int16_t *)&txdataF[0][(l2*frame_parms->ofdm_symbol_size) + re_offset])[1] = (int16_t)(((int32_t)(amp) * x_n_im[l][n])>>15);
|
|
//((int16_t *)txptr[0][re_offset])[0] = (int16_t)((int32_t)amp * x_n_re[(12*l)+n])>>15;
|
|
//((int16_t *)txptr[0][re_offset])[1] = (int16_t)((int32_t)amp * x_n_im[(12*l)+n])>>15;
|
|
//txptr[re_offset] = (x_n_re[(12*l)+n]<<16) + x_n_im[(12*l)+n];
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch0] mapping to RE \t amp=%d \tofdm_symbol_size=%d \tN_RB_DL=%d \tfirst_carrier_offset=%d \ttxptr(%u)=(x_n(l=%d,n=%d)=(%d,%d))\n",
|
|
amp,frame_parms->ofdm_symbol_size,frame_parms->N_RB_DL,frame_parms->first_carrier_offset,(l2*frame_parms->ofdm_symbol_size) + re_offset,
|
|
l2,n,((int16_t *)&txdataF[0][(l2*frame_parms->ofdm_symbol_size) + re_offset])[0],
|
|
((int16_t *)&txdataF[0][(l2*frame_parms->ofdm_symbol_size) + re_offset])[1]);
|
|
#endif
|
|
re_offset++;
|
|
if (re_offset>= frame_parms->ofdm_symbol_size)
|
|
re_offset-=frame_parms->ofdm_symbol_size;
|
|
}
|
|
}
|
|
}
|
|
|
|
void nr_generate_pucch1(PHY_VARS_NR_UE *ue,
|
|
int32_t **txdataF,
|
|
NR_DL_FRAME_PARMS *frame_parms,
|
|
int16_t amp,
|
|
int nr_slot_tx,
|
|
fapi_nr_ul_config_pucch_pdu *pucch_pdu) {
|
|
|
|
uint16_t m0 = pucch_pdu->initial_cyclic_shift;
|
|
uint64_t payload = pucch_pdu->payload;
|
|
uint8_t startingSymbolIndex = pucch_pdu->start_symbol_index;
|
|
uint8_t nrofSymbols = pucch_pdu->nr_of_symbols;
|
|
uint16_t startingPRB = pucch_pdu->prb_start + pucch_pdu->bwp_start;
|
|
uint8_t timeDomainOCC = pucch_pdu->time_domain_occ_idx;
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] start function at slot(nr_slot_tx)=%d payload=%lu m0=%d nrofSymbols=%d startingSymbolIndex=%d startingPRB=%d second_hop_prb=%d timeDomainOCC=%d nr_bit=%d\n",
|
|
nr_slot_tx,payload,m0,nrofSymbols,startingSymbolIndex,startingPRB,pucch_pdu->second_hop_prb,timeDomainOCC,pucch_pdu->n_bit);
|
|
#endif
|
|
/*
|
|
* Implement TS 38.211 Subclause 6.3.2.4.1 Sequence modulation
|
|
*
|
|
*/
|
|
// complex-valued symbol d_re, d_im containing complex-valued symbol d(0):
|
|
int16_t d_re=0, d_im=0;
|
|
|
|
if (pucch_pdu->n_bit == 1) { // using BPSK if M_bit=1 according to TC 38.211 Subclause 5.1.2
|
|
d_re = (payload&1)==0 ? (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15) : -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im = (payload&1)==0 ? (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15) : -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (pucch_pdu->n_bit == 2) { // using QPSK if M_bit=2 according to TC 38.211 Subclause 5.1.2
|
|
if (((payload&1)==0) && (((payload>>1)&1)==0)) {
|
|
d_re = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15); // 32767/sqrt(2) = 23170 (ONE_OVER_SQRT2)
|
|
d_im = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((payload&1)==0) && (((payload>>1)&1)==1)) {
|
|
d_re = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((payload&1)==1) && (((payload>>1)&1)==0)) {
|
|
d_re = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((payload&1)==1) && (((payload>>1)&1)==1)) {
|
|
d_re = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
}
|
|
// printf("d_re=%d\td_im=%d\n",(int)d_re,(int)d_im);
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] sequence modulation: payload=%lx \tde_re=%d \tde_im=%d\n",payload,d_re,d_im);
|
|
#endif
|
|
/*
|
|
* Defining cyclic shift hopping TS 38.211 Subclause 6.3.2.2.2
|
|
*/
|
|
// alpha is cyclic shift
|
|
double alpha;
|
|
// lnormal is the OFDM symbol number in the PUCCH transmission where l=0 corresponds to the first OFDM symbol of the PUCCH transmission
|
|
//uint8_t lnormal = 0 ;
|
|
// lprime is the index of the OFDM symbol in the slot that corresponds to the first OFDM symbol of the PUCCH transmission in the slot given by [5, TS 38.213]
|
|
uint8_t lprime = startingSymbolIndex;
|
|
// mcs = 0 except for PUCCH format 0
|
|
uint8_t mcs=0;
|
|
// r_u_v_alpha_delta_re and r_u_v_alpha_delta_im tables containing the sequence y(n) for the PUCCH, when they are multiplied by d(0)
|
|
// r_u_v_alpha_delta_dmrs_re and r_u_v_alpha_delta_dmrs_im tables containing the sequence for the DM-RS.
|
|
int16_t r_u_v_alpha_delta_re[12],r_u_v_alpha_delta_im[12],r_u_v_alpha_delta_dmrs_re[12],r_u_v_alpha_delta_dmrs_im[12];
|
|
/*
|
|
* in TS 38.213 Subclause 9.2.1 it is said that:
|
|
* for PUCCH format 0 or PUCCH format 1, the index of the cyclic shift
|
|
* is indicated by higher layer parameter PUCCH-F0-F1-initial-cyclic-shift
|
|
*/
|
|
/*
|
|
* the complex-valued symbol d_0 shall be multiplied with a sequence r_u_v_alpha_delta(n): y(n) = d_0 * r_u_v_alpha_delta(n)
|
|
*/
|
|
// the value of u,v (delta always 0 for PUCCH) has to be calculated according to TS 38.211 Subclause 6.3.2.2.1
|
|
uint8_t u=0,v=0;//,delta=0;
|
|
// if frequency hopping is disabled, intraSlotFrequencyHopping is not provided
|
|
// n_hop = 0
|
|
// if frequency hopping is enabled, intraSlotFrequencyHopping is provided
|
|
// n_hop = 0 for first hop
|
|
// n_hop = 1 for second hop
|
|
uint8_t n_hop = 0;
|
|
// Intra-slot frequency hopping shall be assumed when the higher-layer parameter intraSlotFrequencyHopping is provided,
|
|
// regardless of whether the frequency-hop distance is zero or not,
|
|
// otherwise no intra-slot frequency hopping shall be assumed
|
|
//uint8_t PUCCH_Frequency_Hopping = 0 ; // from higher layers
|
|
uint8_t intraSlotFrequencyHopping = 0;
|
|
|
|
if (pucch_pdu->freq_hop_flag) {
|
|
intraSlotFrequencyHopping=1;
|
|
}
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] intraSlotFrequencyHopping = %d \n",intraSlotFrequencyHopping);
|
|
#endif
|
|
/*
|
|
* Implementing TS 38.211 Subclause 6.3.2.4.2 Mapping to physical resources
|
|
*/
|
|
//int32_t *txptr;
|
|
uint32_t re_offset=0;
|
|
int i=0;
|
|
#define MAX_SIZE_Z 168 // this value has to be calculated from mprime*12*table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_noHop[pucch_symbol_length]+m*12+n
|
|
int16_t z_re[MAX_SIZE_Z],z_im[MAX_SIZE_Z];
|
|
int16_t z_dmrs_re[MAX_SIZE_Z],z_dmrs_im[MAX_SIZE_Z];
|
|
|
|
for (int l=0; l<nrofSymbols; l++) {
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] for symbol l=%d, lprime=%d\n",
|
|
l,lprime);
|
|
#endif
|
|
// y_n contains the complex value d multiplied by the sequence r_u_v
|
|
int16_t y_n_re[12],y_n_im[12];
|
|
|
|
if ((intraSlotFrequencyHopping == 1) && (l >= (int)floor(nrofSymbols/2))) n_hop = 1; // n_hop = 1 for second hop
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] entering function nr_group_sequence_hopping with n_hop=%d, nr_slot_tx=%d\n",
|
|
n_hop,nr_slot_tx);
|
|
#endif
|
|
pucch_GroupHopping_t pucch_GroupHopping = pucch_pdu->group_hop_flag + (pucch_pdu->sequence_hop_flag<<1);
|
|
nr_group_sequence_hopping(pucch_GroupHopping,pucch_pdu->hopping_id,n_hop,nr_slot_tx,&u,&v); // calculating u and v value
|
|
alpha = nr_cyclic_shift_hopping(pucch_pdu->hopping_id,m0,mcs,l,lprime,nr_slot_tx);
|
|
|
|
for (int n=0; n<12; n++) {
|
|
r_u_v_alpha_delta_re[n] = (int16_t)(((((int32_t)(round(32767*cos(alpha*n))) * table_5_2_2_2_2_Re[u][n])>>15)
|
|
- (((int32_t)(round(32767*sin(alpha*n))) * table_5_2_2_2_2_Im[u][n])>>15))); // Re part of base sequence shifted by alpha
|
|
r_u_v_alpha_delta_im[n] = (int16_t)(((((int32_t)(round(32767*cos(alpha*n))) * table_5_2_2_2_2_Im[u][n])>>15)
|
|
+ (((int32_t)(round(32767*sin(alpha*n))) * table_5_2_2_2_2_Re[u][n])>>15))); // Im part of base sequence shifted by alpha
|
|
r_u_v_alpha_delta_dmrs_re[n] = (int16_t)(((((int32_t)(round(32767*cos(alpha*n))) * table_5_2_2_2_2_Re[u][n])>>15)
|
|
- (((int32_t)(round(32767*sin(alpha*n))) * table_5_2_2_2_2_Im[u][n])>>15))); // Re part of DMRS base sequence shifted by alpha
|
|
r_u_v_alpha_delta_dmrs_im[n] = (int16_t)(((((int32_t)(round(32767*cos(alpha*n))) * table_5_2_2_2_2_Im[u][n])>>15)
|
|
+ (((int32_t)(round(32767*sin(alpha*n))) * table_5_2_2_2_2_Re[u][n])>>15))); // Im part of DMRS base sequence shifted by alpha
|
|
r_u_v_alpha_delta_dmrs_re[n] = (int16_t)(((int32_t)(amp*r_u_v_alpha_delta_dmrs_re[n]))>>15);
|
|
r_u_v_alpha_delta_dmrs_im[n] = (int16_t)(((int32_t)(amp*r_u_v_alpha_delta_dmrs_im[n]))>>15);
|
|
// printf("symbol=%d\tr_u_v_re=%d\tr_u_v_im=%d\n",l,r_u_v_alpha_delta_re[n],r_u_v_alpha_delta_im[n]);
|
|
// PUCCH sequence = DM-RS sequence multiplied by d(0)
|
|
y_n_re[n] = (int16_t)(((((int32_t)(r_u_v_alpha_delta_re[n])*d_re)>>15)
|
|
- (((int32_t)(r_u_v_alpha_delta_im[n])*d_im)>>15))); // Re part of y(n)
|
|
y_n_im[n] = (int16_t)(((((int32_t)(r_u_v_alpha_delta_re[n])*d_im)>>15)
|
|
+ (((int32_t)(r_u_v_alpha_delta_im[n])*d_re)>>15))); // Im part of y(n)
|
|
// printf("symbol=%d\tr_u_v_dmrs_re=%d\tr_u_v_dmrs_im=%d\n",l,r_u_v_alpha_delta_dmrs_re[n],r_u_v_alpha_delta_dmrs_im[n]);
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] sequence generation \tu=%d \tv=%d \talpha=%lf \tr_u_v_alpha_delta[n=%d]=(%d,%d) \ty_n[n=%d]=(%d,%d)\n",
|
|
u,v,alpha,n,r_u_v_alpha_delta_re[n],r_u_v_alpha_delta_im[n],n,y_n_re[n],y_n_im[n]);
|
|
#endif
|
|
}
|
|
|
|
/*
|
|
* The block of complex-valued symbols y(n) shall be block-wise spread with the orthogonal sequence wi(m)
|
|
* (defined in table_6_3_2_4_1_2_Wi_Re and table_6_3_2_4_1_2_Wi_Im)
|
|
* z(mprime*12*table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_noHop[pucch_symbol_length]+m*12+n)=wi(m)*y(n)
|
|
*
|
|
* The block of complex-valued symbols r_u_v_alpha_dmrs_delta(n) for DM-RS shall be block-wise spread with the orthogonal sequence wi(m)
|
|
* (defined in table_6_3_2_4_1_2_Wi_Re and table_6_3_2_4_1_2_Wi_Im)
|
|
* z(mprime*12*table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_noHop[pucch_symbol_length]+m*12+n)=wi(m)*y(n)
|
|
*
|
|
*/
|
|
// the orthogonal sequence index for wi(m) defined in TS 38.213 Subclause 9.2.1
|
|
// the index of the orthogonal cover code is from a set determined as described in [4, TS 38.211]
|
|
// and is indicated by higher layer parameter PUCCH-F1-time-domain-OCC
|
|
// In the PUCCH_Config IE, the PUCCH-format1, timeDomainOCC field
|
|
uint8_t w_index = timeDomainOCC;
|
|
// N_SF_mprime_PUCCH_1 contains N_SF_mprime from table 6.3.2.4.1-1 (depending on number of PUCCH symbols nrofSymbols, mprime and intra-slot hopping enabled/disabled)
|
|
uint8_t N_SF_mprime_PUCCH_1;
|
|
// N_SF_mprime_PUCCH_1 contains N_SF_mprime from table 6.4.1.3.1.1-1 (depending on number of PUCCH symbols nrofSymbols, mprime and intra-slot hopping enabled/disabled)
|
|
uint8_t N_SF_mprime_PUCCH_DMRS_1;
|
|
// N_SF_mprime_PUCCH_1 contains N_SF_mprime from table 6.3.2.4.1-1 (depending on number of PUCCH symbols nrofSymbols, mprime=0 and intra-slot hopping enabled/disabled)
|
|
uint8_t N_SF_mprime0_PUCCH_1;
|
|
// N_SF_mprime_PUCCH_1 contains N_SF_mprime from table 6.4.1.3.1.1-1 (depending on number of PUCCH symbols nrofSymbols, mprime=0 and intra-slot hopping enabled/disabled)
|
|
uint8_t N_SF_mprime0_PUCCH_DMRS_1;
|
|
// mprime is 0 if no intra-slot hopping / mprime is {0,1} if intra-slot hopping
|
|
uint8_t mprime = 0;
|
|
|
|
if (intraSlotFrequencyHopping == 0) { // intra-slot hopping disabled
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] block-wise spread with the orthogonal sequence wi(m) if intraSlotFrequencyHopping = %d, intra-slot hopping disabled\n",
|
|
intraSlotFrequencyHopping);
|
|
#endif
|
|
N_SF_mprime_PUCCH_1 = table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_noHop[nrofSymbols-1]; // only if intra-slot hopping not enabled (PUCCH)
|
|
N_SF_mprime_PUCCH_DMRS_1 = table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_noHop[nrofSymbols-1]; // only if intra-slot hopping not enabled (DM-RS)
|
|
N_SF_mprime0_PUCCH_1 = table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_noHop[nrofSymbols-1]; // only if intra-slot hopping not enabled mprime = 0 (PUCCH)
|
|
N_SF_mprime0_PUCCH_DMRS_1 = table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_noHop[nrofSymbols-1]; // only if intra-slot hopping not enabled mprime = 0 (DM-RS)
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] w_index = %d, N_SF_mprime_PUCCH_1 = %d, N_SF_mprime_PUCCH_DMRS_1 = %d, N_SF_mprime0_PUCCH_1 = %d, N_SF_mprime0_PUCCH_DMRS_1 = %d\n",
|
|
w_index, N_SF_mprime_PUCCH_1,N_SF_mprime_PUCCH_DMRS_1,N_SF_mprime0_PUCCH_1,N_SF_mprime0_PUCCH_DMRS_1);
|
|
#endif
|
|
|
|
for (int m=0; m < N_SF_mprime_PUCCH_1; m++) {
|
|
for (int n=0; n<12 ; n++) {
|
|
z_re[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*y_n_re[n])>>15)
|
|
- (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*y_n_im[n])>>15));
|
|
z_im[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*y_n_im[n])>>15)
|
|
+ (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*y_n_re[n])>>15));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] block-wise spread with wi(m) (mprime=%d, m=%d, n=%d) z[%d] = ((%d * %d - %d * %d), (%d * %d + %d * %d)) = (%d,%d)\n",
|
|
mprime, m, n, (mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n,
|
|
table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m],y_n_re[n],table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m],y_n_im[n],
|
|
table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m],y_n_im[n],table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m],y_n_re[n],
|
|
z_re[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n],z_im[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n]);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
for (int m=0; m < N_SF_mprime_PUCCH_DMRS_1; m++) {
|
|
for (int n=0; n<12 ; n++) {
|
|
z_dmrs_re[(mprime*12*N_SF_mprime0_PUCCH_DMRS_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_DMRS_1][w_index][m])*r_u_v_alpha_delta_dmrs_re[n])>>15)
|
|
- (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_DMRS_1][w_index][m])*r_u_v_alpha_delta_dmrs_im[n])>>15));
|
|
z_dmrs_im[(mprime*12*N_SF_mprime0_PUCCH_DMRS_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_DMRS_1][w_index][m])*r_u_v_alpha_delta_dmrs_im[n])>>15)
|
|
+ (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_DMRS_1][w_index][m])*r_u_v_alpha_delta_dmrs_re[n])>>15));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] block-wise spread with wi(m) (mprime=%d, m=%d, n=%d) z[%d] = ((%d * %d - %d * %d), (%d * %d + %d * %d)) = (%d,%d)\n",
|
|
mprime, m, n, (mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n,
|
|
table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m],r_u_v_alpha_delta_dmrs_re[n],table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m],r_u_v_alpha_delta_dmrs_im[n],
|
|
table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m],r_u_v_alpha_delta_dmrs_im[n],table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m],r_u_v_alpha_delta_dmrs_re[n],
|
|
z_dmrs_re[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n],z_dmrs_im[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n]);
|
|
#endif
|
|
// printf("gNB entering l=%d\tdmrs_re=%d\tdmrs_im=%d\n",l,z_dmrs_re[(mprime*12*N_SF_mprime0_PUCCH_DMRS_1)+(m*12)+n],z_dmrs_re[(mprime*12*N_SF_mprime0_PUCCH_DMRS_1)+(m*12)+n]);
|
|
}
|
|
}
|
|
}
|
|
|
|
if (intraSlotFrequencyHopping == 1) { // intra-slot hopping enabled
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] block-wise spread with the orthogonal sequence wi(m) if intraSlotFrequencyHopping = %d, intra-slot hopping enabled\n",
|
|
intraSlotFrequencyHopping);
|
|
#endif
|
|
N_SF_mprime_PUCCH_1 = table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_m0Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 0 (PUCCH)
|
|
N_SF_mprime_PUCCH_DMRS_1 = table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_m0Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 0 (DM-RS)
|
|
N_SF_mprime0_PUCCH_1 = table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_m0Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 0 (PUCCH)
|
|
N_SF_mprime0_PUCCH_DMRS_1 = table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_m0Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 0 (DM-RS)
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] w_index = %d, N_SF_mprime_PUCCH_1 = %d, N_SF_mprime_PUCCH_DMRS_1 = %d, N_SF_mprime0_PUCCH_1 = %d, N_SF_mprime0_PUCCH_DMRS_1 = %d\n",
|
|
w_index, N_SF_mprime_PUCCH_1,N_SF_mprime_PUCCH_DMRS_1,N_SF_mprime0_PUCCH_1,N_SF_mprime0_PUCCH_DMRS_1);
|
|
#endif
|
|
|
|
for (mprime = 0; mprime<2; mprime++) { // mprime can get values {0,1}
|
|
for (int m=0; m < N_SF_mprime_PUCCH_1; m++) {
|
|
for (int n=0; n<12 ; n++) {
|
|
z_re[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*y_n_re[n])>>15)
|
|
- (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*y_n_im[n])>>15));
|
|
z_im[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*y_n_im[n])>>15)
|
|
+ (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*y_n_re[n])>>15));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] block-wise spread with wi(m) (mprime=%d, m=%d, n=%d) z[%d] = ((%d * %d - %d * %d), (%d * %d + %d * %d)) = (%d,%d)\n",
|
|
mprime, m, n, (mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n,
|
|
table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m],y_n_re[n],table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m],y_n_im[n],
|
|
table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m],y_n_im[n],table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m],y_n_re[n],
|
|
z_re[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n],z_im[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n]);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
for (int m=0; m < N_SF_mprime_PUCCH_DMRS_1; m++) {
|
|
for (int n=0; n<12 ; n++) {
|
|
z_dmrs_re[(mprime*12*N_SF_mprime0_PUCCH_DMRS_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_DMRS_1][w_index][m])*r_u_v_alpha_delta_dmrs_re[n])>>15)
|
|
- (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_DMRS_1][w_index][m])*r_u_v_alpha_delta_dmrs_im[n])>>15));
|
|
z_dmrs_im[(mprime*12*N_SF_mprime0_PUCCH_DMRS_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_DMRS_1][w_index][m])*r_u_v_alpha_delta_dmrs_im[n])>>15)
|
|
+ (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_DMRS_1][w_index][m])*r_u_v_alpha_delta_dmrs_re[n])>>15));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] block-wise spread with wi(m) (mprime=%d, m=%d, n=%d) z[%d] = ((%d * %d - %d * %d), (%d * %d + %d * %d)) = (%d,%d)\n",
|
|
mprime, m, n, (mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n,
|
|
table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m],r_u_v_alpha_delta_dmrs_re[n],table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m],r_u_v_alpha_delta_dmrs_im[n],
|
|
table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m],r_u_v_alpha_delta_dmrs_im[n],table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m],r_u_v_alpha_delta_dmrs_re[n],
|
|
z_dmrs_re[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n],z_dmrs_im[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n]);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
N_SF_mprime_PUCCH_1 = table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_m1Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 1 (PUCCH)
|
|
N_SF_mprime_PUCCH_DMRS_1 = table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_m1Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 1 (DM-RS)
|
|
}
|
|
}
|
|
|
|
if ((intraSlotFrequencyHopping == 1) && (l<floor(nrofSymbols/2))) { // intra-slot hopping enabled, we need to calculate new offset PRB
|
|
startingPRB = startingPRB + pucch_pdu->second_hop_prb;
|
|
}
|
|
|
|
if ((startingPRB < (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 0)) { // if number RBs in bandwidth is even and current PRB is lower band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*startingPRB) + frame_parms->first_carrier_offset;
|
|
}
|
|
|
|
if ((startingPRB >= (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 0)) { // if number RBs in bandwidth is even and current PRB is upper band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*(startingPRB-(frame_parms->N_RB_DL>>1)));
|
|
}
|
|
|
|
if ((startingPRB < (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB is lower band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*startingPRB) + frame_parms->first_carrier_offset;
|
|
}
|
|
|
|
if ((startingPRB > (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB is upper band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*(startingPRB-(frame_parms->N_RB_DL>>1))) + 6;
|
|
}
|
|
|
|
if ((startingPRB == (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB contains DC
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*startingPRB) + frame_parms->first_carrier_offset;
|
|
}
|
|
|
|
//txptr = &txdataF[0][re_offset];
|
|
for (int n=0; n<12; n++) {
|
|
if ((n==6) && (startingPRB == (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) {
|
|
// if number RBs in bandwidth is odd and current PRB contains DC, we need to recalculate the offset when n=6 (for second half PRB)
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size);
|
|
}
|
|
|
|
if (l%2 == 1) { // mapping PUCCH according to TS38.211 subclause 6.4.1.3.1
|
|
((int16_t *)&txdataF[0][re_offset])[0] = z_re[i+n];
|
|
((int16_t *)&txdataF[0][re_offset])[1] = z_im[i+n];
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] mapping PUCCH to RE \t amp=%d \tofdm_symbol_size=%d \tN_RB_DL=%d \tfirst_carrier_offset=%d \tz_pucch[%d]=txptr(%u)=(x_n(l=%d,n=%d)=(%d,%d))\n",
|
|
amp,frame_parms->ofdm_symbol_size,frame_parms->N_RB_DL,frame_parms->first_carrier_offset,i+n,re_offset,
|
|
l,n,((int16_t *)&txdataF[0][re_offset])[0],((int16_t *)&txdataF[0][re_offset])[1]);
|
|
#endif
|
|
}
|
|
|
|
if (l%2 == 0) { // mapping DM-RS signal according to TS38.211 subclause 6.4.1.3.1
|
|
((int16_t *)&txdataF[0][re_offset])[0] = z_dmrs_re[i+n];
|
|
((int16_t *)&txdataF[0][re_offset])[1] = z_dmrs_im[i+n];
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] mapping DM-RS to RE \t amp=%d \tofdm_symbol_size=%d \tN_RB_DL=%d \tfirst_carrier_offset=%d \tz_dm-rs[%d]=txptr(%u)=(x_n(l=%d,n=%d)=(%d,%d))\n",
|
|
amp,frame_parms->ofdm_symbol_size,frame_parms->N_RB_DL,frame_parms->first_carrier_offset,i+n,re_offset,
|
|
l,n,((int16_t *)&txdataF[0][re_offset])[0],((int16_t *)&txdataF[0][re_offset])[1]);
|
|
#endif
|
|
// printf("gNb l=%d\ti=%d\treoffset=%d\tre=%d\tim=%d\n",l,i,re_offset,z_dmrs_re[i+n],z_dmrs_im[i+n]);
|
|
}
|
|
|
|
re_offset++;
|
|
}
|
|
|
|
if (l%2 == 1) i+=12;
|
|
}
|
|
}
|
|
|
|
#if 0
|
|
void nr_generate_pucch1_old(PHY_VARS_NR_UE *ue,
|
|
int32_t **txdataF,
|
|
NR_DL_FRAME_PARMS *frame_parms,
|
|
PUCCH_CONFIG_DEDICATED *pucch_config_dedicated,
|
|
uint64_t payload,
|
|
int16_t amp,
|
|
int nr_slot_tx,
|
|
uint8_t m0,
|
|
uint8_t nrofSymbols,
|
|
uint8_t startingSymbolIndex,
|
|
uint16_t startingPRB,
|
|
uint16_t startingPRB_intraSlotHopping,
|
|
uint8_t timeDomainOCC,
|
|
uint8_t nr_bit) {
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] start function at slot(nr_slot_tx)=%d payload=%d m0=%d nrofSymbols=%d startingSymbolIndex=%d startingPRB=%d startingPRB_intraSlotHopping=%d timeDomainOCC=%d nr_bit=%d\n",
|
|
nr_slot_tx,payload,m0,nrofSymbols,startingSymbolIndex,startingPRB,startingPRB_intraSlotHopping,timeDomainOCC,nr_bit);
|
|
#endif
|
|
/*
|
|
* Implement TS 38.211 Subclause 6.3.2.4.1 Sequence modulation
|
|
*
|
|
*/
|
|
// complex-valued symbol d_re, d_im containing complex-valued symbol d(0):
|
|
int16_t d_re, d_im;
|
|
|
|
if (nr_bit == 1) { // using BPSK if M_bit=1 according to TC 38.211 Subclause 5.1.2
|
|
d_re = (payload&1)==0 ? (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15) : -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im = (payload&1)==0 ? (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15) : -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (nr_bit == 2) { // using QPSK if M_bit=2 according to TC 38.211 Subclause 5.1.2
|
|
if (((payload&1)==0) && (((payload>>1)&1)==0)) {
|
|
d_re = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15); // 32767/sqrt(2) = 23170 (ONE_OVER_SQRT2)
|
|
d_im = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((payload&1)==0) && (((payload>>1)&1)==1)) {
|
|
d_re = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((payload&1)==1) && (((payload>>1)&1)==0)) {
|
|
d_re = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((payload&1)==1) && (((payload>>1)&1)==1)) {
|
|
d_re = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
}
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] sequence modulation: payload=%x \tde_re=%d \tde_im=%d\n",payload,d_re,d_im);
|
|
#endif
|
|
/*
|
|
* Defining cyclic shift hopping TS 38.211 Subclause 6.3.2.2.2
|
|
*/
|
|
// alpha is cyclic shift
|
|
double alpha;
|
|
// lnormal is the OFDM symbol number in the PUCCH transmission where l=0 corresponds to the first OFDM symbol of the PUCCH transmission
|
|
uint8_t lnormal = 0 ;
|
|
// lprime is the index of the OFDM symbol in the slot that corresponds to the first OFDM symbol of the PUCCH transmission in the slot given by [5, TS 38.213]
|
|
uint8_t lprime = startingSymbolIndex;
|
|
// mcs = 0 except for PUCCH format 0
|
|
uint8_t mcs=0;
|
|
// r_u_v_alpha_delta_re and r_u_v_alpha_delta_im tables containing the sequence for the DM-RS.
|
|
// When they are multiplied by d(0), they become the sequence y(n) for the PUCCH
|
|
int16_t r_u_v_alpha_delta_re[12],r_u_v_alpha_delta_im[12];
|
|
/*
|
|
* in TS 38.213 Subclause 9.2.1 it is said that:
|
|
* for PUCCH format 0 or PUCCH format 1, the index of the cyclic shift
|
|
* is indicated by higher layer parameter PUCCH-F0-F1-initial-cyclic-shift
|
|
*/
|
|
/*
|
|
* the complex-valued symbol d_0 shall be multiplied with a sequence r_u_v_alpha_delta(n): y(n) = d_0 * r_u_v_alpha_delta(n)
|
|
*/
|
|
// the value of u,v (delta always 0 for PUCCH) has to be calculated according to TS 38.211 Subclause 6.3.2.2.1
|
|
uint8_t u=0,v=0,delta=0;
|
|
// if frequency hopping is disabled, intraSlotFrequencyHopping is not provided
|
|
// n_hop = 0
|
|
// if frequency hopping is enabled, intraSlotFrequencyHopping is provided
|
|
// n_hop = 0 for first hop
|
|
// n_hop = 1 for second hop
|
|
uint8_t n_hop = 0;
|
|
// Intra-slot frequency hopping shall be assumed when the higher-layer parameter intraSlotFrequencyHopping is provided,
|
|
// regardless of whether the frequency-hop distance is zero or not,
|
|
// otherwise no intra-slot frequency hopping shall be assumed
|
|
//uint8_t PUCCH_Frequency_Hopping = 0 ; // from higher layers
|
|
uint8_t intraSlotFrequencyHopping = 0;
|
|
|
|
if (startingPRB != startingPRB_intraSlotHopping) {
|
|
intraSlotFrequencyHopping=1;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] intraSlotFrequencyHopping=%d \n",intraSlotFrequencyHopping);
|
|
#endif
|
|
// n_hop = 1 for second hop;
|
|
// FIXME
|
|
// When hopping will be activated we have to implement this function differently as PUCH signal generation depends on n_hop value for u,v calculation
|
|
}
|
|
|
|
// y_n contains the complex value d multiplied by the sequence r_u_v
|
|
int16_t y_n_re[12],y_n_im[12];
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] entering function nr_group_sequence_hopping with n_hop=%d, nr_slot_tx=%d\n",
|
|
n_hop,nr_slot_tx);
|
|
#endif
|
|
nr_group_sequence_hopping(ue->pucch_config_common_nr->pucch_GroupHopping,ue->pucch_config_common_nr->hoppingId,n_hop,nr_slot_tx,&u,&v); // calculating u and v value
|
|
alpha = nr_cyclic_shift_hopping(ue->pucch_config_common_nr->hoppingId,m0,mcs,lnormal,lprime,nr_slot_tx);
|
|
|
|
for (int n=0; n<12; n++) {
|
|
r_u_v_alpha_delta_re[n] = (int16_t)(((((int32_t)(round(32767*cos(alpha*n))) * table_5_2_2_2_2_Re[u][n])>>15)
|
|
- (((int32_t)(round(32767*sin(alpha*n))) * table_5_2_2_2_2_Im[u][n])>>15))); // Re part of base sequence shifted by alpha
|
|
r_u_v_alpha_delta_im[n] = (int16_t)(((((int32_t)(round(32767*cos(alpha*n))) * table_5_2_2_2_2_Im[u][n])>>15)
|
|
+ (((int32_t)(round(32767*sin(alpha*n))) * table_5_2_2_2_2_Re[u][n])>>15))); // Im part of base sequence shifted by alpha
|
|
// PUCCH sequence = DM-RS sequence multiplied by d(0)
|
|
y_n_re[n] = (int16_t)(((((int32_t)(r_u_v_alpha_delta_re[n])*d_re)>>15)
|
|
- (((int32_t)(r_u_v_alpha_delta_im[n])*d_im)>>15))); // Re part of y(n)
|
|
y_n_im[n] = (int16_t)(((((int32_t)(r_u_v_alpha_delta_re[n])*d_im)>>15)
|
|
+ (((int32_t)(r_u_v_alpha_delta_im[n])*d_re)>>15))); // Im part of y(n)
|
|
// DM-RS sequence
|
|
r_u_v_alpha_delta_re[n] = (int16_t)(((int32_t)amp*r_u_v_alpha_delta_re[n])>>15);
|
|
r_u_v_alpha_delta_im[n] = (int16_t)(((int32_t)amp*r_u_v_alpha_delta_im[n])>>15);
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] sequence generation \tu=%d \tv=%d \talpha=%lf \tr_u_v_alpha_delta[n=%d]=(%d,%d) \ty_n[n=%d]=(%d,%d)\n",
|
|
u,v,alpha,n,r_u_v_alpha_delta_re[n],r_u_v_alpha_delta_im[n],n,y_n_re[n],y_n_im[n]);
|
|
#endif
|
|
}
|
|
|
|
/*
|
|
* The block of complex-valued symbols y(n) shall be block-wise spread with the orthogonal sequence wi(m)
|
|
* (defined in table_6_3_2_4_1_2_Wi_Re and table_6_3_2_4_1_2_Wi_Im)
|
|
* z(mprime*12*table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_noHop[pucch_symbol_length]+m*12+n)=wi(m)*y(n)
|
|
*
|
|
* The block of complex-valued symbols r_u_v_alpha_delta(n) for DM-RS shall be block-wise spread with the orthogonal sequence wi(m)
|
|
* (defined in table_6_3_2_4_1_2_Wi_Re and table_6_3_2_4_1_2_Wi_Im)
|
|
* z(mprime*12*table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_noHop[pucch_symbol_length]+m*12+n)=wi(m)*y(n)
|
|
*
|
|
* We are not implementing intra-slot hopping at the moment (so mprime=0)FIXME!
|
|
*/
|
|
#define MAX_SIZE_Z 168 // this value has to be calculated from mprime*12*table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_noHop[pucch_symbol_length]+m*12+n
|
|
int16_t z_re[MAX_SIZE_Z],z_im[MAX_SIZE_Z];
|
|
int16_t z_dmrs_re[MAX_SIZE_Z],z_dmrs_im[MAX_SIZE_Z];
|
|
// the orthogonal sequence index for wi(m) defined in TS 38.213 Subclause 9.2.1
|
|
// the index of the orthogonal cover code is from a set determined as described in [4, TS 38.211]
|
|
// and is indicated by higher layer parameter PUCCH-F1-time-domain-OCC
|
|
// In the PUCCH_Config IE, the PUCCH-format1, timeDomainOCC field FIXME!
|
|
uint8_t w_index = timeDomainOCC; // to be filled with the value of timeDomainOCC, higher layers parameters FIXME !!!
|
|
// N_SF_mprime_PUCCH_1 contains N_SF_mprime from table 6.3.2.4.1-1 (depending on number of PUCCH symbols nrofSymbols, mprime and intra-slot hopping enabled/disabled)
|
|
uint8_t N_SF_mprime_PUCCH_1;
|
|
// N_SF_mprime_PUCCH_1 contains N_SF_mprime from table 6.4.1.3.1.1-1 (depending on number of PUCCH symbols nrofSymbols, mprime and intra-slot hopping enabled/disabled)
|
|
uint8_t N_SF_mprime_PUCCH_DMRS_1;
|
|
// N_SF_mprime_PUCCH_1 contains N_SF_mprime from table 6.3.2.4.1-1 (depending on number of PUCCH symbols nrofSymbols, mprime=0 and intra-slot hopping enabled/disabled)
|
|
uint8_t N_SF_mprime0_PUCCH_1;
|
|
// N_SF_mprime_PUCCH_1 contains N_SF_mprime from table 6.4.1.3.1.1-1 (depending on number of PUCCH symbols nrofSymbols, mprime=0 and intra-slot hopping enabled/disabled)
|
|
uint8_t N_SF_mprime0_PUCCH_DMRS_1;
|
|
// mprime is 0 if no intra-slot hopping / mprime is {0,1} if intra-slot hopping
|
|
uint8_t mprime = 0;
|
|
|
|
if (intraSlotFrequencyHopping == 0) { // intra-slot hopping disabled
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] block-wise spread with the orthogonal sequence wi(m) if intraSlotFrequencyHopping = %d\n",
|
|
intraSlotFrequencyHopping);
|
|
#endif
|
|
N_SF_mprime_PUCCH_1 = table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_noHop[nrofSymbols-1]; // only if intra-slot hopping not enabled (PUCCH)
|
|
N_SF_mprime_PUCCH_DMRS_1 = table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_noHop[nrofSymbols-1]; // only if intra-slot hopping not enabled (DM-RS)
|
|
N_SF_mprime0_PUCCH_1 = table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_noHop[nrofSymbols-1]; // only if intra-slot hopping not enabled mprime = 0 (PUCCH)
|
|
N_SF_mprime0_PUCCH_DMRS_1 = table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_noHop[nrofSymbols-1]; // only if intra-slot hopping not enabled mprime = 0 (DM-RS)
|
|
|
|
for (int m=0; m < N_SF_mprime_PUCCH_1; m++) {
|
|
for (int n=0; n<12 ; n++) {
|
|
z_re[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*y_n_re[n])>>15)
|
|
- (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*y_n_im[n])>>15));
|
|
z_im[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*y_n_im[n])>>15)
|
|
+ (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*y_n_re[n])>>15));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t\t z_pucch[%d] \t= ((%d \t* %d \t-%d \t* %d), (%d \t* %d \t+%d \t*%d)) = (%d,%d)\n",
|
|
(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n,
|
|
table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m],y_n_re[n],table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m],y_n_im[n],
|
|
table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m],y_n_im[n],table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m],y_n_re[n],
|
|
z_re[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n],z_im[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n]);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
for (int m=0; m < N_SF_mprime_PUCCH_DMRS_1; m++) {
|
|
for (int n=0; n<12 ; n++) {
|
|
z_dmrs_re[(mprime*12*N_SF_mprime0_PUCCH_DMRS_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*r_u_v_alpha_delta_re[n])>>15)
|
|
- (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*r_u_v_alpha_delta_im[n])>>15));
|
|
z_dmrs_im[(mprime*12*N_SF_mprime0_PUCCH_DMRS_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*r_u_v_alpha_delta_im[n])>>15)
|
|
+ (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*r_u_v_alpha_delta_re[n])>>15));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t\t z_dm-rs[%d] = ((),()) =(%d,%d)\n",
|
|
(mprime*12*N_SF_mprime0_PUCCH_DMRS_1)+(m*12)+n,z_dmrs_re[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n],z_dmrs_im[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n]);
|
|
#endif
|
|
}
|
|
}
|
|
}
|
|
|
|
if (intraSlotFrequencyHopping == 1) { // intra-slot hopping enabled
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] block-wise spread with the orthogonal sequence wi(m) if intraSlotFrequencyHopping = %d\n",
|
|
intraSlotFrequencyHopping);
|
|
#endif
|
|
N_SF_mprime_PUCCH_1 = table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_m0Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 0 (PUCCH)
|
|
N_SF_mprime_PUCCH_DMRS_1 = table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_m0Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 0 (DM-RS)
|
|
N_SF_mprime0_PUCCH_1 = table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_m0Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 0 (PUCCH)
|
|
N_SF_mprime0_PUCCH_DMRS_1 = table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_m0Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 0 (DM-RS)
|
|
|
|
for (mprime = 0; mprime<2; mprime++) { // mprime can get values {0,1}
|
|
for (int m=0; m < N_SF_mprime_PUCCH_1; m++) {
|
|
for (int n=0; n<12 ; n++) {
|
|
z_re[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*y_n_re[n])>>15)
|
|
- (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*y_n_im[n])>>15));
|
|
z_im[(mprime*12*N_SF_mprime0_PUCCH_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*y_n_im[n])>>15)
|
|
+ (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*y_n_re[n])>>15));
|
|
}
|
|
}
|
|
|
|
for (int m=0; m < N_SF_mprime_PUCCH_DMRS_1; m++) {
|
|
for (int n=0; n<12 ; n++) {
|
|
z_dmrs_re[(mprime*12*N_SF_mprime0_PUCCH_DMRS_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*r_u_v_alpha_delta_re[n])>>15)
|
|
- (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*r_u_v_alpha_delta_im[n])>>15));
|
|
z_dmrs_im[(mprime*12*N_SF_mprime0_PUCCH_DMRS_1)+(m*12)+n] = (int16_t)((((int32_t)(table_6_3_2_4_1_2_Wi_Re[N_SF_mprime_PUCCH_1][w_index][m])*r_u_v_alpha_delta_im[n])>>15)
|
|
+ (((int32_t)(table_6_3_2_4_1_2_Wi_Im[N_SF_mprime_PUCCH_1][w_index][m])*r_u_v_alpha_delta_re[n])>>15));
|
|
}
|
|
}
|
|
|
|
N_SF_mprime_PUCCH_1 = table_6_3_2_4_1_1_N_SF_mprime_PUCCH_1_m1Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 1 (PUCCH)
|
|
N_SF_mprime_PUCCH_DMRS_1 = table_6_4_1_3_1_1_1_N_SF_mprime_PUCCH_1_m1Hop[nrofSymbols-1]; // only if intra-slot hopping enabled mprime = 1 (DM-RS)
|
|
}
|
|
}
|
|
|
|
/*
|
|
* Implementing TS 38.211 Subclause 6.3.2.4.2 Mapping to physical resources
|
|
*/
|
|
int32_t *txptr;
|
|
uint32_t re_offset;
|
|
int i=0;
|
|
|
|
for (int l=0; l<nrofSymbols; l++) {
|
|
if ((intraSlotFrequencyHopping == 1) && (l<floor(nrofSymbols/2))) { // intra-slot hopping enabled, we need to calculate new PRB, FIXME!!!
|
|
startingPRB = startingPRB + startingPRB_intraSlotHopping;
|
|
}
|
|
|
|
if ((startingPRB < (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 0)) { // if number RBs in bandwidth is even and current PRB is lower band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*startingPRB) + frame_parms->first_carrier_offset;
|
|
}
|
|
|
|
if ((startingPRB >= (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 0)) { // if number RBs in bandwidth is even and current PRB is upper band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*(startingPRB-(frame_parms->N_RB_DL>>1)));
|
|
}
|
|
|
|
if ((startingPRB < (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB is lower band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*startingPRB) + frame_parms->first_carrier_offset;
|
|
}
|
|
|
|
if ((startingPRB > (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB is upper band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*(startingPRB-(frame_parms->N_RB_DL>>1))) + 6;
|
|
}
|
|
|
|
if ((startingPRB == (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB contains DC
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*startingPRB) + frame_parms->first_carrier_offset;
|
|
}
|
|
|
|
txptr = &txdataF[0][re_offset];
|
|
|
|
for (int n=0; n<12; n++) {
|
|
if ((n==6) && (startingPRB == (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) {
|
|
// if number RBs in bandwidth is odd and current PRB contains DC, we need to recalculate the offset when n=6 (for second half PRB)
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size);
|
|
}
|
|
|
|
if (l%2 == 1) { // mapping PUCCH according to TS38.211 subclause 6.4.1.3.1
|
|
((int16_t *)&txdataF[0][re_offset])[0] = z_re[i+n];
|
|
((int16_t *)&txdataF[0][re_offset])[1] = z_im[i+n];
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] mapping PUCCH to RE \t amp=%d \tofdm_symbol_size=%d \tN_RB_DL=%d \tfirst_carrier_offset=%d \tz_pucch[%d]=txptr(%d)=(x_n(l=%d,n=%d)=(%d,%d))\n",
|
|
amp,frame_parms->ofdm_symbol_size,frame_parms->N_RB_DL,frame_parms->first_carrier_offset,i+n,re_offset,
|
|
l,n,((int16_t *)&txdataF[0][re_offset])[0],((int16_t *)&txdataF[0][re_offset])[1]);
|
|
#endif
|
|
}
|
|
|
|
if (l%2 == 0) { // mapping DM-RS signal according to TS38.211 subclause 6.4.1.3.1
|
|
((int16_t *)&txdataF[0][re_offset])[0] = z_dmrs_re[i+n];
|
|
((int16_t *)&txdataF[0][re_offset])[1] = z_dmrs_im[i+n];
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch1] mapping DM-RS to RE \t amp=%d \tofdm_symbol_size=%d \tN_RB_DL=%d \tfirst_carrier_offset=%d \tz_dm-rs[%d]=txptr(%d)=(x_n(l=%d,n=%d)=(%d,%d))\n",
|
|
amp,frame_parms->ofdm_symbol_size,frame_parms->N_RB_DL,frame_parms->first_carrier_offset,i+n,re_offset,
|
|
l,n,((int16_t *)&txdataF[0][re_offset])[0],((int16_t *)&txdataF[0][re_offset])[1]);
|
|
#endif
|
|
}
|
|
|
|
re_offset++;
|
|
}
|
|
|
|
if (l%2 == 1) i+=12;
|
|
}
|
|
}
|
|
#endif //0
|
|
|
|
inline void nr_pucch2_3_4_scrambling(uint16_t M_bit,uint16_t rnti,uint16_t n_id,uint64_t *B64,uint8_t *btilde) __attribute__((always_inline));
|
|
inline void nr_pucch2_3_4_scrambling(uint16_t M_bit,uint16_t rnti,uint16_t n_id,uint64_t *B64,uint8_t *btilde) {
|
|
uint32_t x1, x2, s=0;
|
|
int i;
|
|
uint8_t c;
|
|
// c_init=nRNTI*2^15+n_id according to TS 38.211 Subclause 6.3.2.6.1
|
|
//x2 = (rnti) + ((uint32_t)(1+nr_slot_tx)<<16)*(1+(fp->Nid_cell<<1));
|
|
x2 = ((rnti)<<15)+n_id;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t\t [nr_pucch2_3_4_scrambling] gold sequence s=%x, M_bit %d\n",s,M_bit);
|
|
#endif
|
|
|
|
uint8_t *btildep=btilde;
|
|
int M_bit2=M_bit > 31 ? 32 : (M_bit&31), M_bit3=M_bit;
|
|
uint32_t B;
|
|
for (int iprime=0;iprime<=(M_bit>>5);iprime++,btildep+=32) {
|
|
s = lte_gold_generic(&x1, &x2, (iprime==0) ? 1 : 0);
|
|
B=((uint32_t*)B64)[iprime];
|
|
for (int n=0;n<M_bit2;n+=8)
|
|
LOG_D(PHY,"PUCCH2 encoded %d : %d,%d,%d,%d,%d,%d,%d,%d\n",n,
|
|
(B>>n)&1,
|
|
(B>>(n+1))&1,
|
|
(B>>(n+2))&1,
|
|
(B>>(n+3))&1,
|
|
(B>>(n+4))&1,
|
|
(B>>(n+5))&1,
|
|
(B>>(n+6))&1,
|
|
(B>>(n+7))&1
|
|
);
|
|
for (i=0; i<M_bit2; i++) {
|
|
c = (uint8_t)((s>>i)&1);
|
|
btildep[i] = (((B>>i)&1) ^ c);
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t\t\t btilde[%d]=%x from unscrambled bit %d and scrambling %d (%x)\n",i+(iprime<<5),btilde[i],((B>>i)&1),c,s>>i);
|
|
#endif
|
|
}
|
|
M_bit3-=32;
|
|
M_bit2=M_bit3 > 31 ? 32 : (M_bit3&31);
|
|
}
|
|
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t\t [nr_pucch2_3_4_scrambling] scrambling M_bit=%d bits\n", M_bit);
|
|
#endif
|
|
}
|
|
void nr_uci_encoding(uint64_t payload,
|
|
uint8_t nr_bit,
|
|
int fmt,
|
|
uint8_t is_pi_over_2_bpsk_enabled,
|
|
uint8_t nrofSymbols,
|
|
uint8_t nrofPRB,
|
|
uint8_t n_SF_PUCCH_s,
|
|
uint8_t intraSlotFrequencyHopping,
|
|
uint8_t add_dmrs,
|
|
uint64_t *b,
|
|
uint16_t *M_bit) {
|
|
/*
|
|
* Implementing TS 38.212 Subclause 6.3.1.2
|
|
*
|
|
*/
|
|
// A is the payload size, to be provided in function call
|
|
uint8_t A = nr_bit;
|
|
// L is the CRC size
|
|
//uint8_t L;
|
|
// E is the rate matching output sequence length as given in TS 38.212 subclause 6.3.1.4.1
|
|
uint16_t E=0,E_init;
|
|
|
|
if (fmt == 2) E = 16*nrofSymbols*nrofPRB;
|
|
|
|
if (fmt == 3) {
|
|
E_init = (is_pi_over_2_bpsk_enabled == 0) ? 24:12;
|
|
|
|
if (nrofSymbols == 4) {
|
|
E = (intraSlotFrequencyHopping == 0)?(E_init*(nrofSymbols-1)*nrofPRB):((E_init*(nrofSymbols-1)*nrofPRB));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("format 3 nrofSymbols =4 and E_init=%d,E=%d\n",E_init,E);
|
|
#endif
|
|
}
|
|
|
|
if (nrofSymbols > 4) {
|
|
E = E_init*(nrofSymbols-2)*nrofPRB;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("format 3 nrofSymbols >4 and E_init=%d,E = %d\n",E_init,E);
|
|
#endif
|
|
}
|
|
|
|
if (nrofSymbols > 9) {
|
|
E = (add_dmrs == 0)?(E_init*(nrofSymbols-2)*nrofPRB):((E_init*(nrofSymbols-4)*nrofPRB));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("format 3 nrofSymbols >9 and E_init=%d,E = %d\n",E_init,E);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
if (fmt == 4) {
|
|
E_init = (is_pi_over_2_bpsk_enabled == 0) ? 24:12;
|
|
|
|
if (nrofSymbols == 4) {
|
|
E = (intraSlotFrequencyHopping == 0)?(E_init*(nrofSymbols-1)/n_SF_PUCCH_s):((E_init*(nrofSymbols-1)/n_SF_PUCCH_s));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("format 4 nrofSymbols =4 and E_init=%d,E=%d\n",E_init,E);
|
|
#endif
|
|
}
|
|
|
|
if (nrofSymbols > 4) {
|
|
E = E_init*(nrofSymbols-2)/n_SF_PUCCH_s;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("format 4 nrofSymbols >4 and E_init=%d,E = %d\n",E_init,E);
|
|
#endif
|
|
}
|
|
|
|
if (nrofSymbols > 9) {
|
|
E = (add_dmrs == 0)?(E_init*(nrofSymbols-2)/n_SF_PUCCH_s):((E_init*(nrofSymbols-4)/n_SF_PUCCH_s));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("format 4 nrofSymbols >9 and E_init=%d,E = %d\n",E_init,E);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
*M_bit = E;
|
|
//int I_seg;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t\t [nr_uci_encoding] start function with fmt=%d, encoding A=%d bits into M_bit=%d (where nrofSymbols=%d,nrofPRB=%d)\n",fmt,A,*M_bit,nrofSymbols,nrofPRB);
|
|
#endif
|
|
|
|
if (A<=11) {
|
|
// procedure in subclause 6.3.1.2.2 (UCI encoded by channel coding of small block lengths -> subclause 6.3.1.3.2)
|
|
// CRC bits are not attached, and coding small block lengths (subclause 5.3.3)
|
|
uint64_t b0= encodeSmallBlock((uint16_t*)&payload,A);
|
|
// repetition for rate-matching up to 16 PRB
|
|
b[0] = b0 | (b0<<32);
|
|
b[1] = b[0];
|
|
b[2] = b[0];
|
|
b[3] = b[0];
|
|
b[4] = b[0];
|
|
b[5] = b[0];
|
|
b[6] = b[0];
|
|
b[7] = b[0];
|
|
AssertFatal(nrofPRB<=16,"Number of PRB >16\n");
|
|
} else if (A>=12) {
|
|
AssertFatal(A<65,"Polar encoding not supported yet for UCI with more than 64 bits\n");
|
|
t_nrPolar_params *currentPtr = nr_polar_params(NR_POLAR_UCI_PUCCH_MESSAGE_TYPE,
|
|
A,
|
|
nrofPRB,
|
|
1,
|
|
NULL);
|
|
polar_encoder_fast(&payload, b, 0,0,currentPtr);
|
|
}
|
|
|
|
}
|
|
//#if 0
|
|
void nr_generate_pucch2(PHY_VARS_NR_UE *ue,
|
|
int32_t **txdataF,
|
|
NR_DL_FRAME_PARMS *frame_parms,
|
|
int16_t amp,
|
|
int nr_slot_tx,
|
|
fapi_nr_ul_config_pucch_pdu *pucch_pdu) {
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch2] start function at slot(nr_slot_tx)=%d with payload=%lu and nr_bit=%d\n",nr_slot_tx, pucch_pdu->payload, pucch_pdu->n_bit);
|
|
#endif
|
|
// b is the block of bits transmitted on the physical channel after payload coding
|
|
uint64_t b[16]; // limit to 1024-bit encoded length
|
|
// M_bit is the number of bits of block b (payload after encoding)
|
|
uint16_t M_bit;
|
|
nr_uci_encoding(pucch_pdu->payload,
|
|
pucch_pdu->n_bit,
|
|
2,0,
|
|
pucch_pdu->nr_of_symbols,
|
|
pucch_pdu->prb_size,
|
|
1,0,0,&b[0],&M_bit);
|
|
/*
|
|
* Implementing TS 38.211
|
|
* Subclauses 6.3.2.5.1 Scrambling (PUCCH format 2)
|
|
* The block of bits b(0),..., b(M_bit-1 ), where M_bit is the number of bits transmitted on the physical channel,
|
|
* shall be scrambled prior to modulation,
|
|
* resulting in a block of scrambled bits btilde(0),...,btilde(M_bit-1) according to
|
|
* btilde(i)=(b(i)+c(i))mod 2
|
|
* where the scrambling sequence c(i) is given by clause 5.2.1.
|
|
* The scrambling sequence generator shall be initialized with c_init=nRNTI*2^15+n_id
|
|
* n_id = {0,1,...,1023} equals the higher-layer parameter Data-scrambling-Identity if configured
|
|
* n_id = N_ID_cell if higher layer parameter not configured
|
|
*/
|
|
uint8_t *btilde = malloc(sizeof(int8_t)*M_bit);
|
|
// rnti is given by the C-RNTI
|
|
uint16_t rnti=pucch_pdu->rnti;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch2] rnti = %d ,\n",rnti);
|
|
#endif
|
|
/*
|
|
* Implementing TS 38.211 Subclause 6.3.2.5.1 scrambling format 2
|
|
*/
|
|
nr_pucch2_3_4_scrambling(M_bit,rnti,pucch_pdu->data_scrambling_id,&b[0],btilde);
|
|
/*
|
|
* Implementing TS 38.211 Subclause 6.3.2.5.2 modulation format 2
|
|
* btilde shall be modulated as described in subclause 5.1 using QPSK
|
|
* resulting in a block of complex-valued modulation symbols d(0),...,d(m_symbol) where m_symbol=M_bit/2
|
|
*/
|
|
//#define ONE_OVER_SQRT2_S 23171 // 32767/sqrt(2) = 23170 (ONE_OVER_SQRT2)
|
|
// complex-valued symbol d(0)
|
|
int16_t *d_re = malloc(sizeof(int16_t)*M_bit);
|
|
int16_t *d_im = malloc(sizeof(int16_t)*M_bit);
|
|
uint16_t m_symbol = (M_bit%2==0) ? M_bit/2 : floor(M_bit/2)+1;
|
|
|
|
for (int i=0; i < m_symbol; i++) { // QPSK modulation subclause 5.1.3
|
|
if (((btilde[2*i]&1)==0) && ((btilde[(2*i)+1]&1)==0)) {
|
|
d_re[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((btilde[2*i]&1)==0) && ((btilde[(2*i)+1]&1)==1)) {
|
|
d_re[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((btilde[2*i]&1)==1) && ((btilde[(2*i)+1]&1)==0)) {
|
|
d_re[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((btilde[2*i]&1)==1) && ((btilde[(2*i)+1]&1)==1)) {
|
|
d_re[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch2] modulation of bit pair btilde(%d,%d), m_symbol=%d, d(%d)=(%d,%d)\n",(btilde[2*i]&1),(btilde[(2*i)+1]&1),m_symbol,i,d_re[i],d_im[i]);
|
|
#endif
|
|
}
|
|
|
|
/*
|
|
* Implementing TS 38.211 Subclause 6.3.2.5.3 Mapping to physical resources
|
|
*/
|
|
//int32_t *txptr;
|
|
uint32_t re_offset=0;
|
|
uint32_t x1, x2, s=0;
|
|
int i=0;
|
|
int m=0;
|
|
uint8_t startingSymbolIndex = pucch_pdu->start_symbol_index;
|
|
uint16_t startingPRB = pucch_pdu->prb_start + pucch_pdu->bwp_start;
|
|
|
|
for (int l=0; l<pucch_pdu->nr_of_symbols; l++) {
|
|
// c_init calculation according to TS38.211 subclause
|
|
x2 = (((1<<17)*((14*nr_slot_tx) + (l+startingSymbolIndex) + 1)*((2*pucch_pdu->dmrs_scrambling_id) + 1)) + (2*pucch_pdu->dmrs_scrambling_id))%(1U<<31);
|
|
|
|
int reset = 1;
|
|
for (int ii=0; ii<=(startingPRB>>2); ii++) {
|
|
s = lte_gold_generic(&x1, &x2, reset);
|
|
reset = 0;
|
|
}
|
|
m = 0;
|
|
for (int rb=0; rb<pucch_pdu->prb_size; rb++) {
|
|
//startingPRB = startingPRB + rb;
|
|
if (((rb+startingPRB) < (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 0)) { // if number RBs in bandwidth is even and current PRB is lower band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*(rb+startingPRB)) + frame_parms->first_carrier_offset;
|
|
}
|
|
|
|
if (((rb+startingPRB) >= (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 0)) { // if number RBs in bandwidth is even and current PRB is upper band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*((rb+startingPRB)-(frame_parms->N_RB_DL>>1)));
|
|
}
|
|
|
|
if (((rb+startingPRB) < (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB is lower band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*(rb+startingPRB)) + frame_parms->first_carrier_offset;
|
|
}
|
|
|
|
if (((rb+startingPRB) > (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB is upper band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*((rb+startingPRB)-(frame_parms->N_RB_DL>>1))) + 6;
|
|
}
|
|
|
|
if (((rb+startingPRB) == (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB contains DC
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*(rb+startingPRB)) + frame_parms->first_carrier_offset;
|
|
}
|
|
|
|
//txptr = &txdataF[0][re_offset];
|
|
int k=0;
|
|
int kk=0;
|
|
|
|
for (int n=0; n<12; n++) {
|
|
if ((n==6) && ((rb+startingPRB) == (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) {
|
|
// if number RBs in bandwidth is odd and current PRB contains DC, we need to recalculate the offset when n=6 (for second half PRB)
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size);
|
|
}
|
|
|
|
if (n%3 != 1) { // mapping PUCCH according to TS38.211 subclause 6.3.2.5.3
|
|
((int16_t *)&txdataF[0][re_offset])[0] = d_re[i+k];
|
|
((int16_t *)&txdataF[0][re_offset])[1] = d_im[i+k];
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch2] (n=%d,i=%d) mapping PUCCH to RE \t amp=%d \tofdm_symbol_size=%d \tN_RB_DL=%d \tfirst_carrier_offset=%d \tz_pucch[%d]=txptr(%u)=(x_n(l=%d,n=%d)=(%d,%d))\n",
|
|
n,i,amp,frame_parms->ofdm_symbol_size,frame_parms->N_RB_DL,frame_parms->first_carrier_offset,i+k,re_offset,
|
|
l,n,((int16_t *)&txdataF[0][re_offset])[0],((int16_t *)&txdataF[0][re_offset])[1]);
|
|
#endif
|
|
k++;
|
|
}
|
|
|
|
if (n%3 == 1) { // mapping DM-RS signal according to TS38.211 subclause 6.4.1.3.2
|
|
((int16_t *)&txdataF[0][re_offset])[0] = (int16_t)((int32_t)(amp*ONE_OVER_SQRT2*(1-(2*((uint8_t)((s>>(2*m))&1)))))>>15);
|
|
((int16_t *)&txdataF[0][re_offset])[1] = (int16_t)((int32_t)(amp*ONE_OVER_SQRT2*(1-(2*((uint8_t)((s>>((2*m)+1))&1)))))>>15);
|
|
m++;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch2] (n=%d,i=%d) mapping DM-RS to RE \t amp=%d \tofdm_symbol_size=%d \tN_RB_DL=%d \tfirst_carrier_offset=%d \tz_dm-rs[%d]=txptr(%u)=(x_n(l=%d,n=%d)=(%d,%d))\n",
|
|
n,i,amp,frame_parms->ofdm_symbol_size,frame_parms->N_RB_DL,frame_parms->first_carrier_offset,i+kk,re_offset,
|
|
l,n,((int16_t *)&txdataF[0][re_offset])[0],((int16_t *)&txdataF[0][re_offset])[1]);
|
|
#endif
|
|
kk++;
|
|
}
|
|
|
|
re_offset++;
|
|
}
|
|
|
|
i+=8;
|
|
|
|
if ((m&((1<<4)-1))==0) {
|
|
s = lte_gold_generic(&x1, &x2, 0);
|
|
m = 0;
|
|
}
|
|
}
|
|
}
|
|
free(d_re);
|
|
free(d_im);
|
|
free(btilde);
|
|
}
|
|
//#if 0
|
|
void nr_generate_pucch3_4(PHY_VARS_NR_UE *ue,
|
|
int32_t **txdataF,
|
|
NR_DL_FRAME_PARMS *frame_parms,
|
|
int16_t amp,
|
|
int nr_slot_tx,
|
|
fapi_nr_ul_config_pucch_pdu *pucch_pdu) {
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] start function at slot(nr_slot_tx)=%d with payload=%lu and nr_bit=%d\n", nr_slot_tx, pucch_pdu->payload, pucch_pdu->n_bit);
|
|
#endif
|
|
// b is the block of bits transmitted on the physical channel after payload coding
|
|
uint64_t b[16];
|
|
// M_bit is the number of bits of block b (payload after encoding)
|
|
uint16_t M_bit;
|
|
// parameter PUCCH-F4-preDFT-OCC-length set of {2,4} -> to use table -1 or -2
|
|
// in format 4, n_SF_PUCCH_s = {2,4}, provided by higher layer parameter PUCCH-F4-preDFT-OCC-length (in format 3 n_SF_PUCCH_s=1)
|
|
uint8_t n_SF_PUCCH_s;
|
|
if (pucch_pdu->format_type == 3)
|
|
n_SF_PUCCH_s = 1;
|
|
else
|
|
n_SF_PUCCH_s = pucch_pdu->pre_dft_occ_len;
|
|
uint8_t is_pi_over_2_bpsk_enabled = pucch_pdu->pi_2bpsk;
|
|
// Intra-slot frequency hopping shall be assumed when the higher-layer parameter intraSlotFrequencyHopping is provided,
|
|
// regardless of whether the frequency-hop distance is zero or not,
|
|
// otherwise no intra-slot frequency hopping shall be assumed
|
|
//uint8_t PUCCH_Frequency_Hopping = 0 ; // from higher layers
|
|
uint8_t intraSlotFrequencyHopping = 0;
|
|
|
|
if (pucch_pdu->freq_hop_flag) {
|
|
intraSlotFrequencyHopping=1;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] intraSlotFrequencyHopping=%d \n",intraSlotFrequencyHopping);
|
|
#endif
|
|
}
|
|
|
|
uint8_t nrofSymbols = pucch_pdu->nr_of_symbols;
|
|
uint16_t nrofPRB = pucch_pdu->prb_size;
|
|
uint16_t startingPRB = pucch_pdu->prb_start + pucch_pdu->bwp_start;
|
|
uint8_t add_dmrs = pucch_pdu->add_dmrs_flag;
|
|
|
|
nr_uci_encoding(pucch_pdu->payload,
|
|
pucch_pdu->n_bit,
|
|
pucch_pdu->format_type,
|
|
is_pi_over_2_bpsk_enabled,
|
|
nrofSymbols,nrofPRB,
|
|
n_SF_PUCCH_s,
|
|
intraSlotFrequencyHopping,
|
|
add_dmrs,
|
|
&b[0],&M_bit);
|
|
/*
|
|
* Implementing TS 38.211
|
|
* Subclauses 6.3.2.6.1 Scrambling (PUCCH formats 3 and 4)
|
|
* The block of bits b(0),..., b(M_bit-1 ), where M_bit is the number of bits transmitted on the physical channel,
|
|
* shall be scrambled prior to modulation,
|
|
* resulting in a block of scrambled bits btilde(0),...,btilde(M_bit-1) according to
|
|
* btilde(i)=(b(i)+c(i))mod 2
|
|
* where the scrambling sequence c(i) is given by clause 5.2.1.
|
|
* The scrambling sequence generator shall be initialized with c_init=nRNTI*2^15+n_id
|
|
* n_id = {0,1,...,1023} equals the higher-layer parameter Data-scrambling-Identity if configured
|
|
* n_id = N_ID_cell if higher layer parameter not configured
|
|
*/
|
|
uint8_t *btilde = malloc(sizeof(int8_t)*M_bit);
|
|
// rnti is given by the C-RNTI
|
|
uint16_t rnti=pucch_pdu->rnti, n_id=0;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] rnti = %d ,\n",rnti);
|
|
#endif
|
|
/*
|
|
* Implementing TS 38.211 Subclause 6.3.2.6.1 scrambling formats 3 and 4
|
|
*/
|
|
nr_pucch2_3_4_scrambling(M_bit,rnti,n_id,&b[0],btilde);
|
|
/*
|
|
* Implementing TS 38.211 Subclause 6.3.2.6.2 modulation formats 3 and 4
|
|
*
|
|
* Subclause 5.1.1 PI/2-BPSK
|
|
* Subclause 5.1.3 QPSK
|
|
*/
|
|
// complex-valued symbol d(0)
|
|
int16_t *d_re = malloc(sizeof(int16_t)*M_bit);
|
|
int16_t *d_im = malloc(sizeof(int16_t)*M_bit);
|
|
uint16_t m_symbol = (M_bit%2==0) ? M_bit/2 : floor(M_bit/2)+1;
|
|
|
|
if (is_pi_over_2_bpsk_enabled == 0) {
|
|
// using QPSK if PUCCH format 3,4 and pi/2-BPSK is not configured, according to subclause 6.3.2.6.2
|
|
for (int i=0; i < m_symbol; i++) { // QPSK modulation subclause 5.1.3
|
|
if (((btilde[2*i]&1)==0) && ((btilde[(2*i)+1]&1)==0)) {
|
|
d_re[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((btilde[2*i]&1)==0) && ((btilde[(2*i)+1]&1)==1)) {
|
|
d_re[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((btilde[2*i]&1)==1) && ((btilde[(2*i)+1]&1)==0)) {
|
|
d_re[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((btilde[2*i]&1)==1) && ((btilde[(2*i)+1]&1)==1)) {
|
|
d_re[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] modulation QPSK of bit pair btilde(%d,%d), m_symbol=%d, d(%d)=(%d,%d)\n",(btilde[2*i]&1),(btilde[(2*i)+1]&1),m_symbol,i,d_re[i],d_im[i]);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
if (is_pi_over_2_bpsk_enabled == 1) {
|
|
// using PI/2-BPSK if PUCCH format 3,4 and pi/2-BPSK is configured, according to subclause 6.3.2.6.2
|
|
m_symbol = M_bit;
|
|
|
|
for (int i=0; i<m_symbol; i++) { // PI/2-BPSK modulation subclause 5.1.1
|
|
if (((btilde[i]&1)==0) && (i%2 == 0)) {
|
|
d_re[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((btilde[i]&1)==0) && (i%2 == 1)) {
|
|
d_re[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((btilde[i]&1)==1) && (i%2 == 0)) {
|
|
d_re[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
if (((btilde[i]&1)==1) && (i%2 == 1)) {
|
|
d_re[i] = (int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
d_im[i] = -(int16_t)(((int32_t)amp*ONE_OVER_SQRT2)>>15);
|
|
}
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] modulation PI/2-BPSK of bit btilde(%d), m_symbol=%d, d(%d)=(%d,%d)\n",(btilde[i]&1),m_symbol,i,d_re[i],d_im[i]);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
/*
|
|
* Implementing Block-wise spreading subclause 6.3.2.6.3
|
|
*/
|
|
// number of PRBs per PUCCH, provided by higher layers parameters PUCCH-F2-number-of-PRBs or PUCCH-F3-number-of-PRBs (for format 4, it is equal to 1)
|
|
// for PUCCH 3 -> nrofPRBs = (2^alpa2 * 3^alpha3 * 5^alpha5)
|
|
// for PUCCH 4 -> nrofPRBs = 1
|
|
// uint8_t nrofPRBs;
|
|
// number of symbols, provided by higher layers parameters PUCCH-F0-F2-number-of-symbols or PUCCH-F1-F3-F4-number-of-symbols
|
|
// uint8_t nrofSymbols;
|
|
// complex-valued symbol d(0)
|
|
int16_t *y_n_re = malloc(sizeof(int16_t)*4*M_bit); // 4 is the maximum number n_SF_PUCCH_s, so is the maximunm size of y_n
|
|
int16_t *y_n_im = malloc(sizeof(int16_t)*4*M_bit);
|
|
// Re part orthogonal sequences w_n(k) for PUCCH format 4 when N_SF_PUCCH4 = 2 (Table 6.3.2.6.3-1)
|
|
// k={0,..11} n={0,1,2,3}
|
|
// parameter PUCCH-F4-preDFT-OCC-index set of {0,1,2,3} -> n
|
|
uint16_t table_6_3_2_6_3_1_Wn_Re[2][12] = {{1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},
|
|
{1, 1, 1, 1, 1, 1,-1,-1,-1,-1,-1,-1}
|
|
};
|
|
// Im part orthogonal sequences w_n(k) for PUCCH format 4 when N_SF_PUCCH4 = 2 (Table 6.3.2.6.3-1)
|
|
// k={0,..11} n={0,1}
|
|
uint16_t table_6_3_2_6_3_1_Wn_Im[2][12] = {{0,0,0,0,0,0,0,0,0,0,0,0},
|
|
{0,0,0,0,0,0,0,0,0,0,0,0}
|
|
};
|
|
// Re part orthogonal sequences w_n(k) for PUCCH format 4 when N_SF_PUCCH4 = 4 (Table 6.3.2.6.3-2)
|
|
// k={0,..11} n={0,1,2.3}
|
|
uint16_t table_6_3_2_6_3_2_Wn_Re[4][12] = {{1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},
|
|
{1, 1, 1, 0, 0, 0,-1,-1,-1, 0, 0, 0},
|
|
{1, 1, 1,-1,-1,-1, 1, 1, 1,-1,-1,-1},
|
|
{1, 1, 1, 0, 0, 0,-1,-1,-1, 0, 0, 0}
|
|
};
|
|
// Im part orthogonal sequences w_n(k) for PUCCH format 4 when N_SF_PUCCH4 = 4 (Table 6.3.2.6.3-2)
|
|
// k={0,..11} n={0,1,2,3}
|
|
uint16_t table_6_3_2_6_3_2_Wn_Im[4][12] = {{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
|
|
{0, 0, 0,-1,-1,-1, 0, 0, 0, 1, 1, 1},
|
|
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
|
|
{0, 0, 0, 1, 1, 1, 0, 0, 0,-1,-1,-1}
|
|
};
|
|
|
|
uint8_t occ_Index = pucch_pdu->pre_dft_occ_idx; // higher layer parameter occ-Index
|
|
|
|
//occ_Index = 1; //only for testing purposes; to be removed FIXME!!!
|
|
if (pucch_pdu->format_type == 3) { // no block-wise spreading for format 3
|
|
|
|
for (int l=0; l < floor(m_symbol/(12*nrofPRB)); l++) {
|
|
for (int k=0; k < (12*nrofPRB); k++) {
|
|
y_n_re[l*(12*nrofPRB)+k] = d_re[l*(12*nrofPRB)+k];
|
|
y_n_im[l*(12*nrofPRB)+k] = d_im[l*(12*nrofPRB)+k];
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] block-wise spreading for format 3 (no block-wise spreading): (l,k)=(%d,%d)\ty_n(%d) = \t(d_re=%d, d_im=%d)\n",
|
|
l,k,l*(12*nrofPRB)+k,d_re[l*(12*nrofPRB)+k],d_im[l*(12*nrofPRB)+k]);
|
|
#endif
|
|
}
|
|
}
|
|
}
|
|
|
|
if (pucch_pdu->format_type == 4) {
|
|
nrofPRB = 1;
|
|
|
|
for (int l=0; l < floor((n_SF_PUCCH_s*m_symbol)/(12*nrofPRB)); l++) {
|
|
for (int k=0; k < (12*nrofPRB); k++) {
|
|
if (n_SF_PUCCH_s == 2) {
|
|
y_n_re[l*(12*nrofPRB)+k] = (uint16_t)(((uint32_t)d_re[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)] * table_6_3_2_6_3_1_Wn_Re[occ_Index][k])
|
|
- ((uint32_t)d_im[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)] * table_6_3_2_6_3_1_Wn_Im[occ_Index][k]));
|
|
y_n_im[l*(12*nrofPRB)+k] = (uint16_t)(((uint32_t)d_im[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)] * table_6_3_2_6_3_1_Wn_Re[occ_Index][k])
|
|
+ ((uint32_t)d_re[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)] * table_6_3_2_6_3_1_Wn_Im[occ_Index][k]));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] block-wise spreading for format 4 (n_SF_PUCCH_s 2) (occ_Index=%d): (l,k)=(%d,%d)\ty_n(%d) = \t(d_re=%d, d_im=%d)\n",
|
|
occ_Index,l,k,l*(12*nrofPRB)+k,y_n_re[l*(12*nrofPRB)+k],y_n_im[l*(12*nrofPRB)+k]);
|
|
// printf("\t\t d_re[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)] = %d\n",d_re[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)]);
|
|
// printf("\t\t d_im[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)] = %d\n",d_im[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)]);
|
|
// printf("\t\t table_6_3_2_6_3_1_Wn_Re[%d][%d] = %d\n",occ_Index,k,table_6_3_2_6_3_1_Wn_Re[occ_Index][k]);
|
|
// printf("\t\t table_6_3_2_6_3_1_Wn_Im[%d][%d] = %d\n",occ_Index,k,table_6_3_2_6_3_1_Wn_Im[occ_Index][k]);
|
|
#endif
|
|
}
|
|
|
|
if (n_SF_PUCCH_s == 4) {
|
|
y_n_re[l*(12*nrofPRB)+k] = (uint16_t)(((uint32_t)d_re[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)] * table_6_3_2_6_3_2_Wn_Re[occ_Index][k])
|
|
- ((uint32_t)d_im[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)] * table_6_3_2_6_3_2_Wn_Im[occ_Index][k]));
|
|
y_n_im[l*(12*nrofPRB)+k] = (uint16_t)(((uint32_t)d_im[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)] * table_6_3_2_6_3_2_Wn_Re[occ_Index][k])
|
|
+ ((uint32_t)d_re[l*(12*nrofPRB/n_SF_PUCCH_s)+k%(12*nrofPRB/n_SF_PUCCH_s)] * table_6_3_2_6_3_2_Wn_Im[occ_Index][k]));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] block-wise spreading for format 4 (n_SF_PUCCH_s 4) (occ_Index=%d): (l,k)=(%d,%d)\ty_n(%d) = \t(d_re=%d, d_im=%d)\n",
|
|
occ_Index,l,k,l*(12*nrofPRB)+k,y_n_re[l*(12*nrofPRB)+k],y_n_im[l*(12*nrofPRB)+k]);
|
|
#endif
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
/*
|
|
* Implementing Transform pre-coding subclause 6.3.2.6.4
|
|
*/
|
|
int16_t *z_re = malloc(sizeof(int16_t)*4*M_bit); // 4 is the maximum number n_SF_PUCCH_s
|
|
int16_t *z_im = malloc(sizeof(int16_t)*4*M_bit);
|
|
#define M_PI 3.14159265358979323846 // pi
|
|
|
|
//int16_t inv_sqrt_nrofPRBs = (int16_t)round(32767/sqrt(12*nrofPRB));
|
|
for (int l=0; l<floor((n_SF_PUCCH_s*m_symbol)/(12*nrofPRB)); l++) {
|
|
for (int k=0; k<(12*nrofPRB); k++) {
|
|
z_re[l*(12*nrofPRB)+k] = 0;
|
|
z_im[l*(12*nrofPRB)+k] = 0;
|
|
|
|
// int16_t z_re_tmp[240] = {0};
|
|
// int16_t z_im_tmp[240] = {0};
|
|
for (int m=0; m<(12*nrofPRB); m++) {
|
|
//z_re[l*(12*nrofPRB)+k] = y_n_re[l*(12*nrofPRB)+m] * (int16_t)(round(32767*cos((2*M_PI*m*k)/(12*nrofPRB))));
|
|
// z_re_tmp[m] = (int16_t)(((int32_t)round(32767/sqrt(12*nrofPRB))*(int16_t)((((int32_t)y_n_re[l*(12*nrofPRB)+m] * (int16_t)round(32767 * cos(2*M_PI*m*k/(12*nrofPRB))))>>15)
|
|
// + (((int32_t)y_n_im[l*(12*nrofPRB)+m] * (int16_t)round(32767 * sin(2*M_PI*m*k/(12*nrofPRB))))>>15)))>>15);
|
|
// z_im_tmp[m] = (int16_t)(((int32_t)round(32767/sqrt(12*nrofPRB))*(int16_t)((((int32_t)y_n_im[l*(12*nrofPRB)+m] * (int16_t)round(32767 * cos(2*M_PI*m*k/(12*nrofPRB))))>>15)
|
|
// - (((int32_t)y_n_re[l*(12*nrofPRB)+m] * (int16_t)round(32767 * sin(2*M_PI*m*k/(12*nrofPRB))))>>15)))>>15);
|
|
z_re[l*(12*nrofPRB)+k] = z_re[l*(12*nrofPRB)+k]
|
|
+ (int16_t)(((int32_t)round(32767/sqrt(12*nrofPRB))*(int16_t)((((int32_t)y_n_re[l*(12*nrofPRB)+m] * (int16_t)round(32767 * cos(2*M_PI*m*k/(12*nrofPRB))))>>15)
|
|
+ (((int32_t)y_n_im[l*(12*nrofPRB)+m] * (int16_t)round(32767 * sin(2*M_PI*m*k/(12*nrofPRB))))>>15)))>>15);
|
|
z_im[l*(12*nrofPRB)+k] = z_im[l*(12*nrofPRB)+k]
|
|
+ (int16_t)(((int32_t)round(32767/sqrt(12*nrofPRB))*(int16_t)((((int32_t)y_n_im[l*(12*nrofPRB)+m] * (int16_t)round(32767 * cos(2*M_PI*m*k/(12*nrofPRB))))>>15)
|
|
- (((int32_t)y_n_re[l*(12*nrofPRB)+m] * (int16_t)round(32767 * sin(2*M_PI*m*k/(12*nrofPRB))))>>15)))>>15);
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
// printf("\t\t z_re_tmp[%d] = %d\n",m,z_re_tmp[m]);
|
|
// printf("\t\t z_im_tmp[%d] = %d\n",m,z_im_tmp[m]);
|
|
// printf("\t [nr_generate_pucch3_4] transform precoding for formats 3 and 4: (l,k,m)=(%d,%d,%d)\tz(%d) = \t(%d, %d)\n",
|
|
// l,k,m,l*(12*nrofPRB)+k,z_re[l*(12*nrofPRB)+k],z_im[l*(12*nrofPRB)+k]);
|
|
#endif
|
|
}
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] transform precoding for formats 3 and 4: (l,k)=(%d,%d)\tz(%d) = \t(%d, %d)\n",
|
|
l,k,l*(12*nrofPRB)+k,z_re[l*(12*nrofPRB)+k],z_im[l*(12*nrofPRB)+k]);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
/*
|
|
* Implementing TS 38.211 Subclauses 6.3.2.5.3 and 6.3.2.6.5 Mapping to physical resources
|
|
*/
|
|
// the value of u,v (delta always 0 for PUCCH) has to be calculated according to TS 38.211 Subclause 6.3.2.2.1
|
|
uint8_t u=0,v=0;//,delta=0;
|
|
// if frequency hopping is disabled, intraSlotFrequencyHopping is not provided
|
|
// n_hop = 0
|
|
// if frequency hopping is enabled, intraSlotFrequencyHopping is provided
|
|
// n_hop = 0 for first hop
|
|
// n_hop = 1 for second hop
|
|
uint8_t n_hop = 0;
|
|
// lnormal is the OFDM symbol number in the PUCCH transmission where l=0 corresponds to the first OFDM symbol of the PUCCH transmission
|
|
//uint8_t lnormal = 0 ;
|
|
// lprime is the index of the OFDM symbol in the slot that corresponds to the first OFDM symbol of the PUCCH transmission in the slot given by [5, TS 38.213]
|
|
//uint8_t lprime = startingSymbolIndex;
|
|
// m0 is the cyclic shift index calculated depending on the Orthogonal sequence index n, according to table 6.4.1.3.3.1-1 from TS 38.211 subclause 6.4.1.3.3.1
|
|
uint8_t m0;
|
|
uint8_t mcs=0;
|
|
|
|
if (pucch_pdu->format_type == 3) m0 = 0;
|
|
|
|
if (pucch_pdu->format_type == 4) {
|
|
if (n_SF_PUCCH_s == 2) {
|
|
m0 = (occ_Index == 0) ? 0 : 6;
|
|
}
|
|
|
|
if (n_SF_PUCCH_s == 4) {
|
|
m0 = (occ_Index == 3) ? 9 : ((occ_Index == 2) ? 3 : ((occ_Index == 1) ? 6 : 0));
|
|
}
|
|
}
|
|
|
|
double alpha;
|
|
uint8_t N_ZC = 12*nrofPRB;
|
|
int16_t *r_u_v_base_re = malloc(sizeof(int16_t)*12*nrofPRB);
|
|
int16_t *r_u_v_base_im = malloc(sizeof(int16_t)*12*nrofPRB);
|
|
//int16_t *r_u_v_alpha_delta_re = malloc(sizeof(int16_t)*12*nrofPRB);
|
|
//int16_t *r_u_v_alpha_delta_im = malloc(sizeof(int16_t)*12*nrofPRB);
|
|
// Next we proceed to mapping to physical resources according to TS 38.211, subclause 6.3.2.6.5 dor PUCCH formats 3 and 4 and subclause 6.4.1.3.3.2 for DM-RS
|
|
//int32_t *txptr;
|
|
uint32_t re_offset=0;
|
|
//uint32_t x1, x2, s=0;
|
|
// intraSlotFrequencyHopping
|
|
// uint8_t intraSlotFrequencyHopping = 0;
|
|
uint8_t table_6_4_1_3_3_2_1_dmrs_positions[11][14] = {
|
|
{(intraSlotFrequencyHopping==0)?0:1,(intraSlotFrequencyHopping==0)?1:0,(intraSlotFrequencyHopping==0)?0:1,0,0,0,0,0,0,0,0,0,0,0}, // PUCCH length = 4
|
|
{1,0,0,1,0,0,0,0,0,0,0,0,0,0}, // PUCCH length = 5
|
|
{0,1,0,0,1,0,0,0,0,0,0,0,0,0}, // PUCCH length = 6
|
|
{0,1,0,0,1,0,0,0,0,0,0,0,0,0}, // PUCCH length = 7
|
|
{0,1,0,0,0,1,0,0,0,0,0,0,0,0}, // PUCCH length = 8
|
|
{0,1,0,0,0,0,1,0,0,0,0,0,0,0}, // PUCCH length = 9
|
|
{0,(add_dmrs==0?0:1),(add_dmrs==0?1:0),(add_dmrs==0?0:1),0,0,(add_dmrs==0?0:1),(add_dmrs==0?1:0),(add_dmrs==0?0:1),0,0,0,0,0}, // PUCCH length = 10
|
|
{0,(add_dmrs==0?0:1),(add_dmrs==0?1:0),(add_dmrs==0?0:1),0,0,(add_dmrs==0?0:1),(add_dmrs==0?1:0),0,(add_dmrs==0?0:1),0,0,0,0}, // PUCCH length = 11
|
|
{0,(add_dmrs==0?0:1),(add_dmrs==0?1:0),0,(add_dmrs==0?0:1),0,0,(add_dmrs==0?0:1),(add_dmrs==0?1:0),0,(add_dmrs==0?0:1),0,0,0}, // PUCCH length = 12
|
|
{0,(add_dmrs==0?0:1),(add_dmrs==0?1:0),0,(add_dmrs==0?0:1),0,0,(add_dmrs==0?0:1),0,(add_dmrs==0?1:0),0,(add_dmrs==0?0:1),0,0}, // PUCCH length = 13
|
|
{0,(add_dmrs==0?0:1),0,(add_dmrs==0?1:0),0,(add_dmrs==0?0:1),0,0,(add_dmrs==0?0:1),0,(add_dmrs==0?1:0),0,(add_dmrs==0?0:1),0} // PUCCH length = 14
|
|
};
|
|
uint16_t k=0;
|
|
|
|
for (int l=0; l<nrofSymbols; l++) {
|
|
if ((intraSlotFrequencyHopping == 1) && (l >= (int)floor(nrofSymbols/2))) n_hop = 1; // n_hop = 1 for second hop
|
|
|
|
pucch_GroupHopping_t pucch_GroupHopping = pucch_pdu->group_hop_flag + (pucch_pdu->sequence_hop_flag<<1);
|
|
nr_group_sequence_hopping(pucch_GroupHopping,pucch_pdu->hopping_id,n_hop,nr_slot_tx,&u,&v); // calculating u and v value
|
|
|
|
// Next we proceed to calculate base sequence for DM-RS signal, according to TS 38.211 subclause 6.4.1.33
|
|
if (nrofPRB >= 3) { // TS 38.211 subclause 5.2.2.1 (Base sequences of length 36 or larger) applies
|
|
int i = 4;
|
|
|
|
while (list_of_prime_numbers[i] < (12*nrofPRB)) i++;
|
|
|
|
N_ZC = list_of_prime_numbers[i+1]; // N_ZC is given by the largest prime number such that N_ZC < (12*nrofPRB)
|
|
double q_base = (N_ZC*(u+1))/31;
|
|
int8_t q = (uint8_t)floor(q_base + (1/2));
|
|
q = ((uint8_t)floor(2*q_base)%2 == 0 ? q+v : q-v);
|
|
|
|
for (int n=0; n<(12*nrofPRB); n++) {
|
|
r_u_v_base_re[n] = (int16_t)(((int32_t)amp*(int16_t)(32767*cos(M_PI*q*(n%N_ZC)*((n%N_ZC)+1)/N_ZC)))>>15);
|
|
r_u_v_base_im[n] = -(int16_t)(((int32_t)amp*(int16_t)(32767*sin(M_PI*q*(n%N_ZC)*((n%N_ZC)+1)/N_ZC)))>>15);
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] generation DM-RS base sequence when nrofPRB=%d >= 3: r_u_v_base[n=%d]=(%d,%d)\n",
|
|
nrofPRB,n,r_u_v_base_re[n],r_u_v_base_im[n]);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
if (nrofPRB == 2) { // TS 38.211 subclause 5.2.2.2 (Base sequences of length less than 36 using table 5.2.2.2-4) applies
|
|
for (int n=0; n<(12*nrofPRB); n++) {
|
|
r_u_v_base_re[n] = (int16_t)(((int32_t)amp*table_5_2_2_2_4_Re[u][n])>>15);
|
|
r_u_v_base_im[n] = (int16_t)(((int32_t)amp*table_5_2_2_2_4_Im[u][n])>>15);
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] generation DM-RS base sequence when nrofPRB=%d == 2: r_u_v_base[n=%d]=(%d,%d)\n",
|
|
nrofPRB,n,r_u_v_base_re[n],r_u_v_base_im[n]);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
if (nrofPRB == 1) { // TS 38.211 subclause 5.2.2.2 (Base sequences of length less than 36 using table 5.2.2.2-2) applies
|
|
for (int n=0; n<(12*nrofPRB); n++) {
|
|
r_u_v_base_re[n] = (int16_t)(((int32_t)amp*table_5_2_2_2_2_Re[u][n])>>15);
|
|
r_u_v_base_im[n] = (int16_t)(((int32_t)amp*table_5_2_2_2_2_Im[u][n])>>15);
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] generation DM-RS base sequence when nrofPRB=%d == 1: r_u_v_base[n=%d]=(%d,%d)\n",
|
|
nrofPRB,n,r_u_v_base_re[n],r_u_v_base_im[n]);
|
|
#endif
|
|
}
|
|
}
|
|
|
|
uint8_t startingSymbolIndex = pucch_pdu->start_symbol_index;
|
|
uint16_t j=0;
|
|
alpha = nr_cyclic_shift_hopping(pucch_pdu->hopping_id,m0,mcs,l,startingSymbolIndex,nr_slot_tx);
|
|
|
|
for (int rb=0; rb<nrofPRB; rb++) {
|
|
if ((intraSlotFrequencyHopping == 1) && (l<floor(nrofSymbols/2))) { // intra-slot hopping enabled, we need to calculate new offset PRB
|
|
startingPRB = startingPRB + pucch_pdu->second_hop_prb;
|
|
}
|
|
|
|
//startingPRB = startingPRB + rb;
|
|
if (((rb+startingPRB) < (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 0)) { // if number RBs in bandwidth is even and current PRB is lower band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*(rb+startingPRB)) + frame_parms->first_carrier_offset;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("1 ");
|
|
#endif
|
|
}
|
|
|
|
if (((rb+startingPRB) >= (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 0)) { // if number RBs in bandwidth is even and current PRB is upper band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*((rb+startingPRB)-(frame_parms->N_RB_DL>>1)));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("2 ");
|
|
#endif
|
|
}
|
|
|
|
if (((rb+startingPRB) < (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB is lower band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*(rb+startingPRB)) + frame_parms->first_carrier_offset;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("3 ");
|
|
#endif
|
|
}
|
|
|
|
if (((rb+startingPRB) > (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB is upper band
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*((rb+startingPRB)-(frame_parms->N_RB_DL>>1))) + 6;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("4 ");
|
|
#endif
|
|
}
|
|
|
|
if (((rb+startingPRB) == (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) { // if number RBs in bandwidth is odd and current PRB contains DC
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size) + (12*(rb+startingPRB)) + frame_parms->first_carrier_offset;
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("5 ");
|
|
#endif
|
|
}
|
|
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("re_offset=%u,(rb+startingPRB)=%d\n",re_offset,(rb+startingPRB));
|
|
#endif
|
|
|
|
//txptr = &txdataF[0][re_offset];
|
|
for (int n=0; n<12; n++) {
|
|
if ((n==6) && ((rb+startingPRB) == (frame_parms->N_RB_DL>>1)) && ((frame_parms->N_RB_DL & 1) == 1)) {
|
|
// if number RBs in bandwidth is odd and current PRB contains DC, we need to recalculate the offset when n=6 (for second half PRB)
|
|
re_offset = ((l+startingSymbolIndex)*frame_parms->ofdm_symbol_size);
|
|
}
|
|
|
|
if (table_6_4_1_3_3_2_1_dmrs_positions[nrofSymbols-4][l] == 0) { // mapping PUCCH according to TS38.211 subclause 6.3.2.5.3
|
|
((int16_t *)&txdataF[0][re_offset])[0] = z_re[n+k];
|
|
((int16_t *)&txdataF[0][re_offset])[1] = z_im[n+k];
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] (l=%d,rb=%d,n=%d,k=%d) mapping PUCCH to RE \t amp=%d \tofdm_symbol_size=%d \tN_RB_DL=%d \tfirst_carrier_offset=%d \tz_pucch[%d]=txptr(%u)=(z(l=%d,n=%d)=(%d,%d))\n",
|
|
l,rb,n,k,amp,frame_parms->ofdm_symbol_size,frame_parms->N_RB_DL,frame_parms->first_carrier_offset,n+k,re_offset,
|
|
l,n,((int16_t *)&txdataF[0][re_offset])[0],((int16_t *)&txdataF[0][re_offset])[1]);
|
|
#endif
|
|
}
|
|
|
|
if (table_6_4_1_3_3_2_1_dmrs_positions[nrofSymbols-4][l] == 1) { // mapping DM-RS signal according to TS38.211 subclause 6.4.1.3.2
|
|
((int16_t *)&txdataF[0][re_offset])[0] = (int16_t)((((int32_t)(32767*cos(alpha*((n+j)%N_ZC)))*r_u_v_base_re[n+j])>>15)
|
|
- (((int32_t)(32767*sin(alpha*((n+j)%N_ZC)))*r_u_v_base_im[n+j])>>15));
|
|
((int16_t *)&txdataF[0][re_offset])[1] = (int16_t)((((int32_t)(32767*cos(alpha*((n+j)%N_ZC)))*r_u_v_base_im[n+j])>>15)
|
|
+ (((int32_t)(32767*sin(alpha*((n+j)%N_ZC)))*r_u_v_base_re[n+j])>>15));
|
|
#ifdef DEBUG_NR_PUCCH_TX
|
|
printf("\t [nr_generate_pucch3_4] (l=%d,rb=%d,n=%d,j=%d) mapping DM-RS to RE \t amp=%d \tofdm_symbol_size=%d \tN_RB_DL=%d \tfirst_carrier_offset=%d \tz_dm-rs[%d]=txptr(%u)=(r_u_v(l=%d,n=%d)=(%d,%d))\n",
|
|
l,rb,n,j,amp,frame_parms->ofdm_symbol_size,frame_parms->N_RB_DL,frame_parms->first_carrier_offset,n+j,re_offset,
|
|
l,n,((int16_t *)&txdataF[0][re_offset])[0],((int16_t *)&txdataF[0][re_offset])[1]);
|
|
#endif
|
|
}
|
|
|
|
re_offset++;
|
|
}
|
|
|
|
if (table_6_4_1_3_3_2_1_dmrs_positions[nrofSymbols-4][l] == 0) k+=12;
|
|
|
|
if (table_6_4_1_3_3_2_1_dmrs_positions[nrofSymbols-4][l] == 1) j+=12;
|
|
}
|
|
}
|
|
free(z_re);
|
|
free(z_im);
|
|
free(btilde);
|
|
}
|
|
|