Armadillo и C++ - Инициализация матрицы с помощью 2D-массива

Question

Мне нужно инициализировать матрицу Armadillo массивом двойников. Я нашел этот конструктор в оригинальной документации:

mat(*aux_mem, n_rows, n_cols, copy_aux_mem = true, strict = true) 

, и я нашел вопрос об этом здесь на SO:

armadillo C++: matrix initialization from array

Проблема заключается в том, оказывается, что конструктор работает только для инициализации одномерными массивами, а не 2D. Если я пытаюсь его использовать так:

double **matrix = new double *[block_size]; 
for(int i = 0; i < block_size; i++) { 
    matrix[i] = new double[block_size]; 
} 

arma::mat arma_matrix(&matrix[0][0], matrix_size, matrix_size, true, true); 

cout << "am: " << arma_matrix[1][0] << endl; 

Я получаю ошибку:

fined_grain:103/init_function: In function ‘void place_user_fn_103(ca::Context&, ca::TokenList<double>&)’: 
fined_grain:103/init_function:61:42: error: invalid types ‘double[int]’ for array subscript 

Итак, что идеальный способ для инициализации матрицы Arma с 2D-массив? Я предпочитаю самое быстрое решение, потому что мне нужно использовать большие матрицы.

Answer 1

Я быстро взглянул на armadillo library documentation, и я вижу следующие проблемы.

1. Аргумент вы передаете к arma_matrix синтаксически правильно, но неправильно. Вы должны использовать:

   double *matrix = new double [block_size*block_size]; 
   arma::mat arma_matrix(matrix, block_size, block_size, true, true); 
   

2. синтаксис для доступа к элементам является:

   cout << "am: " << arma_matrix(1, 0) << endl; 
   

   Вы также можете использовать:

   int row = 1; 
   int col = 0; 
   cout << "am: " << arma_matrix[row*block_size+col] << endl;