anchor.diff:  call find_permutation_of_group_and_dset(p_with_d,g,dmember,perm)
   allocate(labeling(nAt))
   labeling = aTyp(perm)
  -write(18,'("Input atom labels: ",200(i2,1x))') aTyp
  -write(18,'("Group order:       ",200(i2,1x))') perm
  -write(18,'("Labeling:          ",200(i2,1x))') labeling
  +!write(18,'("Input atom labels: ",200(i2,1x))') aTyp
  +!write(18,'("Group order:       ",200(i2,1x))') perm
  +!write(18,'("Labeling:          ",200(i2,1x))') labeling
   
   ENDSUBROUTINE find_labeling_from_atom_basis
   
  
positive.code:    p_with_d(1:3,pLindx:pLindx+n-1) = p
     p_with_d(4,  pLindx:pLindx+n-1) = iD
  enddo
  do i = 1, 4
     write(17,'(400(i2,1x))') p_with_d(i,:)
  enddo
  
  call matrix_inverse(matmul(L,Ainv),T)
  write(18,'("Map from group:",/,3(3(f7.3,1x),/),/)') transpose(T)
  
  ! If we want to do this, we need the supercell vectors
  !do iAt = 1, nAt
  
negative.code: real(dp) :: Ainv(3,3), T(3,3), gtemp(3)
  logical err, mapped
  
  !open(18,file="debug_find_labeling.out")
  !write(18,'(3/,"Takes a list of vectors, maps them into the group, determines the labeling",2/)')
  
  call matrix_inverse(A,Ainv,err)
  if(err) stop "Coplanar vectors in find_labeling_from_atom_basis"
  