Efficiently solve $n \times n$ linear system

53 Views Asked by At

You already have the LU factorization with partial pivoting of $A \in R^{n \times n}$ (i.e. $PA=LU$). Explain how to efficiently solve the $ n \times n $ linear system $$A^T Ax=b$$ without computing another PLU factorization.

Where do I even start??

2

There are 2 best solutions below

0
On

Solve $A^T y = b$ first, then $Ax = y$.

Well, this was intended as a hint.

More explicitly:

First solve $A^T P^T y = b$ for $y$, which is equivalent to $U^T L^T y = b$.

Then solve $A x = P^T y$ for $x$, premulitplying by $P$ gives $PA x = y$, which is equivalent to $LU x = y$.

0
On

Hint: Given that $PA=LU$, take the transpose of each side. That tells you something about $A^T$.