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

fix the bug: the randomize function cannot properly handle variables with prior

Browse source
This commit is contained in:
Zhenwen Dai 2014-08-28 17:57:29 +01:00
parent 7bdbe21cd1
commit 91d1cd3131
1 changed files with 8 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

10
GPy/core/parameterization/parameter_core.py
View file

@ -778,9 +778,15 @@ class OptimizationHandlable(Indexable):
"""
# first take care of all parameters (from N(0,1))
x = rand_gen(size=self._size_transformed(), *args, **kwargs)
# now draw from prior where possible
[np.put(x, ind, p.rvs(ind.size)) for p, ind in self.priors.iteritems() if not p is None]
self.updates = False # Switch off the updates
self.optimizer_array = x # makes sure all of the tied parameters get the same init (since there's only one prior object...)
# now draw from prior where possible
x = self.param_array.copy()
[np.put(x, ind, p.rvs(ind.size)) for p, ind in self.priors.iteritems() if not p is None]
unfixlist = np.ones((self.size,),dtype=np.bool)
unfixlist[self.constraints[__fixed__]] = False
self.param_array[unfixlist] = x[unfixlist]
self.updates = True
#===========================================================================
# For shared memory arrays. This does nothing in Param, but sets the memory