[Tests][PkFEM] Extend Pk 2D interpolation test to also test for k = {2,3};

Due to that functionality extension, I also renamed the test to reflect
the fact that multiple values of k are tested.
[CMake] Fix headercheck issues with CMake.

Ignore Eigen backend headers if Eigen not found.
[Tests] Fix array initializer in testordering
[Tests] Improve readability of closing #endif by adding comment
[CMake] Add "GRIDDIM 2" to Alberta flags.

Fixes configure warnings.
[Headercheck] Fix headercheck by adding missing includes.
[FEM] Add a new exception type for lookups with invalid GeometryType
[PkFEM] Make sure we are really passed a simplex cell in find()
[Communication] Change DataHandleProvider and EntityIndexCache to keep per-leaf size information

DataHandleProvider and EntityIndexCache cooperate using an informal
protocol to calculate the per-entity DOF count. The DataHandleProvider
obtains this information by visiting the leaves of the Ordering tree and
accumulating the size information at those leaf Orderings.

While the accumulated total is sufficient in the standard PDELab case,
MultiDomain requires the per-leaf information for certain combinations
of subdomain and processor boundaries, in particular if a subdomain
exists on the cells on one side of the processor boundary but not on the
other and the user has chosen a discretization with DOFs on entities
with codim > 0. In this case, the per-entity size information on the
different ranks can become inconsistent.

This patch changes the accumulation process to store the per-leaf sizes
in an array that is exported via the new method offsets() of the
EntityIndexCache.
[Communication] Add PolymorphicBufferWrapper to allow sending heterogeneous data in a single grid communication

PolymorphicBufferWrapper takes a grid MessageBuffer for data of type
char and allows you to write/read data with arbitrary POD types.
[Communication] Add support for communicating entities with partially inconsistent DOF counts

This patch is required for supporting parallel computations with
MultiDomain, where it can be possible for entities of codim > 0 to have
different numbers of DOFs attached to them on different MPI ranks.

The problem arises if an entity belongs to a subdomain, but there is no
associated cell that also belongs to the subdomain on one MPI rank. In
that case, the GridFunctionSpace will not now about this membership and
thus consider the subspace assoicated with the subdomain to be empty,
causing a crash when too much / too few data is sent / received.

In order to fix this problem, there is now a new method sendLeafSizes()
on the DataHandleProvider. If that method returns true, the generic data
handle infrastructure will automatically communicate per-LeafOrdering
sizes and only enforce sizes to be consistent for those leaf
orderings. By default, this method always returns false, so most
existing applications should not even notice the difference between the
old and the new implementation.

Please note that this change requires wrapping the message buffer passed
in by the grid, as we now have to send heterogeneous data (e.g. sizes
and doubles) over the same data handle. This is not supported by the
standard data handles, so we use a DataHandle for char and marshal
everything to a simple byte stream. However, these changes only apply to
per-DOF communication. Per-entity communication does not require the
additional overhead, as there is no ambiguity as to the subdomain
membership of entities at the grid level.
[Cleanup] Remove unused typedefs and variables in constraints-free implementation of scatter_jacobian()
[release,boilerplate] Update links in the documentation to the new location of tutorial.

This is thorrowly needed as there is a lack of source code documentation.
[autotools] Distribute dune/pdelab/test(/fem)/grids with tarball.

The symbolic link in fem is not included, instead the folder grids
is included a second time. Maybe there is a better solution.
Merge branch 'feature/fix-parallel-multidomain' into releases/2.0

Extend PDELab with required support code to allow parallel computations
in MultiDomain again.

* feature/fix-parallel-multidomain:
  [Communication] Add support for communicating entities with partially inconsistent DOF counts
  [Communication] Add PolymorphicBufferWrapper to allow sending heterogeneous data in a single grid communication
  [Communication] Change DataHandleProvider and EntityIndexCache to keep per-leaf size information
[VariableMonomFEM] Check whether bound entities have correct GeometryType

VariableMonomFEM only supports a single GeometryType for grid cells,
even though the fact that you can pass in an arbitrary mapper seems to
suggest otherwise; the returned finite elements are only correct for
the GeometryType passed to the FEM's constructor or the default
GeometryType, which is a cube.

Unfortunately, the FEM did not check whether the grid cell passed to
find() actually had the correct GeometryType, leading to weird, subtle
bugs in the ordering infrastructure later on. This was noticed by
Christoph Grüninger, who ran testordering with assertions turned on,
promptly tripping over a check to make sure that the current
GeometryType is actually in use with the Ordering. Thanks for reporting,
Christoph!

This was really nasty to debug - it's fixed now by checking the
GeometryType in find(), but it's still possible to create a wrong FEM
for a given grid. After the release, we should change the FEM to simply
take the GridView in its constructor and automate the whole process,
maybe even adding support for different GeometryTypes. But for now, the
quick fix will have to do.
[Documentation][Cleanup] Re-word updated example notice in pdelab.hh to avoid compiler warning

The updated hint on where to find examples for the boilerplate layer in
boilerplate/pdelab.hh unfortunately contained the character sequence
"/*" as part of the path to the examples, triggering a compiler warning
about nested comments.

Reworded to avoid the problematic path name.
[Testordering] Initialize VariableMonomFEM with the correct GeometryType

VariableMonomFEM for a simplex ALUGrid was initialized with the default
constructor, which assumes a hexahedral grid, creating subtle
inconsistencies in the ordering data structures. These only caused a
test failure when compiled with assertions, which is why the problem
went unnoticed for that long.

Fixed by explicitly calling the correct constructor. Thanks to Christoph
Grüninger for reporting the original issue.
[CMake] Link FEM tests against DUNE libraries
[ConvectionDiffusionDG] Silence unused variable warnings