Most formulations of Varadhan's formula
$$\lim _{t \to 0_+} 4t \log p_t(x,y) = -d(x,y)^2$$
that I have encountered do not specify where $(x,y)$ lives, so until today I imagined that $(x,y) \in M \times M$. It seems that this is not the case, though, since Berger in "A Panoramic View of Riemannian Geometry" explicitly states (Theorem 168, page 424) that the formula is valid only for $x$ within the injectivity radius at $y$, and that outside this subset wild things happen, citing "Short Time Behaviour of the Heat kernel and its Asymptotic Derivatives" by Malliavin and Stroock. So I found the aforementioned article and, surprise!, the authors claim that Varadhan's formula holds for all $(x,y) \in M \times M$ (formula 0.3, second page of the article)!
Now, I am completely confused. Could you please clarify for me what is the correct, rigorous, formal statement of Varadhan's formula on Riemannian manifolds? (To make things worse, if I recall correctly, Varadhan's result was originally obtained for open subsets of $\mathbb R^n$.)
(Let me add that the heat equation differs by a factor of $\frac 1 2$ in front of the Laplacian between the two references above, but the versions of Varadhan's formula are identical, which makes me believe that Berger is at fault here.)
On a funnier note, is differential geometry the branch with the sloppiest texts in mathematics, or can one get even lower? I'm having troubles with a mistake in texts by Chavel and by Donnely, and I'm amazed about how these guys claim (false) things without bothering to prove them (which would detect their mistakes before publication). I am even more amazed at the referees who let these pass...