163 lines
6.7 KiB
C
163 lines
6.7 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/MODULATION/beamforming.c
|
|
* \brief
|
|
* \author X. JIANG, F. Kaltenberger, R. KNOPP
|
|
* \date 2016
|
|
* \version 0.1
|
|
* \company Eurecom
|
|
* \email: xiwen.jiang@eurecom.fr,florian.kaltenberger@eurecom.fr,raymond.knopp@eurecom.fr
|
|
* \note
|
|
* \warning
|
|
*/
|
|
|
|
/*! \function beam_precoding_one_eNB
|
|
* \brief Create and Implementation of beamforming in one eNB
|
|
* \author TY Hsu, SY Yeh(fdragon), TH Wang(Judy)
|
|
* \date 2018
|
|
* \version 0.1
|
|
* \company ISIP@NCTU and Eurecom
|
|
* \email: tyhsu@cs.nctu.edu.tw,fdragon.cs96g@g2.nctu.edu.tw,Tsu-Han.Wang@eurecom.fr
|
|
* \note
|
|
* \warning
|
|
*/
|
|
|
|
|
|
#include "PHY/defs_common.h"
|
|
#include "PHY/defs_eNB.h"
|
|
#include "PHY/phy_extern.h"
|
|
#include "PHY/CODING/coding_defs.h"
|
|
#include "PHY/CODING/coding_extern.h"
|
|
#include "PHY/CODING/lte_interleaver_inline.h"
|
|
#include "PHY/LTE_TRANSPORT/transport_eNB.h"
|
|
#include "modulation_eNB.h"
|
|
#include "common/utils/LOG/vcd_signal_dumper.h"
|
|
|
|
int beam_precoding(int32_t **txdataF,
|
|
int32_t **txdataF_BF,
|
|
LTE_DL_FRAME_PARMS *frame_parms,
|
|
int32_t ***beam_weights,
|
|
int symbol,
|
|
int aa,
|
|
int p)
|
|
{
|
|
//LOG_I(PHY,"Starting precoding for symbol %d, physical antenna %d, logical port %d\n",symbol,aa,p);
|
|
int rb_offset_neg = frame_parms->ofdm_symbol_size - (6*frame_parms->N_RB_DL);
|
|
/* if(rb<(frame_parms->N_RB_DL/2)) {
|
|
rb_offset = rb*12;
|
|
}else if(rb>=(frame_parms->N_RB_DL/2)) {
|
|
rb_offset = 12*(rb-(frame_parms->N_RB_DL/2));
|
|
//rb_offset = 1+(12*(rb-(frame_parms->N_RB_DL>>1)));
|
|
}*/
|
|
//printf("aa = %d, p = %d, symbol = %d, rb = %d, rb_offset = %d, orisma = %d\n",aa,p,symbol,rb,rb_offset,rb_offset+symbol*12);
|
|
//printf("aa = %d, p = %d, symbol = %d\n",aa,p,symbol);
|
|
|
|
|
|
//VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_MULTADD,1);
|
|
|
|
multadd_cpx_vector((int16_t*)&txdataF[p][rb_offset_neg+(symbol*frame_parms->ofdm_symbol_size)],
|
|
(int16_t*)&beam_weights[p][aa][rb_offset_neg],
|
|
(int16_t*)&txdataF_BF[aa][rb_offset_neg+(symbol*frame_parms->ofdm_symbol_size)],
|
|
0,
|
|
6*frame_parms->N_RB_DL,
|
|
15);
|
|
multadd_cpx_vector((int16_t*)&txdataF[p][(symbol*frame_parms->ofdm_symbol_size)],
|
|
(int16_t*)&beam_weights[p][aa][0],
|
|
(int16_t*)&txdataF_BF[aa][(symbol*frame_parms->ofdm_symbol_size)],
|
|
0,
|
|
7*frame_parms->N_RB_DL, // to allow for extra RE at the end, 12 useless multipy-adds (first one at DC and 11 at end)
|
|
15);
|
|
/*
|
|
if (last_prb == 1 && rb>=(frame_parms->N_RB_DL/2)) {
|
|
// do extra complex multiplication for remaining RE
|
|
uint16_t re = frame_parms->ofdm_symbol_size-1;
|
|
((int16_t*)&txdataF_BF[aa][(rb_offset+1)+re])[0] += (int16_t)((((int16_t*)&txdataF[p][(rb_offset+1)+(symbol*frame_parms->ofdm_symbol_size)+re])[0]*((int16_t*)&beam_weights[p][aa][(rb_offset+1)+re])[0])>>15);
|
|
((int16_t*)&txdataF_BF[aa][(rb_offset+1)+re])[0] -= (int16_t)((((int16_t*)&txdataF[p][(rb_offset+1)+(symbol*frame_parms->ofdm_symbol_size)+re])[1]*((int16_t*)&beam_weights[p][aa][(rb_offset+1)+re])[1])>>15);
|
|
((int16_t*)&txdataF_BF[aa][(rb_offset+1)+re])[1] += (int16_t)((((int16_t*)&txdataF[p][(rb_offset+1)+(symbol*frame_parms->ofdm_symbol_size)+re])[0]*((int16_t*)&beam_weights[p][aa][(rb_offset+1)+re])[1])>>15);
|
|
((int16_t*)&txdataF_BF[aa][(rb_offset+1)+re])[1] += (int16_t)((((int16_t*)&txdataF[p][(rb_offset+1)+(symbol*frame_parms->ofdm_symbol_size)+re])[1]*((int16_t*)&beam_weights[p][aa][(rb_offset+1)+re])[0])>>15);
|
|
}
|
|
*/
|
|
//VCD_SIGNAL_DUMPER_DUMP_FUNCTION_BY_NAME(VCD_SIGNAL_DUMPER_FUNCTIONS_MULTADD,0);
|
|
|
|
return 0;
|
|
}
|
|
|
|
|
|
int beam_precoding_one_eNB(int32_t **txdataF,
|
|
int32_t **txdataF_BF,
|
|
int32_t ***beam_weights,
|
|
int subframe,
|
|
int nb_antenna_ports,
|
|
int nb_tx, // total physical antenna
|
|
LTE_DL_FRAME_PARMS *frame_parms
|
|
)
|
|
{
|
|
int p, re, symbol, aa; // loop index
|
|
int re_offset;
|
|
int ofdm_symbol_size = frame_parms->ofdm_symbol_size;
|
|
int symbols_per_tti = frame_parms->symbols_per_tti;
|
|
int nb_antenna_ports_eNB = frame_parms->nb_antenna_ports_eNB; // logic antenna ports
|
|
|
|
|
|
re_offset = ofdm_symbol_size*symbols_per_tti*subframe;
|
|
|
|
|
|
// txdataF_BF[aa][re] = sum(beam_weghts[p][aa][re]*txdataF[p][re]), p=0~nb_antenna_ports-1
|
|
// real part and image part need to compute separately
|
|
|
|
for (aa=0; aa<nb_tx; aa++) {
|
|
memset(txdataF_BF[aa],0,sizeof(int32_t)*(ofdm_symbol_size*symbols_per_tti));
|
|
for(p=0;p<nb_antenna_ports;p++){
|
|
if (p<nb_antenna_ports_eNB || p==5){
|
|
for (symbol=0; symbol<symbols_per_tti; symbol++){
|
|
|
|
multadd_cpx_vector((int16_t*)&txdataF[p][symbol*ofdm_symbol_size+re_offset],
|
|
(int16_t*)beam_weights[p][aa],
|
|
(int16_t*)&txdataF_BF[aa][symbol*ofdm_symbol_size],
|
|
0,
|
|
ofdm_symbol_size,
|
|
15);
|
|
|
|
|
|
/*
|
|
for (re=0; re<ofdm_symbol_size; re++){
|
|
// direct
|
|
((int16_t*)&txdataF_BF[aa][re])[0] += (int16_t)((((int16_t*)&txdataF[p][re+symbol*ofdm_symbol_size+re_offset])[0]*((int16_t*)&beam_weights[p][aa][re])[0])>>15);
|
|
((int16_t*)&txdataF_BF[aa][re])[0] -= (int16_t)((((int16_t*)&txdataF[p][re+symbol*ofdm_symbol_size+re_offset])[1]*((int16_t*)&beam_weights[p][aa][re])[1])>>15);
|
|
((int16_t*)&txdataF_BF[aa][re])[1] += (int16_t)((((int16_t*)&txdataF[p][re+symbol*ofdm_symbol_size+re_offset])[0]*((int16_t*)&beam_weights[p][aa][re])[1])>>15);
|
|
((int16_t*)&txdataF_BF[aa][re])[1] += (int16_t)((((int16_t*)&txdataF[p][re+symbol*ofdm_symbol_size+re_offset])[1]*((int16_t*)&beam_weights[p][aa][re])[0])>>15);
|
|
}
|
|
*/
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
return 0;
|
|
}
|