anchor.diff:      end if
   
       open(19,file="debug_get_rotation_perms_lists.out")
  -
       ! Number of HNFs (superlattices); Index of the superlattices;
       ! Number of d-vectors in d set
       nH = size(HNF,3); n = determinant(HNF(:,:,1)); nD = size(RPList(1)%v,2)
  
positive.code:     dRPList%nL = nOp  ! Number of operations that fix the parent
      dRPList%RotIndx => null() !lattice (but may permute the d-vectors)
      
  
      do iOp = 1, nOp ! Try each operation in turn and see how the
                      ! d-vectors are permuted for each
         rd = matmul(rot(:,:,iOp),d)+spread(shift(:,iOp),2,nD) ! Rotate each d and add the shift
  
negative.code: 
      ! Beginning of main routine for enumeration
      open(23,file="VERSION.enum")
      write(23,'(A)') "v2.0.4-49-g067d-dirty"
      close(23)
  
      ![TODO] Get rid of all the junk that crept in (writing files, making inactives table, etc. These should all be in routines so that this main routine is still readable)
  