KMM.m
function [laKMM, laMM, BiGraph, A, OBJ, Ah, laKMMh] = KMM_mmconv(X, c, m, k)
% [laKMM, laMM, BiGraph..., Anc, ~, ~, ~]= KMM(X', c, m,k) : K-Multiple-Means
% Input:
% - X: the data matrix of size nFea...4;
m=floor(sqrt(n*c));k=5; fig = 1;
tic
[laKMM,~,~,A,~,Ah,laKMMh ]= KMM(X', c, m,k) ;
toc;...cm = colormap(jet(c+2));
for i=1:c
plot(X(laKMM==rl(i),1),X(laKMM==rl(i),2),'*', 'color',...',0.3*cm(c+2,:),'MarkerSize',5); hold on;
end
%
rl = randperm(c);
for j=1:(size(Ah,2)/m)
A2=Ah(: