Неявная матрица DAE Mass Python

Question

У меня есть проблема M * y '= f (y), которая будет решена в Python, где M - матрица масс, y' - производная, а y - вектор, такой, что y1 , y2 и т. д. относится к разным точкам в r. Кто-нибудь использовал массовую матрицу на аналогичной проблеме в Python?

Проблема представляет собой двумерную задачу в r- и z-направлении. R-направление дискретизировано, чтобы уменьшить проблему до 1D-задачи. Матричная матрица является диагональной матрицей с единицами и нулями на диагонали.

Answer 1

Поскольку матрица масс является сингулярной, это «дифференциально-алгебраическое уравнение». Вы можете найти готовые решения для DAE, такие как IDA-решатель из библиотеки SUNDIALS. SUNDIALS имеет привязки python в пакете scikit.odes.