The Benefits of Using a Consistent Tangent Operator for Viscoelastoplastic Computations in Geodynamics