anchor.diff:  
   kpr = pureDEGMV(g, kp, 'T')
   write(ilog,*) "k-point in: reciprocal space       reciprocal Cartesian space"
  -write(ilog,"(3F10.6,A,3F10.6)") , (kp(j), j=1,3), "     ", (kpr(j), j=1,3)
  -
  +write(ilog,"(3F10.6,A,3F10.6)") (kp(j), j=1,3), "     ", (kpr(j), j=1,3)
   
   ! write KPOINTS file ###########################################
   
  
positive.code:     read(iunt,fmt=*) dk
      read(iunt,fmt=*) band
      read(iunt,fmt=*) prg
      if(prg .eq. 'C') then
          read(iunt,'(3e20.12)') ((f(i,j),j=1,3),i=1,3)
      else if(prg .eq. 'V') then
          read(iunt,*) ((f(i,j),j=1,3),i=1,3)
      end if
  close(iunt)
  
  write(ilog,"(A,I5,F10.6)") "band, dk: ", band, dk
  
  
negative.code: ! y
  do i = 1, size(A)
      i1 = A(i)
      call set_next_eigeval(iunt, g, kpr, 0, i1, 0, w, dk, prg)
  end do
  
  ! z
  do i = 1, size(A)
      i1 = A(i)
      call set_next_eigeval(iunt, g, kpr, 0, 0, i1, w, dk, prg)
  end do
  
  ! xy
  