Работа с Rcpp::List матриц с плотными или разреженными типами матриц ⇐ C++

[Anonymous]

Мне нужно выполнить вычисления для каждой записи матрицы в Rcpp::List. Список содержит произвольную комбинацию плотных и разреженных матриц. Этот список передается в C++ из R. Мой вопрос касается как подходов к реализации, так и эффективности. Моя текущая мысль состоит в том, чтобы перебрать каждый элемент списка, проверить его тип (плотный или разреженный) и сгенерировать Arma::mat или Arma::sp_mat для элемента и использовать его для расчеты, желательно без копирования. Меня смущает следующее:

• Как проверить, является ли тип матрицы плотным или разреженным, из Rcpp::List[d], где d является индексом?
• Я могу использовать расширенный конструктор для Arma::mat через NumericMatrix, который позволяет избежать копирования (как показано ниже), но я не знаю, как сделать то же самое с арма::sp_mat. Эти матрицы большие, и я хочу избежать копирования.

Код: Выделить всё

Rcpp::NumericMatrix Ad = RcppListObject[d];
arma::mat ad=arma::mat( Ad.begin(), Ad.nrow(), Ad.ncol(), false );
//any needed calculations based on the d^{th} entry uses `ad`.

• Тесно связанный со 2, мне нужно будет выполнить множество вычислений с d^{th матрица в списке в разных точках алгоритма. Нужно ли мне каждый раз выполнять два шага, отмеченные в пункте 2, прежде чем я смогу выполнять с его помощью вычисления, или есть более эффективный путь? (Например, может ли Arma::field помочь и обработать комбинацию матриц Arma::mat и Arma::sp_mat без копирования? )

Большое спасибо!

Подробнее здесь: https://stackoverflow.com/questions/792 ... trix-types