Я провожу некоторые численные моделирования и требую решения больших линейных систем. Производные операторы редки (5-точечный шаблон для оператора Лапласа), и поэтому я использую sparse.linalg.spsolve
, чтобы получить решение. Стандартный scipy
универсальный решатель linalg.solve
дает возможность предположить, что матрица обладает определенными свойствами, например, симметричным положительно определенным (мой случай). Эта информация может быть использована для выбора лучшего решателя для данного случая.
spsolve
, однако, по-видимому, отсутствует эта опция, и мне было интересно, если есть способ воспользоваться моими свойствами матрицы. Я видел, что другие сталкивались с подобными проблемами , но пока не получили ответа.
Заранее спасибо