asqtad_cps_wrapper_qprop.h

Go to the documentation of this file.

// -*- C++ -*-
// $Id: asqtad_cps_wrapper_qprop.h,v 3.7 2007/12/08 11:26:53 mcneile Exp $
/*! \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 

---