symtensor.cc

Go to the documentation of this file.

/*! \file
 *  \brief Compute symmetric tensors
 */
 
#include "chromabase.h"
#include "util/ferm/symtensor.h"
 
namespace Chroma 
{
 
  //! Anonymous namespace
  /*! \ingroup ferm */
  namespace
  {
    multi3d<int> sym_tensor3d;
    bool initP = false;
 
 
    //! Init gamma matrix multiplication table factors
    /*!
     * \ingroup ferm
     */
    void symTensor3dInit()
    {
      START_CODE();
 
      sym_tensor3d.resize(3,3,3);
      sym_tensor3d = 0;
 
      // d = |\epsilon^{i,j,k}|
      // Permutations: +(0,1,2)+(1,2,0)+(2,0,1)+(1,0,2)+(0,2,1)+(2,1,0)
 
      sym_tensor3d(0,1,2) = 1;
      sym_tensor3d(1,2,0) = 1;
      sym_tensor3d(2,0,1) = 1;
 
      sym_tensor3d(1,0,2) = 1;
      sym_tensor3d(0,2,1) = 1;
      sym_tensor3d(2,1,0) = 1;
  
      initP = true;
 
      END_CODE();
    }
  } // end namespace
 
 
  //! Return 3d symmetric tensor
  /*!
   * \ingroup ferm
   *
   * \return  \f$s_{ijk} = |\epsilon_{ijk}|\f$
   */
  int symTensor3d(int i, int j, int k)
  {
    if (! initP)
      symTensor3dInit();
 
    if (i < 0 || i > 2 || j < 0 || j > 2 || k < 0 || k > 2)
    {
      QDPIO::cerr << __func__ << ": indices out of bounds: i,j,k=" 
                  << i << " " << j << " " << k << std::endl;
      QDP_abort(1);
    }
 
    return sym_tensor3d(i,j,k);
  }
 
}  // end namespace Chroma