anchor.diff:      end do
       deallocate(multi0)
   
  +    allocate(phi(0:n_multi), stat=err)
  +
     end subroutine compute_multiplicities
   
     subroutine finalize()
  
positive.code:   end subroutine weight_std
  
    subroutine hnsb(estimate,err_estimate, err)
      use dirichlet_mod, only: h_dir
      real(real64), intent(out) :: estimate,err_estimate
      integer(int32), intent(out) :: err
      real(real64)              :: rslt,nrm
  
negative.code: 
      if (allocated(multi_z)) then
         deallocate(multi_z)
      end if
  
      if (allocated(phi)) then
         deallocate(phi)
      end if
  
    end subroutine finalize
  