From c8a824425bbb5b543d41340d91a78ecc9273011e Mon Sep 17 00:00:00 2001
From: Fabian Joswig <fabian.joswig@ed.ac.uk>
Date: Tue, 5 Apr 2022 10:58:22 +0100
Subject: [PATCH] Error message added if another conserved current than vector
 is requested for Wilson type fermions.

---
 .../WilsonFermionImplementation.h             | 20 +++++++++++++++----
 1 file changed, 16 insertions(+), 4 deletions(-)

diff --git a/Grid/qcd/action/fermion/implementation/WilsonFermionImplementation.h b/Grid/qcd/action/fermion/implementation/WilsonFermionImplementation.h
index 40e774b3..69a7fb16 100644
--- a/Grid/qcd/action/fermion/implementation/WilsonFermionImplementation.h
+++ b/Grid/qcd/action/fermion/implementation/WilsonFermionImplementation.h
@@ -599,16 +599,22 @@ void WilsonFermion<Impl>::ContractConservedCurrent(PropagatorField &q_in_1,
                                                    Current curr_type,
                                                    unsigned int mu)
 {
+  if(curr_type != Current::Vector)
+  {
+    std::cout << GridLogError << "Only the conserved vector current is implemented so far." << std::endl;
+    exit(1);
+  }
+
   Gamma g5(Gamma::Algebra::Gamma5);
   conformable(_grid, q_in_1.Grid());
   conformable(_grid, q_in_2.Grid());
   conformable(_grid, q_out.Grid());
   auto UGrid= this->GaugeGrid();
 
-  PropagatorField tmp_shifted(UGrid); 
-  PropagatorField g5Lg5(UGrid); 
-  PropagatorField R(UGrid); 
-  PropagatorField gmuR(UGrid); 
+  PropagatorField tmp_shifted(UGrid);
+  PropagatorField g5Lg5(UGrid);
+  PropagatorField R(UGrid);
+  PropagatorField gmuR(UGrid);
 
     Gamma::Algebra Gmu [] = {
     Gamma::Algebra::GammaX,
@@ -647,6 +653,12 @@ void WilsonFermion<Impl>::SeqConservedCurrent(PropagatorField &q_in,
                                               unsigned int tmax,
 					      ComplexField &lattice_cmplx)
 {
+  if(curr_type != Current::Vector)
+  {
+    std::cout << GridLogError << "Only the conserved vector current is implemented so far." << std::endl;
+    exit(1);
+  }
+
   conformable(_grid, q_in.Grid());
   conformable(_grid, q_out.Grid());
   assert(0);