We present numerical evolutions of three equal mass black holes using the moving puncture approach. We calculate puncture initial data for three black holes solving the constraint equations by means of a high order multigrid elliptic solver.
Galaviz Vilchis, J., Brugmann, B., & Cao, Z. (2010). Numerical evolution of multiple black holes with accurate initial data. Physical Review D, 82(2), 024005-1 - 024005-14. https://doi.org/10.1103/PhysRevD.82.024005