weak_gauge_init.cc

Go to the documentation of this file.

/*! \file
 *  \brief Create a weak config
 */
 
#include "util/gauge/gauge_init_factory.h"
#include "util/gauge/gauge_init_aggregate.h"
 
#include "util/gauge/weak_gauge_init.h"
#include "util/gauge/weak_field.h"
 
namespace Chroma
{
 
  //! Hooks to register the class
  namespace WeakGaugeInitEnv
  {
    //! Callback function
    GaugeInit* createSource(XMLReader& xml_in,
                            const std::string& path)
    {
      return new GaugeIniter(Params(xml_in, path));
    }
 
    //! Name to be used
    const std::string name = "WEAK_FIELD";
 
    //! Local registration flag
    static bool registered = false;
 
    //! Register all the factories
    bool registerAll() 
    {
      bool success = true; 
      if (! registered)
      {
        success &= Chroma::TheGaugeInitFactory::Instance().registerObject(name, createSource);
        registered = true;
      }
      return success;
    }
 
 
    // Parameters for running code
    Params::Params(XMLReader& xml, const std::string& path)
    {
      XMLReader paramtop(xml, path);
    }
 
 
    //! Parameters for running code
    void Params::writeXML(XMLWriter& xml, const std::string& path) const
    {
      push(xml, path);
    
      int version = 1;
      write(xml, "cfg_type", WeakGaugeInitEnv::name);
 
      pop(xml);
    }
 
 
    // Initialize the gauge field
    void
    GaugeIniter::operator()(XMLReader& gauge_file_xml,
                            XMLReader& gauge_xml,
                            multi1d<LatticeColorMatrix>& u) const
    {
      u.resize(Nd);
      QDPIO::cout << "Starting up a weak field config" << std::endl;
      weakField(u);
 
      XMLBufferWriter file_xml, record_xml;
      push(file_xml, "gauge");
      write(file_xml, "id", int(0));
      pop(file_xml);
      push(record_xml, "weak_field");
      pop(record_xml);
 
      gauge_file_xml.open(file_xml);
      gauge_xml.open(record_xml);
    }
  }
}