anchor.diff:                  if (dets(i)%Ms == dets(j) % Ms) then
                       excitation = get_excitation(dets(i)%f, dets(j)%f)
                       if (excitation%nexcit <= 2) then
  -                        root_det = decode_det(dets(i)%f)
  -                        hamil(i,j) = get_hmatel(root_det, excitation)
  +                        hamil(i,j) = get_hmatel(dets(i)%f, excitation)
                       else
                           hamil(i,j) = 0.0_dp
                       end if
  
positive.code:             end do
  
          case(1)
  
              ! One electron operator
              hmatel = hmatel + get_one_e_int(excitation%from_orb(1), excitation%to_orb(1)) 
  
negative.code: 
          integer :: ierr, i, j
          type(excit) :: excitation
  
          allocate(hamil(ndets,ndets), stat=ierr)
  
  