Блог

У меня есть k треугольные матрицы размера d x d . Для каждого из них я хочу решить для n векторов (размерности d ). Насколько я знаю, numpy.linalg.solve это позволяет мне, используя широковещательную передачу, выполнять все решения одновременно. Однако он использует Fortran GESV и не может использовать треугольную структуру. С другой стороны, scipy.linalg.solve_triangular кажется, что он способен обрабатывать треугольные матрицы, но не поддерживает несколько одновременных решений.

Есть ли способ достичь обоих?