A background to the constitutive modeling of elastic-inelastic material response is provided to highlight the uniqueness of the Eulerian formulation of general nonlinear fully anisotropic thermoelastic-inelastic materials proposed in Rubin (1994) [1]. This model introduced Eulerian evolution equations for a triad of microstructural vectors that characterize elastic deformations and anisotropic orientations. Components of tensors which transform like the Cauchy stress referred to these vectors are insensitive to superposed rigid body motions so they can be used to formulate general elastically and inelastically anisotropic constitutive equations. This paper develops a strongly objective, robust numerical algorithm for integrating the evolution equations for the microstructural vectors. This algorithm can easily be implemented into computer codes to simplify the use of general anisotropic constitutive equations for thermoelastic-inelastic material response.