У меня есть вопрос, как решить линейную систему с помощью броненосца. А*х=б. Но моя матрица А довольно большая. Это arma::cx_mat A(40000,40000). Когда я запускаю x=arma::solve(A,b), броненосец попытается выделить новую память для хранения A, а затем использовать новую память для декомпозиции LU и не будет использовать память, в которой хранится исходный A. определение решения можно найти в файле fn_solve.hpp. Но вы уже видите, что моя пятерка довольно большая. Выделение новой памяти быстро израсходует всю память. Для моей проблемы мне больше не нужен оригинальный A. Итак, могу ли я использовать какой-либо метод, чтобы позволить броненосцу использовать исходную память для A для декомпозиции LU?
Большое тебе спасибо
