  We incorporate realistic, tabulated equations of state into fully
relativistic simulations of magnetized neutron stars along with a neutrino
leakage scheme which accounts for cooling via neutrino emission. Both these
improvements utilize open-source code (GR1D) and tables from
http://stellarcollapse.org. Our implementation makes use of a novel method for
the calculation of the optical depth which simplifies its use with distributed
adaptive mesh refinement. We present various tests with and without
magnetization and preliminary results both from single stars and from the
merger of a binary system.
