In these notes, we consider the problem of finding the logarithm or the square root of a real matrix. It is known that for every real n x n matrix, A, if no real eigenvalue of A is negative or zero, then A has a real logarithm, that is, there is a real matrix, X, such that e^X = A. Furthermore, if the eigenvalues, xi, of X satisfy the property -pi < Im(xi) < pi, then X is unique. It is also known that under the same condition every real n x n matrix, A, has a real square root, that is, there is a real matrix, X, such that X^2 = A. Moreover, if the eigenvalues, rho e^{i theta}, of X satisfy the condition -pi/2 < theta < pi/2, then X is unique. These theorems are the theoretical basis for various numerical methods for exponentiating a matrix or for computing its logarithm using a method known as scaling and squaring (resp. inverse scaling and squaring). Such methods play an important role in the log-Euclidean framework due to Arsigny, Fillard, Pennec and Ayache and its applications to medical imaging. Actually, there is a necessary and sufficient condition for a real matrix to have a real logarithm (or a real square root) but it is fairly subtle as it involves the parity of the number of Jordan blocks associated with negative eigenvalues. As far as I know, with the exception of Highams recent book, proofs of these results are scattered in the literature and it is not easy to locate them. Moreover, Highams excellent book assumes a certain level of background in linear algebra that readers interested in the topics of this paper may not possess so we feel that a more elementary presentation might be a valuable supplement to Higham. In these notes, I present a unified exposition of these results and give more direct proofs of some of them using the Real Jordan Form.
In the computation of Feynman integrals which evaluate to multiple polylogarithms one encounters quite often square roots. To express the Feynman integral in terms of multiple polylogarithms, one seeks a transformation of variables, which rationalizes the square roots. In this paper, we give an algorithm for rationalizing roots. The algorithm is applicable whenever the algebraic hypersurface associated with the root has a point of multiplicity $(d-1)$, where $d$ is the degree of the algebraic hypersurface. We show that one can use the algorithm iteratively to rationalize multiple roots simultaneously. Several examples from high energy physics are discussed.
In this paper we give a new proof of the ELSV formula. First, we refine an argument of Okounkov and Pandharipande in order to prove (quasi-)polynomiality of Hurwitz numbers without using the ELSV formula (the only way to do that before used the ELSV formula). Then, using this polynomiality we give a new proof of the Bouchard-Mari~no conjecture. After that, using the correspondence between the Givental group action and the topological recursion coming from matrix models, we prove the equivalence of the Bouchard-Mari~no conjecture and the ELSV formula (it is a refinement of an argument by Eynard).