如何在blas中移植tptrs?

时间:2013-06-17 23:04:54

标签: c matrix fortran blas

如何在blas中移植tptrs?

我想解决:

XA = B

但似乎tptrs只能让我解决:

AX = B

或者,使用'transpose'标志,在tptrs:

A'X = B

,重新排列是:

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

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

1 个答案:

答案 0 :(得分:2)

TPTRS不是BLAS例程;这是一个LAPACK程序。

如果A与B和X相比相对较小,那么将其打包成“普通”三角矩阵并使用BLAS例程TRSM是一个很好的选择,它采用“侧”参数,允许您指定XA = B. A是mxm,B是nxm,解包增加了m ^ 2个操作,与O(nm ^ 2)操作相比,这将是一个很小的开销来进行求解。