An implicit numerical algorithm for solving the general relativistic hydrodynamical equations around accreting compact objects

The results of several test calculations that verify the accuracy and robustness of the algorithm are shown. Extending the algorithm to enable solving the non-ideal MHD equations in the general relativistic regime is the subject of our ongoing research.