Estou tentando resolverAx = b onde a matriz A pode ser grande perto de1M x 1M em tamanho, é escasso e simétrico, mas pode não ser definido positivamente.

O problema é que pode levar muito tempo para calcular a decomposição usando o métodoobjeto sparseLU em eigen e será uma ideia armazenar a matriz sparseLU em vez da matriz original, de modo que sempre que executarmos operações semelhantes usando a mesma matriz A, possamos acelerar as coisas, não precisando recalcular o valor

Uma rápida pesquisa no stackoverflow e no google retornouesta, esta eesta para matriz esparsa para serialização da matriz eigen. No entanto, não tenho certeza se o mesmo código pode ser aplicado ao objeto sparseLU.

Talvez eu deva reformular minha pergunta:

Como posso armazenar a matriz decomposta em um arquivo?

Todos os métodos atuais se concentram no armazenamento da matriz original, mas eu quero armazenar a matriz decomposta. Existe alguma maneira de fazer isso? Obrigado.