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

Question

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

我要解决:

XA = B

但似乎 tptrs 只能让我解决:

AX = B

或者，在 tptrs 中使用 'transpose' 标志:

A'X = B

其中，重新排列是:

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

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

Answer 1

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

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