seoprec_constdet_one_flavor_rat_monomial_w.cc

Go to the documentation of this file.

/*! @file
 * @brief One-flavor collection of even-odd preconditioned 4D ferm monomials
 */
 
#include "update/molecdyn/monomial/seoprec_constdet_one_flavor_rat_monomial_w.h"
#include "update/molecdyn/monomial/monomial_factory.h"
 
#include "update/molecdyn/monomial/rat_approx_factory.h"
#include "update/molecdyn/monomial/rat_approx_aggregate.h"
 
#include "actions/ferm/fermacts/fermacts_aggregate_w.h"
#include "actions/ferm/fermacts/fermact_factory_w.h"
 
namespace Chroma 
{ 
 
namespace SymEvenOddPrecConstDetOneFlavorWilsonTypeFermRatMonomialEnv
{
namespace
{
//! Callback
Monomial< multi1d<LatticeColorMatrix>,
multi1d<LatticeColorMatrix> >* createMonomial(XMLReader& xml, const std::string& path)
{
        return new SymEvenOddPrecConstDetOneFlavorWilsonTypeFermRatMonomial(
                        OneFlavorWilsonTypeFermRatMonomialParams(xml, path));
}
 
//! Local registration flag
bool registered = false;
}
 
const std::string name("ONE_FLAVOR_SEOPREC_CONSTDET_FERM_RAT_MONOMIAL");
 
//! Register all the factories
bool registerAll()
{
        bool success = true;
        if (! registered)
        {
                success &= WilsonTypeFermActs4DEnv::registerAll();
                success &= RationalApproxAggregateEnv::registerAll();
                success &= TheMonomialFactory::Instance().registerObject(name, createMonomial);
                registered = true;
        }
        return success;
}
} //end namespace EvenOddPrec OneFlavorWilsonFermRatMonomialEnv
 
 
 
// Constructor
SymEvenOddPrecConstDetOneFlavorWilsonTypeFermRatMonomial::SymEvenOddPrecConstDetOneFlavorWilsonTypeFermRatMonomial(
                const OneFlavorWilsonTypeFermRatMonomialParams& param)
{
        START_CODE();
 
        QDPIO::cout << "Constructor: " << __func__ << std::endl;
 
        actionInvParam = param.numer.action.invParam;
        forceInvParam  = param.numer.force.invParam;
        num_pf         = param.num_pf;
 
        //*********************************************************************
        // Fermion action
        {
                std::istringstream is(param.numer.fermact.xml);
                XMLReader fermact_reader(is);
                QDPIO::cout << "Construct fermion action= " << param.numer.fermact.id << std::endl;
 
                const WilsonTypeFermAct<T,P,Q>* tmp_act = TheWilsonTypeFermActFactory::Instance().createObject(
                                param.numer.fermact.id,
                                fermact_reader,
                                param.numer.fermact.path);
 
                const SymEvenOddPrecWilsonTypeFermAct<T,P,Q>* downcast =
                                dynamic_cast<const SymEvenOddPrecWilsonTypeFermAct<T,P,Q>*>(tmp_act);
 
                // Check success of the downcast
                if( downcast == 0x0 ) {
                        QDPIO::cerr << "Unable to downcast FermAct to SymEvenOddPrecWilsonTypeFermAct in SymEvenOddPrecOneFlavorWilsonTypeFermRatMonomial()" << std::endl;
                        QDP_abort(1);
                }
 
                fermact = downcast;
        }
 
 
        //*********************************************************************
        // Action rational approx
        {
                std::istringstream is(param.numer.action.ratApprox.xml);
                XMLReader approx_reader(is);
                QDPIO::cout << "Construct action rational approx= " << param.numer.action.ratApprox.id << std::endl;
 
                Handle<RationalApprox> approx(TheRationalApproxFactory::Instance().createObject(
                                param.numer.action.ratApprox.id,
                                approx_reader,
                                param.numer.action.ratApprox.path));
 
                (*approx)(spfe, sipfe);
        }
 
        //*********************************************************************
        // Force rational approx
        {
                std::istringstream is(param.numer.force.ratApprox.xml);
                XMLReader approx_reader(is);
                QDPIO::cout << "Construct force rational approx= " << param.numer.force.ratApprox.id << std::endl;
 
                Handle<RationalApprox> approx(TheRationalApproxFactory::Instance().createObject(
                                param.numer.force.ratApprox.id,
                                approx_reader,
                                param.numer.force.ratApprox.path));
 
                RemezCoeff_t  fipfe;  // discard
                (*approx)(fpfe, fipfe);
        }
        //*********************************************************************
 
        QDPIO::cout << "Finished constructing: " << __func__ << std::endl;
 
        END_CODE();
}
 
} //end namespace Chroma