Reaction calculation 4190

resolve the merge conflict...
Ref #4190 
Leave the computeqpelasticitytensor in its old form. I will move out the jacobian_mult calculation from here.