Here we report a developed high performance and simplified version of the code denominated rio, which can be easily extended, for the generalized BSSN formulation. We implement a code which is regular at the center of symmetry, without use a special procedure for regularization, as usual. We get exponential convergence for constraints. The numerical algorithm is based on the Galerkin-Collocation method developed successfully for diverse physical scenarios by the Numerical Relativity Group at uerj. For the sake of clarity in presentation, we consider here the most simple case to display the most salient features of the procedure. Thus, we focus on the definite tests of the new numerical framework. The timing and performance of the code show that we can reach a better accuracy close to the machine precision, for the Hamiltonian and momentum constraints. rio will be an open source code; currently it is under continuous development to consider more general and realistic problems.

• Received 19 May 2021
• Revised 20 August 2021
• Accepted 28 September 2021

DOI:https://doi.org/10.1103/PhysRevD.104.084065