Умножение матриц с использованием NumericMatrix и NumericVector в Rcpp
Мне интересно, есть ли способ расчета умножения матриц с использованием классов NumericMatrix и NumericVector. Я задаюсь вопросом, есть ли какой-нибудь простой способ помочь мне избежать следующего цикла для проведения этого вычисления. Я просто хочу рассчитать X% *% бета.
// assume X and beta are initialized and X is of dimension (nsites, p),
// beta is a NumericVector with p elements.
for(int j = 0; j < nsites; j++)
{
temp = 0;
for(int l = 0; l < p; l++) temp = temp + X(j,l) * beta[l];
}
Заранее большое спасибо!