def test_construction():
m = TensorProductMean(lambda x: x**2)
x = np.random.randn(10, 1)
c = Cache()
yield m, x
yield m, x, c
yield m, Observed(x)
yield m, Observed(x), c
def test_compare_noisy_kernel_and_additive_component_kernel():
def NoisyKernelCopy(k1, k2):
return AdditiveComponentKernel({
Component('wiggly'): k1,
Component('noise'): k2
}, latent=[Component('wiggly')])
x = np.random.randn(10, 2)
k1 = EQ()
k2 = RQ(1e-2)
kn1 = NoisyKernel(k1, k2)
kn2 = NoisyKernelCopy(k1, k2)
yield ok, allclose(kn1(Latent(x)), k1(x)), 'noisy kernel 1'
yield ok, allclose(kn1(Latent(x), Observed(x)), k1(x)), 'noisy kernel 2'
yield ok, allclose(kn1(Observed(x), Latent(x)), k1(x)), 'noisy kernel 3'
yield ok, allclose(kn1(Observed(x)), (k1 + k2)(x)), 'noisy kernel 4'
yield ok, allclose(kn2(Latent(x)), k1(x)), 'noisy copy 1'
yield ok, allclose(kn2(Latent(x), Observed(x)), k1(x)), 'noisy copy 2'
yield ok, allclose(kn2(Observed(x), Latent(x)), k1(x)), 'noisy copy 3'
yield ok, allclose(kn2(Observed(x)), (k1 + k2)(x)), 'noisy copy 4'
yield ok, allclose(kn2(Latent(x), Component('noise')(x)),
ZeroKernel()(x)), 'additive 1'
yield ok, allclose(kn2(Component('noise')(x), Latent(x)),
ZeroKernel()(x)), 'additive 2'
yield ok, allclose(kn2(Observed(x), Component('noise')(x)),
k2(x)), 'additive 3'
yield ok, allclose(kn2(Component('noise')(x), Observed(x)),
k2(x)), 'additive 4'
yield eq, str(NoisyKernel(EQ(), RQ(1))), \
'NoisyKernel(EQ(), RQ(1))'
def test_construction():
k = EQ()
x = np.random.randn(10, 1)
k(x)
k(x, x)
k(Observed(x))
k(Observed(x), Observed(x))
k(x, Observed(x))
k(Observed(x), x)
k.elwise(x)
k.elwise(x, x)
k.elwise(Observed(x))
k.elwise(Observed(x), Observed(x))
k.elwise(x, Observed(x))
k.elwise(Observed(x), x)
def test_construction():
m = TensorProductMean(lambda x: x**2)
x = np.random.randn(10, 1)
m(x)
m(Observed(x))
def test_construction():
k = EQ()
x = np.random.randn(10, 1)
c = Cache()
yield k, x
yield k, x, c
yield k, x, x
yield k, x, x, c
yield k, Observed(x)
yield k, Observed(x), c
yield k, Observed(x), Observed(x)
yield k, Observed(x), Observed(x), c
yield k, x, Observed(x)
yield k, x, Observed(x), c
yield k, Observed(x), x
yield k, Observed(x), x, c
yield k.elwise, x
yield k.elwise, x, c
yield k.elwise, x, x
yield k.elwise, x, x, c
yield k.elwise, Observed(x)
yield k.elwise, Observed(x), c
yield k.elwise, Observed(x), Observed(x)
yield k.elwise, Observed(x), Observed(x), c
yield k.elwise, x, Observed(x)
yield k.elwise, x, Observed(x), c
yield k.elwise, Observed(x), x
yield k.elwise, Observed(x), x, c