portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-24 12:45:56 +01:00
Code Releases Activity

Small change

Browse Source
This commit is contained in:
Peter Boyle 2021-05-12 10:57:04 -04:00
parent 69a2c8769a
commit 3575278b57
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
Grid/qcd/action/fermion/SchurFactoredFermionOperator.h
View File

@ -67,11 +67,13 @@ class SchurFactoredFermionOperator : public Impl
{ {
INHERIT_IMPL_TYPES(Impl); INHERIT_IMPL_TYPES(Impl);
public:
FermionOperator<Impl> & DirichletFermOp; FermionOperator<Impl> & DirichletFermOp;
FermionOperator<Impl> & FermOp; FermionOperator<Impl> & FermOp;
OperatorFunction<FermionField> &Solver; OperatorFunction<FermionField> &Solver;
Coordinate Block; Coordinate Block;
public:
SchurFactoredFermionOperator(FermionOperator<Impl> & _FermOp, SchurFactoredFermionOperator(FermionOperator<Impl> & _FermOp,
FermionOperator<Impl> & _DirichletFermOp, FermionOperator<Impl> & _DirichletFermOp,
OperatorFunction<FermionField> &_Solver, OperatorFunction<FermionField> &_Solver,
@ -84,6 +86,8 @@ public:
// FIXME -- could check that the DirichletFermOp block matches this. // FIXME -- could check that the DirichletFermOp block matches this.
// Pass in Dirichlet FermOp because we really need two dirac operators // Pass in Dirichlet FermOp because we really need two dirac operators
// as double stored gauge fields differ. // as double stored gauge fields differ.
assert(_FermOp.FermionGrid() = _DirichletFermOp.FermionGrid()); // May not be true in future if change communicator scheme
}; };
enum Domain { Omega=0, OmegaBar=1 }; enum Domain { Omega=0, OmegaBar=1 };