portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-09-20 09:15:38 +01:00
Code Releases Activity

Chulwoo's Zmobius test

Browse Source
This commit is contained in:
paboyle 2017-03-30 13:30:25 +09:00
parent 2b115929dc
commit cb9a297a0a
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
tests/solver/Test_zmobius_cg_prec.cc
View File

@ -81,10 +81,12 @@ int main(int argc, char** argv) {
RealD M5 = 1.8; RealD M5 = 1.8;
std::vector < std::complex<double> > omegas; std::vector < std::complex<double> > omegas;
for(int i=0;i<Ls;i++){ for(int i=0;i<Ls;i++){
std::complex<double> temp (0.25+0.00*i, 0.0+0.00*i); double imag = 0.;
omegas.push_back(temp); if (i==0) imag=1.;
if (i==Ls-1) imag=-1.;
std::complex<double> temp (0.25+0.01*i, imag*0.01);
omegas.push_back(temp);
} }
// DomainWallFermionR Ddwf(Umu, *FGrid, *FrbGrid, *UGrid, *UrbGrid, mass, M5);
ZMobiusFermionR Ddwf(Umu, *FGrid, *FrbGrid, *UGrid, *UrbGrid, mass, M5, omegas,1.,0.); ZMobiusFermionR Ddwf(Umu, *FGrid, *FrbGrid, *UGrid, *UrbGrid, mass, M5, omegas,1.,0.);
LatticeFermion src_o(FrbGrid); LatticeFermion src_o(FrbGrid);