c - 如何在 blas 中为 tptrs 进行转置？

Question

如何在 blas 中为 tptrs 进行转置？

我想解决：

XA = B

但似乎 tptrs 只能让我解决：

AX = B

或者，在 tptrs 中使用“转置”标志：

A'X = B

其中，重新排列是：

(A'X)' = B'
X'A = B'

所以，我可以用它来解决XA = B，但我必须先自己手动转置B，然后再转置答案。我是否错过了一些避免转置的技巧？

Answer 1

TPTRS 不是 BLAS 例程；这是一个LAPACK例程。

如果 A 与 B 和 X 相比相对较小，那么将其解压缩为“正常”三角矩阵并使用 BLAS 例程 TRSM 是一个不错的选择，该例程采用“side”参数允许您指定 XA = B。如果 A 是 mxm并且 B 是 nxm，解包增加了 m^2 操作，与解决问题的 O(nm^2) 操作相比，这将是少量的开销。