Add Schur and Hessenberg functions and a matrix square root. And some tests.