c     f77 -o full2H full2H.f; full2H > 2Hkpt.full.dat
      PROGRAM FULL2H
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      PARAMETER (NX=13,SX=0.5D0)
      DIMENSION b(3,3)
      
      b(1,1) = 1.d0
      b(2,1) = 1./sqrt(3.d0)
      b(3,1) = 0.d0
      
      b(1,2) = 0.d0
      b(2,2) = 2/sqrt(3.d0)
      b(3,2) = 0.d0
      
      b(1,3) = 0.d0
      b(2,3) = 0.d0
      b(3,3) = sqrt(3.d0/8.d0)
      
      write(6,*) nx*nx*nx
      do i=1,nx
      do j=1,nx 
      do k=1,nx
      c1 = (i-1+sx)/nx
      c2 = (j-1+sx)/nx
      c3 = (k-1+sx)/nx
      fkx = 2 * (c1*b(1,1) + c2*b(1,2) + c3*b(1,3))
      fky = 2 * (c1*b(2,1) + c2*b(2,2) + c3*b(2,3))
      fkz = 2 * (c1*b(3,1) + c2*b(3,2) + c3*b(3,3))
      write (6,10) fkx, fky, fkz, 1.d0/nx/nx/nx
10    format(4f18.14)
      enddo
      enddo
      enddo
      stop
      end