of the entries of the global residual vector that do not need
 initial non-equilibrium compensation

Use const references where necessary.

Rename a the method from getNumberOfComponents() to getNumberOfGlobalComponents()
in class LocalToGlobalIndexMap.

Interface and time loop changes. Processes are in the
next commit.

Main reason for modification is the getFlux function, which for
particular processes (ComponentTransport) relies on all solutions being
available.
This modification ripples through all of the code, especially through
the boundary conditions.
Modified `applyNaturalBC`, `preTimestep`, `integrate`, and `assemble`
functions.

preTimestep now accepts a vector of solution together with
a process_id.
In postTimestep only the vector of solutions has been changed.