Numerical implementation of a $$J_2$$- and $$J_3$$-dependent plasticity model based on a spectral decomposition of the stress deviator