anchor.diff:          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))
  
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:         use proc_pointers, only: sc0_ptr
          use hamiltonian, only: get_hmatel
          use hamiltonian_data, only: hmatel_t
          use parallel, only: parent
  
          type(sys_t), intent(in) :: sys
          type(qmc_in_t), intent(in) :: qmc_in
  