umfpack - x = (AtA)^-1 b である UMFPACK を使用してシステムを解く

Question

疎行列 A があり、システム (AtA)X = B を解く必要があります。現在、At に A を掛けて、AtA を使用してシステムを解いています。より高速な代替手段はありますか (A を UMFPACK に渡して因数分解をより迅速に検出させるなど)?

ちなみにAは長方形です。