anchor.diff:            dap = 0
             do iD = 1, nD ! Loop over each row in the (d,g) table
                ! LA^-1(v_i+(RAL^-1)G)
  +             ! v_i is the lattice vector that must be added to return a rotated d-vector to the unit cell
                rgp = matmul(Tinv,(-spread(RPList(iH)%v(:,iD,iOp),2,n)+matmul(matmul(Op(iH)%rot(:,:,iOp),T),g)))
                rag = transpose(matmul(transpose(arrow_basis),op(iH)%rot(:,:,iOp)))
                if (.not. equal(rgp,nint(rgp),eps)) stop "Transform left big fractional parts"
  
positive.code:                 write(11,'("failed at ",2i3)')  jP, iP
                  write(11,'("testperm ",20i2,/)') testperm
                  write(11,'("permlist ",8i2)') transpose(rpl(iL)%perm)
                  !write(*,'(20i2,/)') testperm
                  print *
                  !do itest = 1,nP
                  !   write(*,'(20i2)') rpl(iL)%perm(itest,:)
                  !enddo
                  do_rotperms_form_groups = .false.
                  !exit lists
                  exit perms
               endif
            enddo
  
negative.code:     endif
      if (full) then; write(14,'("full list of labelings (including incomplete labelings) is used")')
      else; write(14,'("partial list of labelings (complete labelings only) is used")'); endif
      write(14,'("(Non)Equivalency list:" ,40(I2,1x))') equivalencies(:)
  
      !write(14,'("Symmetry of the primary lattice is of order ",i2)')
  
  