Как получить правильные матрицы L и U из LU-разложения разреженной матрицы A, не используя scipy.sparse.linalg.splu ()? - PullRequest
0 голосов
/ 08 декабря 2018

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

...