Решение системы с помощью UMFPACK, где x = (AtA)^-1 b
У меня есть разреженная матрица A, и мне нужно решить систему (AtA)X = B. Сейчас я умножаю At на A и использую AtA для решения системы. Есть ли более быстрая альтернатива (например, передача A в UMFPACK и ускорение поиска факторизации)?
Кстати, прямоугольник.