The Harm3D code solves the relativistic MHD equations on arbitrary dynamical spacetimes in arbitrary coordinate systems.

LazEv is the in-house numerical relativity code for binary black hole simulations. It's based on the "moving puncture approach" which produced a breakthrough in numerical relativity in 2005.