如何在 blas 中为 trtrs(或 tptrs)进行转置?

计算科学 线性代数 线性求解器 布拉斯 矩阵
2021-11-30 15:00:33

如何在 blas 中为 trtrs(或 tptrs)进行转置?

我想解决:

XA = B

但似乎 trtrs 只能让我解决:

AX = B

或者,在 trtrs 中使用“转置”标志:

A'X = B

其中,重新排列是:

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

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

请注意,我已通读解决xA=b为了x使用 LAPACK 和 BLAS,但据我所知,它似乎没有解释如何避免两个手动转置?或者也许我以某种方式错过了一些东西?

1个回答

既然 xTRRTS 所做的基本上都是调用 xTRSM(来自 BLAS),为什么不直接调用呢?xTRSM 包含一个 SIDE 参数,可让您准确指定所需内容。