%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Copyright: 2010 Integrated Sytems Laboratory, ETH Zurich
%%            http://www.iis.ee.ethz.ch/~sha3
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function x = roundr(x)

for i=1:16,
    x(i+16) = mod(x(i) + x(i+16), 2^32);
end

for i=1:16,
    x(i) = bin2dec(dec2bin(circshift(double( dec2bin(x(i),32)-48 ),[0 -7]))');
end

for i=1:8,
    temp = x(i);
    x(i) = x(i+8);
    x(i+8) = temp;
end

for i=1:16,
    x(i) = bitxor(x(i), x(i+16));
end

for i=17:4:32,
    temp = x(i);
    x(i) = x(i+2);
    x(i+2) = temp;
    
    temp = x(i+1);
    x(i+1) = x(i+3);
    x(i+3) = temp;
end

%--------------------------------------------------------------------------

for i=1:16,
    x(i+16) = mod(x(i) + x(i+16), 2^32);
end

for i=1:16,
    x(i) = bin2dec(dec2bin(circshift(double(dec2bin(x(i),32)-48),[0 -11]))');
end

for i=1:8:16,
    temp = x(i);
    x(i) = x(i+4);
    x(i+4) = temp;
    
    temp = x(i+1);
    x(i+1) = x(i+5);
    x(i+5) = temp;
    
    temp = x(i+2);
    x(i+2) = x(i+6);
    x(i+6) = temp;
    
    temp = x(i+3);
    x(i+3) = x(i+7);
    x(i+7) = temp;
end

for i=1:16,
    x(i) = bitxor(x(i), x(i+16));
end

for i=17:2:32,
    temp = x(i);
    x(i) = x(i+1);
    x(i+1) = temp;
end