code: def test_search5():
      # Test that the search/product work well together
      search = rft.IntrinsicVolumes([3, 4, 5])
      product = rft.IntrinsicVolumes([1, 2])
      prodsearch = product * search
      x = np.linspace(0, 5, 101)
      g1 = rft.Gaussian()
      g2 = rft.Gaussian(product=product)
      z = 0
      for i in range(prodsearch.mu.shape[0]):
          z += g1.density(x, i) * prodsearch.mu[i]
      y = g2(x, search=search)
      assert_almost_equal(y, z)
  
imports: import numpy
  import scipy.stats
  from scipy.special import gammaln, hermitenorm
  import pytest
  from numpy.testing import assert_almost_equal
  from ..None import rft
label_desc: Use image node_id 18 for calling a global function with example usage: test_search5() without return types