可逆方阵 A 的逆记为,A−1,需满足 AA−1=I。
在 BLAS 的各种实现中,一般都不会直接给出 matrix inverse 的直接实现,其实矩阵(方阵)的逆是可以通过 gemm()和gesvd()操作得到。
gemm()
gesvd()
实值可逆方阵 A,其 SVD 分解如下:
其中:
V,U 均为正交矩阵,
S 为对角矩阵;
所以有:
也即: