global N5 mass_dwf
N5=8, mass_dwf=0.1
load u1
load u2
load u3
load u4
Initialise_Dirac_W;
x=rand(12*N*N5,1)+i*rand(12*N*N5,1);
y=Mult_DWF(x,N5);
y=Mult_DWF_H(x,N5);