key_hadron_2pt_corr.cc

Go to the documentation of this file.

/*! \file
 * \brief Hadron 2pt correlators
 */
 
#include "util/ferm/key_hadron_2pt_corr.h"
 
namespace Chroma
{
#if 0
  namespace
  {
    //! Error output
    std::ostream& operator<<(std::ostream& os, const Array<int>& d)
    {
      if (d.size() > 0)
      {
        os << d[0];
 
        for(int i=1; i < d.size(); ++i)
          os << " " << d[i];
      }
 
      return os;
    }
 
    //! Error output
    std::ostream& operator<<(std::ostream& os, const Array<ComplexD>& d)
    {
      if (d.size() > 0)
      {
        for(int i=0; i < d.size(); ++i)
        {
          if (i > 0)
            os << "  ";
 
          os << "(" << real(d[i]) << "," << imag(d[i]) << ")";
        }
      }
 
      return os;
    }
 
    //! Error output
    template<typename T>
    std::ostream& operator<<(std::ostream& os, const std::vector<T>& d)
    {
      if (d.size() > 0)
      {
        os << d[0];
 
        for(int i=1; i < d.size(); ++i)
          os << " " << d[i];
      }
 
      return os;
    }
  }
 
 
  //----------------------------------------------------------------------------
  //! Used for error output
  std::ostream& operator<<(std::ostream& os, const KeyHadron2PtCorr_t& d)
  {
    os << "KeyHadron2PtCorr_t:"
       << " src_name= " << d.src_name
       << " src_smear= " << d.src_smear
       << " src_lorentz= " << d.src_lorentz
       << " src_spin= " << d.src_spin
       << " snk_name= " << d.snk_name
       << " snk_smear= " << d.snk_smear
       << " snk_lorentz= " << d.snk_lorentz
       << " snk_spin= " << d.snk_spin
       << " mass= " << d.mass
       << " mom= " << d.mom
       << " ensemble= " << d.ensemble
       << " num_vecs= " << d.num_vecs
       << std::endl;
 
    return os;
  }
#endif
 
 
  //----------------------------------------------------------------------------
  // Read a key
  void read(XMLReader& xml, const std::string& path, KeyHadron2PtCorr_t& param)
  {
    XMLReader paramtop(xml, path);
 
    read(paramtop, "num_vecs", param.num_vecs);
    read(paramtop, "src_name", param.src_name);
    read(paramtop, "src_smear", param.src_smear);
    read(paramtop, "src_lorentz", param.src_lorentz);
    read(paramtop, "src_spin", param.src_spin);
    read(paramtop, "snk_name", param.snk_name);
    read(paramtop, "snk_smear", param.snk_smear);
    read(paramtop, "snk_lorentz", param.snk_lorentz);
    read(paramtop, "snk_spin", param.snk_spin);
    read(paramtop, "mom", param.mom);
    read(paramtop, "mass", param.mass);
    read(paramtop, "ensemble", param.ensemble);
  }
 
  //! KeyHadron2PtCorr writer
  void write(XMLWriter& xml, const std::string& path, const KeyHadron2PtCorr_t& param)
  {
    push(xml, path);
 
    write(xml, "num_vecs", param.num_vecs);
    write(xml, "src_name", param.src_name);
    write(xml, "src_smear", param.src_smear);
    write(xml, "src_lorentz", param.src_lorentz);
    write(xml, "src_spin", param.src_spin);
    write(xml, "snk_name", param.snk_name);
    write(xml, "snk_smear", param.snk_smear);
    write(xml, "snk_lorentz", param.snk_lorentz);
    write(xml, "snk_spin", param.snk_spin);
    write(xml, "mom", param.mom);
    write(xml, "mass", param.mass);
    write(xml, "ensemble", param.ensemble);
 
    pop(xml);
  }
 
 
  //----------------------------------------------------------------------------
  //! KeyHadron2PtCorr reader
  void read(BinaryReader& bin, KeyHadron2PtCorr_t& param)
  {
    read(bin, param.src_name, 128);
    read(bin, param.src_smear, 128);
    read(bin, param.src_lorentz);
    read(bin, param.src_spin);
    read(bin, param.snk_name, 128);
    read(bin, param.snk_smear, 128);
    read(bin, param.snk_lorentz);
    read(bin, param.snk_spin);
    read(bin, param.mom);
    read(bin, param.num_vecs);
    read(bin, param.mass, 128);
    read(bin, param.ensemble, 1024);
  }
 
  //! Hadron2PtCorr write
  void write(BinaryWriter& bin, const KeyHadron2PtCorr_t& param)
  {
    write(bin, param.src_name);
    write(bin, param.src_smear);
    write(bin, param.src_lorentz);
    write(bin, param.src_spin);
    write(bin, param.snk_name);
    write(bin, param.snk_smear);
    write(bin, param.snk_lorentz);
    write(bin, param.snk_spin);
    write(bin, param.mom);
    write(bin, param.num_vecs);
    write(bin, param.mass);
    write(bin, param.ensemble);
  }
 
} // namespace Chroma