text: #   Copyright 2020 The PyMC Developers
  #
  #   Licensed under the Apache License, Version 2.0 (the "License");
  #   you may not use this file except in compliance with the License.
  #   You may obtain a copy of the License at
  #
  #       http://www.apache.org/licenses/LICENSE-2.0
  #
  #   Unless required by applicable law or agreed to in writing, software
  #   distributed under the License is distributed on an "AS IS" BASIS,
  #   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
  #   See the License for the specific language governing permissions and
  #   limitations under the License.
  
  from theano import function
  import theano.tensor as tt
  
  import pymc3.distributions.special as ps
  import scipy.special as ss
  import numpy as np
  
  from .checks import close_to
  
  
  def test_functions():
      xvals = list(map(np.atleast_1d, [0.01, 0.1, 2, 100, 10000]))
  
      x = tt.dvector("x")
      x.tag.test_value = xvals[0]
  
      p = tt.iscalar("p")
      p.tag.test_value = 1
  
      gammaln = function([x], ps.gammaln(x))
      psi = function([x], ps.psi(x))
      function([x, p], ps.multigammaln(x, p))
      for x in xvals:
          check_vals(gammaln, ss.gammaln, x)
      for x in xvals[1:]:
          check_vals(psi, ss.psi, x)
  
  
  """
  scipy.special.multigammaln gives bad values if you pass a non scalar to a
  In [14]:
  
      import scipy.special
      scipy.special.multigammaln([2.1], 3)
      Out[14]:
          array([ 1.76253257,  1.60450306,  1.66722239])
  """
  
  
  def t_multigamma():
      xvals = list(map(np.atleast_1d, [0, 0.1, 2, 100]))
  
      x = tt.dvector("x")
      x.tag.test_value = xvals[0]
  
      p = tt.iscalar("p")
      p.tag.test_value = 1
  
      multigammaln = function([x, p], ps.multigammaln(x, p))
  
      def ssmultigammaln(a, b):
          return ss.multigammaln(a[0], b)
  
      for p in [0, 1, 2, 3, 4, 100]:
          for x in xvals:
              check_vals(multigammaln, ssmultigammaln, x, p)
  
  
  def check_vals(fn1, fn2, *args):
      v = fn1(*args)
      close_to(v, fn2(*args), 1e-6)
  