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

Added first draft of polynomial kernel.

Browse source
This commit is contained in:
Neil Lawrence 2013-08-17 09:07:09 +02:00
parent a570198029
commit 3f6c01898d
4 changed files with 27 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

1
GPy/kern/parts/__init__.py
View file

@ -12,6 +12,7 @@ import mlp
import periodic_exponential
import periodic_Matern32
import periodic_Matern52
import poly
import prod_orthogonal
import prod
import rational_quadratic

2
GPy/kern/parts/kernpart.py
View file

@ -56,5 +56,5 @@ class Kernpart(object):
raise NotImplementedError
def dpsi2_dmuS(self,dL_dpsi2,Z,mu,S,target_mu,target_S):
raise NotImplementedError
def dK_dX(self,X,X2,target):
def dK_dX(self, dL_dK, X, X2, target):
raise NotImplementedError

9
GPy/kern/parts/mlp.py
View file

@ -111,10 +111,10 @@ class MLP(Kernpart):
gX = np.zeros((X2.shape[0], X.shape[1], X.shape[0]))
for i in range(X.shape[0]):
gX[:, :, i] = self._dK_dX_point(X, X2, i)
gX[:, :, i] = self._dK_dX_point(dL_dK, X, X2, target, i)
def _dK_dX_point(self, X, X2, i):
def _dK_dX_point(self, dL_dK, X, X2, target, i):
"""Gradient with respect to one point of X"""
inner_prod = self._K_inner_prod[i, :].T
@ -127,10 +127,11 @@ class MLP(Kernpart):
#arg = numer/denom
gX = np.zeros(X2.shape)
denom3 = denom*denom*denom
gX = np.zeros((X2.shape[0], X2.shape[1]))
for j in range(X2.shape[1]):
gX[:, j]=X2[:, j]/denom - vec2*X[i, j]*numer/denom3
gX[:, j] =X2[:, j]/denom - vec2*X[i, j]*numer/denom3
gX[:, j] = four_over_tau*self.weight_variance*self.variance*gX[:, j]/np.sqrt(1-arg*arg)
target[i, :]
def _K_computations(self, X, X2):