It would be nice to have the return-mapping in stress spaces with constraints other than plane-stress. This has an easy implementation by creating an abstract method that defines the C and Q matrices, and their eigen-decompositions. Then various constrained stress spaces can be implemented by just defining the matrices and decompositions.