eKuf = tf.transpose(expectation(pXnew, (kernel, feature))) // M x N (psi1)
Kuu = feature.Kuu(kernel, jitter=default_jitter()) // M x M
Luu = tf.cholesky(Kuu) // M x M
if not white:
q_mu = tf.matrix_triangular_solve(Luu, q_mu, lower=True)
After Change
q_sqrt_r = tf.matrix_band_part(q_sqrt, -1, 0) // D x M x M
eKuf = tf.transpose(expectation(pXnew, (kernel, feature))) // M x N (psi1)
Luu = tf.cholesky(Kuu(feature, kernel, jitter=default_jitter()))
if not white:
q_mu = tf.matrix_triangular_solve(Luu, q_mu, lower=True)
Luu_tiled = tf.tile(Luu[None, :, :], [num_func, 1, 1]) // remove line once issue 216 is fixed
