Мне нужно решить систему линейных уравнений Lx = b, где x всегда является вектором (массив 3x1), L является массивом Nx3, а b является вектором Nx1. N обычно колеблется от 4 до примерно 10. У меня нет проблем с решением этого с помощью

scipy.linalg.lstsq (L, б)

Тем не менее, мне нужно делать это много раз (что-то вроде 200x200 = 40000 раз), так как x на самом деле что-то связанное с каждым пикселем в изображении. Таким образом, x фактически хранится в массиве PxQx3, где P и Q примерно 200-300, а последнее число «3» относится к вектору x. Прямо сейчас я просто перебираю каждый столбец и строку и решаю уравнение один за другим. Как эффективно решить эти 40000 различных переопределенных систем линейных уравнений, возможно, используя некоторые методы векторизации или другие специальные методы?

Спасибо