From fb444d893354a5c82c0ab3030204c14c8759e6ff Mon Sep 17 00:00:00 2001 From: Mike Croucher Date: Mon, 7 Sep 2015 16:23:10 +0100 Subject: [PATCH] Used scipy.log1p since it gives more consistent results cross-platform --- GPy/likelihoods/link_functions.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/GPy/likelihoods/link_functions.py b/GPy/likelihoods/link_functions.py index 3d753395..30ad32ad 100644 --- a/GPy/likelihoods/link_functions.py +++ b/GPy/likelihoods/link_functions.py @@ -2,6 +2,7 @@ # Licensed under the BSD 3-clause license (see LICENSE.txt) import numpy as np +import scipy from ..util.univariate_Gaussian import std_norm_cdf, std_norm_pdf import scipy as sp from ..util.misc import safe_exp, safe_square, safe_cube, safe_quad, safe_three_times @@ -67,7 +68,7 @@ class Probit(GPTransformation): .. math:: g(f) = \\Phi^{-1} (mu) - + """ def transf(self,f): return std_norm_cdf(f) @@ -140,7 +141,7 @@ class Log_ex_1(GPTransformation): """ def transf(self,f): - return np.log1p(safe_exp(f)) + return scipy.log1p(safe_exp(f)) def dtransf_df(self,f): ef = safe_exp(f)