From 7c7ea35ffb7439d1028ab7937844849e1322d66d Mon Sep 17 00:00:00 2001 From: Guido Cossu Date: Tue, 5 Apr 2016 16:28:09 +0900 Subject: [PATCH] Putting the Traceless Antihermitian part outside the deriv in pseudofermion actions --- lib/qcd/action/pseudofermion/OneFlavourEvenOddRational.h | 2 +- .../action/pseudofermion/OneFlavourEvenOddRationalRatio.h | 2 +- lib/qcd/action/pseudofermion/OneFlavourRational.h | 2 +- lib/qcd/action/pseudofermion/OneFlavourRationalRatio.h | 2 +- lib/qcd/action/pseudofermion/TwoFlavourEvenOdd.h | 2 +- lib/qcd/action/pseudofermion/TwoFlavourEvenOddRatio.h | 5 +++-- 6 files changed, 8 insertions(+), 7 deletions(-) diff --git a/lib/qcd/action/pseudofermion/OneFlavourEvenOddRational.h b/lib/qcd/action/pseudofermion/OneFlavourEvenOddRational.h index d7d08516..1b836b4c 100644 --- a/lib/qcd/action/pseudofermion/OneFlavourEvenOddRational.h +++ b/lib/qcd/action/pseudofermion/OneFlavourEvenOddRational.h @@ -201,7 +201,7 @@ namespace Grid{ } - dSdU = Ta(dSdU); + //dSdU = Ta(dSdU); }; }; diff --git a/lib/qcd/action/pseudofermion/OneFlavourEvenOddRationalRatio.h b/lib/qcd/action/pseudofermion/OneFlavourEvenOddRationalRatio.h index 8c8bd7a7..735cb284 100644 --- a/lib/qcd/action/pseudofermion/OneFlavourEvenOddRationalRatio.h +++ b/lib/qcd/action/pseudofermion/OneFlavourEvenOddRationalRatio.h @@ -256,7 +256,7 @@ namespace Grid{ } - dSdU = Ta(dSdU); + //dSdU = Ta(dSdU); }; }; diff --git a/lib/qcd/action/pseudofermion/OneFlavourRational.h b/lib/qcd/action/pseudofermion/OneFlavourRational.h index cb3a96d7..65a99afa 100644 --- a/lib/qcd/action/pseudofermion/OneFlavourRational.h +++ b/lib/qcd/action/pseudofermion/OneFlavourRational.h @@ -186,7 +186,7 @@ namespace Grid{ } - dSdU = Ta(dSdU); + //dSdU = Ta(dSdU); }; }; diff --git a/lib/qcd/action/pseudofermion/OneFlavourRationalRatio.h b/lib/qcd/action/pseudofermion/OneFlavourRationalRatio.h index e6533d15..2207e119 100644 --- a/lib/qcd/action/pseudofermion/OneFlavourRationalRatio.h +++ b/lib/qcd/action/pseudofermion/OneFlavourRationalRatio.h @@ -242,7 +242,7 @@ namespace Grid{ } - dSdU = Ta(dSdU); + //dSdU = Ta(dSdU); }; }; diff --git a/lib/qcd/action/pseudofermion/TwoFlavourEvenOdd.h b/lib/qcd/action/pseudofermion/TwoFlavourEvenOdd.h index e3ac2b2e..601bf9ec 100644 --- a/lib/qcd/action/pseudofermion/TwoFlavourEvenOdd.h +++ b/lib/qcd/action/pseudofermion/TwoFlavourEvenOdd.h @@ -173,7 +173,7 @@ namespace Grid{ FermOp.MeeDeriv(tmp , X, Y,DaggerYes); dSdU=dSdU+tmp; */ - dSdU = Ta(dSdU); + //dSdU = Ta(dSdU); }; diff --git a/lib/qcd/action/pseudofermion/TwoFlavourEvenOddRatio.h b/lib/qcd/action/pseudofermion/TwoFlavourEvenOddRatio.h index b1d73a7c..b2ae4d21 100644 --- a/lib/qcd/action/pseudofermion/TwoFlavourEvenOddRatio.h +++ b/lib/qcd/action/pseudofermion/TwoFlavourEvenOddRatio.h @@ -188,8 +188,9 @@ namespace Grid{ assert(NumOp.ConstEE() == 1); assert(DenOp.ConstEE() == 1); - dSdU = -Ta(dSdU); - + //dSdU = -Ta(dSdU); + dSDu = -dSdU; + }; }; }