From 9d9692d4398e78374948af4b939c9b06ba6abc6c Mon Sep 17 00:00:00 2001 From: Vera Guelpers Date: Thu, 3 May 2018 16:40:16 +0100 Subject: [PATCH] Fix double vs float in boundary phases --- lib/qcd/action/fermion/DomainWallFermion.h | 4 ++-- lib/qcd/action/fermion/FermionOperator.h | 4 ++-- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/lib/qcd/action/fermion/DomainWallFermion.h b/lib/qcd/action/fermion/DomainWallFermion.h index 99f64865..56179f26 100644 --- a/lib/qcd/action/fermion/DomainWallFermion.h +++ b/lib/qcd/action/fermion/DomainWallFermion.h @@ -63,7 +63,7 @@ namespace Grid { LatticeCoordinate(coor, nu + shift); ph = ph + twist[nu]*coor*((1./(in._grid->_fdimensions[nu+shift]))); } - in_buf = exp((RealD)(2.0*M_PI)*ci*ph*(-1.0))*in; + in_buf = exp((Real)(2.0*M_PI)*ci*ph*(-1.0))*in; if(fiveD){//FFT only on temporal and spatial dimensions std::vector mask(Nd+1,1); mask[0] = 0; @@ -78,7 +78,7 @@ namespace Grid { } //phase for boundary condition - out = out * exp((RealD)(2.0*M_PI)*ci*ph); + out = out * exp((Real)(2.0*M_PI)*ci*ph); }; virtual void FreePropagator(const FermionField &in,FermionField &out,RealD mass,std::vector twist) { diff --git a/lib/qcd/action/fermion/FermionOperator.h b/lib/qcd/action/fermion/FermionOperator.h index 1ef99b85..f39e7511 100644 --- a/lib/qcd/action/fermion/FermionOperator.h +++ b/lib/qcd/action/fermion/FermionOperator.h @@ -115,14 +115,14 @@ namespace Grid { LatticeCoordinate(coor, nu); ph = ph + twist[nu]*coor*((1./(in._grid->_fdimensions[nu]))); } - in_buf = exp((RealD)(2.0*M_PI)*ci*ph*(-1.0))*in; + in_buf = exp((Real)(2.0*M_PI)*ci*ph*(-1.0))*in; theFFT.FFT_all_dim(in_k,in_buf,FFT::forward); this->MomentumSpacePropagator(prop_k,in_k,mass,twist); theFFT.FFT_all_dim(out,prop_k,FFT::backward); //phase for boundary condition - out = out * exp((RealD)(2.0*M_PI)*ci*ph); + out = out * exp((Real)(2.0*M_PI)*ci*ph); }; virtual void FreePropagator(const FermionField &in,FermionField &out,RealD mass) {