Я использую mpmath для произвольной десятичной точности. Я создаю большие квадратные матрицы (30 х 30 и 100 х 100). В моем коде я выполняю разложение по единственному значению и инверсия матрицы с использованием встроенных пакетов mpmath.
Моя проблема в том, что mpmath работает медленно, даже с серверной частью gmpy
. Мне нужна точность до 50 знаков после запятой (если решение быстрое, я предпочитаю, чтобы оно масштабировалось до большего количества знаков после запятой).
Есть ли решение для ускорения этих задач линейной алгебры в python?
Кто-то задал аналогичный вопрос здесь, но есть 2 отличия:
- Ответы не касались разложения по сингулярным значениям
- Ответы давали методы оценки обратного, но они не пытались показать, что приближение к истинному ответу быстрее, чем метод mpmath. Я попробовал решение, приведенное в этом посте, и обнаружил, что оно медленнее, чем внутренний алгоритм mpmath.