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

Double or single

Browse Source
This commit is contained in:
Quadro 2021-05-25 18:21:48 -04:00
parent 60f814733d
commit aac1736617
1 changed files with 3 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Grid/qcd/action/momentum/DirichletFilter.h
View File

@ -37,11 +37,6 @@ NAMESPACE_BEGIN(Grid);
template<typename MomentaField> template<typename MomentaField>
struct DirichletFilter: public MomentumFilterBase<MomentaField> struct DirichletFilter: public MomentumFilterBase<MomentaField>
{ {
typedef typename MomentaField::vector_type vector_type; //SIMD-vectorized complex type
typedef typename MomentaField::scalar_type scalar_type; //scalar complex type
typedef iScalar<iScalar<iScalar<vector_type> > > ScalarType; //complex phase for each site
Coordinate Block; Coordinate Block;
DirichletFilter(const Coordinate &_Block): Block(_Block){} DirichletFilter(const Coordinate &_Block): Block(_Block){}
@ -54,7 +49,8 @@ struct DirichletFilter: public MomentumFilterBase<MomentaField>
// Zero strictly links crossing between domains // Zero strictly links crossing between domains
//////////////////////////////////////////////////// ////////////////////////////////////////////////////
LatticeInteger coor(grid); LatticeInteger coor(grid);
LatticeColourMatrix zz(grid); zz = Zero(); typedef decltype(PeekIndex<LorentzIndex>(P,0)) MatrixType;
MatrixType zz(grid); zz = Zero();
Coordinate Global=grid->GlobalDimensions(); Coordinate Global=grid->GlobalDimensions();
for(int mu=0;mu<Nd;mu++) { for(int mu=0;mu<Nd;mu++) {
if ( (Block[mu] <= Global[mu]) && (Block[mu]>1) ) { if ( (Block[mu] <= Global[mu]) && (Block[mu]>1) ) {