transf.cc

Go to the documentation of this file.

/*! \file
 *  \brief Insertion/Extraction utilities for vectors/fermions/propagators
 */
 
#include "chromabase.h"
#include "util/ferm/transf.h"
 
namespace Chroma
{
  //! Convert (insert) a LatticeColorVector into a LatticeFermion
  /*!
   * \ingroup ferm
   *
   * \param a      Source color std::vector
   * \param b      Destination fermion
   * \param spin_index   Spin index
   */
  void CvToFerm(const LatticeColorVectorF& a, LatticeFermionF& b, 
                int spin_index)
  {
    pokeSpin(b, a, spin_index);
  }
 
  //! Convert (insert) a LatticeColorVector into a LatticeFermion
  /*!
   * \ingroup ferm
   *
   * \param a      Source color std::vector
   * \param b      Destination fermion
   * \param spin_index   Spin index
   */
  void CvToFerm(const LatticeColorVectorD& a, LatticeFermionD& b, 
                int spin_index)
  {
    pokeSpin(b, a, spin_index);
  }
 
 
  //! Convert (insert) a LatticeColorVector into a LatticeStaggeredFermion
  /*!
   * \ingroup ferm
   *
   * \param a      Source color std::vector
   * \param b      Destination fermion
   */
  void CvToFerm(const LatticeColorVectorF& a, LatticeStaggeredFermionF& b)
  {
    pokeSpin(b, a, 0);
  }
 
  //! Convert (insert) a LatticeColorVector into a LatticeStaggeredFermion
  /*!
   * \ingroup ferm
   *
   * \param a      Source color std::vector
   * \param b      Destination fermion
   */
  void CvToFerm(const LatticeColorVectorD& a, LatticeStaggeredFermionD& b)
  {
    pokeSpin(b, a, 0);
  }
 
  //! Convert a LatticeStaggeredFermion into a LatticeColorVector (extract)
  /*!
   * \ingroup ferm
   *
   * \param a      Source Fermion
   * \param b      Destination ColorVector
   */
  void FermToCv(const LatticeStaggeredFermionD& a, LatticeColorVectorD& b)
    {
      b = peekSpin(a, 0);
    }
 
  //! Convert  a LatticeStaggeredFermion into a LatticeColorVector (extract)
  /*!
   * \ingroup ferm
   *
   * \param a      Source Fermion
   * \param b      Destination ColorVector
   */
  void FermToCv(const LatticeStaggeredFermionF& a, LatticeColorVectorF& b)
    {
          b = peekSpin(a,0);
    }
 
 
 
  //! Insert a LatticeFermion into a LatticePropagator
  /*!
   * \ingroup ferm
   *
   * \param a      Source fermion
   * \param b      Destination propagator
   * \param color_index  Color index
   * \param spin_index   Spin index
   */
  void FermToProp(const LatticeFermionF& a, LatticePropagatorF& b, 
                  int color_index, int spin_index)
  {
    for(int j = 0; j < Ns; ++j)
    {
      LatticeColorMatrixF bb = peekSpin(b, j, spin_index);
      LatticeColorVectorF aa = peekSpin(a, j);
 
      for(int i = 0; i < Nc; ++i)
        pokeColor(bb, peekColor(aa, i), i, color_index);
 
      pokeSpin(b, bb, j, spin_index);
    }
  }
 
  //! Insert a LatticeFermion into a LatticePropagator
  /*!
   * \ingroup ferm
   *
   * \param a      Source fermion
   * \param b      Destination propagator
   * \param color_index  Color index
   * \param spin_index   Spin index
   */
  void FermToProp(const LatticeFermionD& a, LatticePropagatorD& b, 
                  int color_index, int spin_index)
  {
    for(int j = 0; j < Ns; ++j)
    {
      LatticeColorMatrixD bb = peekSpin(b, j, spin_index);
      LatticeColorVectorD aa = peekSpin(a, j);
 
      for(int i = 0; i < Nc; ++i)
        pokeColor(bb, peekColor(aa, i), i, color_index);
 
      pokeSpin(b, bb, j, spin_index);
    }
  }
 
#ifndef QDP_IS_QDPJIT_NO_NVPTX
  //-------------------------------------------------------------------------------
  //! Insert a LatticeFermion into a LatticeColorVectorSpinMatrix
  /*!
   * \ingroup ferm
   *
   * \param a      Source fermion
   * \param b      Destination propagator
   * \param color_index  Color index
   * \param spin_index   Spin index
   */
  void FermToProp(const LatticeFermionF& a, LatticeColorVectorSpinMatrixF& b, 
                  int spin_index)
  {
    for(int j = 0; j < Ns; ++j)
      pokeSpin(b, peekSpin(a, j), j, spin_index);
  }
 
  
  //-------------------------------------------------------------------------------
  //! Insert a LatticeFermion into a LatticeColorVectorSpinMatrix
  /*!
   * \ingroup ferm
   *
   * \param a      Source fermion
   * \param b      Destination propagator
   * \param color_index  Color index
   * \param spin_index   Spin index
   */
  void FermToProp(const LatticeFermionD& a, LatticeColorVectorSpinMatrixD& b, 
                  int spin_index)
  {
    for(int j = 0; j < Ns; ++j)
      pokeSpin(b, peekSpin(a, j), j, spin_index);
  }
#endif
  
