portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

Cleaning up previous debug lines

Browse Source
This commit is contained in:
Guido Cossu 2016-12-13 07:53:43 +00:00
parent 5c74b6028b
commit 2fb92dbc6e
4 changed files with 4 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
lib/algorithms/LinearOperator.h
View File

@ -234,17 +234,11 @@ namespace Grid {
virtual RealD MpcDag (const Field &in, Field &out){
Field tmp(in._grid);
std::cout << "norm in :" << norm2(in) << std::endl;
_Mat.MeooeDag(in,tmp);
std::cout << "norm tmp :" << norm2(tmp) << std::endl;
_Mat.MooeeInvDag(tmp,out);
std::cout << "norm out :" << norm2(out) << std::endl;
_Mat.MeooeDag(out,tmp);
std::cout << "norm tmp :" << norm2(tmp) << std::endl;
_Mat.MooeeDag(in,out);
std::cout << "norm out :" << norm2(out) << std::endl;
return axpy_norm(out,-1.0,tmp,out);
}
};

4
lib/qcd/action/fermion/CayleyFermion5D.cc
View File

@ -214,9 +214,7 @@ void CayleyFermion5D<Impl>::MeooeDag5D (const FermionField &psi, FermionField
std::vector<Coeff_t> lower=cs;
upper[Ls-1]=-mass*upper[Ls-1];
lower[0] =-mass*lower[0];
std::cout << "MeooeDag5D: psi: " << norm2(psi) << std::endl;
M5Ddag(psi,psi,Din,lower,diag,upper);
std::cout << "MeooeDag5D: Din: " << norm2(Din) << std::endl;
}
template<class Impl>
@ -282,9 +280,7 @@ void CayleyFermion5D<Impl>::MeooeDag (const FermionField &psi, FermionField &
} else {
this->DhopOE(psi,tmp,DaggerYes);
}
std::cout << "MeooeDag: tmp: " << norm2(tmp) << std::endl;
MeooeDag5D(tmp,chi);
std::cout << "MeooeDag: chi: " << norm2(chi) << std::endl;
}
template<class Impl>

6
tests/hmc/Test_hmc_EODWFRatio.cc
View File

@ -60,9 +60,9 @@ public:
// Gauge action
WilsonGaugeActionR Waction(5.6);
Real mass=0.04;
Real pv =1.0;
RealD M5=1.5;
Real mass= 0.04;
Real pv = 1.0;
RealD M5 = 1.5;
FermionAction DenOp(U,*FGrid,*FrbGrid,*UGrid,*UrbGrid,mass,M5);
FermionAction NumOp(U,*FGrid,*FrbGrid,*UGrid,*UrbGrid,pv,M5);

1
tests/solver/Test_dwf_cg_prec_LsVec.cc
View File

@ -64,6 +64,7 @@ int main(int argc, char** argv) {
GridParallelRNG RNG4(UGrid);
RNG4.SeedFixedIntegers(seeds4);
std::cout << GridLogMessage << "Generating source field" << std::endl;
LatticeFermion src(FGrid);
random(RNG5, src);
LatticeFermion result(FGrid);