portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

Bug fix. Guess not initialised in refresh step; didn't hit before due to luck in not having a vector

Browse Source 
created with NAN data.
This commit is contained in:
paboyle 2015-12-21 15:34:35 +00:00
parent 2f41691c11
commit 36e6f9ac7b
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
lib/qcd/action/pseudofermion/TwoFlavourRatio.h
View File

@ -58,6 +58,7 @@ namespace Grid{
MdagMLinearOperator<FermionOperator<Impl> ,FermionField> MdagMOp(NumOp); MdagMLinearOperator<FermionOperator<Impl> ,FermionField> MdagMOp(NumOp);
DenOp.Mdag(eta,Phi); // Mdag eta DenOp.Mdag(eta,Phi); // Mdag eta
tmp = zero;
ActionSolver(MdagMOp,Phi,tmp); // (VdagV)^-1 Mdag eta = V^-1 Vdag^-1 Mdag eta ActionSolver(MdagMOp,Phi,tmp); // (VdagV)^-1 Mdag eta = V^-1 Vdag^-1 Mdag eta
NumOp.M(tmp,Phi); // Vdag^-1 Mdag eta NumOp.M(tmp,Phi); // Vdag^-1 Mdag eta
@ -78,8 +79,8 @@ namespace Grid{
MdagMLinearOperator<FermionOperator<Impl> ,FermionField> MdagMOp(DenOp); MdagMLinearOperator<FermionOperator<Impl> ,FermionField> MdagMOp(DenOp);
X=zero;
NumOp.Mdag(Phi,Y); // Y= Vdag phi NumOp.Mdag(Phi,Y); // Y= Vdag phi
X=zero;
ActionSolver(MdagMOp,Y,X); // X= (MdagM)^-1 Vdag phi ActionSolver(MdagMOp,Y,X); // X= (MdagM)^-1 Vdag phi
DenOp.M(X,Y); // Y= Mdag^-1 Vdag phi DenOp.M(X,Y); // Y= Mdag^-1 Vdag phi
@ -105,12 +106,12 @@ namespace Grid{
GaugeField force(NumOp.GaugeGrid()); GaugeField force(NumOp.GaugeGrid());
X=zero;
//Y=Vdag phi //Y=Vdag phi
//X = (Mdag M)^-1 V^dag phi //X = (Mdag M)^-1 V^dag phi
//Y = (Mdag)^-1 V^dag phi //Y = (Mdag)^-1 V^dag phi
NumOp.Mdag(Phi,Y); // Y= Vdag phi NumOp.Mdag(Phi,Y); // Y= Vdag phi
X=zero;
DerivativeSolver(MdagMOp,Y,X); // X= (MdagM)^-1 Vdag phi DerivativeSolver(MdagMOp,Y,X); // X= (MdagM)^-1 Vdag phi
DenOp.M(X,Y); // Y= Mdag^-1 Vdag phi DenOp.M(X,Y); // Y= Mdag^-1 Vdag phi