Putting IntegrationMethod after GlobalDim enables automatic
template argument deduction for that type.

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.

Fixes multiple copyright messages in the doxygen docu.

ag '\\file' -L --ignore=ThirdParty *Lib |xargs  sed -i 's/\(.*\)\\copyright/\1\\file\n\1\\copyright/'

These are almost the same and the difference is in the
handling of the shape functions.

Clang format run on near-the-header lines.

Happy New Year

Mainly the shape matrices are affected.

Changing the initialization order in the local assemblers
is required.

Creating integration method object for every assembly
is not cheap because of the std::pow call. Since the
method does not change, it can be stored once per
assembler.