The situation is not at all different from that in Newtonian physics. If you have Newton's law for gravity, i.e.
$$ \vec F = -\frac{G m_1 m_2}{\vert\vec x_1 - \vec x_2\vert^3} (\vec x_1 - \vec x_2) $$
in order to exactly describe any motion, in principle you need to know some initial configuration in which you know the exact position of any object in any galaxy. But in practical calculations it is enough to know the state of some subsystem which interacts only weakly with the rest. Like, for example, just consider the motion of one planet around the sun. Or more accurately the motion of all planets in the solar system, but ignoring other stars in the galaxy.
The situation is exactly the same in GR. You cannot measure all positions of each object (and if you could the equations would be far to complicated to ever be solved), but you can consider subsystems.
Common solutions are for example the Schwarzschild solution for a spherical object, or the Friedmann–Lemaître–Robertson–Walker solution, which models cosmology.