From aac1736617df886856926eb82ee7c40ddd3ec086 Mon Sep 17 00:00:00 2001
From: Quadro <Peter Boyle>
Date: Tue, 25 May 2021 18:21:48 -0400
Subject: [PATCH] Double or single

---
 Grid/qcd/action/momentum/DirichletFilter.h | 10 +++-------
 1 file changed, 3 insertions(+), 7 deletions(-)

diff --git a/Grid/qcd/action/momentum/DirichletFilter.h b/Grid/qcd/action/momentum/DirichletFilter.h
index 1014d64e..85010ed1 100644
--- a/Grid/qcd/action/momentum/DirichletFilter.h
+++ b/Grid/qcd/action/momentum/DirichletFilter.h
@@ -37,11 +37,6 @@ NAMESPACE_BEGIN(Grid);
 template<typename 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;
   
   DirichletFilter(const Coordinate &_Block): Block(_Block){}
@@ -53,8 +48,9 @@ struct DirichletFilter: public MomentumFilterBase<MomentaField>
     ////////////////////////////////////////////////////
     // Zero strictly links crossing between domains
     ////////////////////////////////////////////////////
-    LatticeInteger coor(grid); 
-    LatticeColourMatrix zz(grid); zz = Zero();
+    LatticeInteger coor(grid);
+    typedef decltype(PeekIndex<LorentzIndex>(P,0)) MatrixType;
+    MatrixType zz(grid); zz = Zero();
     Coordinate Global=grid->GlobalDimensions();
     for(int mu=0;mu<Nd;mu++) {
       if ( (Block[mu] <= Global[mu]) && (Block[mu]>1) ) {