apunkt/GPy
1
0
Fork 0
mirror of https://github.com/SheffieldML/GPy.git synced 2026-04-30 15:26:23 +02:00
Code Issues Projects Releases Packages Wiki Activity Actions

warning messages optional for re-constraining

Browse source
This commit is contained in:
Max Zwiessele 2013-10-14 17:59:04 +01:00
parent e0d7de9b3b
commit 3a97eaf684
2 changed files with 20 additions and 18 deletions
Download patch file Download diff file Expand all files Collapse all files

4
GPy/core/model.py
View file

@ -250,7 +250,7 @@ class Model(Parameterized):
else:
self._set_params_transformed(initial_parameters)
def ensure_default_constraints(self):
def ensure_default_constraints(self, warning=True):
"""
Ensure that any variables which should clearly be positive
have been constrained somehow. The method performs a regular
@ -268,7 +268,7 @@ class Model(Parameterized):
if not (i in currently_constrained):
to_make_positive.append(i)
if len(to_make_positive):
self.constrain_positive(np.asarray(to_make_positive))
self.constrain_positive(np.asarray(to_make_positive), warning=warning)
def objective_function(self, x):
"""

34
GPy/core/parameterized.py
View file

@ -193,17 +193,17 @@ class Parameterized(object):
else:
self.fixed_indices, self.fixed_values = [], []
def constrain_negative(self, regexp):
def constrain_negative(self, regexp, warning=True):
""" Set negative constraints. """
self.constrain(regexp, transformations.NegativeLogexp())
self.constrain(regexp, transformations.NegativeLogexp(), warning)
def constrain_positive(self, regexp):
def constrain_positive(self, regexp, warning=True):
""" Set positive constraints. """
self.constrain(regexp, transformations.Logexp())
self.constrain(regexp, transformations.Logexp(), warning)
def constrain_bounded(self, regexp, lower, upper):
def constrain_bounded(self, regexp, lower, upper, warning=True):
""" Set bounded constraints. """
self.constrain(regexp, transformations.Logistic(lower, upper))
self.constrain(regexp, transformations.Logistic(lower, upper), warning)
def all_constrained_indices(self):
if len(self.constrained_indices) or len(self.fixed_indices):
@ -211,17 +211,18 @@ class Parameterized(object):
else:
return np.empty(shape=(0,))
def constrain(self, regexp, transform):
def constrain(self, regexp, transform, warning=True):
assert isinstance(transform, transformations.Transformation)
matches = self.grep_param_names(regexp)
overlap = set(matches).intersection(set(self.all_constrained_indices()))
if overlap:
self.unconstrain(np.asarray(list(overlap)))
print 'Warning: re-constraining these parameters'
pn = self._get_param_names()
for i in overlap:
print pn[i]
if warning:
print 'Warning: re-constraining these parameters'
pn = self._get_param_names()
for i in overlap:
print pn[i]
self.constrained_indices.append(matches)
self.constraints.append(transform)
@ -229,7 +230,7 @@ class Parameterized(object):
x[matches] = transform.initialize(x[matches])
self._set_params(x)
def constrain_fixed(self, regexp, value=None):
def constrain_fixed(self, regexp, value=None, warning=True):
"""
:param regexp: which parameters need to be fixed.
@ -249,10 +250,11 @@ class Parameterized(object):
overlap = set(matches).intersection(set(self.all_constrained_indices()))
if overlap:
self.unconstrain(np.asarray(list(overlap)))
print 'Warning: re-constraining these parameters'
pn = self._get_param_names()
for i in overlap:
print pn[i]
if warning:
print 'Warning: re-constraining these parameters'
pn = self._get_param_names()
for i in overlap:
print pn[i]
self.fixed_indices.append(matches)
if value != None: