23 if (in.
Ndim() != 5 || out.
Ndim() != 5) {
24 errorQuda(
"Staggered dslash requires 5-d fermion fields");
28 errorQuda(
"Input and output spinor precisions don't match in dslash_quda");
32 errorQuda(
"ColorSpinorFields are not single parity, in = %d, out = %d",
54 const double &k)
const
117 errorQuda(
"Preconditioned solution requires a preconditioned solve_type");
133 errorQuda(
"The optimized Kahler-Dirac operator is not built through createCoarseOp");
160 if (&
dirac !=
this) {
200 errorQuda(
"MdagM is no longer defined for DiracStaggeredPC. Use M instead.\n");
const ColorSpinorField & Odd() const
QudaSiteSubset SiteSubset() const
const ColorSpinorField & Even() const
bool newTmp(ColorSpinorField **, const ColorSpinorField &) const
void deleteTmp(ColorSpinorField **, const bool &reset) const
Dirac & operator=(const Dirac &dirac)
virtual void checkFullSpinor(const ColorSpinorField &, const ColorSpinorField &) const
check full spinors are compatible (check geometry ?)
int commDim[QUDA_MAX_DIM]
virtual void DslashXpay(ColorSpinorField &out, const ColorSpinorField &in, const QudaParity parity, const ColorSpinorField &x, const double &k) const
Xpay version of Dslash.
virtual void reconstruct(ColorSpinorField &x, const ColorSpinorField &b, const QudaSolutionType) const
DiracStaggered & operator=(const DiracStaggered &dirac)
virtual ~DiracStaggered()
DiracStaggered(const DiracParam ¶m)
virtual void MdagM(ColorSpinorField &out, const ColorSpinorField &in) const
Apply MdagM operator which may be optimized.
virtual void checkParitySpinor(const ColorSpinorField &, const ColorSpinorField &) const
Check parity spinors are usable (check geometry ?)
void createCoarseOp(GaugeField &Y, GaugeField &X, const Transfer &T, double kappa, double mass, double mu=0., double mu_factor=0.) const
Create the coarse staggered operator.
virtual void prepare(ColorSpinorField *&src, ColorSpinorField *&sol, ColorSpinorField &x, ColorSpinorField &b, const QudaSolutionType) const
virtual void Dslash(ColorSpinorField &out, const ColorSpinorField &in, const QudaParity parity) const
apply 'dslash' operator for the DiracOp. This may be e.g. AD
virtual void M(ColorSpinorField &out, const ColorSpinorField &in) const
Apply M for the dirac op. E.g. the Schur Complement operator.
virtual void reconstruct(ColorSpinorField &x, const ColorSpinorField &b, const QudaSolutionType) const
DiracStaggeredPC & operator=(const DiracStaggeredPC &dirac)
virtual ~DiracStaggeredPC()
virtual void MdagM(ColorSpinorField &out, const ColorSpinorField &in) const
Apply MdagM operator which may be optimized.
virtual void M(ColorSpinorField &out, const ColorSpinorField &in) const
Apply M for the dirac op. E.g. the Schur Complement operator.
DiracStaggeredPC(const DiracParam ¶m)
virtual void prepare(ColorSpinorField *&src, ColorSpinorField *&sol, ColorSpinorField &x, ColorSpinorField &b, const QudaSolutionType) const
QudaPrecision Precision() const
QudaTransferType getTransferType() const
quda::mgarray< double > mu_factor
@ QUDA_PARITY_SITE_SUBSET
@ QUDA_TRANSFER_OPTIMIZED_KD
enum QudaSolutionType_s QudaSolutionType
@ QUDA_MATPCDAG_MATPC_SOLUTION
enum QudaParity_s QudaParity
void ax(double a, ColorSpinorField &x)
void ApplyStaggered(ColorSpinorField &out, const ColorSpinorField &in, const GaugeField &U, double a, const ColorSpinorField &x, int parity, bool dagger, const int *comm_override, TimeProfile &profile)
Apply the staggered dslash operator to a color-spinor field.
void StaggeredCoarseOp(GaugeField &Y, GaugeField &X, const Transfer &T, const cudaGaugeField &gauge, const cudaGaugeField *XinvKD, double mass, QudaDiracType dirac, QudaMatPCType matpc)
Coarse operator construction from a fine-grid operator (Staggered)