portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-10 06:00:45 +01:00
Code Releases Activity

Zmobius test was wrong! (only mobius) checking in again

Browse Source
This commit is contained in:
Chulwoo Jung 2017-03-16 23:04:28 -04:00
parent 0b61f75c9e
commit e9712bc7fb
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);