Fixed bug in hisq force unitarization

lattice · May 7, 2016 · aa9b0be · aa9b0be
1 parent a242e6d
commit aa9b0be
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/lib/unitarize_force_quda.cu b/lib/unitarize_force_quda.cu
@@ -499,7 +499,7 @@ namespace { // anonymous
 
 	getUnitarizeForceSite<double>(result, v, oprod, arg);
 
-	arg.force.save((Float*)(oprod.data), idx, dir, parity);
+	arg.force.save((Float*)(result.data), idx, dir, parity);
       } // 4*4528 flops per site
       return;
     } // getUnitarizeForceField
@@ -517,7 +517,7 @@ namespace { // anonymous
 
 	    getUnitarizeForceSite<double>(result, v, oprod, arg);
 
-	    arg.force.save((Float*)(oprod.data), i, dir, parity);
+	    arg.force.save((Float*)(result.data), i, dir, parity);
 	  }
 	}
       }