umfpack - 使用 UMFPACK 求解系统，其中 x = (AtA)^-1 b

Question

我有一个稀疏矩阵 A，需要求解系统 (AtA)X = B。现在我将 At 乘以 A 并使用 AtA 来求解系统。是否有更快的替代方案（例如将 A 传递给 UMFPACK 并让它更快地找到分解）？

顺便说一下，A是矩形的。