[[Imported from SVN: r10003]]

[[Imported from SVN: r10002]]

[[Imported from SVN: r10001]]

[[Imported from SVN: r10000]]

[[Imported from SVN: r9999]]

[[Imported from SVN: r9998]]

[[Imported from SVN: r9997]]

Described in:
 Christof Melcher, "Chiral skyrmions in the plane", Proc. of the Royal Society, online DOI DOI: 10.1098/rspa.2014.0394

[[Imported from SVN: r9996]]

[[Imported from SVN: r9995]]

[[Imported from SVN: r9994]]

    
This removes a compiler deprecation warning.  By using the usual preprocessor magic
we retain compatibility with the dune 2.3 release.

[[Imported from SVN: r9993]]

[[Imported from SVN: r9992]]

There is no need to have our own implementation for this, we had it purely for
historical reasons.

Also, this change fixes a subtle bug:  Our own implementation used AlienElementLocalBasis,
which returns zero if the LocalBasis is evaluated outside of the reference element.
However, this happened sometimes due to rounding error on grids with very small quadrilateral
elements.

[[Imported from SVN: r9991]]

[[Imported from SVN: r9990]]

[[Imported from SVN: r9989]]

[[Imported from SVN: r9988]]

[[Imported from SVN: r9987]]

[[Imported from SVN: r9986]]

[[Imported from SVN: r9985]]

[[Imported from SVN: r9984]]

[[Imported from SVN: r9983]]

[[Imported from SVN: r9982]]

[[Imported from SVN: r9981]]

[[Imported from SVN: r9980]]

[[Imported from SVN: r9979]]

[[Imported from SVN: r9978]]

[[Imported from SVN: r9977]]

This already saves a lot of disk space, because the class issues less whitespace
than our own code.  Plus, switching to binary data will be very easy now.

[[Imported from SVN: r9976]]

[[Imported from SVN: r9975]]

[[Imported from SVN: r9974]]

[[Imported from SVN: r9973]]

[[Imported from SVN: r9972]]

[[Imported from SVN: r9971]]

[[Imported from SVN: r9970]]

[[Imported from SVN: r9969]]

[[Imported from SVN: r9968]]

[[Imported from SVN: r9967]]

[[Imported from SVN: r9966]]

[[Imported from SVN: r9965]]

Previously, the number of multigrid transfer operators was always equal to
the number of grid levels minus one.  That worked out correctly when doing
first-order FE spaces.  However, for higher-order spaces, the up-most transfer
operator transfers from the higher-order space to a first-order space
_on the same grid level_.  But since the total number of transfer operators
was grid levels minus one, this implied that the coarsest grid leve 
was never used, and the coarse grid solver would operate on the second-coarsest
grid.

This patch adds the extra transfer operator in this case.  That means we now
call the coarse grid solver on smaller problems, which makes a noticeable
run-time difference.

[[Imported from SVN: r9964]]