Fast numerics for the spin orbit equation with realistic tidal dissipation and constant eccentricity