code: def test_dkln2():
      dim, offset = 3, 4.0
      m1 = np.zeros(dim)
      P1 = np.eye(dim)
      m2 = offset * np.ones(dim)
      P2 = np.eye(dim)
      assert dkl_gaussian(m1, P1, m2, P2) == 0.5 * dim * offset**2
  
file_path: nipy/nipy/algorithms/clustering/tests/test_bgmm.py
imports: import numpy
  import numpy.random
  from ..bgmm import BGMM, VBGMM, dirichlet_eval, dkl_gaussian, multinomial
label_desc: Use image node_id 4 for calling a global function with example usage: test_dkln2() without return types