anchor.diff:  
               do ireport = 1, qmc_in%nreport
   
  -                call init_dmqmc_report_loop(bloom_stats, dmqmc_estimates, qs%spawn_store%rspawn)
  +                call init_dmqmc_report_loop(dmqmc_in%calc_excit_dist, bloom_stats, dmqmc_estimates, qs%spawn_store%rspawn)
                   tot_nparticles_old = qs%psip_list%tot_nparticles
   
                   do icycle = 1, qmc_in%ncycles
  
positive.code:         call bloom_stats_init_report_loop(bloom_stats)
  
          rspawn = 0.0_p
  
          if (allocated(dmqmc_estimates%excit_dist)) dmqmc_estimates%excit_dist = 0.0_p
          dmqmc_estimates%trace = 0.0_p
          dmqmc_estimates%numerators = 0.0_p
  
  
negative.code:         use bloom_handler, only: bloom_stats_t, bloom_stats_init_report_loop
          use dmqmc_data, only: dmqmc_estimates_t
  
          logical, intent(in) :: calc_excit_dist
          type(bloom_stats_t), intent(inout) :: bloom_stats
          type(dmqmc_estimates_t), intent(inout) :: dmqmc_estimates
          real(p), intent(out) :: rspawn
  