如何在blas中移植tptrs?
我想解决:
XA = B
但似乎tptrs只能让我解决:
AX = B
或者,使用'transpose'标志,在tptrs:
A'X = B
,重新排列是:
(A'X)' = B'
X'A = B'
所以,我可以用它来解决XA = B,但我必须先自己手动转置B,然后再转换答案。我错过了一些技巧以避免不得不进行转置吗?
答案 0 :(得分:2)
TPTRS不是BLAS例程;这是一个LAPACK程序。
如果A与B和X相比相对较小,那么将其打包成“普通”三角矩阵并使用BLAS例程TRSM是一个很好的选择,它采用“侧”参数,允许您指定XA = B. A是mxm,B是nxm,解包增加了m ^ 2个操作,与O(nm ^ 2)操作相比,这将是一个很小的开销来进行求解。