Einstein's theory of gravitation, general relativity, is a purely geometric theory.
In a recent question I wanted to know what the relation of Brownian motion to the Helmholtz equation is and got a very thorough answer from George Lowther.
He pointed out that there is, roughly speaking, a very general relation of semi-elliptic second order differential operators of the form
$$Af = \frac12 a^{ij}f_{,ij} + b^i f_{,i} - cf = 0$$
to a "killed" Brownian motion. (I used some summation convention and $,i = \frac{\partial}{\partial x^i}$.)
Now, the Einstein field equations
$$R_{\mu\nu}-\frac12 g_{\mu\nu}R = \frac{8\pi G}{c^4}T_{\mu\nu}$$
are coupled hyperbolic-elliptic partial differential equations (I dropped the cosmological constant here). Can we somehow adopt the relation of a random process to this kind of equation, or
Is there a way to interprete the Einstein equations stochastically?