Match x and p in computeCloverForce and computeCloverSigmaOprod.

lattice · Jan 4, 2024 · 66a5bb5 · 66a5bb5
1 parent fd50676
commit 66a5bb5
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 7 deletions.
diff --git a/lib/clover_force.cpp b/lib/clover_force.cpp
@@ -67,7 +67,7 @@ namespace quda
     GaugeField oprod(param);
 
     // derivative of the wilson operator it correspond to deriv_Sb(OE,...) plus  deriv_Sb(EO,...) in tmLQCD
-    computeCloverForce(force, gauge, x, p, coeff);
+    computeCloverForce(force, gauge, inv_param.dagger == QUDA_DAG_YES ? p : x, inv_param.dagger == QUDA_DAG_YES ? x : p, coeff);
     // derivative of the determinant of the sw term, second term of (A12) in hep-lat/0112051,  sw_deriv(EE, mnl->mu) in tmLQCD
     if (!detratio) computeCloverSigmaTrace(oprod, clover, sigma_coeff, other_parity);
 

diff --git a/lib/clover_outer_product.cu b/lib/clover_outer_product.cu
@@ -142,10 +142,10 @@ namespace quda {
 
       for (auto i = 0u; i < x.size(); i++) {
         for (int parity = 0; parity < 2; parity++) {
-          const ColorSpinorField &inA = (parity & 1) ? x[i].Odd() : x[i].Even();
-          const ColorSpinorField &inB = (parity & 1) ? p[i].Even() : p[i].Odd();
-          const ColorSpinorField &inC = (parity & 1) ? p[i].Odd() : p[i].Even();
-          const ColorSpinorField &inD = (parity & 1) ? x[i].Even() : x[i].Odd();
+          const ColorSpinorField &inA = (parity & 1) ? p[i].Odd() : p[i].Even();
+          const ColorSpinorField &inB = (parity & 1) ? x[i].Even() : x[i].Odd();
+          const ColorSpinorField &inC = (parity & 1) ? x[i].Odd() : x[i].Even();
+          const ColorSpinorField &inD = (parity & 1) ? p[i].Even() : p[i].Odd();
 
           getProfile().TPSTART(QUDA_PROFILE_COMMS);
           exchangeGhost(inB, parity, dag);

diff --git a/lib/interface_quda.cpp b/lib/interface_quda.cpp
@@ -912,8 +912,6 @@ void loadCloverQuda(void *h_clover, void *h_clovinv, QudaInvertParam *inv_param)
   popVerbosity();
 }
 
-void freeSloppyCloverQuda();
-
 void loadSloppyCloverQuda(const QudaPrecision *prec)
 {
   freeSloppyCloverQuda();