def kernel_call_ynone(X, length_scale=length_scale, periodicity=periodicity): dists = squareform_pdist(X, metric='euclidean') arg = numpy_pi * dists / periodicity sin_of_arg = numpy.sin(arg) K = numpy.exp(-2 * (sin_of_arg / length_scale)**2) return K
def kernel_call_ynone(X, length_scale=1.2, periodicity=1.1, pi=3.141592653589793): dists = squareform_pdist(X, metric='euclidean') arg = dists / periodicity * pi sin_of_arg = numpy.sin(arg) K = numpy.exp((sin_of_arg / length_scale)**2 * (-2)) return K
def kernel_rational_quadratic_none(X, length_scale=1.0, alpha=2.0): dists = squareform_pdist(X, metric='sqeuclidean') cst = py_pow(length_scale, 2) cst = py_mul(cst, alpha, 2) t_cst = py_make_float_array(cst) tmp = dists / t_cst t_one = py_make_float_array(1) base = tmp + t_one t_alpha = py_make_float_array(py_opp(alpha)) K = numpy.power(base, t_alpha) return K
def kernel_call_ynone(X, length_scale=1.2, periodicity=1.1, pi=3.141592653589793): dists = squareform_pdist(X, metric='euclidean') t_pi = py_make_float_array(pi) t_periodicity = py_make_float_array(periodicity) arg = dists / t_periodicity * t_pi sin_of_arg = numpy.sin(arg) t_2 = py_make_float_array(2) t__2 = py_make_float_array(-2) t_length_scale = py_make_float_array(length_scale) K = numpy.exp((sin_of_arg / t_length_scale) ** t_2 * t__2) return K