maki/openairinterface5g
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

code cleanup for abstraction

Browse Source 
git-svn-id: http://svn.eurecom.fr/openair4G/trunk@5112 818b1a75-f10b-46b9-bf7c-635c3b92a50f
This commit is contained in:
Florian Kaltenberger
2014-03-03 17:53:36 +00:00
parent 6969100a68
commit d78d39e6a0
1 changed files with 101 additions and 117 deletions
Download Patch File Download Diff File Expand all files Collapse all files
openair1/SIMULATION/LTE_PHY/dlsim.c
+92 -108
View File
@@ -185,7 +185,7 @@ int main(int argc, char **argv) {
int s,Kr,Kr_bytes; int s,Kr,Kr_bytes;
double sigma2, sigma2_dB=10,SNR,snr0=-2.0,snr1,rate,saving_bler=1; double sigma2, sigma2_dB=10,SNR,snr0=-2.0,snr1,rate;
double snr_step=1,input_snr_step=1, snr_int=30; double snr_step=1,input_snr_step=1, snr_int=30;
LTE_DL_FRAME_PARMS *frame_parms; LTE_DL_FRAME_PARMS *frame_parms;
@@ -1213,7 +1213,6 @@ int main(int argc, char **argv) {
for (ch_realization=0;ch_realization<n_ch_rlz;ch_realization++){ for (ch_realization=0;ch_realization<n_ch_rlz;ch_realization++){
if(abstx){ if(abstx){
printf("**********************Channel Realization Index = %d **************************\n", ch_realization); printf("**********************Channel Realization Index = %d **************************\n", ch_realization);
saving_bler=1;
} }
for (SNR=snr0;SNR<snr1;SNR+=snr_step) { for (SNR=snr0;SNR<snr1;SNR+=snr_step) {
@@ -1664,86 +1663,83 @@ int main(int argc, char **argv) {
if(abstx){ if(abstx){
if(saving_bler==0) if (trials==0 && round==0) {
if (trials==0 && round==0) { // calculate freq domain representation to compute SINR
// calculate freq domain representation to compute SINR freq_channel(eNB2UE[0], NB_RB,2*NB_RB + 1);
freq_channel(eNB2UE[0], NB_RB,2*NB_RB + 1); // snr=pow(10.0,.1*SNR);
// snr=pow(10.0,.1*SNR); fprintf(csv_fd,"%f,",SNR);
fprintf(csv_fd,"%f,",SNR);
for (u=0;u<2*NB_RB;u++){ for (u=0;u<2*NB_RB;u++){
for (aarx=0;aarx<eNB2UE[0]->nb_rx;aarx++) { for (aarx=0;aarx<eNB2UE[0]->nb_rx;aarx++) {
for (aatx=0;aatx<eNB2UE[0]->nb_tx;aatx++) { for (aatx=0;aatx<eNB2UE[0]->nb_tx;aatx++) {
channelx = eNB2UE[0]->chF[aarx+(aatx*eNB2UE[0]->nb_rx)][u].x; channelx = eNB2UE[0]->chF[aarx+(aatx*eNB2UE[0]->nb_rx)][u].x;
channely = eNB2UE[0]->chF[aarx+(aatx*eNB2UE[0]->nb_rx)][u].y; channely = eNB2UE[0]->chF[aarx+(aatx*eNB2UE[0]->nb_rx)][u].y;
fprintf(csv_fd,"%e+i*(%e),",channelx,channely); fprintf(csv_fd,"%e+i*(%e),",channelx,channely);
}
}
} }
}
}
if(num_rounds>1){ if(num_rounds>1){
freq_channel(eNB2UE[1], NB_RB,2*NB_RB + 1); freq_channel(eNB2UE[1], NB_RB,2*NB_RB + 1);
for (u=0;u<2*NB_RB;u++){ for (u=0;u<2*NB_RB;u++){
for (aarx=0;aarx<eNB2UE[1]->nb_rx;aarx++) { for (aarx=0;aarx<eNB2UE[1]->nb_rx;aarx++) {
for (aatx=0;aatx<eNB2UE[1]->nb_tx;aatx++) { for (aatx=0;aatx<eNB2UE[1]->nb_tx;aatx++) {
channelx = eNB2UE[1]->chF[aarx+(aatx*eNB2UE[1]->nb_rx)][u].x; channelx = eNB2UE[1]->chF[aarx+(aatx*eNB2UE[1]->nb_rx)][u].x;
channely = eNB2UE[1]->chF[aarx+(aatx*eNB2UE[1]->nb_rx)][u].y; channely = eNB2UE[1]->chF[aarx+(aatx*eNB2UE[1]->nb_rx)][u].y;
fprintf(csv_fd,"%e+i*(%e),",channelx,channely); fprintf(csv_fd,"%e+i*(%e),",channelx,channely);
}
}
}
freq_channel(eNB2UE[2], NB_RB,2*NB_RB + 1);
for (u=0;u<2*NB_RB;u++){
for (aarx=0;aarx<eNB2UE[2]->nb_rx;aarx++) {
for (aatx=0;aatx<eNB2UE[2]->nb_tx;aatx++) {
channelx = eNB2UE[2]->chF[aarx+(aatx*eNB2UE[2]->nb_rx)][u].x;
channely = eNB2UE[2]->chF[aarx+(aatx*eNB2UE[2]->nb_rx)][u].y;
fprintf(csv_fd,"%e+i*(%e),",channelx,channely);
}
}
}
freq_channel(eNB2UE[3], NB_RB,2*NB_RB + 1);
for (u=0;u<2*NB_RB;u++){
for (aarx=0;aarx<eNB2UE[3]->nb_rx;aarx++) {
for (aatx=0;aatx<eNB2UE[3]->nb_tx;aatx++) {
channelx = eNB2UE[3]->chF[aarx+(aatx*eNB2UE[3]->nb_rx)][u].x;
channely = eNB2UE[3]->chF[aarx+(aatx*eNB2UE[3]->nb_rx)][u].y;
fprintf(csv_fd,"%e+i*(%e),",channelx,channely);
}
}
} }
} }
} }
} freq_channel(eNB2UE[2], NB_RB,2*NB_RB + 1);
for (u=0;u<2*NB_RB;u++){
for (aarx=0;aarx<eNB2UE[2]->nb_rx;aarx++) {
for (aatx=0;aatx<eNB2UE[2]->nb_tx;aatx++) {
channelx = eNB2UE[2]->chF[aarx+(aatx*eNB2UE[2]->nb_rx)][u].x;
channely = eNB2UE[2]->chF[aarx+(aatx*eNB2UE[2]->nb_rx)][u].y;
fprintf(csv_fd,"%e+i*(%e),",channelx,channely);
}
}
}
freq_channel(eNB2UE[3], NB_RB,2*NB_RB + 1);
//AWGN for (u=0;u<2*NB_RB;u++){
// This is the SNR on the PDSCH for OFDM symbols without pilots -> rho_A for (aarx=0;aarx<eNB2UE[3]->nb_rx;aarx++) {
sigma2_dB = 10*log10((double)tx_lev) +10*log10((double)PHY_vars_eNB->lte_frame_parms.ofdm_symbol_size/(double)(NB_RB*12)) - SNR - get_pa_dB(PHY_vars_eNB->pdsch_config_dedicated); for (aatx=0;aatx<eNB2UE[3]->nb_tx;aatx++) {
sigma2 = pow(10,sigma2_dB/10); channelx = eNB2UE[3]->chF[aarx+(aatx*eNB2UE[3]->nb_rx)][u].x;
if (n_frames==1) channely = eNB2UE[3]->chF[aarx+(aatx*eNB2UE[3]->nb_rx)][u].y;
printf("Sigma2 %f (sigma2_dB %f,%f,%f )\n",sigma2,sigma2_dB,10*log10((double)PHY_vars_eNB->lte_frame_parms.ofdm_symbol_size/(double)(NB_RB*12)),get_pa_dB(PHY_vars_eNB->pdsch_config_dedicated)); fprintf(csv_fd,"%e+i*(%e),",channelx,channely);
}
for (i=0; i<2*frame_parms->samples_per_tti; i++) { }
for (aa=0;aa<PHY_vars_eNB->lte_frame_parms.nb_antennas_rx;aa++) { }
//printf("s_re[0][%d]=> %f , r_re[0][%d]=> %f\n",i,s_re[aa][i],i,r_re[aa][i]);
((short*) PHY_vars_UE->lte_ue_common_vars.rxdata[aa])[(2*subframe*PHY_vars_UE->lte_frame_parms.samples_per_tti)+2*i] =
(short) (r_re[aa][i] + sqrt(sigma2/2)*gaussdouble(0.0,1.0));
((short*) PHY_vars_UE->lte_ue_common_vars.rxdata[aa])[(2*subframe*PHY_vars_UE->lte_frame_parms.samples_per_tti)+2*i+1] =
(short) (r_im[aa][i] + (iqim*r_re[aa][i]) + sqrt(sigma2/2)*gaussdouble(0.0,1.0));
} }
} }
}
// lte_sync_time_init(PHY_vars_eNB->lte_frame_parms,lte_ue_common_vars); //AWGN
// lte_sync_time(lte_ue_common_vars->rxdata, PHY_vars_eNB->lte_frame_parms); // This is the SNR on the PDSCH for OFDM symbols without pilots -> rho_A
// lte_sync_time_free(); sigma2_dB = 10*log10((double)tx_lev) +10*log10((double)PHY_vars_eNB->lte_frame_parms.ofdm_symbol_size/(double)(NB_RB*12)) - SNR - get_pa_dB(PHY_vars_eNB->pdsch_config_dedicated);
sigma2 = pow(10,sigma2_dB/10);
if (n_frames==1)
printf("Sigma2 %f (sigma2_dB %f,%f,%f )\n",sigma2,sigma2_dB,10*log10((double)PHY_vars_eNB->lte_frame_parms.ofdm_symbol_size/(double)(NB_RB*12)),get_pa_dB(PHY_vars_eNB->pdsch_config_dedicated));
/* for (i=0; i<2*frame_parms->samples_per_tti; i++) {
for (aa=0;aa<PHY_vars_eNB->lte_frame_parms.nb_antennas_rx;aa++) {
//printf("s_re[0][%d]=> %f , r_re[0][%d]=> %f\n",i,s_re[aa][i],i,r_re[aa][i]);
((short*) PHY_vars_UE->lte_ue_common_vars.rxdata[aa])[(2*subframe*PHY_vars_UE->lte_frame_parms.samples_per_tti)+2*i] =
(short) (r_re[aa][i] + sqrt(sigma2/2)*gaussdouble(0.0,1.0));
((short*) PHY_vars_UE->lte_ue_common_vars.rxdata[aa])[(2*subframe*PHY_vars_UE->lte_frame_parms.samples_per_tti)+2*i+1] =
(short) (r_im[aa][i] + (iqim*r_re[aa][i]) + sqrt(sigma2/2)*gaussdouble(0.0,1.0));
}
}
// lte_sync_time_init(PHY_vars_eNB->lte_frame_parms,lte_ue_common_vars);
// lte_sync_time(lte_ue_common_vars->rxdata, PHY_vars_eNB->lte_frame_parms);
// lte_sync_time_free();
/*
// optional: read rx_frame from file // optional: read rx_frame from file
if ((rx_frame_file = fopen("rx_frame.dat","r")) == NULL) if ((rx_frame_file = fopen("rx_frame.dat","r")) == NULL)
{ {
@@ -1757,27 +1753,27 @@ int main(int argc, char **argv) {
printf("Read %d bytes\n",result); printf("Read %d bytes\n",result);
fclose(rx_frame_file); fclose(rx_frame_file);
*/ */
if (n_frames==1) { if (n_frames==1) {
printf("RX level in null symbol %d\n",dB_fixed(signal_energy(&PHY_vars_UE->lte_ue_common_vars.rxdata[0][160+OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES],OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2))); printf("RX level in null symbol %d\n",dB_fixed(signal_energy(&PHY_vars_UE->lte_ue_common_vars.rxdata[0][160+OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES],OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2)));
printf("RX level in data symbol %d\n",dB_fixed(signal_energy(&PHY_vars_UE->lte_ue_common_vars.rxdata[0][160+(2*OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES)],OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2))); printf("RX level in data symbol %d\n",dB_fixed(signal_energy(&PHY_vars_UE->lte_ue_common_vars.rxdata[0][160+(2*OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES)],OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2)));
printf("rx_level Null symbol %f\n",10*log10(signal_energy_fp(r_re,r_im,1,OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2,256+(OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES)))); printf("rx_level Null symbol %f\n",10*log10(signal_energy_fp(r_re,r_im,1,OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2,256+(OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES))));
printf("rx_level data symbol %f\n",10*log10(signal_energy_fp(r_re,r_im,1,OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2,256+(2*OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES)))); printf("rx_level data symbol %f\n",10*log10(signal_energy_fp(r_re,r_im,1,OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES/2,256+(2*OFDM_SYMBOL_SIZE_COMPLEX_SAMPLES))));
} }
if (PHY_vars_eNB->lte_frame_parms.Ncp == 0) { // normal prefix if (PHY_vars_eNB->lte_frame_parms.Ncp == 0) { // normal prefix
pilot1 = 4; pilot1 = 4;
pilot2 = 7; pilot2 = 7;
pilot3 = 11; pilot3 = 11;
} }
else { // extended prefix else { // extended prefix
pilot1 = 3; pilot1 = 3;
pilot2 = 6; pilot2 = 6;
pilot3 = 9; pilot3 = 9;
} }
// Inner receiver scheduling for 3 slots // Inner receiver scheduling for 3 slots
for (Ns=(2*subframe);Ns<((2*subframe)+3);Ns++) { for (Ns=(2*subframe);Ns<((2*subframe)+3);Ns++) {
for (l=0;l<pilot2;l++) { for (l=0;l<pilot2;l++) {
if (n_frames==1) if (n_frames==1)
@@ -2131,14 +2127,13 @@ int main(int argc, char **argv) {
//saving PMI in case of Transmission Mode > 5 //saving PMI in case of Transmission Mode > 5
if(abstx){ if(abstx){
if(saving_bler==0) if (trials==0 && round==0 && transmission_mode>=5){
if (trials==0 && round==0 && transmission_mode>=5){ for (iii=0; iii<NB_RB; iii++){
for (iii=0; iii<NB_RB; iii++){ //fprintf(csv_fd, "%d, %d", (PHY_vars_UE->lte_ue_pdsch_vars[eNB_id]->pmi_ext[iii]),(PHY_vars_UE->lte_ue_pdsch_vars[eNB_id_i]->pmi_ext[iii]));
//fprintf(csv_fd, "%d, %d", (PHY_vars_UE->lte_ue_pdsch_vars[eNB_id]->pmi_ext[iii]),(PHY_vars_UE->lte_ue_pdsch_vars[eNB_id_i]->pmi_ext[iii])); fprintf(csv_fd,"%x,%x,",(PHY_vars_UE->lte_ue_pdsch_vars[eNB_id]->pmi_ext[iii]),(PHY_vars_UE->lte_ue_pdsch_vars[eNB_id]->pmi_ext[iii]));
fprintf(csv_fd,"%x,%x,",(PHY_vars_UE->lte_ue_pdsch_vars[eNB_id]->pmi_ext[iii]),(PHY_vars_UE->lte_ue_pdsch_vars[eNB_id]->pmi_ext[iii])); msg(" %x",(PHY_vars_UE->lte_ue_pdsch_vars[eNB_id]->pmi_ext[iii]));
msg(" %x",(PHY_vars_UE->lte_ue_pdsch_vars[eNB_id]->pmi_ext[iii]));
}
} }
}
} }
PHY_vars_UE->dlsch_ue[0][0]->rnti = (common_flag==0) ? n_rnti: SI_RNTI; PHY_vars_UE->dlsch_ue[0][0]->rnti = (common_flag==0) ? n_rnti: SI_RNTI;
@@ -2395,22 +2390,11 @@ int main(int argc, char **argv) {
blerr[3] = (double)errs[3]/(round_trials[3]); blerr[3] = (double)errs[3]/(round_trials[3]);
fprintf(csv_fd,"%e,%e,%e,%e;\n",blerr[0],blerr[1],blerr[2],blerr[3]); fprintf(csv_fd,"%e,%e,%e,%e;\n",blerr[0],blerr[1],blerr[2],blerr[3]);
} }
else else {
{ fprintf(csv_fd,"%e;\n",blerr[0]);
}
if(saving_bler==0)
fprintf(csv_fd,"%e;\n",blerr[0]);
}
} //ABStraction } //ABStraction
if(num_rounds==1){
bler= (double)errs[0]/(round_trials[0]);
if (bler<1)
{snr_step = input_snr_step; saving_bler = 0;}
else
{snr_step = 1; saving_bler = 1;}
}
if (((double)errs[0]/(round_trials[0]))<1e-2) if (((double)errs[0]/(round_trials[0]))<1e-2)
break; break;
}// SNR }// SNR