Hi DruckerPrager3D::_calcStressElastoplastic has the following line const PylithScalar d = sqrt(ae * ae * devStressInitialProd + 2.0 * ae * This eqn in the manual, Eq.(5.110), has a factor of 2 in the 1st term. Thanks Birendra