From 8925654c078e762fd62e354ba82bc1824eeae18f Mon Sep 17 00:00:00 2001
From: Robert Schmidt <robert.schmidt@openairinterface.org>
Date: Mon, 16 Feb 2026 14:57:45 +0100
Subject: [PATCH] Remove unused 4G MMSE code

Since the parent commit, USE_MMSE is not defined. Thus, the
corresponding code is commented out. Since there are no linker problems,
it is obvious that this code is never called, and we can remove it.
---
 CMakeLists.txt                                |   1 -
 .../PHY/LTE_UE_TRANSPORT/dlsch_demodulation.c | 434 ------------------
 .../linear_preprocessing_rec.c                | 119 -----
 .../linear_preprocessing_rec.h                |  12 -
 4 files changed, 566 deletions(-)
 delete mode 100644 openair1/PHY/LTE_UE_TRANSPORT/linear_preprocessing_rec.c
 delete mode 100644 openair1/PHY/LTE_UE_TRANSPORT/linear_preprocessing_rec.h

diff --git a/CMakeLists.txt b/CMakeLists.txt
index 4fd22fb1c2..98dba37442 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -926,7 +926,6 @@ set(PHY_SRC_UE
   ${OPENAIR1_DIR}/PHY/LTE_TRANSPORT/sss_gen.c
   ${OPENAIR1_DIR}/PHY/LTE_UE_TRANSPORT/dlsch_demodulation.c
   ${OPENAIR1_DIR}/PHY/LTE_UE_TRANSPORT/dlsch_llr_computation.c
-  #${OPENAIR1_DIR}/PHY/LTE_UE_TRANSPORT/linear_preprocessing_rec.c
   ${OPENAIR1_DIR}/PHY/LTE_UE_TRANSPORT/dlsch_decoding.c
   ${OPENAIR1_DIR}/PHY/LTE_UE_TRANSPORT/dci_tools_ue.c
   ${OPENAIR1_DIR}/PHY/LTE_UE_TRANSPORT/uci_tools_ue.c
diff --git a/openair1/PHY/LTE_UE_TRANSPORT/dlsch_demodulation.c b/openair1/PHY/LTE_UE_TRANSPORT/dlsch_demodulation.c
index cf83b40a7f..b8157652c2 100644
--- a/openair1/PHY/LTE_UE_TRANSPORT/dlsch_demodulation.c
+++ b/openair1/PHY/LTE_UE_TRANSPORT/dlsch_demodulation.c
@@ -41,10 +41,6 @@
 #include <stdlib.h>
 #include <string.h>
 #include <linux/version.h>
-#include "linear_preprocessing_rec.h"
-
-//#define DEBUG_MMSE
-//#define USE_MMSE
 
 /* dynamic shift for LLR computation for TM3/4
  * set as command line argument, see lte-softmodem.c
@@ -348,436 +344,6 @@ static void prec2A_TM4_128(int pmi, simde__m128i *ch0, simde__m128i *ch1)
   //print_shorts("prec2A_TM4 ch1 (end):",ch1);
 }
 
-#ifdef USE_MMSE
-static void precode_channel_est(int32_t **dl_ch_estimates_ext,
-                                LTE_DL_FRAME_PARMS *frame_parms,
-                                LTE_UE_PDSCH *pdsch_vars,
-                                unsigned char symbol,
-                                unsigned short nb_rb,
-                                MIMO_mode_t mimo_mode)
-{
-  unsigned short rb;
-  simde__m128i *dl_ch0_128,*dl_ch1_128;
-  unsigned char aarx=0,symbol_mod,pilots=0;
-  unsigned char *pmi_ext = pdsch_vars->pmi_ext;
-  symbol_mod = (symbol>=(7-frame_parms->Ncp)) ? symbol-(7-frame_parms->Ncp) : symbol;
-
-  if ((symbol_mod == 0) || (symbol_mod == (4-frame_parms->Ncp)))
-    pilots=1;
-
-  for (aarx=0; aarx<frame_parms->nb_antennas_rx; aarx++) {
-    dl_ch0_128          = (simde__m128i *)&dl_ch_estimates_ext[aarx][symbol*frame_parms->N_RB_DL*12]; // this is h11
-    dl_ch1_128          = (simde__m128i *)&dl_ch_estimates_ext[2+aarx][symbol*frame_parms->N_RB_DL*12]; // this is h12
-
-    for (rb=0; rb<nb_rb; rb++) {
-      if (mimo_mode==LARGE_CDD) {
-        prec2A_TM3_128(&dl_ch0_128[0],&dl_ch1_128[0]);
-        prec2A_TM3_128(&dl_ch0_128[1],&dl_ch1_128[1]);
-
-        if (pilots==0) {
-          prec2A_TM3_128(&dl_ch0_128[2],&dl_ch1_128[2]);
-        }
-      } else if (mimo_mode==DUALSTREAM_UNIFORM_PRECODING1) {
-        prec2A_TM4_128(0,&dl_ch0_128[0],&dl_ch1_128[0]);
-        prec2A_TM4_128(0,&dl_ch0_128[1],&dl_ch1_128[1]);
-
-        if (pilots==0) {
-          prec2A_TM4_128(0,&dl_ch0_128[2],&dl_ch1_128[2]);
-        }
-      } else if (mimo_mode==DUALSTREAM_UNIFORM_PRECODINGj) {
-        prec2A_TM4_128(1,&dl_ch0_128[0],&dl_ch1_128[0]);
-        prec2A_TM4_128(1,&dl_ch0_128[1],&dl_ch1_128[1]);
-
-        if (pilots==0) {
-          prec2A_TM4_128(1,&dl_ch0_128[2],&dl_ch1_128[2]);
-        }
-      } else if (mimo_mode==DUALSTREAM_PUSCH_PRECODING) {
-        prec2A_TM4_128(pmi_ext[rb],&dl_ch0_128[0],&dl_ch1_128[0]);
-        prec2A_TM4_128(pmi_ext[rb],&dl_ch0_128[1],&dl_ch1_128[1]);
-
-        if (pilots==0) {
-          prec2A_TM4_128(pmi_ext[rb],&dl_ch0_128[2],&dl_ch1_128[2]);
-        }
-      } else {
-        LOG_E(PHY,"Unknown MIMO mode\n");
-        return;
-      }
-
-      if (pilots==0) {
-        dl_ch0_128+=3;
-        dl_ch1_128+=3;
-      } else {
-        dl_ch0_128+=2;
-        dl_ch1_128+=2;
-      }
-    }
-  }
-}
-
-static void mult_mmse_chan_est(float complex **Wmmse,
-                        float complex **dl_ch_estimates_ext_f,
-                        int n_tx,
-                        int n_rx,
-                        int length,
-                        int start_point) {
-  short re;
-  int aarx, aatx;
-  float complex *chan_est_re =  malloc(n_tx*n_rx*sizeof(float complex));
-  float complex *chan_est_mmse_re =  malloc(n_tx*n_rx*sizeof(float complex));
-  float complex *Wmmse_re =  malloc(n_tx*n_rx*sizeof(float complex));
-
-  for (re=0; re<length; re++) {
-    for (aatx=0; aatx<n_tx; aatx++) {
-      for (aarx=0; aarx<n_rx; aarx++) {
-        chan_est_re[aatx*n_rx + aarx] = dl_ch_estimates_ext_f[aatx*n_rx + aarx][re];
-        Wmmse_re[aatx*n_rx + aarx] = Wmmse[aatx*n_rx + aarx][re];
-#ifdef DEBUG_MMSE
-
-        if (re==0)
-          printf("mult_mmse_chan_est: chan_est_re[%d] = (%f, %f)\n", aatx*n_rx + aarx, creal(chan_est_re[aatx*n_rx + aarx]), cimag(chan_est_re[aatx*n_rx + aarx]));
-
-#endif
-      }
-    }
-
-    mutl_matrix_matrix_col_based(Wmmse_re, chan_est_re, n_rx, n_tx, n_rx, n_tx, chan_est_mmse_re);
-
-    for (aatx=0; aatx<n_tx; aatx++) {
-      for (aarx=0; aarx<n_rx; aarx++) {
-        dl_ch_estimates_ext_f[aatx*n_rx + aarx][re] = chan_est_mmse_re[aatx*n_rx + aarx];
-#ifdef DEBUG_MMSE
-
-        if (re==0)
-          printf("mult_mmse_chan_est: dl_ch_estimates_ext_f[%d][%d] = (%f, %f)\n", aatx*n_rx + aarx, re, creal(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re]), cimag(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re]));
-
-#endif
-      }
-    }
-  }
-
-  free(Wmmse_re);
-  free(chan_est_re);
-  free(chan_est_mmse_re);
-}
-static void mult_mmse_rxdataF(float complex **Wmmse,
-                              float complex **rxdataF_ext_f,
-                              int n_tx,
-                              int n_rx,
-                              int length,
-                              int start_point) {
-  short re;
-  int aarx, aatx;
-  float complex rxdata_re[n_rx];
-  float complex rxdata_mmse_re[n_rx];
-  float complex Wmmse_re[n_tx][n_rx];
-
-  for (re=0; re<length; re++) {
-    for (aarx=0; aarx<n_rx; aarx++) {
-      rxdata_re[aarx] = rxdataF_ext_f[aarx][re];
-#ifdef DEBUG_MMSE
-
-      if (re==0)
-        printf("mult_mmse_rxdataF before: rxdata_re[%d] = (%f, %f)\n", aarx, creal(rxdata_re[aarx]), cimag(rxdata_re[aarx]));
-
-#endif
-    }
-
-    for (aatx=0; aatx<n_tx; aatx++) {
-      for (aarx=0; aarx<n_rx; aarx++) {
-        Wmmse_re[aatx][aarx] = Wmmse[aatx * n_rx + aarx][re];
-      }
-    }
-
-    mutl_matrix_matrix_col_based((float complex *)Wmmse_re, rxdata_re, n_rx, n_tx, n_rx, 1, rxdata_mmse_re);
-
-    for (aarx=0; aarx<n_rx; aarx++) {
-      rxdataF_ext_f[aarx][re] = rxdata_mmse_re[aarx];
-#ifdef DEBUG_MMSE
-
-      if (re==0)
-        printf("mult_mmse_rxdataF after: rxdataF_ext_f[%d] = (%f, %f)\n", aarx, creal(rxdataF_ext_f[aarx][re]), cimag(rxdataF_ext_f[aarx][re]));
-
-#endif
-    }
-  }
-}
-static void float_to_rxdataF(int32_t **rxdataF_ext, float complex **rxdataF_f, int n_tx, int n_rx, int length, int start_point)
-{
-  for (int aarx = 0; aarx < n_rx; aarx++) {
-    for (int re = 0; re < length; re++) {
-      int16_t imag, real;
-      if (cimag(rxdataF_f[aarx][re])<-1)
-        imag = -INT16_MAX;
-      else if (cimag(rxdataF_f[aarx][re])>=1)
-        imag = INT16_MAX;
-      else
-        imag = cimag(rxdataF_f[aarx][re])*32768;
-
-      if (creal(rxdataF_f[aarx][re])<-1)
-        real = -INT16_MAX;
-      else if (creal(rxdataF_f[aarx][re])>=1)
-        real = INT16_MAX;
-      else
-        real = creal(rxdataF_f[aarx][re])*32768;
-
-      rxdataF_ext[aarx][start_point + re] = (((int32_t)imag)<<16)|(((int32_t)real) & 0xffff);
-#ifdef DEBUG_MMSE
-
-      if (re==0) {
-        printf(" float_to_rxdataF: real = %f, imag = %f\n",creal(rxdataF_f[aarx][re]), cimag(rxdataF_f[aarx][re]));
-        printf("float_to_rxdataF: real fixed = %d, imag fixed = %d\n", real, imag);
-        printf("float_to_rxdataF: ant %d, re = %d, rxdataF_ext = %d \n", aarx, re,  rxdataF_ext[aarx][start_point + re]);
-      }
-
-#endif
-    }
-  }
-}
-
-static void float_to_chan_est(int32_t **dl_ch_estimates_ext,
-                       float complex **dl_ch_estimates_ext_f,
-                       int n_tx,
-                       int n_rx,
-                       int length,
-                       int start_point)
-{
-  for (int aatx = 0; aatx < n_tx; aatx++) {
-    for (int aarx = 0; aarx < n_rx; aarx++) {
-      for (int re = 0; re < length; re++) {
-        int16_t imag, real;
-        if (cimag(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re])<-1)
-          imag = -INT16_MAX;
-        else if (cimag(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re])>=1)
-          imag = INT16_MAX;
-        else
-          imag = cimag(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re])*32768;
-
-        if (creal(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re])<-1)
-          real = -INT16_MAX;
-        else if (creal(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re])>=1)
-          real = INT16_MAX;
-        else
-          real = creal(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re])*32768;
-
-        dl_ch_estimates_ext[aatx*n_rx + aarx][start_point + re] = (((int32_t)imag)<<16)|((int32_t)real & 0xffff);
-#ifdef DEBUG_MMSE
-
-        if (re==0) {
-          printf(" float_to_chan_est: chan est real = %f, chan est imag = %f\n",creal(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re]), cimag(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re]));
-          printf("float_to_chan_est: real fixed = %d, imag fixed = %d\n", real, imag);
-          printf("float_to_chan_est: ant %d, re = %d, dl_ch_estimates_ext = %d \n", aatx*n_rx + aarx, re,  dl_ch_estimates_ext[aatx*n_rx + aarx][start_point + re]);
-        }
-
-#endif
-      }
-    }
-  }
-}
-
-static void chan_est_to_float(int32_t **dl_ch_estimates_ext,
-                              float complex **dl_ch_estimates_ext_f,
-                              int n_tx,
-                              int n_rx,
-                              int length,
-                              int start_point)
-{
-  for (int aatx = 0; aatx < n_tx; aatx++) {
-    for (int aarx = 0; aarx < n_rx; aarx++) {
-      for (int re = 0; re < length; re++) {
-        int16_t imag = (int16_t)(dl_ch_estimates_ext[aatx * n_rx + aarx][start_point + re] >> 16);
-        int16_t real = (int16_t)(dl_ch_estimates_ext[aatx * n_rx + aarx][start_point + re] & 0xffff);
-        dl_ch_estimates_ext_f[aatx*n_rx + aarx][re] = (float)(real/(32768.0)) + I*(float)(imag/(32768.0));
-#ifdef DEBUG_MMSE
-
-        if (re==0) {
-          printf("ant %d, re = %d, real = %d, imag = %d \n", aatx*n_rx + aarx, re, real, imag);
-          printf("ant %d, re = %d, real = %f, imag = %f \n", aatx*n_rx + aarx, re, creal(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re]), cimag(dl_ch_estimates_ext_f[aatx*n_rx + aarx][re]));
-        }
-
-#endif
-      }
-    }
-  }
-}
-
-static void rxdataF_to_float(int32_t **rxdataF_ext, float complex **rxdataF_f, int n_rx, int length, int start_point)
-{
-  for (int aarx = 0; aarx < n_rx; aarx++) {
-    for (int re = 0; re < length; re++) {
-      int16_t imag = (int16_t)(rxdataF_ext[aarx][start_point + re] >> 16);
-      int16_t real = (int16_t)(rxdataF_ext[aarx][start_point + re] & 0xffff);
-      rxdataF_f[aarx][re] = (float)(real/(32768.0)) + I*(float)(imag/(32768.0));
-#ifdef DEBUG_MMSE
-
-      if (re==0) {
-        printf("rxdataF_to_float: aarx = %d, real= %d, imag = %d\n", aarx, real, imag);
-        //printf("rxdataF_to_float: rxdataF_ext[%d][%d] = %d\n", aarx, start_point + re, rxdataF_ext[aarx][start_point + re]);
-        //printf("rxdataF_to_float: ant %d, re = %d, rxdataF_f real = %f, rxdataF_f imag = %f \n", aarx, re, creal(rxdataF_f[aarx][re]), cimag(rxdataF_f[aarx][re]));
-      }
-
-#endif
-    }
-  }
-}
-
-static void mmse_processing_core(int32_t **rxdataF_ext, 
-                          int32_t **dl_ch_estimates_ext,
-                          int noise_power,
-                          int n_tx,
-                          int n_rx,
-                          int length,
-                          int start_point) {
-  int aatx, aarx, re;
-  float imag;
-  float real;
-  float complex **W_MMSE= malloc(n_tx*n_rx*sizeof(float complex *));
-
-  for (int j=0; j<n_tx*n_rx; j++) {
-    W_MMSE[j] = malloc(sizeof(float complex)*length);
-  }
-
-  float complex *H=  malloc(n_tx*n_rx*sizeof(float complex));
-  float complex *W_MMSE_re=  malloc(n_tx*n_rx*sizeof(float complex));
-  float complex **dl_ch_estimates_ext_flcpx = malloc(n_tx*n_rx*sizeof(float complex *));
-
-  for (int j=0; j<n_tx*n_rx; j++) {
-    dl_ch_estimates_ext_flcpx[j] = malloc(sizeof(float complex)*length);
-  }
-
-  float complex **rxdataF_ext_flcpx = malloc(n_rx*sizeof(float complex *));
-
-  for (int j=0; j<n_rx; j++) {
-    rxdataF_ext_flcpx[j] = malloc(sizeof(float complex)*length);
-  }
-
-  chan_est_to_float(dl_ch_estimates_ext,
-                    dl_ch_estimates_ext_flcpx,
-                    n_tx,
-                    n_rx,
-                    length,
-                    start_point);
-
-  for (re=0; re<length; re++) {
-    for (aatx=0; aatx<n_tx; aatx++) {
-      for (aarx=0; aarx<n_rx; aarx++) {
-        imag = cimag(dl_ch_estimates_ext_flcpx[aatx*n_rx + aarx][re]);
-        real = creal(dl_ch_estimates_ext_flcpx[aatx*n_rx + aarx][re]);
-        H[aatx*n_rx + aarx] = real+ I*imag;
-      }
-    }
-
-    compute_MMSE(H, n_tx, noise_power, W_MMSE_re);
-
-    for (aatx=0; aatx<n_tx; aatx++) {
-      for (aarx=0; aarx<n_rx; aarx++) {
-        W_MMSE[aatx*n_rx + aarx][re] = W_MMSE_re[aatx*n_rx + aarx];
-      }
-    }
-  }
-
-  rxdataF_to_float(rxdataF_ext,
-                   rxdataF_ext_flcpx,
-                   n_rx,
-                   length,
-                   start_point);
-  mult_mmse_rxdataF(W_MMSE,
-                    rxdataF_ext_flcpx,
-                    n_tx,
-                    n_rx,
-                    length,
-                    start_point);
-  mult_mmse_chan_est(W_MMSE,
-                     dl_ch_estimates_ext_flcpx,
-                     n_tx,
-                     n_rx,
-                     length,
-                     start_point);
-  float_to_rxdataF(rxdataF_ext,
-                   rxdataF_ext_flcpx,
-                   n_tx,
-                   n_rx,
-                   length,
-                   start_point);
-  float_to_chan_est(dl_ch_estimates_ext,
-                    dl_ch_estimates_ext_flcpx,
-                    n_tx,
-                    n_rx,
-                    length,
-                    start_point);
-  free(W_MMSE);
-  free(H);
-  free(W_MMSE_re);
-  free(dl_ch_estimates_ext_flcpx);
-  free(rxdataF_ext_flcpx);
-}
-
-static void mmse_processing_oai(LTE_UE_PDSCH *pdsch_vars,
-                                LTE_DL_FRAME_PARMS *frame_parms,
-                                PHY_MEASUREMENTS *measurements,
-                                unsigned char first_symbol_flag,
-                                MIMO_mode_t mimo_mode,
-                                unsigned short mmse_flag,
-                                int noise_power,
-                                unsigned char symbol,
-                                unsigned short nb_rb)
-{
-  int **rxdataF_ext           = pdsch_vars->rxdataF_ext;
-  int **dl_ch_estimates_ext   = pdsch_vars->dl_ch_estimates_ext;
-  unsigned char *pmi_ext      = pdsch_vars->pmi_ext;
-  int avg_00[frame_parms->nb_antenna_ports_eNB*frame_parms->nb_antennas_rx];
-  int avg_01[frame_parms->nb_antenna_ports_eNB*frame_parms->nb_antennas_rx];
-  int symbol_mod, length, start_point, nre;
-  symbol_mod = (symbol>=(7-frame_parms->Ncp)) ? symbol-(7-frame_parms->Ncp) : symbol;
-
-  if (((symbol_mod == 0) || (symbol_mod == (frame_parms->Ncp-1)))&&(frame_parms->nb_antenna_ports_eNB!=1))
-    nre=8;
-  else if (((symbol_mod == 0) || (symbol_mod == (frame_parms->Ncp-1)))&&(frame_parms->nb_antenna_ports_eNB==1))
-    nre=10;
-  else
-    nre=12;
-
-  length  = nre*nb_rb;
-  start_point = symbol*nb_rb*12;
-  mmse_processing_core(rxdataF_ext,
-                       dl_ch_estimates_ext,
-                       noise_power,
-                       frame_parms->nb_antenna_ports_eNB,
-                       frame_parms->nb_antennas_rx,
-                       length,
-                       start_point);
-
-  /*dlsch_channel_aver_band(dl_ch_estimates_ext,
-                          frame_parms,
-                          chan_avg,
-                          symbol,
-                          nb_rb);
-
-
-   for (aatx=0; aatx<frame_parms->nb_antenna_ports_eNB; aatx++)
-     for (aarx=0; aarx<frame_parms->nb_antennas_rx; aarx++) {
-       H[aatx*frame_parms->nb_antennas_rx + aarx] = (float)(chan_avg[aatx*frame_parms->nb_antennas_rx + aarx].r/(32768.0)) + I*(float)(chan_avg[aatx*frame_parms->nb_antennas_rx + aarx].i/(32768.0));
-      // printf("H [%d] = (%f, %f) \n", aatx*frame_parms->nb_antennas_rx + aarx, creal(H[aatx*frame_parms->nb_antennas_rx + aarx]), cimag(H[aatx*frame_parms->nb_antennas_rx + aarx]));
-   }*/
-
-  if (first_symbol_flag == 1) {
-    dlsch_channel_level_TM34(dl_ch_estimates_ext,
-                             frame_parms,
-                             pmi_ext,
-                             avg_00,
-                             avg_01,
-                             symbol,
-                             nb_rb,
-                             mmse_flag,
-                             mimo_mode);
-    avg_00[0] = (log2_approx(avg_00[0])/2) + dlsch_demod_shift+4;// + 2 ;//+ 4;
-    avg_01[0] = (log2_approx(avg_01[0])/2) + dlsch_demod_shift+4;// + 2 ;//+ 4;
-    pdsch_vars->log2_maxh0 = cmax(avg_00[0],0);
-    pdsch_vars->log2_maxh1 = cmax(avg_01[0],0);
-  }
-}
-
-#endif
-
 int rx_pdsch(PHY_VARS_UE *ue,
              PDSCH_t type,
              unsigned char eNB_id,
diff --git a/openair1/PHY/LTE_UE_TRANSPORT/linear_preprocessing_rec.c b/openair1/PHY/LTE_UE_TRANSPORT/linear_preprocessing_rec.c
deleted file mode 100644
index c1c1672088..0000000000
--- a/openair1/PHY/LTE_UE_TRANSPORT/linear_preprocessing_rec.c
+++ /dev/null
@@ -1,119 +0,0 @@
-/* These functions compute linear preprocessing for
-the UE using LAPACKE and CBLAS modules of
-LAPACK libraries.
-MMSE and MMSE whitening filters are available.
-Functions are using RowMajor storage of the
-matrices, like in conventional C. Traditional
-Fortran functions of LAPACK employ ColumnMajor
-data storage. */
-
-#include<stdio.h>
-#include<math.h>
-#include<complex.h>
-#include <stdlib.h>
-#include <linux/version.h>
-#include <cblas.h>
-#include <string.h>
-#include <linux/version.h>
-#include <lapacke_utils.h>
-#include <lapacke.h>
-//#define DEBUG_PREPROC
-#include "common/utils/utils.h"
-
-static void conjugate_transpose(int N, float complex* A, float complex* Result)
-{
-  // Computes C := alpha*op(A)*op(B) + beta*C,
-  float complex alpha = 1.0 + I * 0;
-  float complex beta = 0.0+I*0;
-  float complex B[N][N];
-  for (int i = 0; i < N; i++)
-    B[i][i] = 1.0 + I * 0;
-  cblas_cgemm(CblasRowMajor, CblasConjTrans, CblasNoTrans, N, N, N, &alpha, A, N, B, N, &beta, Result, N);
-}
-
-static void H_hermH_plus_sigma2I(int N, float complex* A, float sigma2, float complex Result[N][N])
-{
-  // C := alpha*op(A)*op(B) + beta*C,
-  float complex alpha = 1.0 + I * 0;
-  float complex beta = 1.0 + I * 0;
-  for (int i = 0; i < N; i++)
-    Result[i][i] = sigma2 * (1.0 + I * 0);
-  cblas_cgemm(CblasRowMajor, CblasConjTrans, CblasNoTrans, N, N, N, &alpha, A, N, A, N, &beta, Result, N);
-}
-
-static void lin_eq_solver(int N, float complex* A, float complex* B, float complex* Result)
-{
-  int IPIV[N][N];
-
-  // Compute LU-factorization
-  LAPACKE_cgetrf(LAPACK_ROW_MAJOR, N, N, A, N, (int*)IPIV);
-
-  // Solve AX=B
-  LAPACKE_cgetrs(LAPACK_ROW_MAJOR, 'N', N, N, A, N, (int*)IPIV, B, N);
-
-  // cgetrs( "N", N, 4, A, lda, IPIV, B, ldb, INFO )
-
-  memcpy(Result, B, N * N * sizeof(float complex));
-}
-
-void mutl_matrix_matrix_col_based(float complex* M0,
-                                  float complex* M1,
-                                  int rows_M0,
-                                  int col_M0,
-                                  int rows_M1,
-                                  int col_M1,
-                                  float complex* Result)
-{
-  float complex alpha = 1.0;
-  float complex beta = 0.0;
-  cblas_cgemm(CblasColMajor,
-              CblasNoTrans,
-              CblasNoTrans,
-              rows_M0,
-              col_M1,
-              col_M0,
-              &alpha,
-              M0,
-              col_M0,
-              M1,
-              rows_M1,
-              &beta,
-              Result,
-              rows_M1);
-
-#ifdef DEBUG_PREPROC
-  for (i = 0; i < rows_M0 * col_M1; ++i)
-    printf(" result[%d] = (%f + i%f)\n", i , creal(Result[i]), cimag(Result[i]));
-#endif
-}
-
-/*FILTERS */
-void compute_MMSE(float complex* H, int order_H, float sigma2, float complex* W_MMSE)
-{
-  int N = order_H;
-  float complex H_hermH_sigmaI[N][N];
-  memset(H_hermH_sigmaI, 0, sizeof(H_hermH_sigmaI));
-  H_hermH_plus_sigma2I(N, H, sigma2, H_hermH_sigmaI);
-
-#ifdef DEBUG_PREPROC
-  int i =0;
-  for(i=0;i<N*N;i++)
-    printf(" H_hermH_sigmaI[%d] = (%f + i%f)\n", i , creal(H_hermH_sigmaI[i]), cimag(H_hermH_sigmaI[i]));
-#endif
-
-  float complex H_herm[N][N];
-  memset(H_herm, 0, sizeof(H_herm));
-  conjugate_transpose(N, H, (float complex*)H_herm); // equals H_herm
-
-#ifdef DEBUG_PREPROC
-  for(i=0;i<N*N;i++)
-    printf(" H_herm[%d] = (%f + i%f)\n", i , creal(H_herm[i]), cimag(H_herm[i]));
-#endif
-
-  lin_eq_solver(N, (float complex*)H_hermH_sigmaI, (float complex*)H_herm, W_MMSE);
-
-#ifdef DEBUG_PREPROC
-  for(i=0;i<N*N;i++)
-    printf(" W_MMSE[%d] = (%f + i%f)\n", i , creal(W_MMSE[i]), cimag(W_MMSE[i]));
-#endif
-}
diff --git a/openair1/PHY/LTE_UE_TRANSPORT/linear_preprocessing_rec.h b/openair1/PHY/LTE_UE_TRANSPORT/linear_preprocessing_rec.h
deleted file mode 100644
index a8b955acfd..0000000000
--- a/openair1/PHY/LTE_UE_TRANSPORT/linear_preprocessing_rec.h
+++ /dev/null
@@ -1,12 +0,0 @@
-
-#include<stdio.h>
-#include<math.h>
-#include<complex.h>
-#include <stdlib.h>
-#include "PHY/defs_UE.h"
-
-/* FUNCTIONS FOR LINEAR PREPROCESSING: MMSE, WHITENNING, etc*/
-/* mutl_matrix_matrix_col_based performs multiplications matrix is column-oriented H[0], H[2]; H[1], H[3]*/
-void mutl_matrix_matrix_col_based(float complex* M0, float complex* M1, int rows_M0, int col_M0, int rows_M1, int col_M1, float complex* Result );
-
-void compute_MMSE(float complex* H, int order_H, float sigma2, float complex* W_MMSE);