If then . The idea is to write P in the form , where Q is orthogonal and R is a positive definite matrix that commutes with D. Then , from which .
The factorisation P = QR is called the polar decomposition of P. The positive definite factor R is the square root of . So start by calculating . This has square root . (Fortunately, this commutes with D, as you can easily check.) Then .