%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Copyright: 2010 Integrated Sytems Laboratory, ETH Zurich
%%            http://www.iis.ee.ethz.ch/~sha3
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%JH KAT Test
M0=zeros(1,64);
M=zeros(1,64);
H0=zeros(1,128);
H0(1,1)=1;

%iv=jh_256(M0,H0);
H=iv;

M=[M,zeros(1,64)];
M(1,1)=208;
M(1,128)=3;
n=2;

 H(1,129:256) = jh_256(M(1,1:64),H(1,1:128));
  if n>1
    for i=1:n-1
         H(1,(1+(1+i)*128):(i+2)*128) = jh_256(M(1,(1+i*64):((1+i)*64)),H(1,1+(i*128):((i+1)*128)));
    end
    hash = output256(H(1,1+(n*128):(128+n*128)));
  else
    hash = output256(H(1,129:256));
  end

dec2hex(hash)