The residuum consists of heat, gas, and liquid mass flow rates, and nodal forces.
While the residuum for the Newton scheme is readily available, some adjustments to the central and forward differences assemblers were necessary.
The removal of almost everywhere unused matrices M and K from the assemble Jacobian set of functions is not included here and will be another MR.
Feature description was added to the changelog
Tests covering your feature were added?
Any new feature or behavior change was documented?