Pyter said:
how the (at least initial) observations of the masses enters the picture, as in the solar system scenario I've brought up?
Assuming you've read post #70, here goes.
To analyze the solar system, we make a number of restrictive assumptions to reduce the number of unknowns. In the simplest case, where we consider the Sun to be the only massive body and the planets to be test objects, and we assume the Sun to be non-rotating and therefore spherically symmetric, we have already reduced the number of unknown functions in the metric to two, and since we are assuming vacuum outside the Sun itself (and we're not trying to solve for the Sun's interior), we have
no unknown functions at all in the SET, since all of its components are just zero. And we already know that in this case, there is just
one unique solution of the EFE, the Schwarzschild solution, in which the metric components are just functions of ##r##, and there is really only
one such function that has to be determined by solving the EFE, which turns out to be ##1 - 2M / r## (with ##M## the Sun's mass), since the same function appears in both ##g_{tt}## and ##g_{rr}## (just in the denominator of the latter). (The theorem that proves this is called Birkhoff's theorem, and the proof is discussed in the Insights article on it that I referenced.)
Notice that nowhere in any of this did any "initial conditions" appear. That's because this solution is
static, so there is no need for any "initial conditions" for the solution; it's the same at all times. The initial conditions for the motion of the planets is only needed to solve the geodesic equation for the orbits of each planet, not to solve the EFE itself. (I already pointed this out to @cianfa in an earlier post.)