From f8376765c680e8a5768a3483162c637828b52600 Mon Sep 17 00:00:00 2001 From: James Hensman Date: Wed, 12 Aug 2015 15:19:04 +0100 Subject: [PATCH] fixed strange bug. In python 3, numbers startin 0 are octal. --- GPy/inference/latent_function_inference/inferenceX.py | 2 +- GPy/util/univariate_Gaussian.py | 1 + 2 files changed, 2 insertions(+), 1 deletion(-) diff --git a/GPy/inference/latent_function_inference/inferenceX.py b/GPy/inference/latent_function_inference/inferenceX.py index 19013b06..c1c9fe2b 100644 --- a/GPy/inference/latent_function_inference/inferenceX.py +++ b/GPy/inference/latent_function_inference/inferenceX.py @@ -66,7 +66,7 @@ class InferenceX(Model): from ...models.ss_mrd import IBPPrior_SSMRD if isinstance(model.variational_prior, IBPPrior) or isinstance(model.variational_prior, IBPPrior_SSMRD): from ...core.parameterization.variational import SpikeAndSlabPrior - self.variational_prior = SpikeAndSlabPrior(pi=05,learnPi=False, group_spike=False) + self.variational_prior = SpikeAndSlabPrior(pi=0.5, learnPi=False, group_spike=False) else: self.variational_prior = model.variational_prior.copy() else: diff --git a/GPy/util/univariate_Gaussian.py b/GPy/util/univariate_Gaussian.py index e84f071f..97d912c2 100644 --- a/GPy/util/univariate_Gaussian.py +++ b/GPy/util/univariate_Gaussian.py @@ -8,6 +8,7 @@ from scipy.special import ndtr as std_norm_cdf #define a standard normal pdf _sqrt_2pi = np.sqrt(2*np.pi) def std_norm_pdf(x): + x = np.clip(x,-1e150,1e150) return np.exp(-np.square(x)/2)/_sqrt_2pi def inv_std_norm_cdf(x):