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

See also: http://en.wikipedia.org/wiki/Polar_decomposition