generateRM.cuh File Reference

Functions
__global__ void	MatrixElementsInSector (Complex_t *dmatTot, Integer_t const LDT, Complex_t const *dmatLoc, Integer_t const LDL, Integer_t const numLoc, Integer_t const n, Integer_t const dloc, Integer_t const momentum, Integer_t const SecDim, Integer_t const *representatives, Integer_t const *periodicity)
__device__ Complex_t< double >	MAGMA_CEXP (Real_t theta)

Function Documentation

MAGMA_CEXP()

__device__ Complex_t< double > MAGMA_CEXP

(

Real_t

theta

)

                                                      {
  return MAGMA_Z_MAKE(cos(theta),sin(theta));
}

MatrixElementsInSector()

__global__ void MatrixElementsInSector	(	Complex_t *	dmatTot,
		Integer_t const	LDT,
		Complex_t const *	dmatLoc,
		Integer_t const	LDL,
		Integer_t const	numLoc,
		Integer_t const	n,
		Integer_t const	dloc,
		Integer_t const	momentum,
		Integer_t const	SecDim,
		Integer_t const *	representatives,
		Integer_t const *	periodicity
	)

                                                                                                                                                                                                                                                                                                                  {
  int idx = blockIdx.x*blockDim.x +threadIdx.x;
  int idy = blockIdx.y*blockDim.y +threadIdx.y;
  if( (idx>idy) || (idx>=SecDim) || (idy>=SecDim) ) return;
  dmatTot[idx+LDT*idy] = dComplexZero<>;
  Real_t theta;
  Integer_t Id1, Id2, locId1, locId2;
  Integer_t const dimMatLoc = (Integer_t)pow((float)dloc,(float)numLoc);
 
  for(int trans1 = 0;trans1 < periodicity[idx]; ++trans1) {
       Id1 = transSpin(representatives[idx],trans1,dloc,n);
    locId1 = Id1%dimMatLoc;
    for(int trans2 = 0;trans2 < periodicity[idy]; ++trans2) {
         Id2 = transSpin(representatives[idy],trans2,dloc,n);
      locId2 = Id2%dimMatLoc;
      if( Id1/dimMatLoc != Id2/dimMatLoc ) continue;
      theta = 2.0*M_PI*momentum*(trans2-trans1)/(Real_t)n;
      dmatTot[idx+LDT*idy] += dmatLoc[locId1+LDL*locId2]*MAGMA_CEXP(theta);
    }
  }
  dmatTot[idx+LDT*idy] /= sqrt( (Real_t)periodicity[idx]*periodicity[idy] );
  dmatTot[idy+LDT*idx]  = conj(dmatTot[idx+LDT*idy]);
}