Before, the calculation was terminated after exactly three iterations, now the stopping criterion is fulfilled if the change in the norm of the polar factor is small enough. Also: add the algorithm by Higham and Noferini, done by Robin Fränzel.