As demonstrated in the affixed file, there are large differences between the expected deformation gradient (DG) and the DG printed by VUMAT in 2D simple shear example. The differences are especially obvious for the first several increments. The reason is unknown. I log the phenomenon here hoping for some explanations some day.