Oh, I see. I was looking at `sqrt(M)`, which for some reason was different than `M**Rational(1, 2)` (I guess Pow doesn't know about matrices).

I'm still not sure if sqrtdenest is the right function to do the job, though, but definitely simplify would do it once it works.  I think eventually this would be done by a simplification function that uses minpoly.