asqtad_cps_wrapper_qprop.h

Go to the documentation of this file.

// -*- C++ -*-
/*! \file
 *  \brief Propagator solver for an even-odd non-preconditioned fermion operator
 *
 *  Solve for the propagator of an even-odd non-preconditioned fermion operator
 */
 
#ifndef ASQTAD_CPS_WRAPPER_QPROP_H
#define ASQTAD_CPS_WRAPPER_QPROP_H
 
 
#include "chromabase.h"
#include "handle.h"
//#include "invtype.h"
#include "syssolver.h"
#include "linearop.h"
#include "fermact.h"
#include "actions/ferm/invert/syssolver_cg_params.h"
#include "state.h"
#include "actions/ferm/invert/syssolver_cg_params.h"
#include "stagtype_fermact_s.h"
#include "actions/ferm/fermstates/asqtad_state.h"
 
namespace Chroma
{
 
  //! QPROP for wrapping CPS asqtad QPROP
  /*! \ingroup qprop
   *
   */
  class AsqtadCPSWrapperQprop : public SystemSolver<LatticeStaggeredFermion>
  {
  public:
    // Typedefs to save typing
    typedef LatticeStaggeredFermion      T;
    typedef multi1d<LatticeColorMatrix>  P;
    typedef multi1d<LatticeColorMatrix>  Q;
 
    //! Constructor
    /*!
      // Keeping the same interface as for the ordinary staggered 
      // qprop...
      //
      // But the M_ and A_ linop handles are no longer kept
      // (are ignored) -- is there a nice way around this ? 
      // Perhaps not
     */
    AsqtadCPSWrapperQprop(
                          const EvenOddStaggeredTypeFermAct<T,P,Q>& S_,
                          Handle< FermState<T,P,Q> > state,
                          const SysSolverCGParams& invParam_);
    
    //! Destructor is automatic
    ~AsqtadCPSWrapperQprop();
 
    //! Return the subset on which the operator acts
    const Subset& subset() const {return all;}
 
    //! Solver the linear system
    /*!
     * \param psi      quark propagator ( Modify )
     * \param chi      source ( Read )
     * \return number of CG iterations
     */
        SystemSolverResults_t operator() (LatticeStaggeredFermion& psi, const LatticeStaggeredFermion& chi) const;
 
 
  private:
    // Hide default constructor
    AsqtadCPSWrapperQprop() {}
 
    Real Mass;
    SysSolverCGParams invParam;
    //    Handle< FermState<T,P,Q> > state;
    Handle<AsqtadConnectStateBase> state;
 
    Handle< EvenOddLinearOperator<T,P,Q> > M;
 
  };
 
} // End namespace
 
#endif