anchor.diff:          cluster%nexcitors = -1
           deexcit_count = 0
           do i = 0, max_size-min_size-1
  -            psize = psize + 1.0_p/2_int_64**(i+1)
  +            psize = psize + 1.0_p/(max_size-min_size+1)!2_int_64**(i+1)
               if (rand < psize) then
                   ! Found size!
                   cluster%nexcitors = i+min_size
  -                cluster%pselect = cluster%pselect/2_int_64**(i+1)
  +                cluster%pselect = cluster%pselect/(max_size-min_size+1)!2_int_64**(i+1)
                   exit
               end if
           end do
  
positive.code: 
          min_cluster_size = 0
          selection_data%nD0_select = 0 ! instead of this number of deterministic selections, these are chosen stochastically
          selection_data%nstochastic_clusters = nattempts!*10**3
          selection_data%nsingle_excitors = 0
  
          selection_data%nclusters = selection_data%nD0_select + selection_data%nsingle_excitors &
  
negative.code:                 proj_energy_cycle = cmplx(0.0, 0.0, p)
                  D0_population_cycle = cmplx(0.0, 0.0, p)
                  D0_population_ucc_cycle = 0.0_p
                  !$omp do schedule(dynamic,200) reduction(+:D0_population_cycle,proj_energy_cycle, D0_population_ucc_cycle, nattempts_spawn,ndeath)
                  do iattempt = 1, selection_data%nclusters
                      ! For OpenMP scalability, have this test inside a single loop rather
  