使用UMFPACK求解系统，其中x =（AtA）^ - 1 b

Question

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

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