[PL/CT] solve reaction equation
In the current master, the update of the nodal concentration is through extrapolation from integration points to mesh nodes in the reaction stage. In this way, the concentration value at the nodes where the Dirichlet boundary condition is imposed cannot be kept as prescribed. Instead, the nodal concentration is updated by solving the discretized reaction equation in such a way that the Dirichlet boundary condition is forcibly satisfied.
Feature description was added to the changelog
Tests covering your feature were added?
Any new feature or behavior change was documented?