diff --git a/lib/qcd/action/fermion/CayleyFermion5Dcache.cc b/lib/qcd/action/fermion/CayleyFermion5Dcache.cc
index 7393b760..560f5dcb 100644
--- a/lib/qcd/action/fermion/CayleyFermion5Dcache.cc
+++ b/lib/qcd/action/fermion/CayleyFermion5Dcache.cc
@@ -135,6 +135,7 @@ void CayleyFermion5D<Impl>::MooeeInv    (const FermionField &psi_i, FermionField
 {
   chi_i.Checkerboard()=psi_i.Checkerboard();
   GridBase *grid=psi_i.Grid();
+
   auto psi = psi_i.View();
   auto chi = chi_i.View();
 
@@ -183,23 +184,12 @@ void CayleyFermion5D<Impl>::MooeeInvDag (const FermionField &psi_i, FermionField
   chi_i.Checkerboard()=psi_i.Checkerboard();
   GridBase *grid=psi_i.Grid();
   int Ls=this->Ls;
+
   auto psi = psi_i.View();
   auto chi = chi_i.View();
 
   assert(psi.Checkerboard() == psi.Checkerboard());
 
-  std::vector<Coeff_t> ueec(Ls);
-  std::vector<Coeff_t> deec(Ls);
-  std::vector<Coeff_t> leec(Ls);
-  std::vector<Coeff_t> ueemc(Ls);
-  std::vector<Coeff_t> leemc(Ls);
-  for(int s=0;s<ueec.size();s++){
-    ueec[s] = conjugate(uee[s]);
-    deec[s] = conjugate(dee[s]);
-    leec[s] = conjugate(lee[s]);
-    ueemc[s]= conjugate(ueem[s]);
-    leemc[s]= conjugate(leem[s]);
-  }
   MooeeInvCalls++;
   MooeeInvTime-=usecond();
 
@@ -211,25 +201,25 @@ void CayleyFermion5D<Impl>::MooeeInvDag (const FermionField &psi_i, FermionField
     chi[ss]=psi[ss];
     for (int s=1;s<Ls;s++){
       spProj5m(tmp,chi[ss+s-1]);
-      chi[ss+s] = psi[ss+s]-ueec[s-1]*tmp;
+      chi[ss+s] = psi[ss+s]-conjugate(uee[s-1])*tmp;
     }
     // U_m^{-\dagger} 
     for (int s=0;s<Ls-1;s++){
       spProj5p(tmp,chi[ss+s]);
-      chi[ss+Ls-1] = chi[ss+Ls-1] - ueemc[s]*tmp;
+      chi[ss+Ls-1] = chi[ss+Ls-1] - conjugate(ueem[s])*tmp;
     }
 
     // L_m^{-\dagger} D^{-dagger}
     for (int s=0;s<Ls-1;s++){
       spProj5m(tmp,chi[ss+Ls-1]);
-      chi[ss+s] = (1.0/deec[s])*chi[ss+s]-(leemc[s]/deec[Ls-1])*tmp;
+      chi[ss+s] = conjugate(1.0/dee[s])*chi[ss+s]-conjugate(leem[s]/dee[Ls-1])*tmp;
     }	
-    chi[ss+Ls-1]= (1.0/deec[Ls-1])*chi[ss+Ls-1];
+    chi[ss+Ls-1]= conjugate(1.0/dee[Ls-1])*chi[ss+Ls-1];
   
     // Apply L^{-dagger}
     for (int s=Ls-2;s>=0;s--){
       spProj5p(tmp,chi[ss+s+1]);
-      chi[ss+s] = chi[ss+s] - leec[s]*tmp;
+      chi[ss+s] = chi[ss+s] - conjugate(lee[s])*tmp;
     }
   });