fermbcs_reader_s.cc

Go to the documentation of this file.

/*! \file
 *  \brief Fermionic BC reader
 */
 
#include "actions/ferm/fermbcs/fermbcs_reader_s.h"
#include "actions/ferm/fermbcs/fermbc_factory_s.h"
#include "actions/ferm/fermbcs/simple_fermbc_s.h"
 
namespace Chroma
{
 
  //! Name and registration
  namespace StaggeredTypeFermBCEnv
  {
    // Helper function for the FermionAction readers
    /*
     * This structure should not be replicated. This routine helps maintain
     * backwards compatibility with the FermionAction readers by looking for
     * either the "boundary" tag or the FermionBC group
     */
    Handle< FermBC<LatticeStaggeredFermion,
                   multi1d<LatticeColorMatrix>,
                   multi1d<LatticeColorMatrix> > > reader(XMLReader& xml_in, 
                                                          const std::string& path)
    {
      XMLReader top(xml_in, path);
 
      std::string fermbc;
      std::string fermbc_path;
      if (top.count("FermionBC") != 0)
      {
        fermbc_path = "FermionBC";
        read(top, fermbc_path + "/FermBC", fermbc);
      }
      else if (top.count("boundary") != 0)
      {
        fermbc_path = ".";
        fermbc = StaggeredTypeSimpleFermBCEnv::name;
      }
      else
      {
        QDPIO::cerr << "Error: neither FermionBC group nor boundary found" << std::endl;
        QDP_abort(1);
      }
 
      Handle< FermBC<LatticeStaggeredFermion,
                     multi1d<LatticeColorMatrix>,
                     multi1d<LatticeColorMatrix> > >
        fbc(TheStaggeredTypeFermBCFactory::Instance().createObject(fermbc,
                                                                   top,
                                                                   fermbc_path));
 
      return fbc;
    }
  }
 
}