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);