Reaction calculation 4190

resolve the merge conflict...
Ref #4190 
@cpritam But it may become invalid when some materials override computeQpElasticityTensor. Could it be possible to rotate it inversely to this rotation in crystal plasticity and rotate it again to intermediate configuration? We don't have to write the getEulerRotation function because we can call .rotate(R.transpose()) to accomplish it.