如何在blas中对tptrs进行转置?
我想解决:

XA = B

但似乎tptrs只让我解决:
AX = B

或者,在tptrs中使用“transpose”标志:
A'X = B

其中,重新排列是:
(A'X)' = B'
X'A = B'

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

最佳答案

TPTRS不是BLAS程序,而是LAPACK程序。
如果A与B和X相比相对较小,则可以将其解压为“正常”三角形矩阵,并使用BLAS例程TRSM,该例程带有“side”参数,允许您指定XA=B。如果A是m X m,B是nxm,则解压将添加m^2操作,这与执行解决。

10-01 14:00