anchor.diff:      integer              :: i_,k_,ni_
       integer              :: err
       integer              :: nmax
  -    integer, allocatable :: multi0(:)
  +    integer, allocatable :: wrk(:)
       real(8)                :: n_empty_bins
       integer              :: n_multi
   
  
positive.code:     ! nmax is the largest number of samples in a bin
      nbins = size(counts)
      nmax = maxval(counts)
      allocate(multi0(nmax),stat=err)
      ! multi0(n) is the number of states with frequency n
      multi0 = 0
      ! take into account the alphabet_size - nbins states with zero frequency
      n_empty_bins = alphabet_size - nbins
      do i_ = 1,nbins
  
negative.code:   integer :: i
    real(8)   :: ni,n_data
    integer              :: mi,nmax,err
    integer, allocatable :: wrk(:)
    logical :: multi = .false.
  
    if (multi) then
  