LazyTB: Extended LTB-solutions with Multiple Interacting Fluids

We present a low entry barrier numerical code LazyTB that solves the Einstein equations in the ADM-formulation for a general spherically symmetric problem with an arbitrary number of fluids with user defined equations of state and optional interactions between the fluids mediated by drag-vectors. We present the ADM equations in terms of carefully chosen dimensionless variables that avoid some numerical instabilities by analytic cancellation of potentially diverging terms. We demonstrate the code by applying it to a number of examples including the standard LTB model, a model with radiation and…