diff --git a/openair1/PHY/NR_UE_TRANSPORT/nr_prach.c b/openair1/PHY/NR_UE_TRANSPORT/nr_prach.c
index 70240b8468d57af06d499892b3df5a61c8d7e6fb..121a1820615d1b0a167b4fbc5c064749f1702fca 100644
--- a/openair1/PHY/NR_UE_TRANSPORT/nr_prach.c
+++ b/openair1/PHY/NR_UE_TRANSPORT/nr_prach.c
@@ -67,7 +67,7 @@ int32_t generate_nr_prach(PHY_VARS_NR_UE *ue, uint8_t gNB_id, uint8_t slot){
uint8_t Mod_id, fd_occasion, preamble_index, restricted_set, not_found;
uint16_t rootSequenceIndex, prach_fmt_id, NCS, *prach_root_sequence_map, preamble_offset = 0;
uint16_t preamble_shift = 0, preamble_index0, n_shift_ra, n_shift_ra_bar, d_start=INT16_MAX, numshift, N_ZC, u, offset, offset2, first_nonzero_root_idx;
- int16_t prach_tmp[98304*2*4] __attribute__((aligned(32)));
+ int16_t prach_tmp[(4688+4*24576)*4*2] __attribute__((aligned(32)));
int16_t Ncp = 0, amp, *prach, *prach2, *prachF, *Xu;
int32_t Xu_re, Xu_im;
@@ -238,897 +238,254 @@ int32_t generate_nr_prach(PHY_VARS_NR_UE *ue, uint8_t gNB_id, uint8_t slot){
preamble_offset,
first_nonzero_root_idx);
- Xu = (int16_t*)ue->X_u[preamble_offset-first_nonzero_root_idx];
-
- #if defined (PRACH_WRITE_OUTPUT_DEBUG)
- LOG_M("X_u.m", "X_u", (int16_t*)ue->X_u[preamble_offset-first_nonzero_root_idx], N_ZC, 1, 1);
- #endif
-
- /********************************************************
- *
- * In function init_prach_tables:
- * to compute quantized roots of unity ru(n) = 32767 * exp j*[ (2 * PI * n) / N_ZC ]
- *
- * In compute_prach_seq:
- * to calculate Xu = DFT xu = xu (inv_u*k) * Xu[0] (This is a Zadoff-Chou sequence property: DFT ZC sequence is another ZC sequence)
- *
- * In generate_prach:
- * to do the cyclic-shifted DFT by multiplying Xu[k] * ru[k*preamble_shift] as:
- * If X[k] = DFT x(n) -> X_shifted[k] = DFT x(n+preamble_shift) = X[k] * exp -j*[ (2*PI*k*preamble_shift) / N_ZC ]
- *
- *********************************************************/
-
- if (mu==1) {
- switch(fp->samples_per_subframe) {
- case 15360:
- // 10, 15 MHz @ 15.36 Ms/s
- if (prach_sequence_length == 0) {
- if (prach_fmt_id == 0 || prach_fmt_id == 1 || prach_fmt_id == 2)
- dftlen = 12288;
- if (prach_fmt_id == 3)
- dftlen = 3072;
- } else { // 839 sequence
- dftlen = 512;
- }
- break;
-
- case 30720:
- // 20, 25, 30 MHz @ 30.72 Ms/s
- if (prach_sequence_length == 0) {
- if (prach_fmt_id == 0 || prach_fmt_id == 1 || prach_fmt_id == 2)
- dftlen = 24576;
- if (prach_fmt_id == 3)
- dftlen = 6144;
- } else { // 839 sequence
- dftlen = 1024;
- }
- break;
-
- case 46080:
- // 40 MHz @ 46.08 Ms/s
- if (prach_sequence_length == 0) {
- if (prach_fmt_id == 0 || prach_fmt_id == 1 || prach_fmt_id == 2)
- dftlen = 36864;
- if (prach_fmt_id == 3)
- dftlen = 9216;
- } else { // 839 sequence
- dftlen = 1536;
- }
- break;
-
- case 61440:
- // 40, 50, 60 MHz @ 61.44 Ms/s
- if (prach_sequence_length == 0) {
- if (prach_fmt_id == 0 || prach_fmt_id == 1 || prach_fmt_id == 2)
- dftlen = 49152;
- if (prach_fmt_id == 3)
- dftlen = 12288;
- } else { // 839 sequence
- dftlen = 2048;
- }
- break;
-
- case 92160:
- // 50, 60, 70, 80, 90 MHz @ 92.16 Ms/s
- if (prach_sequence_length == 0) {
- if (prach_fmt_id == 0 || prach_fmt_id == 1 || prach_fmt_id == 2)
- dftlen = 73728;
- if (prach_fmt_id == 3)
- dftlen = 18432;
- } else { // 839 sequence
- dftlen = 3072;
- }
- break;
-
- case 122880:
- // 70, 80, 90, 100 MHz @ 122.88 Ms/s
- if (prach_sequence_length == 0) {
- if (prach_fmt_id == 0 || prach_fmt_id == 1 || prach_fmt_id == 2)
- dftlen = 98304;
- if (prach_fmt_id == 3)
- dftlen = 24576;
- } else { // 839 sequence
- dftlen = 4096;
- }
- break;
-
- default:
- AssertFatal(1==0,"sample rate %f MHz not supported for numerology %d\n", fp->samples_per_subframe / 1000.0, mu);
- }
- }
- else if (mu==3) {
- if (fp->threequarter_fs)
- AssertFatal(1==0,"3/4 sampling not supported for numerology %d\n",mu);
-
- if (prach_sequence_length == 0)
- AssertFatal(1==0,"long prach not supported for numerology %d\n",mu);
-
- if (fp->N_RB_UL == 32)
- dftlen=512;
- else if (fp->N_RB_UL == 66)
- dftlen=1024;
- else
- AssertFatal(1==0,"N_RB_UL %d not support for numerology %d\n",fp->N_RB_UL,mu);
- }
-
-
- for (offset=0,offset2=0; offset<N_ZC; offset++,offset2+=preamble_shift) {
-
- if (offset2 >= N_ZC)
- offset2 -= N_ZC;
-
- Xu_re = (((int32_t)Xu[offset<<1]*amp)>>15);
- Xu_im = (((int32_t)Xu[1+(offset<<1)]*amp)>>15);
- prachF[k++]= ((Xu_re*nr_ru[offset2<<1]) - (Xu_im*nr_ru[1+(offset2<<1)]))>>15;
- prachF[k++]= ((Xu_im*nr_ru[offset2<<1]) + (Xu_re*nr_ru[1+(offset2<<1)]))>>15;
-
- if (k==dftlen) k=0;
- }
-
- #if defined (PRACH_WRITE_OUTPUT_DEBUG)
- LOG_M("prachF.m", "prachF", &prachF[1804], 1024, 1, 1);
- LOG_M("Xu.m", "Xu", Xu, N_ZC, 1, 1);
- #endif
-
+ // Ncp and dftlen here is given in terms of T_s wich is 30.72MHz sampling
if (prach_sequence_length == 0) {
-
- AssertFatal(prach_fmt_id < 4, "Illegal PRACH format %d for sequence length 839\n", prach_fmt_id);
-
- // Ncp here is given in terms of T_s wich is 30.72MHz sampling
switch (prach_fmt_id) {
case 0:
Ncp = 3168;
+ dftlen = 24576;
break;
+
case 1:
Ncp = 21024;
+ dftlen = 24576;
break;
+
case 2:
Ncp = 4688;
+ dftlen = 24576;
break;
+
case 3:
Ncp = 3168;
+ dftlen = 6144;
break;
- }
+ default:
+ AssertFatal(1==0, "Illegal PRACH format %d for sequence length 839\n", prach_fmt_id);
+ break;
+ }
} else {
-
switch (prach_fmt_id) {
case 4: //A1
- Ncp = 288/(1<<mu);
+ Ncp = 288 >> mu;
break;
+
case 5: //A2
- Ncp = 576/(1<<mu);
+ Ncp = 576 >> mu;
break;
+
case 6: //A3
- Ncp = 864/(1<<mu);
+ Ncp = 864 >> mu;
break;
+
case 7: //B1
- Ncp = 216/(1<<mu);
- break;
+ Ncp = 216 >> mu;
+ break;
+
/*
case 4: //B2
- Ncp = 360/(1<<mu);
+ Ncp = 360 >> mu;
break;
+
case 5: //B3
- Ncp = 504/(1<<mu);
+ Ncp = 504 >> mu;
break;
*/
+
case 8: //B4
- Ncp = 936/(1<<mu);
+ Ncp = 936 >> mu;
break;
+
case 9: //C0
- Ncp = 1240/(1<<mu);
+ Ncp = 1240 >> mu;
break;
+
case 10: //C2
- Ncp = 2048/(1<<mu);
+ Ncp = 2048 >> mu;
break;
+
default:
AssertFatal(1==0,"Unknown PRACH format ID %d\n", prach_fmt_id);
break;
}
+ dftlen = 2048 >> mu;
}
- #ifdef NR_PRACH_DEBUG
- LOG_D(PHY, "PRACH [UE %d] Ncp %d, dftlen %d \n", Mod_id, Ncp, dftlen);
- #endif
-
//actually what we should be checking here is how often the current prach crosses a 0.5ms boundary. I am not quite sure for which paramter set this would be the case, so I will ignore it for now and just check if the prach starts on a 0.5ms boundary
- uint8_t use_extended_prach_prefix = 0;
if(fp->numerology_index == 0) {
if (prachStartSymbol == 0 || prachStartSymbol == 7)
- use_extended_prach_prefix = 1;
+ Ncp += 16;
}
else {
if (slot%(fp->slots_per_subframe/2)==0 && prachStartSymbol == 0)
- use_extended_prach_prefix = 1;
+ Ncp += 16;
}
- if (mu == 3) {
- switch (fp->samples_per_subframe) {
- case 61440: // 32 PRB case, 61.44 Msps
- Ncp<<=1; //to account for 61.44Mbps
- // This is after cyclic prefix
- prach2 = prach+(Ncp<<1); //times 2 for complex samples
- if (prach_sequence_length == 0)
- AssertFatal(1==0,"no long PRACH for this PRACH size %d\n",fp->N_RB_UL);
- else {
- if (use_extended_prach_prefix)
- Ncp+=32; // 16*kappa, kappa=2 for 61.44Msps
- prach2 = prach+(Ncp<<1); //times 2 for complex samples
- if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
- idft(IDFT_512,prachF,prach2,1);
- // here we have |empty | Prach512 |
- if (prach_fmt_id != 9) {
- memmove(prach2+(512<<1),prach2,(512<<2));
- prach_len = (512*2)+Ncp;
- }
- else prach_len = (512*1)+Ncp;
- memmove(prach,prach+(512<<1),(Ncp<<2));
- // here we have |Prefix | Prach512 | Prach512 (if ! 0xc0) |
- } else if (prach_fmt_id == 5) { // 6x512
- idft(IDFT_512,prachF,prach2,1);
- // here we have |empty | Prach512 |
- memmove(prach2+(512<<1),prach2,(512<<2));
- // here we have |empty | Prach512 | Prach512| empty512 | empty512 |
- memmove(prach2+(512<<2),prach2,(512<<3));
- // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 |
- memmove(prach,prach+(512<<1),(Ncp<<2));
- // here we have |Prefix | Prach512 |
- prach_len = (512*4)+Ncp;
- } else if (prach_fmt_id == 6) { // 6x512
- idft(IDFT_512,prachF,prach2,1);
- // here we have |empty | Prach512 |
- memmove(prach2+(512<<1),prach2,(512<<2));
- // here we have |empty | Prach512 | Prach512| empty512 | empty512 | empty512 | empty512
- memmove(prach2+(512<<2),prach2,(512<<3));
- // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 | empty512 | empty512
- memmove(prach2+(512<<3),prach2,(512<<3));
- // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512
- memmove(prach,prach+(512<<1),(Ncp<<2));
- // here we have |Prefix | Prach512 |
- prach_len = (512*6)+Ncp;
- } else if (prach_fmt_id == 8) { // 12x512
- idft(IDFT_512,prachF,prach2,1);
- // here we have |empty | Prach512 |
- memmove(prach2+(512<<1),prach2,(512<<2));
- // here we have |empty | Prach512 | Prach512| empty512 | empty512 | empty512 | empty512
- memmove(prach2+(512<<2),prach2,(512<<3));
- // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 | empty512 | empty512
- memmove(prach2+(512<<3),prach2,(512<<3));
- // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512
- memmove(prach2+(512<<1)*6,prach2,(512<<2)*6);
- // here we have |empty | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512|
- memmove(prach,prach+(512<<1),(Ncp<<2));
- // here we have |Prefix | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512| Prach512 | Prach512 | Prach512 | Prach512|
- prach_len = (512*12)+Ncp;
- }
- }
- break;
+ switch(fp->samples_per_subframe) {
+ case 7680:
+ // 5 MHz @ 7.68 Ms/s
+ Ncp >>= 2;
+ dftlen >>= 2;
+ break;
- case 122880: // 66 PRB case, 122.88 Msps
- Ncp<<=2; //to account for 122.88Mbps
- // This is after cyclic prefix
- prach2 = prach+(Ncp<<1); //times 2 for complex samples
- if (prach_sequence_length == 0)
- AssertFatal(1==0,"no long PRACH for this PRACH size %d\n",fp->N_RB_UL);
- else {
- if (use_extended_prach_prefix)
- Ncp+=64; // 16*kappa, kappa=4 for 122.88Msps
- prach2 = prach+(Ncp<<1); //times 2 for complex samples
- if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
- idft(IDFT_1024,prachF,prach2,1);
- // here we have |empty | Prach1024 |
- if (prach_fmt_id != 9) {
- memmove(prach2+(1024<<1),prach2,(1024<<2));
- prach_len = (1024*2)+Ncp;
- }
- else prach_len = (1024*1)+Ncp;
- memmove(prach,prach+(1024<<1),(Ncp<<2));
- // here we have |Prefix | Prach1024 | Prach1024 (if ! 0xc0) |
- } else if (prach_fmt_id == 5) { // 6x1024
- idft(IDFT_1024,prachF,prach2,1);
- // here we have |empty | Prach1024 |
- memmove(prach2+(1024<<1),prach2,(1024<<2));
- // here we have |empty | Prach1024 | Prach1024| empty1024 | empty1024 |
- memmove(prach2+(1024<<2),prach2,(1024<<3));
- // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 |
- memmove(prach,prach+(1024<<1),(Ncp<<2));
- // here we have |Prefix | Prach1024 |
- prach_len = (1024*4)+Ncp;
- } else if (prach_fmt_id == 6) { // 6x1024
- idft(IDFT_1024,prachF,prach2,1);
- // here we have |empty | Prach1024 |
- memmove(prach2+(1024<<1),prach2,(1024<<2));
- // here we have |empty | Prach1024 | Prach1024| empty1024 | empty1024 | empty1024 | empty1024
- memmove(prach2+(1024<<2),prach2,(1024<<3));
- // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 | empty1024 | empty1024
- memmove(prach2+(1024<<3),prach2,(1024<<3));
- // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024
- memmove(prach,prach+(1024<<1),(Ncp<<2));
- // here we have |Prefix | Prach1024 |
- prach_len = (1024*6)+Ncp;
- } else if (prach_fmt_id == 8) { // 12x1024
- idft(IDFT_1024,prachF,prach2,1);
- // here we have |empty | Prach1024 |
- memmove(prach2+(1024<<1),prach2,(1024<<2));
- // here we have |empty | Prach1024 | Prach1024| empty1024 | empty1024 | empty1024 | empty1024
- memmove(prach2+(1024<<2),prach2,(1024<<3));
- // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 | empty1024 | empty1024
- memmove(prach2+(1024<<3),prach2,(1024<<3));
- // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024
- memmove(prach2+(1024<<1)*6,prach2,(1024<<2)*6);
- // here we have |empty | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024|
- memmove(prach,prach+(1024<<1),(Ncp<<2));
- // here we have |Prefix | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024| Prach1024 | Prach1024 | Prach1024 | Prach1024|
- prach_len = (1024*12)+Ncp;
- }
- }
- break;
+ case 15360:
+ // 10, 15 MHz @ 15.36 Ms/s
+ Ncp >>= 1;
+ dftlen >>= 1;
+ break;
- default:
- AssertFatal(1==0,"sample rate %f MHz not supported for numerology %d\n", fp->samples_per_subframe / 1000.0, mu);
- }
- } else if (mu == 1) {
- switch (fp->samples_per_subframe) {
- case 15360: // full sampling @ 15.36 Ms/s
- Ncp = Ncp/2; // to account for 15.36 Ms/s
- // This is after cyclic prefix
- prach2 = prach+(2*Ncp); // times 2 for complex samples
- if (prach_sequence_length == 0){
- if (prach_fmt_id == 0) { // 24576 samples @ 30.72 Ms/s, 12288 samples @ 15.36 Ms/s
- idft(IDFT_12288,prachF,prach2,1);
- // here we have | empty | Prach12288 |
- memmove(prach,prach+(12288<<1),(Ncp<<2));
- // here we have | Prefix | Prach12288 |
- prach_len = 12288+Ncp;
- } else if (prach_fmt_id == 1) { // 24576 samples @ 30.72 Ms/s, 12288 samples @ 15.36 Ms/s
- idft(IDFT_12288,prachF,prach2,1);
- // here we have | empty | Prach12288 | empty12288 |
- memmove(prach2+(12288<<1),prach2,(12288<<2));
- // here we have | empty | Prach12288 | Prach12288 |
- memmove(prach,prach+(12288<<2),(Ncp<<2));
- // here we have | Prefix | Prach12288 | Prach12288 |
- prach_len = (12288*2)+Ncp;
- } else if (prach_fmt_id == 2) { // 24576 samples @ 30.72 Ms/s, 12288 samples @ 15.36 Ms/s
- idft(IDFT_12288,prachF,prach2,1);
- // here we have | empty | Prach12288 | empty12288 | empty12288 | empty12288 |
- memmove(prach2+(12288<<1),prach2,(12288<<2));
- // here we have | empty | Prach12288 | Prach12288 | empty12288 | empty12288 |
- memmove(prach2+(12288<<2),prach2,(12288<<3));
- // here we have | empty | Prach12288 | Prach12288 | Prach12288 | Prach12288 |
- memmove(prach,prach+(12288<<3),(Ncp<<2));
- // here we have | Prefix | Prach12288 | Prach12288 | Prach12288 | Prach12288 |
- prach_len = (12288*4)+Ncp;
- } else if (prach_fmt_id == 3) { // 6144 samples @ 30.72 Ms/s, 3072 samples @ 15.36 Ms/s
- idft(IDFT_3072,prachF,prach2,1);
- // here we have | empty | Prach3072 | empty3072 | empty3072 | empty3072 |
- memmove(prach2+(3072<<1),prach2,(3072<<2));
- // here we have | empty | Prach3072 | Prach3072 | empty3072 | empty3072 |
- memmove(prach2+(3072<<2),prach2,(3072<<3));
- // here we have | empty | Prach3072 | Prach3072 | Prach3072 | Prach3072 |
- memmove(prach,prach+(3072<<3),(Ncp<<2));
- // here we have | Prefix | Prach3072 | Prach3072 | Prach3072 | Prach3072 |
- prach_len = (3072*4)+Ncp;
- }
- } else { // short PRACH sequence
- if (use_extended_prach_prefix)
- Ncp += 8; // 16*kappa, kappa=0.5 for 15.36 Ms/s
- prach2 = prach+(2*Ncp); // times 2 for complex samples
- if (prach_fmt_id == 9) {
- idft(IDFT_512,prachF,prach2,1);
- // here we have | empty | Prach512 |
- memmove(prach,prach+(512<<1),(Ncp<<2));
- // here we have | Prefix | Prach512 |
- prach_len = (512*1)+Ncp;
- } else if (prach_fmt_id == 4 || prach_fmt_id == 7) {
- idft(IDFT_512,prachF,prach2,1);
- // here we have | empty | Prach512 | empty512 |
- memmove(prach2+(512<<1),prach2,(512<<2));
- // here we have | empty | Prach512 | Prach512 |
- memmove(prach,prach+(512<<1),(Ncp<<2));
- // here we have | Prefix | Prach512 | Prach512 |
- prach_len = (512*2)+Ncp;
- } else if (prach_fmt_id == 5) { // 4x512
- idft(IDFT_512,prachF,prach2,1);
- // here we have | empty | Prach512 | empty512 | empty512 | empty512 |
- memmove(prach2+(512<<1),prach2,(512<<2));
- // here we have | empty | Prach512 | Prach512 | empty512 | empty512 |
- memmove(prach2+(512<<2),prach2,(512<<3));
- // here we have | empty | Prach512 | Prach512 | Prach512 | Prach512 |
- memmove(prach,prach+(512<<1),(Ncp<<2));
- // here we have | Prefix | Prach512 | Prach512 | Prach512 | Prach512 |
- prach_len = (512*4)+Ncp;
- } else if (prach_fmt_id == 6) { // 6x512
- idft(IDFT_512,prachF,prach2,1);
- // here we have | empty | Prach512 | empty512 | empty512 | empty512 | empty512 | empty512 |
- memmove(prach2+(512<<1),prach2,(512<<2));
- // here we have | empty | Prach512 | Prach512 | empty512 | empty512 | empty512 | empty512 |
- memmove(prach2+(512<<2),prach2,(512<<3));
- // here we have | empty | Prach512 | Prach512 | Prach512 | Prach512 | empty512 | empty512 |
- memmove(prach2+(512<<3),prach2,(512<<3));
- // here we have | empty | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 |
- memmove(prach,prach+(512<<1),(Ncp<<2));
- // here we have | Prefix | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 |
- prach_len = (512*6)+Ncp;
- } else if (prach_fmt_id == 8) { // 12x512
- idft(IDFT_512,prachF,prach2,1);
- // here we have | empty | Prach512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 |
- memmove(prach2+(512<<1),prach2,(512<<2));
- // here we have | empty | Prach512 | Prach512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 |
- memmove(prach2+(512<<2),prach2,(512<<3));
- // here we have | empty | Prach512 | Prach512 | Prach512 | Prach512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 |
- memmove(prach2+(512<<3),prach2,(512<<3));
- // here we have | empty | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | empty512 | empty512 | empty512 | empty512 | empty512 | empty512 |
- memmove(prach2+(512<<1)*6,prach2,(512<<2)*6);
- // here we have | empty | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 |
- memmove(prach,prach+(512<<1),(Ncp<<2));
- // here we have | Prefix | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 | Prach512 |
- prach_len = (512*12)+Ncp;
- }
- }
- break;
+ case 30720:
+ // 20, 25, 30 MHz @ 30.72 Ms/s
+ Ncp = Ncp;
+ dftlen = dftlen;
+ break;
- case 30720: // full sampling @ 30.72 Ms/s
- Ncp = Ncp*1; // to account for 30.72 Ms/s
- // This is after cyclic prefix
- prach2 = prach+(2*Ncp); // times 2 for complex samples
- if (prach_sequence_length == 0){
- if (prach_fmt_id == 0) { // 24576 samples @ 30.72 Ms/s
- idft(IDFT_24576,prachF,prach2,1);
- // here we have | empty | Prach24576 |
- memmove(prach,prach+(24576<<1),(Ncp<<2));
- // here we have | Prefix | Prach24576 |
- prach_len = 24576+Ncp;
- } else if (prach_fmt_id == 1) { // 24576 samples @ 30.72 Ms/s
- idft(IDFT_24576,prachF,prach2,1);
- // here we have | empty | Prach24576 | empty24576 |
- memmove(prach2+(24576<<1),prach2,(24576<<2));
- // here we have | empty | Prach24576 | Prach24576 |
- memmove(prach,prach+(24576<<2),(Ncp<<2));
- // here we have | Prefix | Prach24576 | Prach24576 |
- prach_len = (24576*2)+Ncp;
- } else if (prach_fmt_id == 2) { // 24576 samples @ 30.72 Ms/s
- idft(IDFT_24576,prachF,prach2,1);
- // here we have | empty | Prach24576 | empty24576 | empty24576 | empty24576 |
- memmove(prach2+(24576<<1),prach2,(24576<<2));
- // here we have | empty | Prach24576 | Prach24576 | empty24576 | empty24576 |
- memmove(prach2+(24576<<2),prach2,(24576<<3));
- // here we have | empty | Prach24576 | Prach24576 | Prach24576 | Prach24576 |
- memmove(prach,prach+(24576<<3),(Ncp<<2));
- // here we have | Prefix | Prach24576 | Prach24576 | Prach24576 | Prach24576 |
- prach_len = (24576*4)+Ncp;
- } else if (prach_fmt_id == 3) { // 6144 samples @ 30.72 Ms/s
- idft(IDFT_6144,prachF,prach2,1);
- // here we have | empty | Prach6144 | empty6144 | empty6144 | empty6144 |
- memmove(prach2+(6144<<1),prach2,(6144<<2));
- // here we have | empty | Prach6144 | Prach6144 | empty6144 | empty6144 |
- memmove(prach2+(6144<<2),prach2,(6144<<3));
- // here we have | empty | Prach6144 | Prach6144 | Prach6144 | Prach6144 |
- memmove(prach,prach+(6144<<3),(Ncp<<2));
- // here we have | Prefix | Prach6144 | Prach6144 | Prach6144 | Prach6144 |
- prach_len = (6144*4)+Ncp;
- }
- } else { // short PRACH sequence
- if (use_extended_prach_prefix)
- Ncp += 16; // 16*kappa, kappa=1 for 30.72Msps
- prach2 = prach+(2*Ncp); // times 2 for complex samples
- if (prach_fmt_id == 9) {
- idft(IDFT_1024,prachF,prach2,1);
- // here we have | empty | Prach1024 |
- memmove(prach,prach+(1024<<1),(Ncp<<2));
- // here we have | Prefix | Prach1024 |
- prach_len = (1024*1)+Ncp;
- } else if (prach_fmt_id == 4 || prach_fmt_id == 7) {
- idft(IDFT_1024,prachF,prach2,1);
- // here we have | empty | Prach1024 | empty1024 |
- memmove(prach2+(1024<<1),prach2,(1024<<2));
- // here we have | empty | Prach1024 | Prach1024 |
- memmove(prach,prach+(1024<<1),(Ncp<<2));
- // here we have | Prefix | Prach1024 | Prach1024 |
- prach_len = (1024*2)+Ncp;
- } else if (prach_fmt_id == 5) { // 4x1024
- idft(IDFT_1024,prachF,prach2,1);
- // here we have | empty | Prach1024 | empty1024 | empty1024 | empty1024 |
- memmove(prach2+(1024<<1),prach2,(1024<<2));
- // here we have | empty | Prach1024 | Prach1024 | empty1024 | empty1024 |
- memmove(prach2+(1024<<2),prach2,(1024<<3));
- // here we have | empty | Prach1024 | Prach1024 | Prach1024 | Prach1024 |
- memmove(prach,prach+(1024<<1),(Ncp<<2));
- // here we have | Prefix | Prach1024 | Prach1024 | Prach1024 | Prach1024 |
- prach_len = (1024*4)+Ncp;
- } else if (prach_fmt_id == 6) { // 6x1024
- idft(IDFT_1024,prachF,prach2,1);
- // here we have | empty | Prach1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 |
- memmove(prach2+(1024<<1),prach2,(1024<<2));
- // here we have | empty | Prach1024 | Prach1024 | empty1024 | empty1024 | empty1024 | empty1024 |
- memmove(prach2+(1024<<2),prach2,(1024<<3));
- // here we have | empty | Prach1024 | Prach1024 | Prach1024 | Prach1024 | empty1024 | empty1024 |
- memmove(prach2+(1024<<3),prach2,(1024<<3));
- // here we have | empty | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 |
- memmove(prach,prach+(1024<<1),(Ncp<<2));
- // here we have | Prefix | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 |
- prach_len = (1024*6)+Ncp;
- } else if (prach_fmt_id == 8) { // 12x1024
- idft(IDFT_1024,prachF,prach2,1);
- // here we have | empty | Prach1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 |
- memmove(prach2+(1024<<1),prach2,(1024<<2));
- // here we have | empty | Prach1024 | Prach1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 |
- memmove(prach2+(1024<<2),prach2,(1024<<3));
- // here we have | empty | Prach1024 | Prach1024 | Prach1024 | Prach1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 |
- memmove(prach2+(1024<<3),prach2,(1024<<3));
- // here we have | empty | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 | empty1024 |
- memmove(prach2+(1024<<1)*6,prach2,(1024<<2)*6);
- // here we have | empty | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 |
- memmove(prach,prach+(1024<<1),(Ncp<<2));
- // here we have | Prefix | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 | Prach1024 |
- prach_len = (1024*12)+Ncp;
- }
- }
- break;
+ case 46080:
+ // 40 MHz @ 46.08 Ms/s
+ Ncp = (Ncp*3)/2;
+ dftlen = (dftlen*3)/2;
+ break;
- case 61440: // full sampling @ 61.44 Ms/s
- Ncp = Ncp*2; // to account for 61.44 Ms/s
- // This is after cyclic prefix
- prach2 = prach+(Ncp<<1); //times 2 for complex samples
- if (prach_sequence_length == 0){
- if (prach_fmt_id == 0) { //24576 samples @ 30.72 Ms/s, 49152 samples @ 61.44 Ms/s
- idft(IDFT_49152,prachF,prach2,1);
- // here we have |empty | Prach49152|
- memmove(prach,prach+(49152<<1),(Ncp<<2));
- // here we have |Prefix | Prach49152|
- prach_len = 49152+Ncp;
- } else if (prach_fmt_id == 1) { //24576 samples @ 30.72 Ms/s, 49152 samples @ 61.44 Ms/s
- idft(IDFT_49152,prachF,prach2,1);
- memmove(prach2+(49152<<1),prach2,(49152<<2));
- // here we have |empty | Prach49152 | Prach49152|
- memmove(prach,prach+(49152<<2),(Ncp<<2));
- // here we have |Prefix | Prach49152 | Prach49152|
- prach_len = (49152*2)+Ncp;
- } else if (prach_fmt_id == 2) { //24576 samples @ 30.72 Ms/s, 49152 samples @ 61.44 Ms/s
- idft(IDFT_49152,prachF,prach2,1);
- memmove(prach2+(49152<<1),prach2,(49152<<2));
- // here we have |empty | Prach49152 | Prach49152| empty49152 | empty49152
- memmove(prach2+(49152<<2),prach2,(49152<<3));
- // here we have |empty | Prach49152 | Prach49152| Prach49152 | Prach49152
- memmove(prach,prach+(49152<<3),(Ncp<<2));
- // here we have |Prefix | Prach49152 | Prach49152| Prach49152 | Prach49152
- prach_len = (49152*4)+Ncp;
- } else if (prach_fmt_id == 3) { // 6144 samples @ 30.72 Ms/s, 12288 samples @ 61.44 Ms/s
- idft(IDFT_12288,prachF,prach2,1);
- memmove(prach2+(12288<<1),prach2,(12288<<2));
- // here we have |empty | Prach12288 | Prach12288| empty12288 | empty12288
- memmove(prach2+(12288<<2),prach2,(12288<<3));
- // here we have |empty | Prach12288 | Prach12288| Prach12288 | Prach12288
- memmove(prach,prach+(12288<<3),(Ncp<<2));
- // here we have |Prefix | Prach12288 | Prach12288| Prach12288 | Prach12288
- prach_len = (12288*4)+Ncp;
- }
- } else { // short PRACH sequence
- if (use_extended_prach_prefix)
- Ncp+=32; // 16*kappa, kappa=2 for 61.44Msps
- prach2 = prach+(Ncp<<1); //times 2 for complex samples
- if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
- idft(IDFT_2048,prachF,prach2,1);
- // here we have |empty | Prach2048 |
- if (prach_fmt_id != 9) {
- memmove(prach2+(2048<<1),prach2,(2048<<2));
- prach_len = (2048*2)+Ncp;
- }
- else prach_len = (2048*1)+Ncp;
- memmove(prach,prach+(2048<<1),(Ncp<<2));
- // here we have |Prefix | Prach2048 | Prach2048 (if ! 0xc0) |
- } else if (prach_fmt_id == 5) { // 6x2048
- idft(IDFT_2048,prachF,prach2,1);
- // here we have |empty | Prach2048 |
- memmove(prach2+(2048<<1),prach2,(2048<<2));
- // here we have |empty | Prach2048 | Prach2048| empty2048 | empty2048 |
- memmove(prach2+(2048<<2),prach2,(2048<<3));
- // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 |
- memmove(prach,prach+(2048<<1),(Ncp<<2));
- // here we have |Prefix | Prach2048 |
- prach_len = (2048*4)+Ncp;
- } else if (prach_fmt_id == 6) { // 6x2048
- idft(IDFT_2048,prachF,prach2,1);
- // here we have |empty | Prach2048 |
- memmove(prach2+(2048<<1),prach2,(2048<<2));
- // here we have |empty | Prach2048 | Prach2048| empty2048 | empty2048 | empty2048 | empty2048
- memmove(prach2+(2048<<2),prach2,(2048<<3));
- // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 | empty2048 | empty2048
- memmove(prach2+(2048<<3),prach2,(2048<<3));
- // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048
- memmove(prach,prach+(2048<<1),(Ncp<<2));
- // here we have |Prefix | Prach2048 |
- prach_len = (2048*6)+Ncp;
- } else if (prach_fmt_id == 8) { // 12x2048
- idft(IDFT_2048,prachF,prach2,1);
- // here we have |empty | Prach2048 |
- memmove(prach2+(2048<<1),prach2,(2048<<2));
- // here we have |empty | Prach2048 | Prach2048| empty2048 | empty2048 | empty2048 | empty2048
- memmove(prach2+(2048<<2),prach2,(2048<<3));
- // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 | empty2048 | empty2048
- memmove(prach2+(2048<<3),prach2,(2048<<3));
- // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048
- memmove(prach2+(2048<<1)*6,prach2,(2048<<2)*6);
- // here we have |empty | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048 | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048|
- memmove(prach,prach+(2048<<1),(Ncp<<2));
- // here we have |Prefix | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048 | Prach2048 | Prach2048| Prach2048 | Prach2048 | Prach2048 | Prach2048|
- prach_len = (2048*12)+Ncp;
- }
- }
- break;
+ case 61440:
+ // 40, 50, 60 MHz @ 61.44 Ms/s
+ Ncp <<= 1;
+ dftlen <<= 1;
+ break;
- case 46080: // threequarter sampling @ 46.08 Ms/s
- Ncp = (Ncp*3)/2;
- prach2 = prach+(Ncp<<1);
- if (prach_sequence_length == 0){
- if (prach_fmt_id == 0) {
- idft(IDFT_36864,prachF,prach2,1);
- // here we have |empty | Prach73728|
- memmove(prach,prach+(36864<<1),(Ncp<<2));
- // here we have |Prefix | Prach73728|
- prach_len = (36864*1)+Ncp;
- } else if (prach_fmt_id == 1) {
- idft(IDFT_36864,prachF,prach2,1);
- memmove(prach2+(36864<<1),prach2,(36864<<2));
- // here we have |empty | Prach73728 | Prach73728|
- memmove(prach,prach+(36864<<2),(Ncp<<2));
- // here we have |Prefix | Prach73728 | Prach73728|
- prach_len = (36864*2)+Ncp;
- } else if (prach_fmt_id == 2) {
- idft(IDFT_36864,prachF,prach2,1);
- memmove(prach2+(36864<<1),prach2,(36864<<2));
- // here we have |empty | Prach73728 | Prach73728| empty73728 | empty73728
- memmove(prach2+(36864<<2),prach2,(36864<<3));
- // here we have |empty | Prach73728 | Prach73728| Prach73728 | Prach73728
- memmove(prach,prach+(36864<<3),(Ncp<<2));
- // here we have |Prefix | Prach73728 | Prach73728| Prach73728 | Prach73728
- prach_len = (36864*4)+Ncp;
- } else if (prach_fmt_id == 3) {
- idft(IDFT_9216,prachF,prach2,1);
- memmove(prach2+(9216<<1),prach2,(9216<<2));
- // here we have |empty | Prach9216 | Prach9216| empty9216 | empty9216
- memmove(prach2+(9216<<2),prach2,(9216<<3));
- // here we have |empty | Prach9216 | Prach9216| Prach9216 | Prach9216
- memmove(prach,prach+(9216<<3),(Ncp<<2));
- // here we have |Prefix | Prach9216 | Prach9216| Prach9216 | Prach9216
- prach_len = (9216*4)+Ncp;
- }
- } else { // short sequence
- if (use_extended_prach_prefix)
- Ncp+=24; // 16*kappa, kappa=1.5 for 46.08Msps
- prach2 = prach+(Ncp<<1); //times 2 for complex samples
- if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
- idft(IDFT_1536,prachF,prach2,1);
- // here we have |empty | Prach1536 |
- if (prach_fmt_id != 9) {
- memmove(prach2+(1536<<1),prach2,(1536<<2));
- prach_len = (1536*2)+Ncp;
- } else prach_len = (1536*1)+Ncp;
-
- memmove(prach,prach+(1536<<1),(Ncp<<2));
- // here we have |Prefix | Prach1536 | Prach1536 (if ! 0xc0) |
-
- } else if (prach_fmt_id == 5) { // 6x1536
- idft(IDFT_1536,prachF,prach2,1);
- // here we have |empty | Prach1536 |
- memmove(prach2+(1536<<1),prach2,(1536<<2));
- // here we have |empty | Prach1536 | Prach1536| empty1536 | empty1536 |
- memmove(prach2+(1536<<2),prach2,(1536<<3));
- // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 |
- memmove(prach,prach+(1536<<1),(Ncp<<2));
- // here we have |Prefix | Prach1536 |
- prach_len = (1536*4)+Ncp;
- } else if (prach_fmt_id == 6) { // 6x1536
- idft(IDFT_1536,prachF,prach2,1);
- // here we have |empty | Prach1536 |
- memmove(prach2+(1536<<1),prach2,(1536<<2));
- // here we have |empty | Prach1536 | Prach1536| empty1536 | empty1536 | empty1536 | empty1536
- memmove(prach2+(1536<<2),prach2,(1536<<3));
- // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 | empty1536 | empty1536
- memmove(prach2+(1536<<3),prach2,(1536<<3));
- // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536
- memmove(prach,prach+(1536<<1),(Ncp<<2));
- // here we have |Prefix | Prach1536 |
- prach_len = (1536*6)+Ncp;
- } else if (prach_fmt_id == 8) { // 12x1536
- idft(IDFT_1536,prachF,prach2,1);
- // here we have |empty | Prach1536 |
- memmove(prach2+(1536<<1),prach2,(1536<<2));
- // here we have |empty | Prach1536 | Prach1536| empty1536 | empty1536 | empty1536 | empty1536
- memmove(prach2+(1536<<2),prach2,(1536<<3));
- // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 | empty1536 | empty1536
- memmove(prach2+(1536<<3),prach2,(1536<<3));
- // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536
- memmove(prach2+(1536<<1)*6,prach2,(1536<<2)*6);
- // here we have |empty | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536 | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536|
- memmove(prach,prach+(1536<<1),(Ncp<<2));
- // here we have |Prefix | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536 | Prach1536 | Prach1536| Prach1536 | Prach1536 | Prach1536 | Prach1536|
- prach_len = (1536*12)+Ncp;
- }
- }
- break;
+ case 92160:
+ // 50, 60, 70, 80, 90 MHz @ 92.16 Ms/s
+ Ncp *= 3;
+ dftlen *= 3;
+ break;
- case 122880: // full sampling @ 122.88 Ms/s
- Ncp<<=2; //to account for 122.88Mbps
- // This is after cyclic prefix
- prach2 = prach+(Ncp<<1); //times 2 for complex samples
- if (prach_sequence_length == 0){
- if (prach_fmt_id == 0) { //24576 samples @ 30.72 Ms/s, 98304 samples @ 122.88 Ms/s
- idft(IDFT_98304,prachF,prach2,1);
- // here we have |empty | Prach98304|
- memmove(prach,prach+(98304<<1),(Ncp<<2));
- // here we have |Prefix | Prach98304|
- prach_len = (98304*1)+Ncp;
- } else if (prach_fmt_id == 1) {
- idft(IDFT_98304,prachF,prach2,1);
- memmove(prach2+(98304<<1),prach2,(98304<<2));
- // here we have |empty | Prach98304 | Prach98304|
- memmove(prach,prach+(98304<<2),(Ncp<<2));
- // here we have |Prefix | Prach98304 | Prach98304|
- prach_len = (98304*2)+Ncp;
- } else if (prach_fmt_id == 2) {
- idft(IDFT_98304,prachF,prach2,1);
- memmove(prach2+(98304<<1),prach2,(98304<<2));
- // here we have |empty | Prach98304 | Prach98304| empty98304 | empty98304
- memmove(prach2+(98304<<2),prach2,(98304<<3));
- // here we have |empty | Prach98304 | Prach98304| Prach98304 | Prach98304
- memmove(prach,prach+(98304<<3),(Ncp<<2));
- // here we have |Prefix | Prach98304 | Prach98304| Prach98304 | Prach98304
- prach_len = (98304*4)+Ncp;
- } else if (prach_fmt_id == 3) { // 4x6144, Ncp 3168
- idft(IDFT_24576,prachF,prach2,1);
- memmove(prach2+(24576<<1),prach2,(24576<<2));
- // here we have |empty | Prach24576 | Prach24576| empty24576 | empty24576
- memmove(prach2+(24576<<2),prach2,(24576<<3));
- // here we have |empty | Prach24576 | Prach24576| Prach24576 | Prach24576
- memmove(prach,prach+(24576<<3),(Ncp<<2));
- // here we have |Prefix | Prach24576 | Prach24576| Prach24576 | Prach24576
- prach_len = (24576*4)+Ncp;
- }
- } else { // short sequence
- if (use_extended_prach_prefix)
- Ncp+=64; // 16*kappa, kappa=4 for 122.88Msps
- prach2 = prach+(Ncp<<1); //times 2 for complex samples
- if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
- idft(IDFT_4096,prachF,prach2,1);
- // here we have |empty | Prach4096 |
- if (prach_fmt_id != 9) {
- memmove(prach2+(4096<<1),prach2,(4096<<2));
- prach_len = (4096*2)+Ncp;
- } else prach_len = (4096*1)+Ncp;
- memmove(prach,prach+(4096<<1),(Ncp<<2));
- // here we have |Prefix | Prach4096 | Prach4096 (if ! 0xc0) |
- } else if (prach_fmt_id == 5) { // 4x4096
- idft(IDFT_4096,prachF,prach2,1);
- // here we have |empty | Prach4096 |
- memmove(prach2+(4096<<1),prach2,(4096<<2));
- // here we have |empty | Prach4096 | Prach4096| empty4096 | empty4096 |
- memmove(prach2+(4096<<2),prach2,(4096<<3));
- // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 |
- memmove(prach,prach+(4096<<1),(Ncp<<2));
- // here we have |Prefix | Prach4096 |
- prach_len = (4096*4)+Ncp;
- } else if (prach_fmt_id == 6) { // 6x4096
- idft(IDFT_4096,prachF,prach2,1);
- // here we have |empty | Prach4096 |
- memmove(prach2+(4096<<1),prach2,(4096<<2));
- // here we have |empty | Prach4096 | Prach4096| empty4096 | empty4096 | empty4096 | empty4096
- memmove(prach2+(4096<<2),prach2,(4096<<3));
- // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 | empty4096 | empty4096
- memmove(prach2+(4096<<3),prach2,(4096<<3));
- // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096
- memmove(prach,prach+(4096<<1),(Ncp<<2));
- // here we have |Prefix | Prach4096 |
- prach_len = (4096*6)+Ncp;
- } else if (prach_fmt_id == 8) { // 12x4096
- idft(IDFT_4096,prachF,prach2,1);
- // here we have |empty | Prach4096 |
- memmove(prach2+(4096<<1),prach2,(4096<<2));
- // here we have |empty | Prach4096 | Prach4096| empty4096 | empty4096 | empty4096 | empty4096
- memmove(prach2+(4096<<2),prach2,(4096<<3));
- // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 | empty4096 | empty4096
- memmove(prach2+(4096<<3),prach2,(4096<<3));
- // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096
- memmove(prach2+(4096<<1)*6,prach2,(4096<<2)*6);
- // here we have |empty | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096 | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096|
- memmove(prach,prach+(4096<<1),(Ncp<<2));
- // here we have |Prefix | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096 | Prach4096 | Prach4096| Prach4096 | Prach4096 | Prach4096 | Prach4096|
- prach_len = (4096*12)+Ncp;
- }
- }
- break;
+ case 122880:
+ // 70, 80, 90, 100 MHz @ 122.88 Ms/s
+ Ncp <<= 2;
+ dftlen <<= 2;
+ break;
- case 92160: // three quarter sampling @ 92.16 Ms/s
- Ncp = (Ncp*3); //to account for 92.16 Msps
- prach2 = prach+(Ncp<<1); //times 2 for complex samples
- if (prach_sequence_length == 0){
- if (prach_fmt_id == 0) {
- idft(IDFT_73728,prachF,prach2,1);
- // here we have |empty | Prach73728|
- memmove(prach,prach+(73728<<1),(Ncp<<2));
- // here we have |Prefix | Prach73728|
- prach_len = (73728*1)+Ncp;
- } else if (prach_fmt_id == 1) {
- idft(IDFT_73728,prachF,prach2,1);
- memmove(prach2+(73728<<1),prach2,(73728<<2));
- // here we have |empty | Prach73728 | Prach73728|
- memmove(prach,prach+(73728<<2),(Ncp<<2));
- // here we have |Prefix | Prach73728 | Prach73728|
- prach_len = (73728*2)+Ncp;
- } if (prach_fmt_id == 2) {
- idft(IDFT_73728,prachF,prach2,1);
- memmove(prach2+(73728<<1),prach2,(73728<<2));
- // here we have |empty | Prach73728 | Prach73728| empty73728 | empty73728
- memmove(prach2+(73728<<2),prach2,(73728<<3));
- // here we have |empty | Prach73728 | Prach73728| Prach73728 | Prach73728
- memmove(prach,prach+(73728<<3),(Ncp<<2));
- // here we have |Prefix | Prach73728 | Prach73728| Prach73728 | Prach73728
- prach_len = (73728*4)+Ncp;
- } else if (prach_fmt_id == 3) {
- idft(IDFT_18432,prachF,prach2,1);
- memmove(prach2+(18432<<1),prach2,(18432<<2));
- // here we have |empty | Prach18432 | Prach18432| empty18432 | empty18432
- memmove(prach2+(18432<<2),prach2,(18432<<3));
- // here we have |empty | Prach18432 | Prach18432| Prach18432 | Prach18432
- memmove(prach,prach+(18432<<3),(Ncp<<2));
- // here we have |Prefix | Prach18432 | Prach18432| Prach18432 | Prach18432
- prach_len = (18432*4)+Ncp;
- }
- } else { // short sequence
- if (use_extended_prach_prefix)
- Ncp+=48; // 16*kappa, kappa=3 for 92.16Msps
- prach2 = prach+(Ncp<<1); //times 2 for complex samples
- if (prach_fmt_id == 4 || prach_fmt_id == 7 || prach_fmt_id == 9) {
- idft(IDFT_3072,prachF,prach2,1);
- // here we have |empty | Prach3072 |
- if (prach_fmt_id != 9) {
- memmove(prach2+(3072<<1),prach2,(3072<<2));
- prach_len = (3072*2)+Ncp;
- } else prach_len = (3072*1)+Ncp;
- memmove(prach,prach+(3072<<1),(Ncp<<2));
- // here we have |Prefix | Prach3072 | Prach3072 (if ! 0xc0) |
- } else if (prach_fmt_id == 6) { // 6x3072
- idft(IDFT_3072,prachF,prach2,1);
- // here we have |empty | Prach3072 |
- memmove(prach2+(3072<<1),prach2,(3072<<2));
- // here we have |empty | Prach3072 | Prach3072| empty3072 | empty3072 | empty3072 | empty3072
- memmove(prach2+(3072<<2),prach2,(3072<<3));
- // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 | empty3072 | empty3072
- memmove(prach2+(3072<<3),prach2,(3072<<3));
- // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072
- memmove(prach,prach+(3072<<1),(Ncp<<2));
- // here we have |Prefix | Prach3072 |
- prach_len = (3072*6)+Ncp;
- } else if (prach_fmt_id == 5) { // 4x3072
- idft(IDFT_3072,prachF,prach2,1);
- // here we have |empty | Prach3072 |
- memmove(prach2+(3072<<1),prach2,(3072<<2));
- // here we have |empty | Prach3072 | Prach3072| empty3072 | empty3072 |
- memmove(prach2+(3072<<2),prach2,(3072<<3));
- // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 |
- memmove(prach,prach+(3072<<1),(Ncp<<2));
- // here we have |Prefix | Prach3072 |
- prach_len = (3072*4)+Ncp;
- } else if (prach_fmt_id == 6) { // 12x3072
- idft(IDFT_3072,prachF,prach2,1);
- // here we have |empty | Prach3072 |
- memmove(prach2+(3072<<1),prach2,(3072<<2));
- // here we have |empty | Prach3072 | Prach3072| empty3072 | empty3072 | empty3072 | empty3072
- memmove(prach2+(3072<<2),prach2,(3072<<3));
- // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 | empty3072 | empty3072
- memmove(prach2+(3072<<3),prach2,(3072<<3));
- // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072
- memmove(prach2+(3072<<1)*6,prach2,(3072<<2)*6);
- // here we have |empty | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072 | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072|
- memmove(prach,prach+(3072<<1),(Ncp<<2));
- // here we have |Prefix | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072 | Prach3072 | Prach3072| Prach3072 | Prach3072 | Prach3072 | Prach3072|
- prach_len = (3072*12)+Ncp;
- }
- }
- break;
+ default:
+ AssertFatal(1==0,"sample rate %f MHz not supported for numerology %d\n", fp->samples_per_subframe / 1000.0, mu);
+ }
- default:
- AssertFatal(1==0,"sample rate %f MHz not supported for numerology %d\n", fp->samples_per_subframe / 1000.0, mu);
+ #ifdef NR_PRACH_DEBUG
+ LOG_I(PHY, "PRACH [UE %d] Ncp %d, dftlen %d \n", Mod_id, Ncp, dftlen);
+ #endif
+
+ /********************************************************
+ *
+ * In function init_prach_tables:
+ * to compute quantized roots of unity ru(n) = 32767 * exp j*[ (2 * PI * n) / N_ZC ]
+ *
+ * In compute_prach_seq:
+ * to calculate Xu = DFT xu = xu (inv_u*k) * Xu[0] (This is a Zadoff-Chou sequence property: DFT ZC sequence is another ZC sequence)
+ *
+ * In generate_prach:
+ * to do the cyclic-shifted DFT by multiplying Xu[k] * ru[k*preamble_shift] as:
+ * If X[k] = DFT x(n) -> X_shifted[k] = DFT x(n+preamble_shift) = X[k] * exp -j*[ (2*PI*k*preamble_shift) / N_ZC ]
+ *
+ *********************************************************/
+
+ Xu = (int16_t*)ue->X_u[preamble_offset-first_nonzero_root_idx];
+
+ #if defined (PRACH_WRITE_OUTPUT_DEBUG)
+ LOG_M("X_u.m", "X_u", (int16_t*)ue->X_u[preamble_offset-first_nonzero_root_idx], N_ZC, 1, 1);
+ #endif
+
+ for (offset=0,offset2=0; offset<N_ZC; offset++,offset2+=preamble_shift) {
+
+ if (offset2 >= N_ZC)
+ offset2 -= N_ZC;
+
+ Xu_re = (((int32_t)Xu[offset<<1]*amp)>>15);
+ Xu_im = (((int32_t)Xu[1+(offset<<1)]*amp)>>15);
+ prachF[k++]= ((Xu_re*nr_ru[offset2<<1]) - (Xu_im*nr_ru[1+(offset2<<1)]))>>15;
+ prachF[k++]= ((Xu_im*nr_ru[offset2<<1]) + (Xu_re*nr_ru[1+(offset2<<1)]))>>15;
+
+ if (k==dftlen) k=0;
+ }
+
+ #if defined (PRACH_WRITE_OUTPUT_DEBUG)
+ LOG_M("prachF.m", "prachF", &prachF[1804], 1024, 1, 1);
+ LOG_M("Xu.m", "Xu", Xu, N_ZC, 1, 1);
+ #endif
+
+ // This is after cyclic prefix
+ prach2 = prach+(2*Ncp); // times 2 for complex samples
+ const idft_size_idx_t idft_size = get_idft(dftlen);
+ idft(idft_size, prachF, prach, 1);
+ memmove(prach2, prach, (dftlen<<2));
+
+ if (prach_sequence_length == 0) {
+ if (prach_fmt_id == 0) {
+ // here we have | empty | Prach |
+ memcpy(prach, prach+(dftlen<<1), (Ncp<<2));
+ // here we have | Prefix | Prach |
+ prach_len = dftlen+Ncp;
+ } else if (prach_fmt_id == 1) {
+ // here we have | empty | Prach | empty |
+ memcpy(prach2+(dftlen<<1), prach2, (dftlen<<2));
+ // here we have | empty | Prach | Prach |
+ memcpy(prach, prach+(dftlen<<2), (Ncp<<2));
+ // here we have | Prefix | Prach | Prach |
+ prach_len = (dftlen*2)+Ncp;
+ } else if (prach_fmt_id == 2 || prach_fmt_id == 3) {
+ // here we have | empty | Prach | empty | empty | empty |
+ memcpy(prach2+(dftlen<<1), prach2, (dftlen<<2));
+ // here we have | empty | Prach | Prach | empty | empty |
+ memcpy(prach2+(dftlen<<2), prach2, (dftlen<<3));
+ // here we have | empty | Prach | Prach | Prach | Prach |
+ memcpy(prach, prach+(dftlen<<3), (Ncp<<2));
+ // here we have | Prefix | Prach | Prach | Prach | Prach |
+ prach_len = (dftlen*4)+Ncp;
+ }
+ } else { // short PRACH sequence
+ if (prach_fmt_id == 9) {
+ // here we have | empty | Prach |
+ memcpy(prach, prach+(dftlen<<1), (Ncp<<2));
+ // here we have | Prefix | Prach |
+ prach_len = (dftlen*1)+Ncp;
+ } else if (prach_fmt_id == 4 || prach_fmt_id == 7) {
+ // here we have | empty | Prach | empty |
+ memcpy(prach2+(dftlen<<1), prach2, (dftlen<<2));
+ // here we have | empty | Prach | Prach |
+ memcpy(prach, prach+(dftlen<<1), (Ncp<<2));
+ // here we have | Prefix | Prach | Prach |
+ prach_len = (dftlen*2)+Ncp;
+ } else if (prach_fmt_id == 5) { // 4xdftlen
+ // here we have | empty | Prach | empty | empty | empty |
+ memcpy(prach2+(dftlen<<1), prach2, (dftlen<<2));
+ // here we have | empty | Prach | Prach | empty | empty |
+ memcpy(prach2+(dftlen<<2), prach2, (dftlen<<3));
+ // here we have | empty | Prach | Prach | Prach | Prach |
+ memcpy(prach, prach+(dftlen<<1), (Ncp<<2));
+ // here we have | Prefix | Prach | Prach | Prach | Prach |
+ prach_len = (dftlen*4)+Ncp;
+ } else if (prach_fmt_id == 6) { // 6xdftlen
+ // here we have | empty | Prach | empty | empty | empty | empty | empty |
+ memcpy(prach2+(dftlen<<1), prach2, (dftlen<<2));
+ // here we have | empty | Prach | Prach | empty | empty | empty | empty |
+ memcpy(prach2+(dftlen<<2), prach2, (dftlen<<3));
+ // here we have | empty | Prach | Prach | Prach | Prach | empty | empty |
+ memcpy(prach2+(dftlen<<3), prach2, (dftlen<<3));
+ // here we have | empty | Prach | Prach | Prach | Prach | Prach | Prach |
+ memcpy(prach, prach+(dftlen<<1), (Ncp<<2));
+ // here we have | Prefix | Prach | Prach | Prach | Prach | Prach | Prach |
+ prach_len = (dftlen*6)+Ncp;
+ } else if (prach_fmt_id == 8) { // 12xdftlen
+ // here we have | empty | Prach | empty | empty | empty | empty | empty | empty | empty | empty | empty | empty | empty |
+ memcpy(prach2+(dftlen<<1), prach2, (dftlen<<2));
+ // here we have | empty | Prach | Prach | empty | empty | empty | empty | empty | empty | empty | empty | empty | empty |
+ memcpy(prach2+(dftlen<<2), prach2, (dftlen<<3));
+ // here we have | empty | Prach | Prach | Prach | Prach | empty | empty | empty | empty | empty | empty | empty | empty |
+ memcpy(prach2+(dftlen<<3), prach2, (dftlen<<3));
+ // here we have | empty | Prach | Prach | Prach | Prach | Prach | Prach | empty | empty | empty | empty | empty | empty |
+ memcpy(prach2+(dftlen<<1)*6, prach2, (dftlen<<2)*6);
+ // here we have | empty | Prach | Prach | Prach | Prach | Prach | Prach | Prach | Prach | Prach | Prach | Prach | Prach |
+ memcpy(prach, prach+(dftlen<<1), (Ncp<<2));
+ // here we have | Prefix | Prach | Prach | Prach | Prach | Prach | Prach | Prach | Prach | Prach | Prach | Prach | Prach |
+ prach_len = (dftlen*12)+Ncp;
}
}