ENH: linalg: add a function to compute the polar decomposition of a matrix.