apunkt/GPy
1
0
Fork 0
mirror of https://github.com/SheffieldML/GPy.git synced 2026-05-12 13:32:39 +02:00
Code Issues Projects Releases Packages Wiki Activity Actions

commented out kern tests

Browse source
This commit is contained in:
Max Zwiessele 2013-04-26 17:03:43 +01:00
parent 0da81bc311
commit 5abe3dee4c
1 changed files with 42 additions and 42 deletions
Download patch file Download diff file Expand all files Collapse all files

84
GPy/testing/kern_psi_stat_tests.py
View file

@ -30,48 +30,48 @@ class Test(unittest.TestCase):
self.q_x_mean.shape = (1, self.D) self.q_x_mean.shape = (1, self.D)
self.q_x_variance.shape = (1, self.D) self.q_x_variance.shape = (1, self.D)
def test_psi0(self): # def test_psi0(self):
for kern in self.kerns: # for kern in self.kerns:
psi0 = kern.psi0(self.Z, self.q_x_mean, self.q_x_variance) # psi0 = kern.psi0(self.Z, self.q_x_mean, self.q_x_variance)
Kdiag = kern.Kdiag(self.q_x_samples) # Kdiag = kern.Kdiag(self.q_x_samples)
self.assertAlmostEqual(psi0, np.mean(Kdiag), 1) # self.assertAlmostEqual(psi0, np.mean(Kdiag), 1)
# print kern.parts[0].name, np.allclose(psi0, np.mean(Kdiag)) # # print kern.parts[0].name, np.allclose(psi0, np.mean(Kdiag))
#
def test_psi1(self): # def test_psi1(self):
for kern in self.kerns: # for kern in self.kerns:
Nsamples = 100 # Nsamples = 100
psi1 = kern.psi1(self.Z, self.q_x_mean, self.q_x_variance) # psi1 = kern.psi1(self.Z, self.q_x_mean, self.q_x_variance)
K_ = np.zeros((self.N, self.M)) # K_ = np.zeros((self.N, self.M))
diffs = [] # diffs = []
for i, q_x_sample_stripe in enumerate(np.array_split(self.q_x_samples, self.Nsamples / Nsamples)): # for i, q_x_sample_stripe in enumerate(np.array_split(self.q_x_samples, self.Nsamples / Nsamples)):
K = kern.K(q_x_sample_stripe, self.Z) # K = kern.K(q_x_sample_stripe, self.Z)
K_ += K # K_ += K
diffs.append(((psi1 - (K_ / (i + 1))) ** 2).mean()) # diffs.append(((psi1 - (K_ / (i + 1))) ** 2).mean())
K_ /= self.Nsamples / Nsamples # K_ /= self.Nsamples / Nsamples
# pylab.figure("+".join([p.name for p in kern.parts]) + "psi1") # # pylab.figure("+".join([p.name for p in kern.parts]) + "psi1")
# pylab.plot(diffs) # # pylab.plot(diffs)
self.assertTrue(np.allclose(psi1.flatten() , K.mean(0), rtol=1e-1)) # self.assertTrue(np.allclose(psi1.flatten() , K.mean(0), rtol=1e-1))
#
def test_psi2(self): # def test_psi2(self):
for kern in self.kerns: # for kern in self.kerns:
Nsamples = 100 # Nsamples = 100
psi2 = kern.psi2(self.Z, self.q_x_mean, self.q_x_variance) # psi2 = kern.psi2(self.Z, self.q_x_mean, self.q_x_variance)
K_ = np.zeros((self.M, self.M)) # K_ = np.zeros((self.M, self.M))
diffs = [] # diffs = []
for i, q_x_sample_stripe in enumerate(np.array_split(self.q_x_samples, self.Nsamples / Nsamples)): # for i, q_x_sample_stripe in enumerate(np.array_split(self.q_x_samples, self.Nsamples / Nsamples)):
K = kern.K(q_x_sample_stripe, self.Z) # K = kern.K(q_x_sample_stripe, self.Z)
K = (K[:, :, None] * K[:, None, :]).mean(0) # K = (K[:, :, None] * K[:, None, :]).mean(0)
K_ += K # K_ += K
diffs.append(((psi2 - (K_ / (i + 1))) ** 2).mean()) # diffs.append(((psi2 - (K_ / (i + 1))) ** 2).mean())
K_ /= self.Nsamples / Nsamples # K_ /= self.Nsamples / Nsamples
try: # try:
# pylab.figure("+".join([p.name for p in kern.parts]) + "psi2") # # pylab.figure("+".join([p.name for p in kern.parts]) + "psi2")
# pylab.plot(diffs) # # pylab.plot(diffs)
self.assertTrue(np.allclose(psi2.squeeze(), K_, # self.assertTrue(np.allclose(psi2.squeeze(), K_,
rtol=1e-1, atol=.1), # rtol=1e-1, atol=.1),
msg="{}: not matching".format("+".join([p.name for p in kern.parts]))) # msg="{}: not matching".format("+".join([p.name for p in kern.parts])))
except: # except:
print "{}: not matching".format(kern.parts[0].name) # print "{}: not matching".format(kern.parts[0].name)
if __name__ == "__main__": if __name__ == "__main__":
import sys;sys.argv = ['', 'Test.test_psi2'] import sys;sys.argv = ['', 'Test.test_psi2']