function y=Mult_DWF_H(x,N5);
% Multiplies a vector by Hermitian Conjugate of
% the Domain Wall Fermion matrix
global N mass_dwf

x=reshape(x,12*N,N5);
%
y(:,1)=x(:,1)-Mult_Dirac_W_H(x(:,1))+P5minus(x(:,2))-mass_dwf*P5plus(x(:,N5));
for j5=2:N5-1;
  y(:,j5)=x(:,j5)-Mult_Dirac_W_H(x(:,j5))+P5minus(x(:,j5+1))+P5plus(x(:,j5-1));
end
y(:,N5)=x(:,N5)-Mult_Dirac_W_H(x(:,N5))-mass_dwf*P5minus(x(:,1))+P5plus(x(:,N5-1));

x=reshape(x,12*N*N5,1);
y=reshape(y,12*N*N5,1);
%
% Copyright (C) 2006-2007 Artan Borici.
% This program is a free software licensed under the terms of the GNU General Public License