anchor.diff:                      if (dmqmc_in%ipdmqmc .and. dmqmc_in%symmetric) &
                               & weighted_sampling%probs(sys%max_number_excitations+1) = &
                               & 0.5*(qs%target_beta-(iteration-1)*qs%dmqmc_factor*qs%tau)
  +
   
                       do idet = 1, qs%psip_list%nstates ! loop over walkers/dets
   
  
positive.code:         call copy_sys_spin_info(sys_copy, sys)
          call dealloc_det_info_t(cdet1, .false.)
          call dealloc_det_info_t(cdet2, .false.)
  
          call dSFMT_end(rng)
  
  
negative.code: 
              ! DMQMC quasi-newton not functional, so we artificially set this value to 0 which should not affect non-QN calcs.
              qs%estimators%proj_energy_old = 0_p
              
              do ireport = 1, nreport
  
                  call init_dmqmc_report_loop(calc_ref_proj_energy, dmqmc_in%calc_excit_dist, dmqmc_in%calc_mom_dist, &
  