On August 17, 2017, the LIGO/VIRGO collaboration detected the first gravitational wave signal coming from the merger of two neutron stars. This groundbreaking discovery, referred to as GW170817, revealed exciting new insights about how heavy elements, such as gold and platinum, are synthesized through a mechanism known as rapid neutron capture (r-process). In order to fully understand these signals, we need to couple observational surveys with numerical simulations of accretion disks surrounding black holes. Simulating the accretion disk and predicting the composition of the outflows is a task that requires efficient computing codes that include general relativity magnetohydrodynamics (GRMHD), neutrino physics and a model for matter at high densities. In this talk I will present the implementation of a tabulated equation of state that takes care of matter at high densities and a neutrino leakage scheme that considers the impact of neutrinos into HARM3D, a GRMHD parallelized code. I'll present tests for both the tabulated equation of state and the leakage scheme, as well as apply the tools to a magnetized torus.