I have a technical question about numerical relativity, hopefully someone can help. In the usual 3 + 1
decomposition in NR, the four gauge freedoms are expressed via the lapse, and three shift components . In a finite element numerical scheme, each grid point will have a value for these 4 components. Now, looking at evolution equations on the other ...