Add EquationBase.set_concentration_vector method.

- prepare empty method of descendants.
Remove Reaction.compute_one_step method.

- remove method in all descendants
- replace the update_solution method
Add Reaction.get_element_for_dof_index method.

- method finds mesh elemnt for local dof index.