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
