function [u_out,ef,s,tmbs,lambda] = signature_full(u_in)
     
  [N,M]       = size(u_in);     % u_in should be 8 bit binary image
   S          = max(N,M);       % use this to set scales
   nNeighbors = 16;             % we use 16 neighbors
   scale      = floor(S/4);     % use this to choose range of lambda's
   lambda     = 1./(1:scale);   % lambda's chosen

  for i = 1:scale

         u_out_t = Graph_anisoTV_L1_v2(u_in,lambda(i),nNeighbors,2);   % do not change the last input argument
	 u_out(i,:,:) = u_out_t;
         [e,s,t] = energy_full(u_out_t,u_in,lambda(i));
         output(i,:) = [e,s,t]/255;

  end

        ef  = output(:,1);
         s  = output(:,2);
      tmbs  = output(:,3);