Merge pull request #481 from lattice/hotfix/twisted-clover

Fixes a bug in asymmetric preconditioning
lattice · Jun 13, 2016 · bc696cf · bc696cf
2 parents ebf6f1e + 476063b
commit bc696cf
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/lib/dirac_twisted_clover.cpp b/lib/dirac_twisted_clover.cpp
@@ -333,15 +333,15 @@ namespace quda {
 	  twistedCloverDslashCuda(&static_cast<cudaColorSpinorField&>(out), *gauge, cs, cI,
 				  static_cast<cudaColorSpinorField*>(tmp1), QUDA_EVEN_PARITY, dagger,
 				  &static_cast<const cudaColorSpinorField&>(in),
-				  QUDA_DEG_DSLASH_CLOVER_TWIST_XPAY, a, kappa, 0.0, 0.0, commDim, profile);
+				  QUDA_DEG_DSLASH_CLOVER_TWIST_XPAY, a, kappa2, 0.0, 0.0, commDim, profile);
 
           flops += (1320ll+96ll)*in.Volume();	 
         } else if (matpcType == QUDA_MATPC_ODD_ODD_ASYMMETRIC) {
 	  Dslash(*tmp1, in, QUDA_EVEN_PARITY);
 	  twistedCloverDslashCuda(&static_cast<cudaColorSpinorField&>(out), *gauge, cs, cI,
 				  static_cast<const cudaColorSpinorField*>(tmp1), QUDA_ODD_PARITY, dagger,
 				  &static_cast<const cudaColorSpinorField&>(in),
-				  QUDA_DEG_DSLASH_CLOVER_TWIST_XPAY, a, kappa, 0.0, 0.0, commDim, profile);
+				  QUDA_DEG_DSLASH_CLOVER_TWIST_XPAY, a, kappa2, 0.0, 0.0, commDim, profile);
           flops += (1320ll+96ll)*in.Volume();
         }else { // symmetric preconditioning
           errorQuda("Invalid matpcType");