asqtad_qprop.cc

Go to the documentation of this file.

/*! \file
 *  \brief Propagator solver for an even-odd non-preconditioned fermion operator
 *
 *  Solve for the propagator of an even-odd non-preconditioned fermion operator
 */
 
#include "chromabase.h"
 
#include "fermact.h"
#include "actions/ferm/invert/invcg1.h"
 
#include "actions/ferm/fermacts/asqtad_fermact_s.h"
#include "actions/ferm/qprop/asqtad_qprop.h"
#include "actions/ferm/invert/syssolver_cg_params.h"
 
namespace Chroma
{
 
  typedef LatticeStaggeredFermion LF;
  typedef multi1d<LatticeColorMatrix> LCM;
 
  SystemSolver<LF>* 
  AsqtadFermAct::qprop(Handle< FermState<LF,LCM,LCM> > state,
                       const GroupXML_t& invParam) const
  {
    const EvenOddStaggeredTypeFermAct<LF,LCM,LCM>& S_cast = 
      dynamic_cast< const EvenOddStaggeredTypeFermAct<LF,LCM,LCM>& >(*this);
 
    std::istringstream  is(invParam.xml);
    XMLReader  paramtop(is);
    SysSolverCGParams params(paramtop, invParam.path);
 
    return new AsqtadQprop(S_cast,
                           state,
                           params);
  }
  
 
}