#include<config.h>
CPS_START_NAMESPACE
/*!\file
  \brief  Definitions of functions computing 
  the characters of some SU(3) representations.
*/

CPS_END_NAMESPACE
#include <util/data_types.h>
CPS_START_NAMESPACE

static IFloat f1o3 = 1.0 / 3.0 ;
static IFloat f2o3 = 2.0 / 3.0 ;

IFloat reChar6(IFloat *p) ;
IFloat imChar6(IFloat *p) ;

IFloat reChar8(IFloat *p) ;

IFloat reChar10(IFloat *p) ;
IFloat imChar10(IFloat *p) ;

IFloat reChar6(IFloat *p)
{
//  return   *(p   ) * *(p   ) - *(p+ 1) * *(p+ 1)
//         + *(p+ 8) * *(p+ 8) - *(p+ 9) * *(p+ 9)
//         + *(p+16) * *(p+16) - *(p+17) * *(p+17)
//         + *(p   ) * *(p+ 8) - *(p+ 1) * *(p+ 9)
//         + *(p   ) * *(p+16) - *(p+ 1) * *(p+17)
//         + *(p+ 8) * *(p+16) - *(p+ 9) * *(p+17)
//         + *(p+ 2) * *(p+ 6) - *(p+ 3) * *(p+ 7)
//         + *(p+ 4) * *(p+12) - *(p+ 5) * *(p+13)
//         + *(p+10) * *(p+14) - *(p+11) * *(p+15) ;

  return   *(p   ) * ( *(p   ) + *(p+ 8) + *(p+16) )
         - *(p+ 1) * ( *(p+ 1) + *(p+ 9) + *(p+17) )
         + *(p+ 8) * ( *(p+ 8) + *(p+16) )
	 - *(p+ 9) * ( *(p+ 9) + *(p+17) )
         + *(p+16) * *(p+16) - *(p+17) * *(p+17)
         + *(p+ 2) * *(p+ 6) - *(p+ 3) * *(p+ 7)
         + *(p+ 4) * *(p+12) - *(p+ 5) * *(p+13)
         + *(p+10) * *(p+14) - *(p+11) * *(p+15) ;
}


IFloat imChar6(IFloat *p)
{
//  return   *(p   ) * *(p+ 1) + *(p   ) * *(p+ 1)
//         + *(p+ 8) * *(p+ 9) + *(p+ 8) * *(p+ 9)
//         + *(p+16) * *(p+17) + *(p+16) * *(p+17)
//         + *(p   ) * *(p+ 9) + *(p+ 1) * *(p+ 8)
//         + *(p   ) * *(p+17) + *(p+ 1) * *(p+16)
//         + *(p+ 8) * *(p+17) + *(p+ 9) * *(p+16)
//         + *(p+ 2) * *(p+ 7) + *(p+ 3) * *(p+ 6)
//         + *(p+ 4) * *(p+13) + *(p+ 5) * *(p+12)
//         + *(p+10) * *(p+15) + *(p+11) * *(p+14) ;

  return   *(p   ) * ( *(p+ 1) + *(p+ 9) + *(p+17) )
         + *(p+ 1) * ( *(p   ) + *(p+ 8) + *(p+16) )
         + *(p+ 8) * ( *(p+ 9) + *(p+17) )
	 + *(p+ 9) * ( *(p+ 8) + *(p+16) )
         + *(p+16) * *(p+17) + *(p+16) * *(p+17)
         + *(p+ 2) * *(p+ 7) + *(p+ 3) * *(p+ 6)
         + *(p+ 4) * *(p+13) + *(p+ 5) * *(p+12)
         + *(p+10) * *(p+15) + *(p+11) * *(p+14) ;
}

IFloat reChar8(IFloat *p) 
{
  return    2.0 * (   *(p   ) * ( *(p+ 8) + *(p+16) )
                    + *(p+ 1) * ( *(p+ 9) + *(p+17) )
                    + *(p+ 8) * *(p+16) + *(p+ 9) * *(p+17) )
	 + f2o3 * (   *(p   ) * *(p   ) + *(p+ 1) * *(p+ 1)
                    + *(p+ 8) * *(p+ 8) + *(p+ 9) * *(p+ 9)
                    + *(p+16) * *(p+16) + *(p+17) * *(p+17) )
         - f1o3 * (   *(p+ 2) * *(p+ 2) + *(p+ 3) * *(p+ 3)
                    + *(p+ 4) * *(p+ 4) + *(p+ 5) * *(p+ 5)
                    + *(p+ 6) * *(p+ 6) + *(p+ 7) * *(p+ 7)
                    + *(p+10) * *(p+10) + *(p+11) * *(p+11)
                    + *(p+12) * *(p+12) + *(p+13) * *(p+13)
                    + *(p+14) * *(p+14) + *(p+15) * *(p+15) ) ;
}

IFloat reChar10(IFloat *p)
{
  return 0.0 ;
}

IFloat imChar10(IFloat *p)
{
  return 0.0 ;
}

CPS_END_NAMESPACE