Я ищу, чтобы найти пример решателя полос в Python, который принимает ту же форму ввода, что и функция scipy.linalg.solve_banded (или несколько советов о том, как написать свою собственную).Я не могу использовать функцию scipy, потому что я решаю матрицу с помощью числа с плавающей точкой mpmath вместо стандартных чисел с плавающей точкой двойной точности, а функции матрицы scipy не поддерживают числа с плавающей точкой mpmath.
Я не уверен, является ли Гауссово исключение и обратное замещение самым быстрым способом, но это, безусловно, возможно.