anchor.diff:                  call convert_excitor_to_determinant(cdet%f, excitation_level, cluster_population)
   
                   ! Normalisation factor for amplitudes...
  -                amplitude = real(cluster_population,p)/(D0_population**(cluster_size-1)) ! TODO: don't accumulate D0 over report loops
  +                amplitude = real(cluster_population,p)/(D0_normalisation**(cluster_size-1)) ! TODO: don't accumulate D0 over report loops
   
                   ! We chose excitors uniformly.
                   !  -> prob. of each excitor is 1./(number of excitors on processor).
  
positive.code:             ! Must be the reference.
              cdet%f = f0
              excitation_level = 0
              amplitude = D0_population ! TODO: don't accumulate D0 over report loops
              ! Only one cluster of this size to choose => p_clust = 1
          case default
              ! Select cluster from the excitors on the current processor in
  
negative.code: 
          logical :: soft_exit
  
          integer :: pos
          logical :: hit
  
          real :: t1, t2
  
          ! Allocate det_info components.
  