Polar decomposition of square matrices in C++ using boost::ublas