/ani/mrses

function res=mrses_mtx(A,B,k,Niter,Ncycle,distmod,block)

  if (nargin<7)

    block=256;

  end

  if (nargin<6)

    distmod=1;

  end

  if (nargin<5)

    Ncycle=1000;

  end

  if (nargin<4)

    Niter=500;

  end

  if (nargin<3)

    k=5;

  end

  if (nargin<2)

    error('As minimum two matrixes needed for MRSES');

  end

  if (nargin>6)

    error('Too much parameters');

  end

  sa=size(A);sb=size(B);

  if (sa(2)==sb(2))

    genes=sa(2);

  else

    error('Features dimension mismatch');

  end

  ctx = mrses_hw();

  mrses_hw(ctx, 1, k, block, single(A), single(B), distmod);

  optki=[];

  for icycle=0:block:(Ncycle-1)

    block_size = min(block, Ncycle - icycle);

%    optki=mrses_hw(ctx, 15, Niter, block_size);

    optki = [optki, mrses_hw(ctx, 15, Niter, block_size)];

  end

%  n = mrses_hw(ctx, 16);

  mrses_hw(ctx);

  %H=[double(n)./Ncycle;1:genes];

  %res=flipud(sortrows(H'));

  res=optki;