portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-11-04 14:04:32 +00:00
Code Releases Activity

more stochastic QED fixes

Browse Source
This commit is contained in:
Antonin Portelli
2016-10-21 15:20:12 +01:00
parent 462921e549
commit 7c8f79b147
2 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
lib/qcd/action/Actions.h
View File

@@ -57,6 +57,7 @@ Author: paboyle <paboyle@ph.ed.ac.uk>
//////////////////////////////////////////// ////////////////////////////////////////////
// Gauge Actions // Gauge Actions
//////////////////////////////////////////// ////////////////////////////////////////////
#include <Grid/qcd/action/gauge/Photon.h>
#include <Grid/qcd/action/gauge/WilsonGaugeAction.h> #include <Grid/qcd/action/gauge/WilsonGaugeAction.h>
#include <Grid/qcd/action/gauge/PlaqPlusRectangleAction.h> #include <Grid/qcd/action/gauge/PlaqPlusRectangleAction.h>

4
lib/qcd/action/gauge/Photon.h
View File

@@ -124,7 +124,7 @@ namespace QCD{
LatticeCoordinate(coor,d); LatticeCoordinate(coor,d);
spNrm = spNrm + coor*coor; spNrm = spNrm + coor*coor;
} }
out = where(spNrm == 0, 0.*out, out); out = where(spNrm == Integer(0), 0.*out, out);
break; break;
} }
@@ -149,7 +149,7 @@ namespace QCD{
template<class Gimpl> template<class Gimpl>
void Photon<Gimpl>::StochasticField(GaugeField &out, GridParallelRNG &rng) void Photon<Gimpl>::StochasticField(GaugeField &out, GridParallelRNG &rng)
{ {
auto *grid = out._grid; auto *grid = dynamic_cast<GridCartesian *>(out._grid);
const unsigned int nd = grid->_ndimension; const unsigned int nd = grid->_ndimension;
GaugeLinkField sqrtK2Inv(grid), r(grid); GaugeLinkField sqrtK2Inv(grid), r(grid);
GaugeField aTilde(grid); GaugeField aTilde(grid);