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

Merge pull request #226 from nils-asmussen/fix/Gauss

Browse Source 
Fix compiling of MSource::Gauss for single precision
This commit is contained in:
Antonin Portelli 2019-08-14 17:50:38 +01:00 committed by GitHub
commit 55c095f620
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
Hadrons/Modules/MSource/Gauss.hpp
View File

@ -114,6 +114,7 @@ void TGauss<FImpl>::execute(void)
const int dim=env().getNd()-1; const int dim=env().getNd()-1;
const Real fact=-0.5/std::pow(par().width,2); const Real fact=-0.5/std::pow(par().width,2);
const Complex i(0.0, 1.0); const Complex i(0.0, 1.0);
const Real Pi(M_PI);
const SitePropagator idMat=[](){ SitePropagator s; s=1.; return s; }(); const SitePropagator idMat=[](){ SitePropagator s; s=1.; return s; }();
ScalarRho=zero; ScalarRho=zero;
@ -129,7 +130,7 @@ void TGauss<FImpl>::execute(void)
LatticeCoordinate(compHelper, mu); LatticeCoordinate(compHelper, mu);
//spatial dimensions of momentum phase //spatial dimensions of momentum phase
ScalarRho+=(i*(mom_[mu]*2*M_PI/Lmu))*component; ScalarRho+=(i*(mom_[mu]*2*Pi/Lmu))*component;
//Gauss distribution //Gauss distribution
component-=Complex(posMu); component-=Complex(posMu);
@ -150,10 +151,10 @@ void TGauss<FImpl>::execute(void)
//time component of momentum phase //time component of momentum phase
LatticeCoordinate(component, dim); LatticeCoordinate(component, dim);
ScalarRho+=(i*(mom_.at(dim)*2*M_PI/env().getDim(dim)))*component; ScalarRho+=(i*(mom_.at(dim)*2*Pi/env().getDim(dim)))*component;
//compute scalar result //compute scalar result
ScalarRho=exp(ScalarRho)*Complex(std::pow(sqrt(2*M_PI)*par().width,-dim)); ScalarRho=exp(ScalarRho)*Complex(std::pow(sqrt(2*Pi)*par().width,-dim));
//select time slices //select time slices
LatticeCoordinate(compHelper, dim); LatticeCoordinate(compHelper, dim);