anchor.diff:          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
  +        if (calc_excit_dist) dmqmc_estimates%excit_dist = 0.0_p
           dmqmc_estimates%trace = 0.0_p
           dmqmc_estimates%numerators = 0.0_p
   
  
positive.code: 
              do ireport = 1, qmc_in%nreport
  
                  call init_dmqmc_report_loop(bloom_stats, dmqmc_estimates, qs%spawn_store%rspawn)
                  tot_nparticles_old = qs%psip_list%tot_nparticles
  
                  do icycle = 1, qmc_in%ncycles
  
negative.code:             ! Calculate and output all requested estimators based on the reduced
              ! density matrix. This is for ground-state RDMs only.
              if (dmqmc_in%rdm%calc_ground_rdm) call call_ground_rdm_procedures(dmqmc_estimates, beta_cycle, dmqmc_in%rdm)
              ! Calculate and output new weights based on the psip distriubtion in
              ! the previous loop.
              if (dmqmc_in%find_weights) call output_and_alter_weights(dmqmc_in, sys%max_number_excitations, &
                                                                       dmqmc_estimates%excit_dist, weighted_sampling)
  