mirror of
https://gitlab.eurecom.fr/oai/openairinterface5g.git
synced 2026-07-13 04:30:28 +00:00
8107939f08
- all RAN code, CI code, configuration files, dockerfiles, in CSSL v1.0
- all deployment code (openshift, charts, ancillary files like shell
scripts), in MIT
- documentation in CC-BY-4.0
- exceptions might apply and are listed in NOTICE
- there is a new LICENSES folder with all licenses
- CONTRIBUTIONS.md has been updated accordingly
For automated changes based on OAI PL v1.1:
perl -i~ -0pe 's/\/\*.*Licensed to the OpenAirInterface.*openairinterface.org\n#?/\/*\n * SPDX-License-Identifier: LicenseRef-CSSL-1.0\n/s' **/*.{c,h,cpp}
perl -i~ -0pe 's/\/\*.*Licensed to the OpenAirInterface.*openairinterface.org\n#?/\/*\n * SPDX-License-Identifier: LicenseRef-CSSL-1.0\n/s' **/*.ts
perl -i~ -0pe 's/<!--.*Licensed to the OpenAirInterface.*openairinterface.org\n.*-->/<!-- SPDX-License-Identifier: LicenseRef-CSSL-1.0 -->/s' **/*.xml
The rest (cmake, files with missing license, cmake) manually.
120 lines
3.2 KiB
C
120 lines
3.2 KiB
C
/*
|
|
* SPDX-License-Identifier: LicenseRef-CSSL-1.0
|
|
*/
|
|
|
|
//#define DEBUG_DAC 1
|
|
#include <math.h>
|
|
#include <stdio.h>
|
|
#include "PHY/TOOLS/tools_defs.h"
|
|
#include "rf.h"
|
|
|
|
|
|
void dac(double *s_re[2],
|
|
double *s_im[2],
|
|
int32_t **input,
|
|
uint32_t input_offset,
|
|
uint32_t nb_tx_antennas,
|
|
uint32_t length,
|
|
double amp_dBm,
|
|
uint8_t B,
|
|
uint32_t meas_length,
|
|
uint32_t meas_offset)
|
|
{
|
|
|
|
int i;
|
|
int aa;
|
|
double V=0.0,amp;
|
|
|
|
for (i=0; i<length; i++) {
|
|
for (aa=0; aa<nb_tx_antennas; aa++) {
|
|
s_re[aa][i] = ((double)(((short *)input[aa]))[((i+input_offset)<<1)])/(1<<(B-1));
|
|
s_im[aa][i] = ((double)(((short *)input[aa]))[((i+input_offset)<<1)+1])/(1<<(B-1));
|
|
|
|
}
|
|
}
|
|
|
|
for (i=meas_offset; i<meas_offset+meas_length; i++) {
|
|
for (aa=0; aa<nb_tx_antennas; aa++) {
|
|
V= V + (s_re[aa][i]*s_re[aa][i]) + (s_im[aa][i]*s_im[aa][i]);
|
|
}
|
|
}
|
|
|
|
V /= (meas_length);
|
|
#ifdef DEBUG_DAC
|
|
LOG_I(OCM,"DAC: 10*log10(V)=%f (%f)\n",10*log10(V),V);
|
|
#endif
|
|
|
|
if (V) {
|
|
amp = pow(10.0,.1*amp_dBm)/V;
|
|
amp = sqrt(amp);
|
|
} else {
|
|
amp = 1;
|
|
}
|
|
|
|
for (i=0; i<length; i++) {
|
|
for (aa=0; aa<nb_tx_antennas; aa++) {
|
|
s_re[aa][i] *= amp;
|
|
s_im[aa][i] *= amp;
|
|
}
|
|
}
|
|
}
|
|
|
|
double dac_fixed_gain(double *s_re[2],
|
|
double *s_im[2],
|
|
int32_t **input,
|
|
uint32_t input_offset,
|
|
uint32_t nb_tx_antennas,
|
|
uint32_t length,
|
|
uint32_t input_offset_meas,
|
|
uint32_t length_meas,
|
|
uint8_t B,
|
|
double txpwr_dBm,
|
|
uint8_t do_amp_compute,
|
|
double *amp1,
|
|
int NB_RE)
|
|
{
|
|
|
|
int i;
|
|
int aa;
|
|
double amp1_local = 0.0;
|
|
double *amp1p = (amp1 == NULL) ? &1_local : amp1;
|
|
double amp = pow(10.0,.05*txpwr_dBm)/sqrt(nb_tx_antennas); //this is amp per tx antenna
|
|
|
|
if (do_amp_compute==1) {
|
|
*amp1p = 0;
|
|
for (aa=0; aa<nb_tx_antennas; aa++) {
|
|
*amp1p += (double)signal_energy((int32_t*)&input[aa][input_offset_meas],length_meas)/NB_RE;
|
|
}
|
|
*amp1p/=nb_tx_antennas;
|
|
*amp1p=sqrt(*amp1p);
|
|
}
|
|
|
|
AssertFatal(amp1p != NULL && *amp1p != 0.0, "Precondition to avoid UB\n");
|
|
|
|
#ifdef DEBUG_DAC
|
|
LOG_I(OCM,"DAC: amp %f, amp1 %f dB (%d,%d), tx_power target %f (actual %f %f),length %d,pos %d\n",
|
|
20*log10(amp),
|
|
20*log10(*amp1p),
|
|
input_offset,
|
|
input_offset_meas,
|
|
txpwr_dBm,
|
|
10*log10((double)signal_energy((int32_t*)&input[0][input_offset_meas],length_meas)/NB_RE),
|
|
10*log10((double)signal_energy((int32_t*)&input[1][input_offset_meas],length_meas)/NB_RE),
|
|
length,
|
|
input_offset_meas);
|
|
|
|
#endif
|
|
|
|
for (i=0; i<length; i++) {
|
|
for (aa=0; aa<nb_tx_antennas; aa++) {
|
|
s_re[aa][i] = amp*((double)(((short *)input[aa]))[((i+input_offset)<<1)])/(*amp1p);
|
|
s_im[aa][i] = amp*((double)(((short *)input[aa]))[((i+input_offset)<<1)+1])/(*amp1p);
|
|
}
|
|
}
|
|
|
|
#ifdef DEBUG_DAC
|
|
LOG_I(OCM,"ener %e\n",signal_energy_fp(s_re,s_im,nb_tx_antennas,length<length_meas?length:length_meas,0));
|
|
#endif
|
|
return(signal_energy_fp(s_re,s_im,nb_tx_antennas,length<length_meas?length:length_meas,0)/NB_RE);
|
|
}
|