Implementing a new recovery scheme for primitive variables in the general relativistic magnetohydrodynamic code Spritz