Я заметил, что scipy.sparse.linalg.splu()
не позволяет мне разложить разреженную матрицу A на правильную матрицу L и U, которую я могу вызвать отдельно.Команда «просто» позволяет мне разложить матрицу и восстановить ее позже, используя матрицу перестановок.Однако для моего кода мне нужно разложить разреженную матрицу A на разреженную матрицу L и U и затем иметь возможность вызывать матрицы L и U отдельно (без матриц перестановок и т. Д.).Это не работает при использовании команды scipy.sparse.linalg.splu()
.Я мог бы использовать scipy.linalg.lu()
, но я не могу применить это к матрице A в разреженном формате.Существуют ли другие способы получения правильных матриц разложения L и U из разреженной матрицы A?Заранее спасибо.