anchor.diff:  
     if (alpha < 1.0e-10_real64) then
        ! if alpha == 0.0 (no pseudocounts)
  -     call plugin(n, hist, estimate)
  +     call plugin(n, counts, estimate)
        return
     end if
   
  -  nbins = size(hist)
  +  nbins = size(counts)
     if (nbins == 1) then 
        estimate = 0.0_real64
        return
     end if
  -  ndata = sum(hist)
  +  ndata = sum(counts)
     estimate = 0.0_real64
     do i = 1,nbins
  -     ni = hist(i) + alpha
  +     ni = counts(i) + alpha
        estimate = estimate - ni*log(ni)
     end do
     ! correct for the (nc - nbins) bins with frequency alpha
  
positive.code: 
  end subroutine plugin
  
  subroutine pseudo(n,hist,nc,alpha,estimate)
    use iso_fortran_env
    ! pseudocount estimator(s)
    ! estimate the bin frequencies using pseudocounts 
  
negative.code: 
  end subroutine dirichlet
  
  subroutine nsb(n,counts,nc,estimate,err_estimate)
    use iso_fortran_env
    use dirichlet_mod, only: compute_multiplicities,dirichlet_finalize
    use nsb_mod, only: hnsb
    use nsb_mod, only: compute_integration_range
    implicit none
  
    integer(int32), intent(in)  :: n
    integer(int32), intent(in)  :: counts(n)
    integer(int32), intent(in)  :: nc
    real(real64),   intent(out) :: estimate
    real(real64),   intent(out) :: err_estimate
  
    if (size(counts) == 1) then 
       estimate = 0.0_real64
       err_estimate = 0.0_real64
       return
    end if
  
    call compute_multiplicities(nc,counts)
  
    call compute_integration_range()
  
  