如何在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操作,这与执行解决。