[[Imported from SVN: r8326]]

[[Imported from SVN: r8310]]

Use the long form for the quadratic energy stiffness.  I tested it: when kappa=1 is does the same as the short form

[[Imported from SVN: r8308]]

[[Imported from SVN: r8307]]

[[Imported from SVN: r8306]]

[[Imported from SVN: r8305]]

[[Imported from SVN: r8304]]

[[Imported from SVN: r8303]]

to set the shear correction factor.

[[Imported from SVN: r8290]]

[[Imported from SVN: r8277]]

[[Imported from SVN: r8276]]

[[Imported from SVN: r8261]]

[[Imported from SVN: r8260]]

[[Imported from SVN: r8259]]

[[Imported from SVN: r8258]]

[[Imported from SVN: r8257]]

choose the correct quadrature order for the matrices that create the norms to measure the discretization error

[[Imported from SVN: r8243]]

[[Imported from SVN: r8242]]

[[Imported from SVN: r8241]]

[[Imported from SVN: r8240]]

store all local dofs for all elements.

This leads to some redundancy, but it scales to 
higher polynomial orders.

[[Imported from SVN: r8239]]

[[Imported from SVN: r8238]]

[[Imported from SVN: r8237]]

[[Imported from SVN: r8236]]

[[Imported from SVN: r8235]]

[[Imported from SVN: r8234]]

[[Imported from SVN: r8233]]

[[Imported from SVN: r8221]]

[[Imported from SVN: r8220]]

[[Imported from SVN: r8219]]

from a Pk-basis to a P1-basis.

Doesn't work beyond P2 yet, I'll fix that in a second.

Should eventually go into dune-solvers.

[[Imported from SVN: r8218]]

[[Imported from SVN: r8217]]

[[Imported from SVN: r8214]]

[[Imported from SVN: r8213]]

[[Imported from SVN: r8212]]

[[Imported from SVN: r8211]]

[[Imported from SVN: r8210]]

[[Imported from SVN: r8209]]

[[Imported from SVN: r8208]]

Quaternion class to the Rotation class.  Rotation axes and
angles are a genuine rotation concept and don't really
belong into a quaternion implementation.

[[Imported from SVN: r8207]]