anchor.diff:  !***************************************************************************************************
   ! Map a list of real space atomic basis vectors and their labels into g-space and extracts the
   ! labeling 
  -SUBROUTINE find_labeling_from_atom_basis(L,A,aBas,aTyp,SNF,labeling)
  +SUBROUTINE find_labeling_from_atom_basis(L,A,aBas,aTyp,SNF,labeling,eps)
   integer, dimension(3,3), intent(in) :: L
   real(dp), dimension(3,3), intent(in) :: A
   real(dp), dimension(:,:), pointer :: aBas ! (in)
   integer, dimension(:), intent(in) :: aTyp
   integer, dimension(3,3), intent(in) :: SNF
   integer, pointer :: labeling(:) ! out
  +real(dp), intent(in) :: eps
   
   integer, dimension(3,size(aBas,2)) :: g
   integer, pointer :: p(:,:)
  
positive.code:    foundHNF = .false.
     do iHNF = 1, nHNF
        if(.not. all(HNF==HNFin(:,:,iHNF))) cycle
     enddo
     if(.not. foundHNF) then
        write(13,'("HNF doesn''t match for str #:",i8," ---Skipping")') StrN
  
negative.code: do ip = 1, nP
     write(17,'("Perm #",i3,":",1x,200(i2,1x))') ip,LattRotList(1)%perm(ip,:)
  enddo
  call find_labeling_from_atom_basis(L(:,:,1),pLV,aBas,aTyp,SNF,labeling,eps)
  write(17,'("Input atom labels: ",200(i2,1x))') aTyp
  write(17,'("Labeling:          ",200(i2,1x))') labeling
  
  