24 using namespace shogun;
25 using namespace Eigen;
56 "Exact inference method can only use Gaussian likelihood function\n")
58 "Labels must be type of CRegressionLabels\n")
100 float64_t result=(eigen_y-eigen_m).dot(eigen_alpha)/2.0+
101 eigen_L.diagonal().array().log().sum()+
m_L.
num_rows*
178 a=L.triangularView<Upper>().adjoint().solve(eigen_y-eigen_m);
179 a=L.triangularView<Upper>().solve(a);
195 Map<VectorXd> eigen_mu(m_mu.
vector, m_mu.
vlen);
213 MatrixXd eigen_V=eigen_L.triangularView<Upper>().adjoint().solve(
235 eigen_Q=eigen_L.triangularView<Upper>().adjoint().solve(
237 eigen_Q=eigen_L.triangularView<Upper>().solve(eigen_Q);
243 eigen_Q-=eigen_alpha*eigen_alpha.transpose();
249 REQUIRE(!strcmp(param->
m_name,
"scale"),
"Can't compute derivative of "
250 "the nagative log marginal likelihood wrt %s.%s parameter\n",
259 result[0]=(eigen_Q.cwiseProduct(eigen_K)*
m_scale*2.0).
sum()/2.0;
267 REQUIRE(!strcmp(param->
m_name,
"sigma"),
"Can't compute derivative of "
268 "the nagative log marginal likelihood wrt %s.%s parameter\n",
283 result[0]=
CMath::sq(sigma)*eigen_Q.trace();
300 "Length of the parameter %s should not be NULL\n", param->
m_name)
338 "Length of the parameter %s should not be NULL\n", param->
m_name)
356 Map<VectorXd> eigen_dmu(dmu.
vector, dmu.
vlen);
359 result[i]=-eigen_dmu.
dot(eigen_alpha);