Например, для двух 3D-массивов в R: N <- 1000 x <- rnorm(N*3*3); dim(x) <- c(N,3,3) y <- rnorm(N*3*3); dim(y) <- c(N,3,3) Я могу сделать следующее перекрестное произведение циклом: gg <- 0 for (n in 1:dim(x)[1]){...

Когда я пытаюсь умножить arma::mat и NumericVector в Rcpp с помощью оператора*, я получаю следующую ошибку: no match for operator* Вот пример того, что я пытаюсь умножить: NumericVector exampleNumericVector(L-1); arma::mat...

Для генерации стандартных однородных случайных чисел с помощью Rcpp я всегда использовал Rcpp::runif(1, 0, 1)[0] [0] связано с Rcpp::runif возвращением векторов . Недавно я обнаружил, что вы также можете использовать R API и вместо этого...