anchor.diff:    integer                        :: N_holes(2), N_particles(2)
     integer                        :: hole_list(N_int*bit_kind_size,2)
     integer                        :: particle_list(N_int*bit_kind_size,2)
  -  integer(bit_kind), allocatable:: preinteresting_det(:,:,:)
  -  allocate (preinteresting_det(N_int,2,N_det))
   
     call bitstring_to_list_ab(hole    , hole_list    , N_holes    , N_int)
     call bitstring_to_list_ab(particle, particle_list, N_particles, N_int)
  
positive.code:       
        do ii=1,preinteresting(0)
          i = preinteresting(ii)
          mobMask(1,1) = iand(negMask(1,1), preinteresting_det(1,1,ii))
          mobMask(1,2) = iand(negMask(1,2), preinteresting_det(1,2,ii))
          nt = popcnt(mobMask(1, 1)) + popcnt(mobMask(1, 2))
          do j=2,N_int
            mobMask(j,1) = iand(negMask(j,1), preinteresting_det(j,1,ii))
            mobMask(j,2) = iand(negMask(j,2), preinteresting_det(j,2,ii))
            nt = nt+ popcnt(mobMask(j, 1)) + popcnt(mobMask(j, 2))
          end do
          
          if(nt <= 4) then
            interesting(0) += 1
            interesting(interesting(0)) = i
            minilist(1,1,interesting(0)) = preinteresting_det(1,1,ii)
            minilist(1,2,interesting(0)) = preinteresting_det(1,2,ii)
            do j=2,N_int
              minilist(j,1,interesting(0)) = preinteresting_det(j,1,ii)
              minilist(j,2,interesting(0)) = preinteresting_det(j,2,ii)
            enddo
            if(nt <= 2) then
              fullinteresting(0) += 1
              fullinteresting(fullinteresting(0)) = i
              fullminilist(1,1,fullinteresting(0)) = preinteresting_det(1,1,ii)
              fullminilist(1,2,fullinteresting(0)) = preinteresting_det(1,2,ii)
              do j=2,N_int
                fullminilist(j,1,fullinteresting(0)) = preinteresting_det(j,1,ii)
                fullminilist(j,2,fullinteresting(0)) = preinteresting_det(j,2,ii)
              enddo
            end if
          end if
  
negative.code:       if(i <= N_det_selectors) then
          preinteresting(0) += 1
          preinteresting(preinteresting(0)) = i
        else if(nt <= 2) then
          prefullinteresting(0) += 1
          prefullinteresting(prefullinteresting(0)) = i
  