anchor.diff:      ! call it y. Then the renyi-2 entropy is then given by -log_2(x/y).
       real(p), allocatable :: renyi_2(:) ! (nrdms)
   end type dmqmc_inst_rdms_t
  -
   
   !--- Type for a ground state RDM ---
   type dmqmc_ground_rdm_t
  
positive.code:     real(p) :: trace
  end type dmqmc_ground_rdm_t
  
  !--- Type for weighted sampling parameters ---
  type dmqmc_weighted_sampling_t
      ! This holds the factors by which the populations on each excitation level
  
negative.code:     integer(i0), allocatable :: end1(:), end2(:)
  end type rdm_t
  
  !--- Type for all instantaneous RDMs ---
  type dmqmc_inst_rdms_t
      ! The total number of rdms beings calculated.
  