anchor.diff:      logical, optional, intent(in) :: polya
       logical, optional, intent(in) :: origCrossOutAlgorithm
   
  -    integer :: iD, i, ivol, LatDim, Scnt, Tcnt, iBlock, HNFcnt, status, iC, j
  +    integer :: iD, i, ivol, LatDim, Scnt, Tcnt, iBlock, HNFcnt, iC, j
       integer, pointer, dimension(:,:,:) :: HNF => null(),SNF => null(), L => null(), R => null()
       integer, pointer :: SNF_labels(:) =>null(), uqSNF(:,:,:) => null()
       integer, pointer, dimension(:,:,:) :: rdHNF =>null()
  
positive.code:     type(RotPermList), intent(in) :: dperms
      type(RotPermList), pointer :: RPList(:)
      real(dp), pointer :: latts(:,:,:)
  
      integer nD
      real(dp), allocatable:: sgrots(:,:,:), sgshift(:,:)
      real(dp), dimension(3,3) :: test_latticei, test_latticej
      integer i, Nhnf, iuq, irot, j, nRot, Nq, status
      integer, allocatable :: temp_hnf(:,:,:)
  
negative.code:     real(dp), intent(in) :: eps
  
      integer iD, nOp, iOp, status!, iC!, nDfull
      integer, pointer :: tList(:,:)
      real(dp) :: rd(size(d,1),size(d,2)),tRD(size(d,1),size(d,2))
      real(dp) :: inv_pLV(3,3) ! Inverse of the pLV matrix
      real(dp), pointer:: tv(:,:,:)
  