NLopt с данными броненосца

Целевая функция NLopt выглядит следующим образом:

double myfunc(const std::vector<double> &x, std::vector<double> &grad, void *my_func_data)

x - данные, которые оптимизируются, grad - вектор градиентов, а my_func_data содержит дополнительные данные.

Я заинтересован в предоставлении матриц Армадилло A и B для аннулирования * my_func_data.

Я возился с функциями члена Armadillo

mat A(5,5);
mat B(5,5);
double* A_mem = A.memptr();
double* B_mem = B.memptr();

который дает мне указатели на матрицы A и B. Я думал об определении другого указателя на эти указатели:

double** CombineMat;
int* Arow = A.n_rows; int* Acols = A.n_cols; //obtain dimensions of A
int* Brows = B.n_rows; int* Bcols = B.n_cols; // dim(B)
CombineMat[0] = A_mem; CombineMat[1] = Arows; CombineMat[2] = Acols;
CombineMat[3] = B_mem; CombineMat[4] = Brows; CombineMat[5] = Bcols;

а затем передать * CombineMat как my_func_data.

Это способ сделать это? Кажется неуклюжим ...После того, как CombineMat пройден, как преобразовать тип void во что-то пригодное для использования, когда я нахожусь внутри myfunc?

ОТВЕТ

Я ответил на свой вопрос с помощьюВот.

mat A(2,2);
A << 1 << 2 << endr << 3 << 4;

mat B(2,2);
B << 5 << 6 << endr << 7 << 8;

mat C[2];
C[0] = A;
C[1] = B;

opt.set_min_objective(myfunc, &C);

Оказавшись внутри myfunc, данные в C могут быть преобразованы обратно в матрицы Armadillo следующим образом:

mat* pC = (mat*)(my_func_data);
mat A = pC[0];
mat B = pC[1];