Multiple scattering of elastic waves in realistic media makes that average
field intensities or energy densities follow diffusive processes. In such regime the successive *P* to *S* energy conversions by distributed random inhomogeneities give rise to
equipartition which means that in the phase space the available elastic energy is distributed in average with equal amounts among the possible states of *P *and *S* waves. In
such diffusive regime the *P* to *S* energy ratio equilibrates in an universal way independent of the particular details of the scattering. It has been demonstrated that averaging
the cross correlations at any two points of an elastic medium subjected to diffuse elastic wavefields leads to the emergence of the Green function, which is the wave field
that would be observed at one position if an impulsive load is applied at the other. In
this work we study the problem of the retrieval of the 2D tensor elastodynamic Green
function in an infinite elastic space containing a circular cylinder inclusion. We illuminate isotropically the elastic space with plane waves. We assume the spectra for both *P*and *S* waves uniform but such that the energy ratio *E _{S}*/