So, as I am sure many of you know, many matrices do not have unique square roots. Even something so simple as the 2x2 identity matrix has infinite results.
So, my question is: When a matrix is put into Mathematica, and we take the square root, does anyone know how Mathematica determines which of these infinite answers to use?
Or, better yet, in the interest of the Lagrangian in the theory I'm working with, it would be even better if I could force Mathematica to return the positive-definite result specifically. Does anyone know how one might go about this?
Any documentation or insight would be greatly appreciated. Thank you.