clustering_the_solutions.m

function [Solution, VectIndex, Cent, reconstr, mean_savg, number_of_clust_sim] = clustering_the_solutions(number_of_sources,nd,sol,normF, Qyes)

minsil_old = -2;
    
    if Qyes == 1
        steps      = 10;
        quants     = quantile(normF, linspace(0.2, 0.01, steps)); %This use to be linspace(0.2, 0.01, steps)     %start with 20% the data 
    else
        steps = 1;
        quants = max(normF);
    end
    
%% Calculating the silhuettes in the different quantiles 
 
for p= 1:steps %steps
       %disp(p);
       ind  = find(normF <= quants(p));
       sol1 = sol(ind,:);
  
       justSource = sol1(:,4:end);
       Sources3D  = zeros( 3, size(justSource,2)/3, size(justSource, 1) );

    for J = 1:size(sol1,1)
        hold = justSource(J,:);

            for kJ = 1:size(justSource,2)/3


                if kJ == 1
                    col = hold(kJ*3-2:kJ*3);
                else
                    col =[col; hold(kJ*3-2:kJ*3)];

                end

            end
        
        Sources3D(:,:,J) = col';

        if J == 1
            colSources = col;
        else  
            colSources = [colSources ; col];
        end

    end
    


    [VectIndex, Cent] = kmeans(colSources, number_of_sources);
     
   
    [ss,h] = silhouette(colSources,VectIndex);
    savg   = grpstats(ss,VectIndex);
    minsil = savg;
    
    if mean(minsil) < minsil_old/2
        break
    end
    
    if size(ind,1) < 5
        break
    end
    
    if min(minsil) > 0.95
        break
    end
    
    minsil_old = mean(minsil);


end

number_of_clust_sim = p;





idx1     = VectIndex==1;
avg_sol  = mean(sol1(VectIndex(idx1),:));
Solution = zeros(number_of_sources,6);

    for i = 1:number_of_sources
    Solution(i,:)= [Cent(i,:) avg_sol(1:3)];
    end
    
reconstr  = mean(normF( ind));
mean_savg = mean(savg);

file_name1 = sprintf('./Results/Solution_%ddet_%dsources.mat',nd, number_of_sources);
save(file_name1, 'Solution', 'VectIndex', 'Cent', 'ss', 'savg', 'reconstr', 'mean_savg','number_of_clust_sim');   

end