diff --git a/matlab/tfgm/tf_filtering/compute_gamma.m b/matlab/tfgm/tf_filtering/compute_gamma.m
deleted file mode 100644
index b5df28f4a34bc2e2983217bc153daa3593f3dd6c..0000000000000000000000000000000000000000
--- a/matlab/tfgm/tf_filtering/compute_gamma.m
+++ /dev/null
@@ -1,8 +0,0 @@
-function gamma_vec = compute_gamma(s_vec)
-% compute gamma
-
-gamma_vec =@(lambda_coef) lambda_coef*s_vec./(lambda_coef*s_vec + 1- s_vec);
-
-end
-
-
diff --git a/matlab/tfgm/tf_filtering/filtering_out_Pareas.m b/matlab/tfgm/tf_filtering/filtering_out_Pareas.m
deleted file mode 100644
index ece8365bf6a0d67fb1d46d61253a1f4c4a065c43..0000000000000000000000000000000000000000
--- a/matlab/tfgm/tf_filtering/filtering_out_Pareas.m
+++ /dev/null
@@ -1,82 +0,0 @@
-function [x_rec, t_arrf, t_evdn, t_ut_x, rank_q, s_vec_list, u_mat_list,...,
-     lambda_vec_opt] = filtering_out_Pareas(x_mix,  mask_labels, dgt, idgt, e_target, tolerance_arrf, proba_arrf)
-
-%%
-
-r =  compute_r(sig_len, sig_len, proba_arrf);
-
-P = length(mask_labels);
-P= 2;
-t_arrf = zeros(P,1);
-t_evdn = zeros(P,1);
-t_ut_x = zeros(P,1);
-rank_q = zeros(P,1);
-
-s_vec_list =cell(P,1);
-u_mat_list = cell(P,1);
-
-
-
-lambda_vec_opt = zeros(P,1);
-
-%%
-for p = 1:P
-p=2;
-    %%
-    mask_p = mask_labels{p};
-    
-    gabmul_p = gen_gabmul_operator(dgt, idgt, mask_p);
-    
-    figure; imagesc(mask_p);
-
-    fprintf('Random EVD of Gabor multiplier number %.f\n',p)
-    fprintf('#coefs in mask: %.f  environ %f %%\n\n', sum(mask_p(:)),sum(mask_p(:)/(size(mask_p,1)*size(mask_p,2))*100));
-    %%
-    
-    tic;
-    q_mat = adaptative_randomized_range_finder(gabmul_p, sig_len, tolerance_arrf, r);
-    t_arrf(p) = toc;
-    
-    rank_q(p) = size(q_mat,2);
-    fprintf('Q size: %.f %.f\n', size(q_mat,1),size(q_mat,2));
-    %%
-    tic;
-    evdn = EVD_nystrom(gabmul_p, q_mat);
-    t_evdn(p) = toc;
-    
-    u_mat = evdn.U;
-    s_vec = diag(evdn.D);
-    s_vec_list{p} = s_vec;
-    u_mat_list{p} = u_mat;
-    
-    %%
-    
-    fprintf('Running times:\n')
-    fprintf('- adaptive_randomized_range_finder: %f s\n',t_arrf(p));
-    fprintf('- evd_nystrom: %f s \n\n', t_evdn(p));
-
-
- %%
- ut_x = u_mat'*x_mix;
- gamma_p =@(lambda_coef)(lambda_coef*s_vec)./(1- (1-lambda_coef)*s_vec);
-  
- x_p=@(lambda_coef) x_mix-u_mat*(gamma_p(lambda_coef).*ut_x);
-  figure; 
-  plot_spectrogram(x_p(1), dgt_params, signal_params, dgt)
-
-%%
-    obj_fun = @(lambda_coef) abs(e_target(p) - norm(mask_p.*dgt(x_p(lambda_coef))));
-    sol = fminbnd(obj_fun, 0,1);
-    lambda_vec_opt(p) = sol;
-   
-end
-%%
-
-
-
-
-%%
-
-
-end
-