multi_syssolver_mdagm_cg_chrono_clover.cc

Go to the documentation of this file.

/*! \file
 *  \brief Solve a MdagM*psi=chi linear system by CG2
 */
 
#include "actions/ferm/invert/multi_syssolver_mdagm_factory.h"
#include "actions/ferm/invert/multi_syssolver_mdagm_aggregate.h"
 
#include "actions/ferm/invert/multi_syssolver_mdagm_cg_chrono_clover.h"
 
namespace Chroma
{
 
 
  //! CG2 system solver namespace
  namespace MdagMMultiSysSolverCGChronoCloverEnv
  {
    //! Callback function
    MdagMMultiSystemSolver<LatticeFermion>* createFerm(XMLReader& xml_in,
                                                       const std::string& path,
                                                       Handle< FermState< LatticeFermion, multi1d<LatticeColorMatrix>, multi1d<LatticeColorMatrix> > > state, 
                                                       Handle< LinearOperator<LatticeFermion> > A)
    {
      return new MdagMMultiSysSolverCGChronoClover(A, state,MultiSysSolverCGChronoCloverParams(xml_in, path));
    }
 
    //! Name to be used
    const std::string name("CG_CHRONO_CLOVER_INVERTER");
 
    //! Local registration flag
    static bool registered = false;
 
    //! Register all the factories
    bool registerAll() 
    {
      bool success = true; 
      if (! registered)
      {
        success &= Chroma::TheMdagMFermMultiSystemSolverFactory::Instance().registerObject(name, createFerm);
        registered = true;
      }
      return success;
    }
  }
 
 
  MultiSysSolverCGChronoCloverParams::MultiSysSolverCGChronoCloverParams(XMLReader& xml, 
                                                       const std::string& path)
  {
    XMLReader paramtop(xml, path);
    try {
      read(paramtop, "CloverParams", clovParams);
      read(paramtop, "MaxIter", MaxIter);
      read(paramtop, "MaxChrono", MaxChrono);
      read(paramtop, "Delta", Delta);
      read(paramtop, "CutoffRsd", CutoffRsd);
      read(paramtop, "RsdTarget", RsdTarget);
    }
    catch(const std::string e ) {
      QDPIO::cout << "Caught: " << e << std::endl;
      throw;
    }
  }
 
  void read(XMLReader& xml, const std::string& path, 
            MultiSysSolverCGChronoCloverParams& p)
  {
    MultiSysSolverCGChronoCloverParams tmp(xml, path);
    p = tmp;
  }
 
  void write(XMLWriter& xml, const std::string& path, 
             const MultiSysSolverCGChronoCloverParams& p) {
    push(xml, path);
    write(xml, "CloverParams", p.clovParams);
    write(xml, "MaxIter", p.MaxIter);
    write(xml, "MaxChrono", p.MaxChrono);
    write(xml, "Delta", p.Delta);
    write(xml, "CutoffRsd", p.CutoffRsd);
    write(xml, "RsdTarget", p.RsdTarget);
 
    
    pop(xml);
 
  }
 
}