anchor.diff:  
     end function log_weight
   
  -  elemental real(real64) function h_bayes(alpha)
  +  elemental real(real64) function h_dir(alpha)
       ! posterior average of the entropy given data and a specific alpha value
       ! computed from histogram multiplicities
       use gamma_funcs, only: digamma
  
positive.code:   end subroutine weight_std
  
    subroutine hnsb(estimate,err_estimate, err)
      use dirichlet_mod, only: h_bayes
      real(real64), intent(out) :: estimate,err_estimate
      integer(int32), intent(out) :: err
      real(real64)              :: rslt,nrm
  
negative.code: 
      err = 0
      if (ascale < 1.e-20) then
         estimate = h_dir(amax)
         err_estimate = 0.0
      else
         call quad(nrm_func,log_alpha1,log_alpha2, nrm, ierr)
  