Есть ли способ получить / установить факторизацию больших разреженных матриц в собственных решателях?

У меня есть программа, в которой системная матрица в разреженном решателе (использующая библиотеку Eigen) часто одинакова и может быть предварительно рассчитана один раз и загружена в виде файла. В идеале я хотел бы сериализовать весь решатель, но, кажется, нет никакого способа сделать это с Eigen.

Я сейчас пользуюсь CholmodSupernodalLLT, так что я полагаю, один из способов может быть подкласс CholmodSupernodalLLT и найти способ сериализации cholmod_factorОднако я не знаю, как сериализовать этот тип.

Я готов использовать другие решатели (предпочтительно LLT) или даже другие библиотеки, если у них есть такая возможность.