anchor.diff:          integer :: ndets, i
           real(p) :: e_max
           type(hmatel_t) :: offdiagel
  +        integer :: iunit
  +        real :: t1, t2
  +
  +        iunit = 6
   
           qs%cheby_prop%using_chebyshev = qmc_in%chebyshev
           qs%cheby_prop%icheb = 1
  
positive.code:             qs%cheby_prop%spectral_range(1) = 0.0_p
              qs%cheby_prop%spectral_range(2) = e_max - qs%ref%H00
              call update_chebyshev(qs%cheby_prop, 0.0_p)
          else
              allocate(qs%cheby_prop%weights(1))
              qs%cheby_prop%weights(1) = 1.0_p ! Similar to QN, the hmatel is always scaled to save on too many checks (branching)
  
negative.code:         qs%cheby_prop%order = qmc_in%chebyshev_order
          
          if (qs%cheby_prop%using_chebyshev) then
              call cpu_time(t1)
  
              allocate(qs%cheby_prop%zeroes(qs%cheby_prop%order))
              allocate(qs%cheby_prop%weights(qs%cheby_prop%order))
              allocate(occ_list_max(sys%nel))
  