三角矩阵的广播求解

Question

我有 k 个大小为 d x d 的三角矩阵。对于每一个，我想求解 n 向量（维度为 d）。 据我所知，numpy.linalg.solve 允许我使用广播同时执行所有求解。但是，它使用 Fortran 的 GESV 并且不能利用三角形结构。 另一方面，scipy.linalg.solve_triangular 似乎能够处理三角矩阵，但不支持多个同时求解。

有没有办法同时实现？