  //! Insert a LatticeStaggeredFermion into a LatticeStaggeredPropagator
  /*!
   * \ingroup ferm
   *
   * \param a      Source fermion
   * \param b      Destination propagator
   * \param color_index  Color index
   */
  void FermToProp(const LatticeStaggeredFermionF& a, LatticeStaggeredPropagatorF& b, 
                  int color_index)
  {
    LatticeColorMatrixF bb = peekSpin(b, 0, 0);
    LatticeColorVectorF aa = peekSpin(a, 0);
 
    for(int i = 0; i < Nc; ++i)
      pokeColor(bb, peekColor(aa, i), i, color_index);
 
    pokeSpin(b, bb, 0, 0);
  }
 
  //! Insert a LatticeStaggeredFermion into a LatticeStaggeredPropagator
  /*!
   * \ingroup ferm
   *
   * \param a      Source fermion
   * \param b      Destination propagator
   * \param color_index  Color index
   */
  void FermToProp(const LatticeStaggeredFermionD& a, LatticeStaggeredPropagatorD& b, 
                  int color_index)
  {
    LatticeColorMatrixD bb = peekSpin(b, 0, 0);
    LatticeColorVectorD aa = peekSpin(a, 0);
 
    for(int i = 0; i < Nc; ++i)
      pokeColor(bb, peekColor(aa, i), i, color_index);
 
    pokeSpin(b, bb, 0, 0);
  }
 
 
 
 
  //! Extract a LatticeFermion from a LatticePropagator
  /*!
   * \ingroup ferm
   *
   * \param a      Source propagator
   * \param b      Destination fermion
   * \param color_index  Color index
   * \param spin_index   Spin index
   */
  void PropToFerm(const LatticePropagatorF& b, LatticeFermionF& a, 
                  int color_index, int spin_index)
  {
    for(int j = 0; j < Ns; ++j)
    {
      LatticeColorMatrixF bb = peekSpin(b, j, spin_index);
      LatticeColorVectorF aa = peekSpin(a, j);
 
      for(int i = 0; i < Nc; ++i)
        pokeColor(aa, peekColor(bb, i, color_index), i);
 
      pokeSpin(a, aa, j);
    }
  }
 
  //! Extract a LatticeFermion from a LatticePropagator
  /*!
   * \ingroup ferm
   *
   * \param a      Source propagator
   * \param b      Destination fermion
   * \param color_index  Color index
   * \param spin_index   Spin index
   */
  void PropToFerm(const LatticePropagatorD& b, LatticeFermionD& a, 
                  int color_index, int spin_index)
  {
    for(int j = 0; j < Ns; ++j)
    {
      LatticeColorMatrixD bb = peekSpin(b, j, spin_index);
      LatticeColorVectorD aa = peekSpin(a, j);
 
      for(int i = 0; i < Nc; ++i)
        pokeColor(aa, peekColor(bb, i, color_index), i);
 
      pokeSpin(a, aa, j);
    }
  }
 
 
 
  //! Extract a LatticeStaggeredFermion from a LatticeStaggeredPropagator
  /*!
   * \ingroup ferm
   *
   * \param a      Source propagator
   * \param b      Destination fermion
   * \param color_index  Color index
   */
  void PropToFerm(const LatticeStaggeredPropagatorF& b, LatticeStaggeredFermionF& a, 
                  int color_index)
  {
    LatticeColorMatrixF bb = peekSpin(b, 0, 0);
    LatticeColorVectorF aa = peekSpin(a, 0);
 
    for(int i = 0; i < Nc; ++i)
      pokeColor(aa, peekColor(bb, i, color_index), i);
 
    pokeSpin(a, aa, 0);
  }
 
  //! Extract a LatticeStaggeredFermion from a LatticeStaggeredPropagator
  /*!
   * \ingroup ferm
   *
   * \param a      Source propagator
   * \param b      Destination fermion
   * \param color_index  Color index
   */
  void PropToFerm(const LatticeStaggeredPropagatorD& b, LatticeStaggeredFermionD& a, 
                  int color_index)
  {
    LatticeColorMatrixD bb = peekSpin(b, 0, 0);
    LatticeColorVectorD aa = peekSpin(a, 0);
 
    for(int i = 0; i < Nc; ++i)
      pokeColor(aa, peekColor(bb, i, color_index), i);
 
    pokeSpin(a, aa, 0);
  }
 
}